Mercurial > repos > yufei-luo > s_mart
diff SMART/galaxy/GetDifferentialExpression.xml @ 15:440ceca58672
Uploaded
author | m-zytnicki |
---|---|
date | Mon, 22 Apr 2013 11:08:07 -0400 |
parents | 769e306b7933 |
children | 94ab73e8a190 |
line wrap: on
line diff
--- a/SMART/galaxy/GetDifferentialExpression.xml Fri Apr 19 10:13:11 2013 -0400 +++ b/SMART/galaxy/GetDifferentialExpression.xml Mon Apr 22 11:08:07 2013 -0400 @@ -51,7 +51,6 @@ $simple $adjusted - #if $optionSimplePara.simplePara == 'Yes': -S $optionSimplePara.paraValue #end if @@ -63,7 +62,6 @@ #if $optionFDR.FDR == 'Yes': -d $optionFDR.FDRValue #end if - $plot $outputFilePNG </command> <inputs> @@ -154,8 +152,8 @@ </when> </conditional> - <param name="simple" type="boolean" truevalue="-s" falsevalue="" checked="false" label="normalize using the number of reads in each condition"/> - <param name="adjusted" type="boolean" truevalue="-a" falsevalue="" checked="false" label="normalize using the number of reads of 'mean' regions"/> + <param name="simple" type="boolean" truevalue="-s" falsevalue="" checked="false" label="Normalize using the number of reads in each condition"/> + <param name="adjusted" type="boolean" truevalue="-a" falsevalue="" checked="false" label="Normalize using the number of reads of interquartile expression region"/> <conditional name="optionSimplePara"> <param name="simplePara" type="select" label="provide the number of reads" > @@ -170,7 +168,7 @@ </conditional> <conditional name="optionFixedSizeFactor"> - <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)"> + <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)"> <option value="Yes">Yes</option> <option value="No" selected="true">No</option> </param> @@ -193,18 +191,21 @@ </when> </conditional> - <param name="plot" type="boolean" truevalue="-p" falsevalue="" checked="false" label="plot option" help="plot cloud plot"/> - </inputs> <outputs> <data name="outputFileGff" format="gff3" label="[GetDifferentialExpression]out file"/> - <data name="outputFilePNG" format="PNG" label="[GetDifferentialExpression]PNG file"> - <filter>plot</filter> - </data> </outputs> <help> - example: python GetDifferentialExpression.py -i input1 -f gff3 -j input2 -g gff3 -k ref -l gff3 -o output.gff3 +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. + +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. + +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. + +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. + +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. </help> </tool>