1
|
1 <tool id="rgDGE" name="Digital differential expression" version="0.01">
|
|
2 <description>for RNA counts</description>
|
|
3 <command interpreter="python">rgDGE.py --output_html "$html_file.files_path" --input_matrix "$input1" --treatcols "$Treat_cols"
|
|
4 --treatname "$treatment_name" --ctrlcols "$Control_cols" --ctrlname "$control_name" --output_tab "$outtab" --output_html
|
|
5 "$html_file" --output_dir "$html_file.files_path" --method "edgeR"
|
|
6 </command>
|
|
7 <inputs>
|
|
8 <param name="input1" type="data" format="tabular" label="Select an input matrix - rows are contigs, columns are counts for each sample"
|
|
9 help="Use the DGE matrix preparation tool to create these matrices from BAM files and a BED file of contigs"/>
|
|
10 <param name="treatment_name" type="text" value="Treatment" size="50" label="Treatment Name"/>
|
|
11 <param name="Treat_cols" label="Select columns containing treatment." type="data_column" data_ref="input1" numerical="True"
|
|
12 multiple="true" use_header_names="true" size="120" display="checkboxes">
|
|
13 <validator type="no_options" message="Please select at least one column."/>
|
|
14 </param>
|
|
15
|
|
16 <param name="control_name" type="text" value="Control" size="50" label="Control Name"/>
|
|
17 <param name="Control_cols" label="Select columns containing control." type="data_column" data_ref="input1" numerical="True"
|
|
18 multiple="true" use_header_names="true" size="120" display="checkboxes">
|
|
19 <validator type="no_options" message="Please select at least one column."/>
|
|
20 </param>
|
|
21
|
|
22 <param name="method" type="select" label="Method" >
|
|
23 <option value="edgeR" selected="true">edgeR</option>
|
|
24 <option value="DESeq">DESeq - coming soonish!</option>
|
|
25 </param>
|
|
26
|
|
27 </inputs>
|
|
28 <outputs>
|
|
29 <data format="tabular" name="outtab" label="${on_string}_DGE.xls"/>
|
|
30 <data format="html" name="html_file" label="${on_string}_DGE.html"/>
|
|
31 </outputs>
|
|
32 <tests>
|
|
33 <test>
|
|
34 <param name='input1' value='DGEtest.xls' ftype='tabular' />
|
|
35 <param name='treatment_name' value='case' />
|
|
36 <param name='control_name' value='control' />
|
|
37 <param name='Treat_cols' value='c3,c6,c9' />
|
|
38 <param name='Control_cols' value='c2,c5,c8' />
|
|
39 <output name='outtab' file='DGEtest1out.xls' ftype='tabular' compare='diff' />
|
|
40 <output name='html_file' file='DGEtest1out.html' ftype='html' compare='diff' lines_diff='20' />
|
|
41 </test>
|
|
42 </tests>
|
|
43 <help>
|
|
44 .. class:: infomark
|
|
45
|
|
46 **What it does**
|
|
47
|
|
48 Performs digital differential gene expression analysis between two groups (eg a treatment and control) each with biological replicates, from count data arranged in a matrix as a tabular file.
|
|
49
|
|
50 **Input**
|
|
51
|
|
52 A matrix consisting of non-negative integers. The matrix must have a unique header row identifiying the samples, as well as a unique set of row names
|
|
53 as the first column. To create appropriate input files, the "Multiple BAMs to matrix" tool takes multiple bam files containing mapped RNAseq reads mapped to a reference genome, and a BED file containing (eg Refseq) contigs. A count matrix containing contigs as rows and samples as
|
|
54 columns will be created and can then be used in this tool to estimate differential gene expression using methods appropriate for count data.
|
|
55
|
|
56 **Output**
|
|
57
|
|
58 A tabular file containing relative expression levels, statistical estimates of differential expression probability, the R script, log and sessionInfo, and some helpful diagnostic plots.
|
|
59
|
|
60 .. class:: infomark
|
|
61
|
|
62 **Attribution**
|
|
63 DGE.R was originally written by: S.Lunke and A.Kaspi.
|
|
64 Wrapping for Galaxy and general kibbitzing from Ross Lazarus.
|
|
65 This tool wrapts the edgeR_ Bioconductor package so all calculations and plots are controlled by that code. See edgeR_ for all documentation and appropriate attribution.
|
|
66 Recommended reference is Mark D. Robinson, Davis J. McCarthy, Gordon K. Smyth, PMCID: PMC2796818
|
|
67
|
|
68 .. _LGPL: http://www.gnu.org/copyleft/lesser.html
|
|
69 .. _edgeR: http://www.bioconductor.org/packages/2.3/bioc/html/edgeR.html
|
|
70
|
|
71 </help>
|
|
72
|
|
73 </tool>
|
|
74
|
|
75
|