|
2
|
1 <tool id="sartools_edger" name="SARTools edgeR" version="@WRAPPER_VERSION@.0">
|
|
0
|
2
|
|
|
3 <!-- [REQUIRED] Tool description displayed after the tool name -->
|
|
|
4 <description>Compare two or more biological conditions in a RNA-Seq framework with edgeR</description>
|
|
|
5
|
|
|
6 <macros>
|
|
|
7 <import>macros.xml</import>
|
|
|
8 </macros>
|
|
|
9
|
|
|
10 <expand macro="requirements"/>
|
|
|
11 <expand macro="stdio"/>
|
|
|
12
|
|
|
13
|
|
|
14 <!-- [REQUIRED] The command to execute -->
|
|
|
15 <command interpreter="python"><![CDATA[
|
|
|
16
|
|
|
17 abims_sartools_edger_wrapper.py
|
|
|
18 ## parameters
|
|
|
19 @COMMAND_BASIC_PARAMETERS@
|
|
|
20 #if str( $advanced_parameters.adv_param ) == "show":
|
|
|
21 @COMMAND_BATCH_PARAM@
|
|
|
22 --alpha $advanced_parameters.alpha
|
|
|
23 --pAdjustMethod $advanced_parameters.pAdjustMethod
|
|
|
24 --cpmCutoff $advanced_parameters.cpmCutoff
|
|
|
25 --geneSelection $advanced_parameters.geneSelection
|
|
|
26 --normalizationMethod $advanced_parameters.normalizationMethod
|
|
|
27 --colors $advanced_parameters.colors
|
|
3
|
28 --forceCairoGraph $advanced_parameters.forceCairoGraph
|
|
0
|
29 #end if
|
|
|
30 ## ouputs
|
|
|
31 @COMMAND_OUTPUTS@
|
|
|
32
|
|
|
33 ]]></command>
|
|
|
34
|
|
|
35 <!-- [REQUIRED] Input files and tool parameters -->
|
|
|
36 <inputs>
|
|
|
37
|
|
|
38 <expand macro="basic_parameters" />
|
|
|
39
|
|
|
40 <conditional name="advanced_parameters" >
|
|
|
41 <param name="adv_param" type="select" label="Advanced Parameters" help="" >
|
|
|
42 <option value="hide" selected="true">Hide</option>
|
|
|
43 <option value="show">Show</option>
|
|
|
44 </param>
|
|
|
45 <when value="hide" />
|
|
|
46 <when value="show">
|
|
|
47 <expand macro="batch_param" />
|
|
|
48 <expand macro="alpha_param" />
|
|
|
49 <expand macro="padjustmethod_param" />
|
|
3
|
50 <param type="integer" value="1" min="0" label="Counts-per-million cut-off to filter low counts" argument="--cpmCutoff" help="Set to 0 to disable filtering. Default is 1." />
|
|
|
51 <param name="geneSelection" type="select" label="Selection of the features in MDSPlot" argument="--gene.selection" help="Default is 'pairwise'." >
|
|
0
|
52 <option value="pairwise" selected="true">pairwise</option>
|
|
|
53 <option value="common">common</option>
|
|
|
54 </param>
|
|
3
|
55 <param type="select" label="Normalization method in calcNormFactors" argument="--normalizationMethod" help="'TMM' (default), 'RLE' (DESeq method) or 'upperquartile'." >
|
|
0
|
56 <option value="TMM" selected="true">TMM</option>
|
|
|
57 <option value="RLE">RLE</option>
|
|
|
58 <option value="upperquartile">upperquartile</option>
|
|
|
59 </param>
|
|
|
60 <expand macro="colors_param" />
|
|
3
|
61 <expand macro="forceCairoGraph_param" />
|
|
0
|
62 </when>
|
|
|
63 </conditional>
|
|
|
64
|
|
|
65 </inputs>
|
|
|
66
|
|
|
67 <!-- [REQUIRED] Output files -->
|
|
|
68 <outputs>
|
|
|
69
|
|
|
70 <expand macro="outputs" />
|
|
|
71
|
|
|
72 </outputs>
|
|
|
73
|
|
|
74 <!-- [OPTIONAL] Tests to be run manually by the Galaxy admin -->
|
|
|
75 <tests>
|
|
|
76 <!-- [HELP] Test files have to be in the ~/test-data directory -->
|
|
|
77 <test>
|
|
|
78 <!-- Test with 2 conditions, 2 replicates, 10 features -->
|
|
|
79 <param name="targetFile" dbkey="?" value="target_small.txt" />
|
|
|
80 <param name="rawDir" value="raw_small.zip" dbkey="?" ftype="zip"/>
|
|
|
81 <param name="adv_param" value="show"/>
|
|
|
82 <output name="log">
|
|
|
83 <assert_contents>
|
|
|
84 <has_text text="KO vs WT 5 4 9" />
|
|
|
85 <has_text text="HTML report created" />
|
|
|
86 </assert_contents>
|
|
|
87 </output>
|
|
|
88 </test>
|
|
|
89 <!-- <test>
|
|
|
90 --> <!-- NOT WORKING YET: Test with 3 conditions, 3 replicates, 10 features, with batch effect -->
|
|
|
91 <!-- <param name="targetFile" dbkey="?" value="targetT048_small.txt" />
|
|
|
92 <param name="rawDir" value="rawT048_small.zip" dbkey="?" ftype="no_unzip.zip"/>
|
|
|
93 <param name="condRef" value="T0"/>
|
|
|
94 <param name="adv_param" value="show"/>
|
|
|
95 <param name="condition" value="true"/>
|
|
|
96 <output name="tables_html" file="SARTools_edgeR_targetT048_small_tables.html" lines_diff="12">
|
|
|
97 <extra_files type="file" name="T4vsT0.complete.txt" value="SARTools_edgeR_T4vsT0_small.complete.txt"/>
|
|
|
98 <extra_files type="file" name="T8vsT0.complete.txt" value="SARTools_edgeR_T8vsT0_small.complete.txt"/>
|
|
|
99 <extra_files type="file" name="T8vsT4.complete.txt" value="SARTools_edgeR_T8vsT4_small.complete.txt"/>
|
|
|
100 </output>
|
|
|
101 </test>
|
|
|
102 --> <test>
|
|
|
103 <!-- Test with 2 conditions, 2 replicates, 8217 features -->
|
|
|
104 <param name="targetFile" dbkey="?" value="target.txt" />
|
|
|
105 <param name="rawDir" value="raw.zip" dbkey="?" ftype="zip"/>
|
|
|
106 <param name="adv_param" value="show"/>
|
|
|
107 <output name="log">
|
|
|
108 <assert_contents>
|
|
|
109 <has_text text="KO vs WT 2691 2713 5404" />
|
|
|
110 <has_text text="HTML report created" />
|
|
|
111 </assert_contents>
|
|
|
112 </output>
|
|
|
113 </test>
|
|
|
114 <!-- <test>
|
|
|
115 --> <!-- NOT WORKING YET: Test with 3 conditions, 3 replicates, 10160 features, with batch effect -->
|
|
|
116 <!-- <param name="targetFile" dbkey="?" value="targetT048.txt" />
|
|
|
117 <param name="rawDir" value="rawT048.zip" dbkey="?" ftype="no_unzip.zip"/>
|
|
|
118 <param name="condRef" value="T0"/>
|
|
|
119 <param name="adv_param" value="show"/>
|
|
|
120 <param name="condition" value="true"/>
|
|
|
121 <output name="tables_html" file="SARTools_edgeR_targetT048_tables.html" lines_diff="14">
|
|
|
122 <extra_files type="file" name="T4vsT0.complete.txt" value="SARTools_edgeR_T4vsT0.complete.txt"/>
|
|
|
123 <extra_files type="file" name="T8vsT0.complete.txt" value="SARTools_edgeR_T8vsT0.complete.txt"/>
|
|
|
124 <extra_files type="file" name="T8vsT4.complete.txt" value="SARTools_edgeR_T8vsT4.complete.txt"/>
|
|
|
125 </output>
|
|
|
126 </test>
|
|
|
127 --> </tests>
|
|
|
128
|
|
|
129 <!-- [OPTIONAL] Help displayed in Galaxy -->
|
|
|
130 <help><![CDATA[
|
|
|
131
|
|
|
132 @HELP_AUTHORS@
|
|
|
133
|
|
|
134 ==============
|
|
|
135 SARTools edgeR
|
|
|
136 ==============
|
|
|
137
|
|
|
138 -----------
|
|
|
139 Description
|
|
|
140 -----------
|
|
|
141
|
|
|
142 @HELP_DESCRIPTION@
|
|
|
143
|
|
|
144
|
|
|
145 -----------
|
|
|
146 Input files
|
|
|
147 -----------
|
|
|
148
|
|
|
149 @HELP_INPUT_FILES@
|
|
|
150
|
|
|
151
|
|
|
152 ----------
|
|
|
153 Parameters
|
|
|
154 ----------
|
|
|
155
|
|
|
156 @HELP_BASIC_PARAMETERS@
|
|
|
157 * **batch:** adjustment variable to use as a batch effect, must be a column of the target file (NULL if no batch effect needs to be taken into account);
|
|
|
158 * **alpha:** significance threshold applied to the adjusted p-values to select the differentially expressed features (default is 0.05);
|
|
|
159 * **pAdjustMethod:** p-value adjustment method for multiple testing [4, 5] ("BH" by default, "BY" or any value of p.adjust.methods);
|
|
|
160 * **cpmCutoff:** counts-per-million cut-off to filter low counts (default is 1, set to 0 to disable filtering);
|
|
|
161 * **gene.selection:** method of selection of the features for the MultiDimensional Scaling plot ("pairwise" by default or common);
|
|
|
162 * **normalizationMethod:** normalization method in calcNormFactors(): "TMM" (default), "RLE" (DESeq method) or "upperquartile";
|
|
|
163 * **colors:** colors used for the figures (one per biological condition), 8 are given by default.
|
|
3
|
164 * **forceCairoGraph:** TRUE or FALSE (default) to force the use of cairo with options(bitmapType="cairo").
|
|
0
|
165
|
|
|
166
|
|
|
167 ------------
|
|
|
168 Output files
|
|
|
169 ------------
|
|
|
170
|
|
|
171 @HELP_OUTPUT_FILES@
|
|
|
172
|
|
|
173
|
|
|
174 ---------------------------------------------------
|
|
|
175
|
|
|
176 [1] G.-K. Smyth. Limma: linear models for microarray data. In R. Gentleman, V. Carey, S. Dudoit, R. Irizarry, and W. Huber, editors, Bioinformatics and Computational Biology Solutions Using R and Bioconductor, pages 397–420. Springer, New York, 2005.
|
|
|
177
|
|
|
178 [2] S. Anders. HTSeq: Analysing high-throughput sequencing data with Python. http://www-huber.embl.de/users/anders/HTSeq/, 2011.
|
|
|
179
|
|
|
180 [3] S. Anders, P.-T. Pyl, and W. Huber. HTSeq - A Python framework to work with high-throughput sequencing data. bioRxiv preprint, 2014. URL: http://dx.doi.org/10.1101/002824.
|
|
|
181
|
|
|
182 [4] Y. Benjamini and Y. Hochberg. Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society B, 57:289–300, 1995.
|
|
|
183
|
|
|
184 [5] Y. Benjamini and D. Yekutieli. The control of the false discovery rate in multiple testing under dependency. Ann. Statist., 29(4):1165–1188, 2001.
|
|
|
185
|
|
|
186
|
|
|
187 ]]></help>
|
|
|
188
|
|
|
189 <citations>
|
|
|
190 <expand macro="common_citations" />
|
|
|
191 </citations>
|
|
|
192
|
|
|
193 </tool>
|
