Mercurial > repos > yufei-luo > s_mart

comparison SMART/galaxy/GetDifferentialExpression.xml @ 31:0ab839023fe4

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author m-zytnicki
date Tue, 30 Apr 2013 14:33:21 -0400
parents 94ab73e8a190
children
comparison
equal deleted inserted replaced
30:5677346472b5 31:0ab839023fe4
1 <tool id="GetDifferentialExpression" name="get differential expression"> 1 <tool id="GetDifferentialExpression" name="get differential expression">
2 <description>Get the differential expression between 2 conditions using Fisher's exact test, on regions defined by a third file.</description> 2 <description>Get the differential expression between 2 conditions using Fisher's exact test, on regions defined by a third file.</description>
3 <requirements>
4 <requirement type="set_environment">PYTHONPATH</requirement>
5 </requirements>
6 <command interpreter="python"> 3 <command interpreter="python">
7 ../Java/Python/GetDifferentialExpression.py -i $formatType.inputFileName1 4 ../Java/Python/GetDifferentialExpression.py -i $formatType.inputFileName1
8 #if $formatType.FormatInputFileName1 == 'bed': 5 #if $formatType.FormatInputFileName1 == 'bed':
9 -f bed 6 -f bed
10 #elif $formatType.FormatInputFileName1 == 'gff': 7 #elif $formatType.FormatInputFileName1 == 'gff':
52 -o $outputFileGff 49 -o $outputFileGff
53 50
54 $simple 51 $simple
55 $adjusted 52 $adjusted
56 53
54
57 #if $optionSimplePara.simplePara == 'Yes': 55 #if $optionSimplePara.simplePara == 'Yes':
58 -S $optionSimplePara.paraValue 56 -S $optionSimplePara.paraValue
59 #end if 57 #end if
60 58
61 #if $optionFixedSizeFactor.FSF == 'Yes': 59 #if $optionFixedSizeFactor.FSF == 'Yes':
63 #end if 61 #end if
64 62
65 #if $optionFDR.FDR == 'Yes': 63 #if $optionFDR.FDR == 'Yes':
66 -d $optionFDR.FDRValue 64 -d $optionFDR.FDRValue
67 #end if 65 #end if
66 $plot $outputFilePNG
68 </command> 67 </command>
69 68
70 <inputs> 69 <inputs>
71 <conditional name="formatType"> 70 <conditional name="formatType">
72 <param name="FormatInputFileName1" type="select" label="Input File Format 1"> 71 <param name="FormatInputFileName1" type="select" label="Input File Format 1">
153 <when value="gtf"> 152 <when value="gtf">
154 <param name="inputFileNameRef" format="gtf" type="data" label="Input Ref File"/> 153 <param name="inputFileNameRef" format="gtf" type="data" label="Input Ref File"/>
155 </when> 154 </when>
156 </conditional> 155 </conditional>
157 156
158 <param name="simple" type="boolean" truevalue="-s" falsevalue="" checked="false" label="Normalize using the number of reads in each condition"/> 157 <param name="simple" type="boolean" truevalue="-s" falsevalue="" checked="false" label="normalize using the number of reads in each condition"/>
159 <param name="adjusted" type="boolean" truevalue="-a" falsevalue="" checked="false" label="Normalize using the number of reads of interquartile expression region"/> 158 <param name="adjusted" type="boolean" truevalue="-a" falsevalue="" checked="false" label="normalize using the number of reads of 'mean' regions"/>
160 159
161 <conditional name="optionSimplePara"> 160 <conditional name="optionSimplePara">
162 <param name="simplePara" type="select" label="provide the number of reads" > 161 <param name="simplePara" type="select" label="provide the number of reads" >
163 <option value="Yes">Yes</option> 162 <option value="Yes">Yes</option>
164 <option value="No" selected="true">No</option> 163 <option value="No" selected="true">No</option>
169 <when value="No"> 168 <when value="No">
170 </when> 169 </when>
171 </conditional> 170 </conditional>
172 171
173 <conditional name="optionFixedSizeFactor"> 172 <conditional name="optionFixedSizeFactor">
174 <param name="FSF" type="select" label="Give the magnification factor for the normalization using fixed size sliding windows in reference regions (leave empty for no such normalization)"> 173 <param name="FSF" type="select" label="give the magnification factor for the normalization using fixed size sliding windows in reference regions (leave empty for no such normalization)">
175 <option value="Yes">Yes</option> 174 <option value="Yes">Yes</option>
176 <option value="No" selected="true">No</option> 175 <option value="No" selected="true">No</option>
177 </param> 176 </param>
178 <when value="Yes"> 177 <when value="Yes">
179 <param name="FSFValue" type="integer" value="0" /> 178 <param name="FSFValue" type="integer" value="0" />
192 </when> 191 </when>
193 <when value="No"> 192 <when value="No">
194 </when> 193 </when>
195 </conditional> 194 </conditional>
196 195
196 <param name="plot" type="boolean" truevalue="-p" falsevalue="" checked="false" label="plot option" help="plot cloud plot"/>
197
197 </inputs> 198 </inputs>
198 199
199 <outputs> 200 <outputs>
200 <data name="outputFileGff" format="gff3" label="[GetDifferentialExpression]out file"/> 201 <data name="outputFileGff" format="gff3" label="[GetDifferentialExpression]out file"/>
202 <data name="outputFilePNG" format="PNG" label="[GetDifferentialExpression]PNG file">
203 <filter>plot</filter>
204 </data>
201 </outputs> 205 </outputs>
202 206
203 <help> 207 <help>
204 This tool compares two sets of data and find the differential expression. One very important component of the tool is the reference set. Actually, to use the tool, you need the two input sets of data, of course, and the reference set. The reference set is a set of genomic coordinates and, for each interval, it will count the number of feature on each sample and compute the differential expression. For each reference interval, it will output the direction of the regulation (up or down, with respect to the first input set), and a *p*-value from a Fisher exact test. 208 example: python GetDifferentialExpression.py -i input1 -f gff3 -j input2 -g gff3 -k ref -l gff3 -o output.gff3
205
206 This reference set seems boring. Why not computing the differential expression without this set? The answer is: the differential expression of what? I cannot guess it. Actually, you might want to compare the expression of genes, of small RNAs, of transposable elements, of anything... So the reference set can be a list of genes, and in this case, you can compute the differential expression of genes. But you can also compute many other things.
207
208 Suppose that you cluster the data of your two input samples (you can do it with the *clusterize* and the *mergeTranscriptLists* tools). You now have a list of all the regions which are transcribed in at least one of the input samples. This can be your reference set. This reference set is interesting since you can detect the differential expression of data which is outside any annotation.
209
210 Suppose now that you clusterize using a sliding window the two input samples (you can do it with the *clusterizeBySlidingWindows* and the *mergeSlidingWindowsClusters* tools). You can now select all the regions of a given size which contain at least one read in one of the two input samples (do it with *selectByTag* and the tag **nbElements**). Again, this can be an other interesting reference set.
211
212 In most cases, the sizes of the two input samples will be different, so you should probably normalize the data, which is an available option. The ---rather crude--- normalization increases the number of data in the least populated sample and decreases the number of data in the most populated sample to the average number of data.
213 </help> 209 </help>
214 </tool> 210 </tool>