Mercurial > repos > jjohnson > cummerbund

comparison cuffdiff_wrapper.xml @ 2:fdf01b3c1841

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
author Jim Johnson <jj@umn.edu>
date Fri, 08 Nov 2013 14:54:01 -0600
parents da7241f92ecf
children f109453ecfa2
comparison
equal deleted inserted replaced
1:ebb9a992508d 2:fdf01b3c1841
1 <tool id="cuffdiff_cummerbund" name="Cuffdiff" version="0.0.6"> 1 <tool id="cuffdiff_cummerbund" name="Cuffdiff for cummeRbund" version="0.0.7">
2 <!-- Wrapper supports Cuffdiff versions v1.3.0-v2.0 --> 2 <!-- Wrapper supports Cuffdiff versions 2.1.0-2.1.1 -->
3 <description>find significant changes in transcript expression, splicing, and promoter use</description> 3 <description>find significant changes in transcript expression, splicing, and promoter use</description>
4 <requirements> 4 <requirements>
5 <requirement type="package">cufflinks</requirement> 5 <requirement type="package" version="2.1.1">cufflinks</requirement>
6 </requirements> 6 </requirements>
7 <command interpreter="python"> 7 <version_command>cuffdiff 2>&amp;1 | head -n 1</version_command>
8 #set sel_outputs = $output_sel.__str__.split(',') 8 <command>
9 cuffdiff_wrapper.py 9 cuffdiff
10 --no-update-check
10 --FDR=$fdr 11 --FDR=$fdr
11 --num-threads="4" 12 --num-threads="4"
12 --min-alignment-count=$min_alignment_count 13 --min-alignment-count=$min_alignment_count
13 14 --library-norm-method=$library_norm_method
14 #if 'cuffdata' in $sel_outputs or not $output_sel: 15 --dispersion-method=$dispersion_method
15 --cuffdatadir=$cuffdata.extra_files_path 16
16 #end if
17 #if 'cummeRbund_db' in $sel_outputs:
18 --cummeRbund_db=$cummeRbund_db
19 #end if
20
21 #if 'isoforms_fpkm_tracking' in $sel_outputs:
22 --isoforms_fpkm_tracking_output=$isoforms_fpkm_tracking
23 #end if
24 #if 'genes_fpkm_tracking' in $sel_outputs:
25 --genes_fpkm_tracking_output=$genes_fpkm_tracking
26 #end if
27 #if 'cds_fpkm_tracking' in $sel_outputs:
28 --cds_fpkm_tracking_output=$cds_fpkm_tracking
29 #end if
30 #if 'tss_groups_fpkm_tracking' in $sel_outputs:
31 --tss_groups_fpkm_tracking_output=$tss_groups_fpkm_tracking
32 #end if
33 #if 'isoforms_exp_diff' in $sel_outputs:
34 --isoforms_exp_output=$isoforms_exp_diff
35 #end if
36 #if 'genes_exp_diff' in $sel_outputs:
37 --genes_exp_output=$genes_exp_diff
38 #end if
39 #if 'tss_groups_exp_diff' in $sel_outputs:
40 --tss_groups_exp_output=$tss_groups_exp_diff
41 #end if
42 #if 'cds_exp_fpkm_tracking' in $sel_outputs:
43 --cds_exp_fpkm_tracking_output=$cds_exp_fpkm_tracking
44 #end if
45 #if 'splicing_diff' in $sel_outputs:
46 --splicing_diff_output=$splicing_diff
47 #end if
48 #if 'cds_diff' in $sel_outputs:
49 --cds_diff_output=$cds_diff
50 #end if
51 #if 'promoters_diff' in $sel_outputs:
52 --promoters_diff_output=$promoters_diff
53 #end if
54 #if 'cds_read_group_tracking' in $sel_outputs:
55 --cds_read_group_tracking=$cds_read_group_tracking
56 #end if
57 #if 'tss_groups_read_group_tracking' in $sel_outputs:
58 --tss_groups_read_group_tracking=$tss_groups_read_group_tracking
59 #end if
60 #if 'genes_read_group_tracking' in $sel_outputs:
61 --genes_read_group_tracking=$genes_read_group_tracking
62 #end if
63 #if 'isoforms_read_group_tracking' in $sel_outputs:
64 --isoforms_read_group_tracking=$isoforms_read_group_tracking
65 #end if
66
67 ## Set advanced data parameters? 17 ## Set advanced data parameters?
68 #if $additional.sAdditional == "Yes": 18 #if $additional.sAdditional == "Yes":
69 -m $additional.frag_mean_len 19 #if $additional.frag_mean_len:
70 -s $additional.frag_len_std_dev 20 -m $additional.frag_mean_len
71 #end if 21 #end if
72 22 #if $additional.frag_len_std_dev:
73 ## Normalization? 23 -s $additional.frag_len_std_dev
74 #if str($do_normalization) == "Yes": 24 #end if
75 -N 25 #if $additional.max_bundle_frags:
26 --max-bundle-frags="$additional.max_bundle_frags"
27 #end if
76 #end if 28 #end if
77 29
78 ## Multi-read correct? 30 ## Multi-read correct?
79 #if str($multiread_correct) == "Yes": 31 #if str($multiread_correct) == "Yes":
80 -u 32 -u
81 #end if 33 #end if
82 34
83 ## Bias correction? 35 ## Bias correction?
84 #if $bias_correction.do_bias_correction == "Yes": 36 #if $bias_correction.do_bias_correction == "Yes":
85 -b 37 -b
86 #if $bias_correction.seq_source.index_source == "history": 38 #if $bias_correction.seq_source.index_source == "history":
87 --ref_file=$bias_correction.seq_source.ref_file 39 ## Custom genome from history.
40 $bias_correction.seq_source.ref_file
88 #else: 41 #else:
89 --ref_file="None" 42 ## Built-in genome.
90 #end if 43 ${__get_data_table_entry__('sam_fa_indexes', 'value', $gtf_input.dbkey, 'path')}
91 --dbkey=${gtf_input.metadata.dbkey} 44 #end if
92 --index_dir=${GALAXY_DATA_INDEX_DIR} 45 #end if
93 #end if 46
94 47 #set labels = '\'' + '\',\''.join( [ str( $condition.name ) for $condition in $conditions ] ) + '\''
48 --labels $labels
49
95 ## Inputs. 50 ## Inputs.
96 --inputA=$gtf_input 51 $gtf_input
97 #if $group_analysis.do_groups == "No": 52 #for $condition in $conditions:
98 --input1=$aligned_reads1 53 #set samples = ','.join( [ str( $sample.sample ) for $sample in $condition.samples ] )
99 --input2=$aligned_reads2 54 $samples
100 #else: 55 #end for
101 ## Replicates. 56
102 --labels 57 ## If build cummerbund db
103 #for $group in $group_analysis.groups 58 #if $build_cummerbund_db:
104 ${group.group} 59 &amp;&amp; echo 'library(cummeRbund)' > cuffData.r
105 #end for 60 #if $bias_correction.do_bias_correction == "Yes":
106 --files 61 #if $bias_correction.seq_source.index_source == "history":
107 #for $group in $group_analysis.groups 62 ## Custom genome from history.
108 #for $file in $group.files: 63 &amp;&amp; echo 'cuff&lt;-readCufflinks( dbFile = "cuffdata.db", gtfFile = "$gtf_input", genome = "$bias_correction.seq_source.ref_file", rebuild = T)' >> cuffData.r
109 ${file.file} 64 #else:
110 #end for 65 ## Built-in genome.
111 , 66 ${__get_data_table_entry__('sam_fa_indexes', 'value', $gtf_input.dbkey, 'path')}
112 #end for 67 &amp;&amp; echo 'cuff&lt;-readCufflinks( dbFile = "cuffdata.db", gtfFile = "$gtf_input", genome = "${__get_data_table_entry__('sam_fa_indexes', 'value', $gtf_input.dbkey, 'path')}", rebuild = T)' >> cuffData.r
113 #end if 68 #end if
114 69 #else
70 &amp;&amp; echo 'cuff&lt;-readCufflinks( dbFile = "cuffdata.db", rebuild = T)' >> cuffData.r
71 #end if
72 &amp;&amp; Rscript --vanilla cuffData.r
73 &amp;&amp; cp cuffdata.db $cummerbund_db
74 #end if
115 </command> 75 </command>
116 <inputs> 76 <inputs>
117 <param format="gtf,gff3" name="gtf_input" type="data" label="Transcripts" help="A transcript GFF3 or GTF file produced by cufflinks, cuffcompare, or other source."/> 77 <param format="gtf,gff3" name="gtf_input" type="data" label="Transcripts" help="A transcript GFF3 or GTF file produced by cufflinks, cuffcompare, or other source."/>
118 <conditional name="group_analysis"> 78
119 <param name="do_groups" type="select" label="Perform replicate analysis" help="Perform cuffdiff with replicates in each group."> 79 <repeat name="conditions" title="Condition" min="2">
120 <option value="No">No</option> 80 <param name="name" title="Condition name" type="text" label="Name"/>
121 <option value="Yes">Yes</option> 81 <repeat name="samples" title="Replicate" min="1">
122 </param> 82 <param name="sample" label="Add replicate" type="data" format="sam,bam"/>
123 <when value="Yes"> 83 </repeat>
124 <repeat name="groups" title="Group"> 84 </repeat>
125 <param name="group" title="Group name" type="text" label="Group name (no spaces or commas)"/> 85 <param name="time_series" type="boolean" checked="false" truevalue="--time-series" falsevalue="" optional="true" label="treat samples as a time-series">
126 <repeat name="files" title="Replicate"> 86 <help>
127 <param name="file" label="Add file" type="data" format="sam,bam"/> 87 Instructs Cuffdiff to analyze the provided samples as a time series, rather than testing for differences between all pairs of samples.
128 </repeat> 88 Samples should be provided in increasing time order at the command line (e.g first time point SAM, second timepoint SAM, etc.)
129 </repeat> 89 </help>
130 </when> 90 </param>
131 <when value="No"> 91
132 <param format="sam,bam" name="aligned_reads1" type="data" label="SAM or BAM file of aligned RNA-Seq reads" help=""/> 92 <param name="library_norm_method" type="select" label="Library normalization method">
133 <param format="sam,bam" name="aligned_reads2" type="data" label="SAM or BAM file of aligned RNA-Seq reads" help=""/> 93 <option value="geometric" selected="True">geometric</option>
134 </when> 94 <option value="classic-fpkm">classic-fpkm</option>
135 </conditional> 95 <option value="quartile">quartile</option>
96 </param>
97
98 <param name="dispersion_method" type="select" label="Dispersion estimation method" help="If using only one sample per condition, you must use 'blind.'">
99 <option value="pooled" selected="True">pooled</option>
100 <option value="per-condition">per-condition</option>
101 <option value="blind">blind</option>
102 </param>
136 103
137 <param name="fdr" type="float" value="0.05" label="False Discovery Rate" help="The allowed false discovery rate."/> 104 <param name="fdr" type="float" value="0.05" label="False Discovery Rate" help="The allowed false discovery rate."/>
138 105
139 <param name="min_alignment_count" type="integer" value="10" label="Min Alignment Count" help="The minimum number of alignments in a locus for needed to conduct significance testing on changes in that locus observed between samples."/> 106 <param name="min_alignment_count" type="integer" value="10" label="Min Alignment Count" help="The minimum number of alignments in a locus for needed to conduct significance testing on changes in that locus observed between samples."/>
140
141 <param name="do_normalization" type="select" label="Perform quartile normalization" help="Removes top 25% of genes from FPKM denominator to improve accuracy of differential expression calls for low abundance transcripts.">
142 <option value="No">No</option>
143 <option value="Yes">Yes</option>
144 </param>
145 107
146 <param name="multiread_correct" type="select" label="Use multi-read correct" help="Tells Cufflinks to do an initial estimation procedure to more accurately weight reads mapping to multiple locations in the genome."> 108 <param name="multiread_correct" type="select" label="Use multi-read correct" help="Tells Cufflinks to do an initial estimation procedure to more accurately weight reads mapping to multiple locations in the genome.">
147 <option value="No" selected="true">No</option> 109 <option value="No" selected="true">No</option>
148 <option value="Yes">Yes</option> 110 <option value="Yes">Yes</option>
149 </param> 111 </param>
166 </conditional> 128 </conditional>
167 </when> 129 </when>
168 <when value="No"></when> 130 <when value="No"></when>
169 </conditional> 131 </conditional>
170 132
133 <param name="include_read_group_files" type="select" label="Include Read Group Datasets" help="Read group datasets provide information on replicates.">
134 <option value="No" selected="true">No</option>
135 <option value="Yes">Yes</option>
136 </param>
137 <param name="build_cummerbund_db" type="boolean" truevalue="yes" falsevalue="no" checked="false" label="Build cummeRbund database"/>
171 <conditional name="additional"> 138 <conditional name="additional">
172 <param name="sAdditional" type="select" label="Set Additional Parameters? (not recommended)"> 139 <param name="sAdditional" type="select" label="Set Additional Parameters? (not recommended for paired-end reads)">
173 <option value="No">No</option> 140 <option value="No">No</option>
174 <option value="Yes">Yes</option> 141 <option value="Yes">Yes</option>
175 </param> 142 </param>
176 <when value="No"></when> 143 <when value="No"></when>
177 <when value="Yes"> 144 <when value="Yes">
178 <param name="frag_mean_len" type="integer" value="200" label="Average Fragment Length"/> 145 <param name="frag_mean_len" type="integer" value="" optional="true" label="Average Fragment Length Default: 200">
179 <param name="frag_len_std_dev" type="integer" value="80" label="Fragment Length Standard Deviation"/> 146 <help>
147 Note: Cufflinks now learns the fragment length mean for each SAM file,
148 so using this option is no longer recommended with paired-end reads.
149 </help>
150 </param>
151 <param name="frag_len_std_dev" type="integer" value="" optional="true" label="Fragment Length Standard Deviation Default: 80">
152 <help>
153 Note: Cufflinks now learns the fragment length mean for each SAM file,
154 so using this option is no longer recommended with paired-end reads.
155 </help>
156 </param>
157 <param name="max_bundle_frags" type="integer" value="" optional="true" label="--max-bundle-frags">
158 <help>
159 Sets the maximum number of fragments a locus may have before being skipped. Skipped loci are listed in skipped.gtf. Default: 1000000
160 </help>
161 <validator type="in_range" message="Value greater than 0" min="1"/>
162 </param>
180 </when> 163 </when>
181 </conditional> 164 </conditional>
182
183 <param name="output_sel" type="select" multiple="true" display="checkboxes" force_select="true" label="Select outputs for history datasets">
184 <option value="cuffdata">cuffdata - html page with links to cuffdiff outputs</option>
185 <option value="cummeRbund_db">cummeRbund database</option>
186 <option value="run_info">run.info</option>
187 <option value="read_groups_info">read_groups.info</option>
188 <option value="splicing_diff">splicing.diff</option>
189 <option value="promoters_diff">promoters.diff</option>
190 <option value="genes_exp_diff">genes_exp.diff</option>
191 <option value="genes_fpkm_tracking">genes.fpkm_tracking</option>
192 <option value="genes_count_tracking">genes.count_tracking</option>
193 <option value="genes_read_group_tracking">genes.read_group_tracking</option>
194 <option value="isoforms_exp_diff">isoforms.exp_diff</option>
195 <option value="isoforms_fpkm_tracking">isoforms.fpkm_tracking</option>
196 <option value="isoforms_count_tracking">isoforms.count_tracking</option>
197 <option value="isoforms_read_group_tracking">isoforms.read_group_tracking</option>
198 <option value="cds_diff">cds.diff</option>
199 <option value="cds_exp_diff">cds_exp.diff</option>
200 <option value="cds_fpkm_tracking">cds.fpkm_tracking</option>
201 <option value="cds_count_tracking">cds.count_tracking</option>
202 <option value="cds_read_group_tracking">cds.read_group_tracking</option>
203 <option value="tss_groups_exp_diff">tss_groups_exp.diff</option>
204 <option value="tss_groups_fpkm_tracking">tss_groups.fpkm_tracking</option>
205 <option value="tss_groups_count_tracking">tss_groups.count_tracking</option>
206 <option value="tss_groups_read_group_tracking">tss_groups.read_group_tracking</option>
207 </param>
208
209 </inputs> 165 </inputs>
210 166
167 <stdio>
168 <regex match="Error" source="both" level="fatal" description="Error"/>
169 <regex match=".*" source="both" level="log" description="tool progress"/>
170 </stdio>
171
211 <outputs> 172 <outputs>
212 <data format="text" name="run_info" label="${tool.name} on ${on_string}: run.info"> 173 <!-- Optional read group datasets. -->
213 <filter>output_sel and 'run_info' in output_sel</filter> 174 <data format="cuffdatadb" name="cummerbund_db" label="${tool.name} on ${on_string}: cummeRbund sqlite Database" >
214 </data> 175 <filter>build_cummerbund_db</filter>
215 <data format="tabular" name="read_groups_info" label="${tool.name} on ${on_string}: read_groups.info"> 176 </data>
216 <filter>output_sel and 'read_groups_info' in output_sel</filter> 177 <data format="tabular" name="isoforms_read_group" label="${tool.name} on ${on_string}: isoforms read group tracking" from_work_dir="isoforms.read_group_tracking" >
217 </data> 178 <filter>(params['include_read_group_files'] == 'Yes'</filter>
218 <data format="tabular" name="splicing_diff" label="${tool.name} on ${on_string}: splicing differential expression testing"> 179 </data>
219 <filter>output_sel and 'splicing_diff' in output_sel</filter> 180 <data format="tabular" name="genes_read_group" label="${tool.name} on ${on_string}: genes read group tracking" from_work_dir="genes.read_group_tracking" >
220 </data> 181 <filter>(params['include_read_group_files'] == 'Yes'</filter>
221 <data format="tabular" name="promoters_diff" label="${tool.name} on ${on_string}: promoters differential expression testing"> 182 </data>
222 <filter>output_sel and 'promoters_diff' in output_sel</filter> 183 <data format="tabular" name="cds_read_group" label="${tool.name} on ${on_string}: CDs read group tracking" from_work_dir="cds.read_group_tracking" >
223 </data> 184 <filter>(params['include_read_group_files'] == 'Yes'</filter>
224 <data format="tabular" name="cds_diff" label="${tool.name} on ${on_string}: CDS overloading diffential expression testing"> 185 </data>
225 <filter>output_sel and 'cds_diff' in output_sel</filter> 186 <data format="tabular" name="tss_groups_read_group" label="${tool.name} on ${on_string}: TSS groups read group tracking" from_work_dir="tss_groups.read_group_tracking" >
226 </data> 187 <filter>(params['include_read_group_files'] == 'Yes'</filter>
227 <data format="tabular" name="cds_exp_diff" label="${tool.name} on ${on_string}: CDS differential expression testing"> 188 </data>
228 <filter>output_sel and 'cds_exp_diff' in output_sel</filter> 189
229 </data> 190 <!-- Standard datasets. -->
230 <data format="tabular" name="cds_fpkm_tracking" label="${tool.name} on ${on_string}: CDS FPKM tracking"> 191 <data format="tabular" name="splicing_diff" label="${tool.name} on ${on_string}: splicing differential expression testing" from_work_dir="splicing.diff" />
231 <filter>output_sel and 'cds_fpkm_tracking' in output_sel</filter> 192 <data format="tabular" name="promoters_diff" label="${tool.name} on ${on_string}: promoters differential expression testing" from_work_dir="promoters.diff" />
232 </data> 193 <data format="tabular" name="cds_diff" label="${tool.name} on ${on_string}: CDS overloading diffential expression testing" from_work_dir="cds.diff" />
233 <data format="tabular" name="cds_count_tracking" label="${tool.name} on ${on_string}: CDS counts"> 194 <data format="tabular" name="cds_exp_fpkm_tracking" label="${tool.name} on ${on_string}: CDS FPKM differential expression testing" from_work_dir="cds_exp.diff" />
234 <filter>output_sel and 'cds_count_tracking' in output_sel</filter> 195 <data format="tabular" name="cds_fpkm_tracking" label="${tool.name} on ${on_string}: CDS FPKM tracking" from_work_dir="cds.fpkm_tracking" />
235 </data> 196 <data format="tabular" name="tss_groups_exp" label="${tool.name} on ${on_string}: TSS groups differential expression testing" from_work_dir="tss_group_exp.diff" />
236 <data format="tabular" name="cds_read_group_tracking" label="${tool.name} on ${on_string}: CDS Read Group tracking"> 197 <data format="tabular" name="tss_groups_fpkm_tracking" label="${tool.name} on ${on_string}: TSS groups FPKM tracking" from_work_dir="tss_groups.fpkm_tracking" />
237 <filter>output_sel and 'cds_read_group_tracking' in output_sel</filter> 198 <data format="tabular" name="genes_exp" label="${tool.name} on ${on_string}: gene differential expression testing" from_work_dir="gene_exp.diff" />
238 </data> 199 <data format="tabular" name="genes_fpkm_tracking" label="${tool.name} on ${on_string}: gene FPKM tracking" from_work_dir="genes.fpkm_tracking" />
239 <data format="tabular" name="tss_groups_exp_diff" label="${tool.name} on ${on_string}: TSS groups differential expression testing"> 200 <data format="tabular" name="isoforms_exp" label="${tool.name} on ${on_string}: transcript differential expression testing" from_work_dir="isoform_exp.diff" />
240 <filter>output_sel and 'tss_groups_exp_diff' in output_sel</filter> 201 <data format="tabular" name="isoforms_fpkm_tracking" label="${tool.name} on ${on_string}: transcript FPKM tracking" from_work_dir="isoforms.fpkm_tracking" />
241 </data>
242 <data format="tabular" name="tss_groups_fpkm_tracking" label="${tool.name} on ${on_string}: TSS groups FPKM tracking">
243 <filter>output_sel and 'tss_groups_fpkm_tracking' in output_sel</filter>
244 </data>
245 <data format="tabular" name="tss_groups_count_tracking" label="${tool.name} on ${on_string}: TSS groups counts">
246 <filter>output_sel and 'tss_groups_count_tracking' in output_sel</filter>
247 </data>
248 <data format="tabular" name="tss_groups_read_group_tracking" label="${tool.name} on ${on_string}: TSS groups Read Group tracking">
249 <filter>output_sel and 'tss_groups_read_group_tracking' in output_sel</filter>
250 </data>
251 <data format="tabular" name="isoforms_exp_diff" label="${tool.name} on ${on_string}: transcript differential expression testing">
252 <filter>output_sel and 'isoforms_exp_diff' in output_sel</filter>
253 </data>
254 <data format="tabular" name="isoforms_fpkm_tracking" label="${tool.name} on ${on_string}: transcript FPKM tracking">
255 <filter>output_sel and 'isoforms_fpkm_tracking' in output_sel</filter>
256 </data>
257 <data format="tabular" name="isoforms_count_tracking" label="${tool.name} on ${on_string}: transcript counts">
258 <filter>output_sel and 'isoforms_count_tracking' in output_sel</filter>
259 </data>
260 <data format="tabular" name="isoforms_read_group_tracking" label="${tool.name} on ${on_string}: transcript Read Group tracking">
261 <filter>output_sel and 'isoforms_read_group_tracking' in output_sel</filter>
262 </data>
263 <data format="tabular" name="genes_exp_diff" label="${tool.name} on ${on_string}: gene differential expression testing">
264 <filter>output_sel and 'genes_exp_diff' in output_sel</filter>
265 </data>
266 <data format="tabular" name="genes_fpkm_tracking" label="${tool.name} on ${on_string}: gene FPKM tracking">
267 <filter>output_sel and 'genes_fpkm_tracking' in output_sel</filter>
268 </data>
269 <data format="tabular" name="genes_count_tracking" label="${tool.name} on ${on_string}: gene counts">
270 <filter>output_sel and 'genes_count_tracking' in output_sel</filter>
271 </data>
272 <data format="tabular" name="genes_read_group_tracking" label="${tool.name} on ${on_string}: gene Read Group tracking">
273 <filter>output_sel and 'genes_read_group_tracking' in output_sel</filter>
274 </data>
275 <data format="cuffdata" name="cuffdata" label="${tool.name} on ${on_string}: cuffdata" >
276 <filter>not output_sel or output_sel and 'cuffdata' in output_sel</filter>
277 </data>
278 <data format="cuffdatadb" name="cummeRbund_db" label="${tool.name} on ${on_string}: cummeRbund sqlite Database" >
279 <filter>output_sel and 'cummeRbund_db' in output_sel</filter>
280 </data>
281 </outputs> 202 </outputs>
282 <stdio>
283 <exit_code range="1:" level="fatal" description="Cufflinks Err" />
284 </stdio>
285
286 203
287 <tests> 204 <tests>
288 <test> 205 <test>
289 <!-- 206 <!--
290 cuffdiff cuffcompare_out5.gtf cuffdiff_in1.sam cuffdiff_in2.sam 207 cuffdiff cuffcompare_out5.gtf cuffdiff_in1.sam cuffdiff_in2.sam
291 --> 208 -->
209 <!--
210 NOTE: as of version 0.0.6 of the wrapper, tests cannot be run because multiple inputs to a repeat
211 element are not supported.
292 <param name="gtf_input" value="cuffcompare_out5.gtf" ftype="gtf" /> 212 <param name="gtf_input" value="cuffcompare_out5.gtf" ftype="gtf" />
293 <param name="do_groups" value="No" /> 213 <param name="do_groups" value="No" />
294 <param name="aligned_reads1" value="cuffdiff_in1.sam" ftype="sam" /> 214 <param name="aligned_reads1" value="cuffdiff_in1.sam" ftype="sam" />
295 <param name="aligned_reads2" value="cuffdiff_in2.sam" ftype="sam" /> 215 <param name="aligned_reads2" value="cuffdiff_in2.sam" ftype="sam" />
296 <!-- Defaults. -->
297 <param name="fdr" value="0.05" /> 216 <param name="fdr" value="0.05" />
298 <param name="min_alignment_count" value="0" /> 217 <param name="min_alignment_count" value="0" />
299 <param name="do_bias_correction" value="No" /> 218 <param name="do_bias_correction" value="No" />
300 <param name="do_normalization" value="No" /> 219 <param name="do_normalization" value="No" />
301 <param name="multiread_correct" value="No"/> 220 <param name="multiread_correct" value="No"/>
302 <param name="sAdditional" value="No"/> 221 <param name="sAdditional" value="No"/>
303 <!--
304 Line diffs are needed because cuffdiff does not produce deterministic output.
305 TODO: can we find datasets that lead to deterministic behavior?
306 -->
307 <output name="splicing_diff" file="cuffdiff_out9.txt"/> 222 <output name="splicing_diff" file="cuffdiff_out9.txt"/>
308 <output name="promoters_diff" file="cuffdiff_out10.txt"/> 223 <output name="promoters_diff" file="cuffdiff_out10.txt"/>
309 <output name="cds_diff" file="cuffdiff_out11.txt"/> 224 <output name="cds_diff" file="cuffdiff_out11.txt"/>
310 <output name="cds_exp_fpkm_tracking" file="cuffdiff_out4.txt"/> 225 <output name="cds_exp_fpkm_tracking" file="cuffdiff_out4.txt"/>
311 <output name="cds_fpkm_tracking" file="cuffdiff_out8.txt"/> 226 <output name="cds_fpkm_tracking" file="cuffdiff_out8.txt"/>
313 <output name="tss_groups_fpkm_tracking" file="cuffdiff_out7.txt"/> 228 <output name="tss_groups_fpkm_tracking" file="cuffdiff_out7.txt"/>
314 <output name="genes_exp" file="cuffdiff_out2.txt" lines_diff="200"/> 229 <output name="genes_exp" file="cuffdiff_out2.txt" lines_diff="200"/>
315 <output name="genes_fpkm_tracking" file="cuffdiff_out6.txt" lines_diff="200"/> 230 <output name="genes_fpkm_tracking" file="cuffdiff_out6.txt" lines_diff="200"/>
316 <output name="isoforms_exp" file="cuffdiff_out1.txt" lines_diff="200"/> 231 <output name="isoforms_exp" file="cuffdiff_out1.txt" lines_diff="200"/>
317 <output name="isoforms_fpkm_tracking" file="cuffdiff_out5.txt" lines_diff="200"/> 232 <output name="isoforms_fpkm_tracking" file="cuffdiff_out5.txt" lines_diff="200"/>
233 -->
318 </test> 234 </test>
319 </tests> 235 </tests>
320 236
321 <help> 237 <help>
322 **Cuffdiff Overview** 238 **Cuffdiff Overview**