Mercurial > repos > iuc > deseq2
comparison deseq2.xml @ 27:d027d1f4984e draft
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit c8e7020d39d581d7374a13cf94c64998a6481e05"
author | iuc |
---|---|
date | Mon, 28 Jun 2021 05:14:25 +0000 |
parents | 6a3a025714d3 |
children | 7ff33c2d525b |
comparison
equal
deleted
inserted
replaced
26:6a3a025714d3 | 27:d027d1f4984e |
---|---|
1 <tool id="deseq2" name="DESeq2" version="2.11.40.6+galaxy1"> | 1 <tool id="deseq2" name="DESeq2" version="@TOOL_VERSION@+galaxy@SUFFIX_VERSION@"> |
2 <description>Determines differentially expressed features from count tables</description> | 2 <description>Determines differentially expressed features from count tables</description> |
3 <macros> | 3 <macros> |
4 <import>deseq2_macros.xml</import> | 4 <import>deseq2_macros.xml</import> |
5 </macros> | 5 </macros> |
6 <requirements> | 6 <expand macro='requirements'/> |
7 <requirement type="package" version="1.22.1">bioconductor-deseq2</requirement> | 7 <expand macro='edam_ontology' /> |
8 <!-- Optional dependency of tximport, needed to import kallisto results https://github.com/galaxyproject/usegalaxy-playbook/issues/161 --> | 8 <expand macro='xrefs'/> |
9 <requirement type="package" version="2.26.2">bioconductor-rhdf5</requirement> | |
10 <requirement type="package" version="1.10.0">bioconductor-tximport</requirement> | |
11 <requirement type="package" version="1.34.1">bioconductor-genomicfeatures</requirement> | |
12 <requirement type="package" version="1.20.2">r-getopt</requirement> | |
13 <requirement type="package" version="0.8.0">r-ggrepel</requirement> | |
14 <requirement type="package" version="3.0.1">r-gplots</requirement> | |
15 <requirement type="package" version="1.0.10">r-pheatmap</requirement> | |
16 <requirement type="package" version="0.2.20">r-rjson</requirement> | |
17 </requirements> | |
18 <stdio> | 9 <stdio> |
19 <regex match="Execution halted" | 10 <regex match="Execution halted" |
20 source="both" | 11 source="both" |
21 level="fatal" | 12 level="fatal" |
22 description="Execution halted." /> | 13 description="Execution halted." /> |
47 #import json | 38 #import json |
48 #import os | 39 #import os |
49 Rscript '${__tool_directory__}/deseq2.R' | 40 Rscript '${__tool_directory__}/deseq2.R' |
50 --cores \${GALAXY_SLOTS:-1} | 41 --cores \${GALAXY_SLOTS:-1} |
51 -o '$deseq_out' | 42 -o '$deseq_out' |
52 #if $pdf: | 43 #if 'pdf' in $output_options.output_selector: |
53 -p '$plots' | 44 -p '$plots' |
54 #end if | 45 #end if |
55 #if $normCounts: | 46 -A $output_options.alpha_ma |
47 #if 'normCounts' in $output_options.output_selector: | |
56 -n '$counts_out' | 48 -n '$counts_out' |
57 #end if | 49 #end if |
58 #if $normRLog: | 50 #if 'normRLog' in $output_options.output_selector: |
59 -r '$rlog_out' | 51 -r '$rlog_out' |
60 #end if | 52 #end if |
61 #if $normVST: | 53 #if 'normVST' in $output_options.output_selector: |
62 -v '$vst_out' | 54 -v '$vst_out' |
63 #end if | 55 #end if |
64 #set $filename_to_element_identifiers = {} | 56 #set $filename_to_element_identifiers = {} |
65 #set $temp_factor_names = list() | 57 #set $temp_factor_names = list() |
66 #for $factor in $select_data.rep_factorName: | 58 #for $factor in $select_data.rep_factorName: |
88 | 80 |
89 $header | 81 $header |
90 | 82 |
91 -f '#echo json.dumps(temp_factor_names)#' | 83 -f '#echo json.dumps(temp_factor_names)#' |
92 -l '#echo json.dumps(filename_to_element_identifiers)#' | 84 -l '#echo json.dumps(filename_to_element_identifiers)#' |
93 #if $esf: | 85 #if $advanced_options.esf: |
94 -e $esf | 86 -e $advanced_options.esf |
95 #end if | 87 #end if |
96 -t $fit_type | 88 -t $advanced_options.fit_type |
97 #if $batch_factors: | 89 #if $batch_factors: |
98 --batch_factors '$batch_factors' | 90 --batch_factors '$batch_factors' |
99 #end if | 91 #end if |
100 #if $outlier_replace_off: | 92 #if $advanced_options.outlier_replace_off: |
101 -a | 93 -a |
102 #end if | 94 #end if |
103 #if $outlier_filter_off: | 95 #if $advanced_options.outlier_filter_off: |
104 -b | 96 -b |
105 #end if | 97 #end if |
106 #if $auto_mean_filter_off: | 98 #if $advanced_options.auto_mean_filter_off: |
107 -c | 99 -c |
108 #end if | 100 #end if |
109 #if $many_contrasts: | 101 #if 'many_contrasts' in $output_options.output_selector |
110 -m | 102 -m |
111 #end if | 103 #end if |
112 #if $tximport.tximport_selector == 'tximport': | 104 #if $tximport.tximport_selector == 'tximport': |
113 -i | 105 -i |
114 -y $tximport.txtype | 106 -y $tximport.txtype |
166 </when> | 158 </when> |
167 </conditional> | 159 </conditional> |
168 </when> | 160 </when> |
169 <when value="count" /> | 161 <when value="count" /> |
170 </conditional> | 162 </conditional> |
171 <param name="pdf" type="boolean" truevalue="1" falsevalue="0" checked="true" | 163 <section name="advanced_options" title="Advanced options"> |
172 label="Visualising the analysis results" | 164 <param name="esf" type="select" label="Method for estimateSizeFactors" |
173 help="output an additional PDF files" /> | 165 help="Method for estimation: either 'ratio', 'poscounts', or 'iterate'. 'ratio' uses the standard median ratio method introduced in DESeq. |
174 <param name="normCounts" type="boolean" truevalue="1" falsevalue="0" checked="false" | 166 The size factor is the median ratio of the sample over a 'pseudosample': for each gene, the geometric mean of all samples. |
175 label="Output normalized counts table" /> | 167 'poscounts' and 'iterate' offer alternative estimators, which can be used even when all genes contain a sample with a zero (a problem |
176 <param name="normRLog" type="boolean" truevalue="1" falsevalue="0" checked="false" | 168 for the default method, as the geometric mean becomes zero, and the ratio undefined). The 'poscounts' estimator deals with a gene with |
177 label="Output rLog normalized table" /> | 169 some zeros, by calculating a modified geometric mean by taking the n-th root of the product of the non-zero counts. This evolved out of |
178 <param name="normVST" type="boolean" truevalue="1" falsevalue="0" checked="false" | 170 use cases with Paul McMurdie's phyloseq package for metagenomic samples. The 'iterate' estimator iterates between estimating the dispersion |
179 label="Output VST normalized table" /> | 171 with a design of ~1, and finding a size factor vector by numerically optimizing the likelihood of the ~1 model."> |
180 <param name="many_contrasts" type="boolean" truevalue="1" falsevalue="0" checked="false" | 172 <option value="" selected="true">No Selection (use default)</option> |
181 label="Output all levels vs all levels of primary factor (use when you have >2 levels for primary factor)" | 173 <option value="ratio">ratio</option> |
182 help=" DESeq2 performs independent filtering by default using the mean of normalized counts as a filter statistic" /> | 174 <option value="poscounts">poscounts</option> |
183 <param name="esf" type="select" label="(Optional) Method for estimateSizeFactors" | 175 <option value="iterate">iterate</option> |
184 help="Method for estimation: either 'ratio', 'poscounts', or 'iterate'. 'ratio' uses the standard median ratio method introduced in DESeq. | 176 </param> |
185 The size factor is the median ratio of the sample over a 'pseudosample': for each gene, the geometric mean of all samples. | 177 <param name="fit_type" type="select" label="Fit type"> |
186 'poscounts' and 'iterate' offer alternative estimators, which can be used even when all genes contain a sample with a zero (a problem | 178 <option value="1" selected="true">parametric</option> |
187 for the default method, as the geometric mean becomes zero, and the ratio undefined). The 'poscounts' estimator deals with a gene with | 179 <option value="2">local</option> |
188 some zeros, by calculating a modified geometric mean by taking the n-th root of the product of the non-zero counts. This evolved out of | 180 <option value="3">mean</option> |
189 use cases with Paul McMurdie's phyloseq package for metagenomic samples. The 'iterate' estimator iterates between estimating the dispersion | 181 </param> |
190 with a design of ~1, and finding a size factor vector by numerically optimizing the likelihood of the ~1 model."> | 182 <param name="outlier_replace_off" type="boolean" truevalue="1" falsevalue="0" checked="false" |
191 <option value="" selected="true">No Selection (use default)</option> | 183 label="Turn off outliers replacement (only affects with >6 replicates)" |
192 <option value="ratio">ratio</option> | 184 help="When there are more than 6 replicates for a given sample, the DESeq2 will automatically replace |
193 <option value="poscounts">poscounts</option> | 185 counts with large Cook’s distance with the trimmed mean over all samples, scaled up by the size factor |
194 <option value="iterate">iterate</option> | 186 or normalization factor for that sample" /> |
195 </param> | 187 <param name="outlier_filter_off" type="boolean" truevalue="1" falsevalue="0" checked="false" |
196 <param name="fit_type" type="select" label="Fit type"> | 188 label="Turn off outliers filtering (only affects with >2 replicates)" |
197 <option value="1" selected="true">parametric</option> | 189 help="When there are more than 2 replicates for a given sample, the DESeq2 will automatically |
198 <option value="2">local</option> | 190 filter genes which contain a Cook’s distance above a cutoff" /> |
199 <option value="3">mean</option> | 191 <param name="auto_mean_filter_off" type="boolean" truevalue="1" falsevalue="0" checked="false" |
200 </param> | 192 label="Turn off independent filtering" |
201 <param name="outlier_replace_off" type="boolean" truevalue="1" falsevalue="0" checked="false" | 193 help=" DESeq2 performs independent filtering by default using the mean of normalized counts as a filter statistic" /> |
202 label="Turn off outliers replacement (only affects with >6 replicates)" | 194 </section> |
203 help="When there are more than 6 replicates for a given sample, the DESeq2 will automatically replace | 195 <section name="output_options" title="Output options"> |
204 counts with large Cook’s distance with the trimmed mean over all samples, scaled up by the size factor | 196 <param name="output_selector" type="select" multiple="True" display="checkboxes" label="Output selector"> |
205 or normalization factor for that sample" /> | 197 <option value="pdf" selected="True">Generate plots for visualizing the analysis results</option> |
206 <param name="outlier_filter_off" type="boolean" truevalue="1" falsevalue="0" checked="false" | 198 <option value="normCounts">Output rLog normalized table</option> |
207 label="Turn off outliers filtering (only affects with >2 replicates)" | 199 <option value="normVST">Output VST normalized table</option> |
208 help="When there are more than 2 replicates for a given sample, the DESeq2 will automatically | 200 <option value="normRLog">Output rLog normalized table</option> |
209 filter genes which contain a Cook’s distance above a cutoff" /> | 201 <option value="many_contrasts">Output all levels vs all levels of primary factor (use when you have >2 levels for primary factor)</option> |
210 <param name="auto_mean_filter_off" type="boolean" truevalue="1" falsevalue="0" checked="false" | 202 </param> |
211 label="Turn off independent filtering" | 203 <param name="alpha_ma" type="float" min="0" max="0.5" value="0.1" label="Alpha value for MA-plot" help="Default value is 0.1. This option is only meaninful when plots are generated" /> |
212 help=" DESeq2 performs independent filtering by default using the mean of normalized counts as a filter statistic" /> | 204 </section> |
213 </inputs> | 205 </inputs> |
214 <outputs> | 206 <outputs> |
215 <data name="deseq_out" format="tabular" label="DESeq2 result file on ${on_string}"> | 207 <data name="deseq_out" format="tabular" label="DESeq2 result file on ${on_string}"> |
216 <filter>many_contrasts is False</filter> | 208 <filter>'many_contrasts' not in output_options['output_selector']</filter> |
217 <actions> | 209 <actions> |
218 <action name="column_names" type="metadata" default="GeneID,Base mean,log2(FC),StdErr,Wald-Stats,P-value,P-adj" /> | 210 <action name="column_names" type="metadata" default="GeneID,Base mean,log2(FC),StdErr,Wald-Stats,P-value,P-adj" /> |
219 </actions> | 211 </actions> |
220 </data> | 212 </data> |
221 <collection name="split_output" type="list" label="DESeq2 result files on ${on_string}"> | 213 <collection name="split_output" type="list" label="DESeq2 result files on ${on_string}"> |
222 <filter>many_contrasts is True</filter> | 214 <filter>output_options['output_selector'] and 'many_contrasts' in output_options['output_selector']</filter> |
223 <discover_datasets pattern="(?P<designation>.+_vs_.+)" format="tabular" directory="." visible="false"/> | 215 <discover_datasets pattern="(?P<designation>.+_vs_.+)" format="tabular" directory="." visible="false"/> |
224 </collection> | 216 </collection> |
225 <data name="plots" format="pdf" label="DESeq2 plots on ${on_string}"> | 217 <data name="plots" format="pdf" label="DESeq2 plots on ${on_string}"> |
226 <filter>pdf == True</filter> | 218 <filter>output_options['output_selector'] and 'pdf' in output_options['output_selector']</filter> |
227 </data> | 219 </data> |
228 <data name="counts_out" format="tabular" label="Normalized counts file on ${on_string}"> | 220 <data name="counts_out" format="tabular" label="Normalized counts file on ${on_string}"> |
229 <filter>normCounts == True</filter> | 221 <filter>output_options['output_selector'] and 'normCounts' in output_options['output_selector']</filter> |
230 </data> | 222 </data> |
231 <data name="rlog_out" format="tabular" label="rLog-Normalized counts file on ${on_string}"> | 223 <data name="rlog_out" format="tabular" label="rLog-Normalized counts file on ${on_string}"> |
232 <filter>normRLog == True</filter> | 224 <filter>output_options['output_selector'] and 'normRLog' in output_options['output_selector']</filter> |
233 </data> | 225 </data> |
234 <data name="vst_out" format="tabular" label="VST-Normalized counts file on ${on_string}"> | 226 <data name="vst_out" format="tabular" label="VST-Normalized counts file on ${on_string}"> |
235 <filter>normVST == True</filter> | 227 <filter>output_options['output_selector'] and 'normVST' in output_options['output_selector']</filter> |
236 </data> | 228 </data> |
237 </outputs> | 229 </outputs> |
238 <tests> | 230 <tests> |
239 <!--Ensure counts files with header works --> | 231 <!--Ensure counts files with header works --> |
240 <test expect_num_outputs="4"> | 232 <test expect_num_outputs="4"> |
247 <repeat name="rep_factorLevel"> | 239 <repeat name="rep_factorLevel"> |
248 <param name="factorLevel" value="Untreated"/> | 240 <param name="factorLevel" value="Untreated"/> |
249 <param name="countsFile" value="GSM461176_untreat_single.counts,GSM461177_untreat_paired.counts,GSM461178_untreat_paired.counts,GSM461182_untreat_single.counts"/> | 241 <param name="countsFile" value="GSM461176_untreat_single.counts,GSM461177_untreat_paired.counts,GSM461178_untreat_paired.counts,GSM461182_untreat_single.counts"/> |
250 </repeat> | 242 </repeat> |
251 </repeat> | 243 </repeat> |
252 <param name="pdf" value="False"/> | 244 <section name="output_options"> |
253 <param name="normCounts" value="True"/> | 245 <param name="output_selector" value="normCounts,normRLog,normVST"/> |
254 <param name="normRLog" value="True"/> | 246 </section> |
255 <param name="normVST" value="True"/> | |
256 <output name="counts_out"> | 247 <output name="counts_out"> |
257 <assert_contents> | 248 <assert_contents> |
258 <has_text_matching expression="GSM461176_untreat_single.counts\tGSM461177_untreat_paired.counts\tGSM461178_untreat_paired.counts\tGSM461182_untreat_single.counts\tGSM461179_treat_single.counts\tGSM461180_treat_paired.counts\tGSM461181_treat_paired.counts" /> | 249 <has_text_matching expression="GSM461176_untreat_single.counts\tGSM461177_untreat_paired.counts\tGSM461178_untreat_paired.counts\tGSM461182_untreat_single.counts\tGSM461179_treat_single.counts\tGSM461180_treat_paired.counts\tGSM461181_treat_paired.counts" /> |
259 <has_text_matching expression="FBgn0000003\t0\t0\t0\t0\t0\t0\t0" /> | 250 <has_text_matching expression="FBgn0000003\t0\t0\t0\t0\t0\t0\t0" /> |
260 </assert_contents> | 251 </assert_contents> |
289 <param name="factorLevel" value="Untreated"/> | 280 <param name="factorLevel" value="Untreated"/> |
290 <param name="countsFile" value="GSM461176_untreat_single.counts,GSM461177_untreat_paired.counts,GSM461178_untreat_paired.counts,GSM461182_untreat_single.counts"/> | 281 <param name="countsFile" value="GSM461176_untreat_single.counts,GSM461177_untreat_paired.counts,GSM461178_untreat_paired.counts,GSM461182_untreat_single.counts"/> |
291 </repeat> | 282 </repeat> |
292 </repeat> | 283 </repeat> |
293 <param name="batch_factors" value="batch_factors.tab"/> | 284 <param name="batch_factors" value="batch_factors.tab"/> |
294 <param name="pdf" value="False"/> | 285 <section name="output_options"> |
295 <param name="normCounts" value="True"/> | 286 <param name="output_selector" value="normCounts"/> |
287 </section> | |
296 <output name="deseq_out"> | 288 <output name="deseq_out"> |
297 <assert_contents> | 289 <assert_contents> |
298 <has_text_matching expression="FBgn0003360\t1933.*\t-2.9.*\t0.1.*\t-26.*\t1.*-152\t4.*-149" /> | 290 <has_text_matching expression="FBgn0003360\t1933.*\t-2.9.*\t0.1.*\t-26.*\t1.*-152\t4.*-149" /> |
299 </assert_contents> | 291 </assert_contents> |
300 </output> | 292 </output> |
311 <param name="factorLevel" value="Untreated"/> | 303 <param name="factorLevel" value="Untreated"/> |
312 <param name="countsFile" value="GSM461176_untreat_single.counts.noheader,GSM461177_untreat_paired.counts.noheader,GSM461178_untreat_paired.counts.noheader,GSM461182_untreat_single.counts.noheader"/> | 304 <param name="countsFile" value="GSM461176_untreat_single.counts.noheader,GSM461177_untreat_paired.counts.noheader,GSM461178_untreat_paired.counts.noheader,GSM461182_untreat_single.counts.noheader"/> |
313 </repeat> | 305 </repeat> |
314 </repeat> | 306 </repeat> |
315 <param name="header" value="False"/> | 307 <param name="header" value="False"/> |
316 <param name="pdf" value="False"/> | 308 <section name="output_options"> |
317 <param name="normCounts" value="True"/> | 309 <param name="output_selector" value="normCounts,normRLog,normVST"/> |
318 <param name="normRLog" value="True"/> | 310 </section> |
319 <param name="normVST" value="True"/> | |
320 <output name="counts_out"> | 311 <output name="counts_out"> |
321 <assert_contents> | 312 <assert_contents> |
322 <has_text_matching expression="GSM461176_untreat_single.counts.noheader\tGSM461177_untreat_paired.counts.noheader\tGSM461178_untreat_paired.counts.noheader\tGSM461182_untreat_single.counts.noheader\tGSM461179_treat_single.counts.noheader\tGSM461180_treat_paired.counts.noheader\tGSM461181_treat_paired.counts.noheader" /> | 313 <has_text_matching expression="GSM461176_untreat_single.counts.noheader\tGSM461177_untreat_paired.counts.noheader\tGSM461178_untreat_paired.counts.noheader\tGSM461182_untreat_single.counts.noheader\tGSM461179_treat_single.counts.noheader\tGSM461180_treat_paired.counts.noheader\tGSM461181_treat_paired.counts.noheader" /> |
323 <has_text_matching expression="FBgn0000003\t0\t0\t0\t0\t0\t0\t0" /> | 314 <has_text_matching expression="FBgn0000003\t0\t0\t0\t0\t0\t0\t0" /> |
324 </assert_contents> | 315 </assert_contents> |
352 <repeat name="rep_factorLevel"> | 343 <repeat name="rep_factorLevel"> |
353 <param name="factorLevel" value="Untreated"/> | 344 <param name="factorLevel" value="Untreated"/> |
354 <param name="countsFile" value="sailfish/sailfish_quant.sf4.tab,sailfish/sailfish_quant.sf5.tab,sailfish/sailfish_quant.sf6.tab"/> | 345 <param name="countsFile" value="sailfish/sailfish_quant.sf4.tab,sailfish/sailfish_quant.sf5.tab,sailfish/sailfish_quant.sf6.tab"/> |
355 </repeat> | 346 </repeat> |
356 </repeat> | 347 </repeat> |
357 <param name="pdf" value="False"/> | 348 <section name="output_options"> |
349 <param name="output_selector" value=""/> | |
350 </section> | |
358 <param name="tximport_selector" value="tximport"/> | 351 <param name="tximport_selector" value="tximport"/> |
359 <param name="txtype" value="sailfish"/> | 352 <param name="txtype" value="sailfish"/> |
360 <param name="mapping_format_selector" value="tabular"/> | 353 <param name="mapping_format_selector" value="tabular"/> |
361 <param name="tabular_file" value="tx2gene.tab"/> | 354 <param name="tabular_file" value="tx2gene.tab"/> |
362 <output name="deseq_out" > | 355 <output name="deseq_out" > |
376 <repeat name="rep_factorLevel"> | 369 <repeat name="rep_factorLevel"> |
377 <param name="factorLevel" value="Untreated"/> | 370 <param name="factorLevel" value="Untreated"/> |
378 <param name="countsFile" value="sailfish/sailfish_quant.sf4.tab,sailfish/sailfish_quant.sf5.tab,sailfish/sailfish_quant.sf6.tab"/> | 371 <param name="countsFile" value="sailfish/sailfish_quant.sf4.tab,sailfish/sailfish_quant.sf5.tab,sailfish/sailfish_quant.sf6.tab"/> |
379 </repeat> | 372 </repeat> |
380 </repeat> | 373 </repeat> |
381 <param name="pdf" value="False"/> | 374 <section name="output_options"> |
375 <param name="output_selector" value=""/> | |
376 </section> | |
382 <param name="tximport_selector" value="tximport"/> | 377 <param name="tximport_selector" value="tximport"/> |
383 <param name="txtype" value="sailfish"/> | 378 <param name="txtype" value="sailfish"/> |
384 <param name="mapping_format_selector" value="gtf"/> | 379 <param name="mapping_format_selector" value="gtf"/> |
385 <param name="gtf_file" value="GRCh38_latest_genomic.gff"/> | 380 <param name="gtf_file" value="GRCh38_latest_genomic.gff"/> |
386 <output name="deseq_out" > | 381 <output name="deseq_out" > |
400 <repeat name="rep_factorLevel"> | 395 <repeat name="rep_factorLevel"> |
401 <param name="factorLevel" value="Untreated"/> | 396 <param name="factorLevel" value="Untreated"/> |
402 <param name="countsFile" value="sailfish_ensembl/sailfish_quant.sf4.tab,sailfish_ensembl/sailfish_quant.sf5.tab,sailfish_ensembl/sailfish_quant.sf6.tab"/> | 397 <param name="countsFile" value="sailfish_ensembl/sailfish_quant.sf4.tab,sailfish_ensembl/sailfish_quant.sf5.tab,sailfish_ensembl/sailfish_quant.sf6.tab"/> |
403 </repeat> | 398 </repeat> |
404 </repeat> | 399 </repeat> |
405 <param name="pdf" value="False"/> | 400 <section name="output_options"> |
401 <param name="output_selector" value=""/> | |
402 </section> | |
406 <param name="tximport_selector" value="tximport"/> | 403 <param name="tximport_selector" value="tximport"/> |
407 <param name="txtype" value="sailfish"/> | 404 <param name="txtype" value="sailfish"/> |
408 <param name="mapping_format_selector" value="gtf"/> | 405 <param name="mapping_format_selector" value="gtf"/> |
409 <param name="gtf_file" value="Homo_sapiens.GRCh38.94.gtf" ftype="gtf"/> | 406 <param name="gtf_file" value="Homo_sapiens.GRCh38.94.gtf" ftype="gtf"/> |
410 <output name="deseq_out" > | 407 <output name="deseq_out" > |
435 <repeat name="rep_factorLevel"> | 432 <repeat name="rep_factorLevel"> |
436 <param name="factorLevel" value="Untreated"/> | 433 <param name="factorLevel" value="Untreated"/> |
437 <param name="groups" value="primary:untreated"/> | 434 <param name="groups" value="primary:untreated"/> |
438 </repeat> | 435 </repeat> |
439 </repeat> | 436 </repeat> |
440 <param name="pdf" value="False"/> | 437 <section name="output_options"> |
438 <param name="output_selector" value=""/> | |
439 </section> | |
441 <param name="tximport_selector" value="tximport"/> | 440 <param name="tximport_selector" value="tximport"/> |
442 <param name="txtype" value="sailfish"/> | 441 <param name="txtype" value="sailfish"/> |
443 <param name="mapping_format_selector" value="tabular"/> | 442 <param name="mapping_format_selector" value="tabular"/> |
444 <param name="tabular_file" value="tx2gene.tab"/> | 443 <param name="tabular_file" value="tx2gene.tab"/> |
445 <output name="deseq_out" > | 444 <output name="deseq_out" > |
470 <repeat name="rep_factorLevel"> | 469 <repeat name="rep_factorLevel"> |
471 <param name="factorLevel" value="Untreated"/> | 470 <param name="factorLevel" value="Untreated"/> |
472 <param name="groups" value="primary:untreated"/> | 471 <param name="groups" value="primary:untreated"/> |
473 </repeat> | 472 </repeat> |
474 </repeat> | 473 </repeat> |
475 <param name="pdf" value="False"/> | 474 <section name="output_options"> |
475 <param name="output_selector" value="many_contrasts"/> | |
476 </section> | |
476 <param name="tximport_selector" value="tximport"/> | 477 <param name="tximport_selector" value="tximport"/> |
477 <param name="txtype" value="sailfish"/> | 478 <param name="txtype" value="sailfish"/> |
478 <param name="mapping_format_selector" value="tabular"/> | 479 <param name="mapping_format_selector" value="tabular"/> |
479 <param name="tabular_file" value="tx2gene.tab"/> | 480 <param name="tabular_file" value="tx2gene.tab"/> |
480 <param name="many_contrasts" value="true"/> | |
481 <output_collection name="split_output" type="list" count="1"> | 481 <output_collection name="split_output" type="list" count="1"> |
482 <element name="Treatment_Treated_vs_Untreated"> | 482 <element name="Treatment_Treated_vs_Untreated"> |
483 <assert_contents> | 483 <assert_contents> |
484 <has_text_matching expression="HOXC11\t0.557.*\t0.324.*\t0.437.*\t0.741.*\t0.458.*\t0.999.*"/> | 484 <has_text_matching expression="HOXC11\t0.557.*\t0.324.*\t0.437.*\t0.741.*\t0.458.*\t0.999.*"/> |
485 </assert_contents> | 485 </assert_contents> |
486 </element> | 486 </element> |
487 </output_collection> | 487 </output_collection> |
488 </test> | |
489 <!--Test alpha_ma option--> | |
490 <test expect_num_outputs="1"> | |
491 <repeat name="rep_factorName"> | |
492 <param name="factorName" value="Treatment"/> | |
493 <repeat name="rep_factorLevel"> | |
494 <param name="factorLevel" value="Treated"/> | |
495 <param name="countsFile" value="sailfish_ensembl/sailfish_quant.sf1.tab,sailfish_ensembl/sailfish_quant.sf2.tab,sailfish_ensembl/sailfish_quant.sf3.tab"/> | |
496 </repeat> | |
497 <repeat name="rep_factorLevel"> | |
498 <param name="factorLevel" value="Untreated"/> | |
499 <param name="countsFile" value="sailfish_ensembl/sailfish_quant.sf4.tab,sailfish_ensembl/sailfish_quant.sf5.tab,sailfish_ensembl/sailfish_quant.sf6.tab"/> | |
500 </repeat> | |
501 </repeat> | |
502 <section name="output_options"> | |
503 <param name="output_selector" value=""/> | |
504 <param name="alpha_ma" value="0.05"/> | |
505 </section> | |
506 <param name="tximport_selector" value="tximport"/> | |
507 <param name="txtype" value="sailfish"/> | |
508 <param name="mapping_format_selector" value="gtf"/> | |
509 <param name="gtf_file" value="Homo_sapiens.GRCh38.94.gtf" ftype="gtf"/> | |
510 <output name="deseq_out" > | |
511 <assert_contents> | |
512 <has_text_matching expression="ENSG00000168671\t1.8841.*\t-0.1180.*\t0.7429.*\t-0.1589.*\t0.8737.*\t0.9999.*" /> | |
513 </assert_contents> | |
514 </output> | |
488 </test> | 515 </test> |
489 </tests> | 516 </tests> |
490 <help><![CDATA[ | 517 <help><![CDATA[ |
491 .. class:: infomark | 518 .. class:: infomark |
492 | 519 |
582 ====== ========================================================== | 609 ====== ========================================================== |
583 | 610 |
584 .. _DESeq2: http://master.bioconductor.org/packages/release/bioc/html/DESeq2.html | 611 .. _DESeq2: http://master.bioconductor.org/packages/release/bioc/html/DESeq2.html |
585 .. _tximport: https://bioconductor.org/packages/devel/bioc/vignettes/tximport/inst/doc/tximport.html | 612 .. _tximport: https://bioconductor.org/packages/devel/bioc/vignettes/tximport/inst/doc/tximport.html |
586 ]]></help> | 613 ]]></help> |
587 <citations> | 614 <expand macro="citations" /> |
588 <citation type="doi">10.1186/s13059-014-0550-8</citation> | |
589 </citations> | |
590 </tool> | 615 </tool> |