changeset 42:8d30aeacb52a draft

Uploaded
author greg
date Wed, 08 Feb 2017 12:57:32 -0500
parents 844033a58b65
children 181076e82999
files multigps.xml
diffstat 1 files changed, 40 insertions(+), 32 deletions(-) [+]
line wrap: on
line diff
--- a/multigps.xml	Wed Dec 14 11:31:49 2016 -0500
+++ b/multigps.xml	Wed Feb 08 12:57:32 2017 -0500
@@ -1,27 +1,29 @@
-<tool id="multigps" name="MultiGPS" version="0.72.0">
+<tool id="multigps" name="MultiGPS" version="0.73">
     <description>analyzes collections of multi-condition ChIP-seq data</description>
     <macros>
         <import>macros.xml</import>
     </macros>
     <requirements>
-        <requirement type="package" version="0.72">multigps</requirement>
+        <requirement type="package" version="0.73">multigps</requirement>
     </requirements>
     <command detect_errors="aggressive">
         <![CDATA[
+            #set output_dir = $output_html.files_path
+            mkdir -p $output_dir &&
             multigps
             ############################
             ## General options
             ############################
-            --expt "$expt"
-            --format "$expt.ext"
-            #if $ctrl is not None:
-                --ctrl "$ctrl"
+            --expt '$expt'
+            --format $expt.ext
+            #if str($ctrl) != 'None':
+                --ctrl '$ctrl'
             #end if
             #if str($verbose) == 'yes':
                 --verbose
             #end if
-            --threads="\${GALAXY_SLOTS:-4}"
-            --geninfo "$chromInfo"
+            --threads=\${GALAXY_SLOTS:-4}
+            --geninfo '$chromInfo'
             ############################
             ## Advanced options
             ############################
@@ -33,17 +35,17 @@
                 #set rloc = $aoc.reads_limits_options_cond
                 #set sdc = $aoc.scale_data_cond
                 #set umc = $aoc.use_motif_cond
-                #if str($umc.use_motif) == "yes":
+                #if str($umc.use_motif) == 'yes':
                     #set rgc = $umc.reference_genome_cond
-                    #if str($rgc.reference_genome_source) == "cached":
+                    #if str($rgc.reference_genome_source) == 'cached':
                         #set seq_dir = os.path.split($rgc.reference_genome.fields.path)[0]
                     #else:
                         ## MultiGPS requires a directory containing reference files, so symlink the history dataset.
                         #import os
                         #import tempfile
-                        #set seq_dir = tempfile.mkdtemp(prefix="tmp-multigps-seq-dir")
+                        #set seq_dir = tempfile.mkdtemp(prefix='tmp-multigps-seq-dir')
                         #set seq_file = str($rgc.reference_genome)
-                        #set tmp_filename = "%s.fa" % str($rgc.reference_genome.dbkey)
+                        #set tmp_filename = '%s.fa' % str($rgc.reference_genome.dbkey)
                         #set tmp_seq_file = os.path.join($seq_dir, $tmp_filename)
                         ln -f -s $tmp_seq_file $seq_file &&
                     #end if
@@ -91,7 +93,7 @@
                 ## Running MultiGPS
                 ############################
                 #if str($aoc.readdistributionfile) != 'None':
-                    --d "$aoc.readdistributionfile"
+                    --d '$aoc.readdistributionfile'
                 #end if
                 --r $aoc.maxtrainingrounds
                 #if str($aoc.nomodelupdate) == 'no':
@@ -122,7 +124,7 @@
                     --mlconfignotshared
                 #end if
                 #if str($aoc.exclude) != 'None':
-                    --exclude "$aoc.exclude"
+                    --exclude '$aoc.exclude'
                 #end if
                 ############################
                 ## MultiGPS priors
@@ -164,15 +166,17 @@
                     --diffp $rbec.diffp
                 #end if
             #end if
-            --out "$output_html.files_path"
-            > "$output_html"
+            --out '$output_html.files_path'
+            && mv $output_dir/*.html $output_html
+            && mv $output_dir/*.table $all_events_table
+            && mv $output_dir/*.counts $replicates_counts
         ]]>
     </command>
     <inputs>
         <param name="expt" type="data" format="bam,bed,scidx" label="Run MultiGPS on">
             <validator type="unspecified_build" />
         </param>
-        <param name="ctrl" type="data" format="bam,bed,scidx" optional="True" label="Optional file containing file containing reads from a control experiment" help="Must be same forat as the input above" />
+        <param name="ctrl" type="data" format="bam,bed,scidx" optional="True" label="Optional file containing reads from a control experiment" help="Must be same forat as the input above" />
         <param name="verbose" type="select" label="Generate intermediate files and extra output?" help="Verbose flag">
             <option value="no" selected="True">No</option>
             <option value="yes">Yes</option>
@@ -411,23 +415,27 @@
         -->
     </inputs>
     <outputs>
-        <data name="output_html" format="html"/>
+        <data name="output_html" format="html" label="Results (HTML) on ${on_string}"/>
+        <data name="all_events_table" format="tabular" label="All events table on ${on_string}"/>
+        <data name="replicates_counts" format="tabular" label="Replicates counts on ${on_string}"/>
     </outputs>
     <tests>
         <test>
-            <param name="expt" value="sacCer3_1.scidx" ftype="scidx" dbkey="sacCer3"/>
-            <param name="binding_model_smoothing" value="no"/>
-            <param name="gaussmodelsmoothing" value="yes"/>
-            <param name="gausssmoothparam" value="3"/>
-            <param name="use_motif" value="yes"/>
-            <param name="reference_genome_source" value="history"/>
-            <param name="reference_genome" value="phiX.fasta" dbkey="phiX"/>
-            <param name="memenmotifs" value="3"/>
-            <param name="mememinw" value="6"/>
-            <param name="mememaxw" value="16"/>
-            <param name="output_process_log" value="yes"/>
-            <output name="output_process" file="output_process1.txt" ftype="txt" lines_diff="12"/>
-            <output name="output_html" file="output_html1.html" ftype="html" lines_diff="12"/>
+            <param name="expt" value="expt_hg19.scidx" ftype="bam" dbkey="sacCer3" />
+            <param name="verbose" value="no" />
+            <param name="advanced_options_cond" value="hide" />
+            <output name="output_html" file="hg19_output_html1.html" ftype="html" lines_diff="12"/>
+            <output name="all_events_table" file="hg19_all_events_table1.tabular" ftype="tabular"/>
+            <output name="replicates_counts" file="hg19_replicates_counts1.tabular" ftype="tabular"/>
+        </test>
+        <test>
+            <param name="expt" value="expt_hg19.scidx" ftype="bam" dbkey="sacCer3" />
+            <param name="ctrl" value="cntrl_hg19.scidx" ftype="bam" dbkey="sacCer3" />
+            <param name="verbose" value="no" />
+            <param name="advanced_options_cond" value="display" />
+            <output name="output_html" file="hg19_output_html2.html" ftype="html" lines_diff="12"/>
+            <output name="all_events_table" file="hg19_all_events_table2.tabular" ftype="tabular"/>
+            <output name="replicates_counts" file="hg19_replicates_counts2.tabular" ftype="tabular"/>
         </test>
     </tests>
     <help>
@@ -446,7 +454,7 @@
 
 * **Loading data:**
 
- - **Optional file containing file containing reads from a control experiment** - must be same format as input experiment
+ - **Optional file containing reads from a control experiment** - must be same format as input experiment
  - **Fixed per-base limit** - Fixed per-base limit (default: estimated from background model).
  - **Poisson threshold for filtering per base** - Look at neighboring positions to decide what the per-base limit should be.
  - **Use non-unique reads** - Use non-unique reads.