Mercurial > repos > jjohnson > cummerbund
annotate cuffdiff_wrapper.xml @ 6:137aab1d9ac1
Add metadata to datatype: CuffDataDB
author | Jim Johnson <jj@umn.edu> |
---|---|
date | Mon, 18 Nov 2013 16:43:15 -0600 |
parents | f109453ecfa2 |
children | b5562b9a55c7 |
rev | line source |
---|---|
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
1 <tool id="cuffdiff_cummerbund" name="Cuffdiff for cummeRbund" version="0.0.7"> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
2 <!-- Wrapper supports Cuffdiff versions 2.1.0-2.1.1 --> |
0 | 3 <description>find significant changes in transcript expression, splicing, and promoter use</description> |
4 <requirements> | |
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
5 <requirement type="package" version="2.1.1">cufflinks</requirement> |
0 | 6 </requirements> |
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
7 <version_command>cuffdiff 2>&1 | head -n 1</version_command> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
8 <command> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
9 cuffdiff |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
10 --no-update-check |
0 | 11 --FDR=$fdr |
12 --num-threads="4" | |
13 --min-alignment-count=$min_alignment_count | |
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
14 --library-norm-method=$library_norm_method |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
15 --dispersion-method=$dispersion_method |
0 | 16 |
17 ## Set advanced data parameters? | |
18 #if $additional.sAdditional == "Yes": | |
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
19 #if $additional.frag_mean_len: |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
20 -m $additional.frag_mean_len |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
21 #end if |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
22 #if $additional.frag_len_std_dev: |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
23 -s $additional.frag_len_std_dev |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
24 #end if |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
25 #if $additional.max_bundle_frags: |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
26 --max-bundle-frags="$additional.max_bundle_frags" |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
27 #end if |
0 | 28 #end if |
29 | |
30 ## Multi-read correct? | |
31 #if str($multiread_correct) == "Yes": | |
32 -u | |
33 #end if | |
34 | |
35 ## Bias correction? | |
36 #if $bias_correction.do_bias_correction == "Yes": | |
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
37 -b |
0 | 38 #if $bias_correction.seq_source.index_source == "history": |
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
39 ## Custom genome from history. |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
40 $bias_correction.seq_source.ref_file |
0 | 41 #else: |
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
42 ## Built-in genome. |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
43 ${__get_data_table_entry__('sam_fa_indexes', 'value', $gtf_input.dbkey, 'path')} |
0 | 44 #end if |
45 #end if | |
46 | |
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
47 #set labels = '\'' + '\',\''.join( [ str( $condition.name ) for $condition in $conditions ] ) + '\'' |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
48 --labels $labels |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
49 |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
50 ## Inputs. |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
51 $gtf_input |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
52 #for $condition in $conditions: |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
53 #set samples = ','.join( [ str( $sample.sample ) for $sample in $condition.samples ] ) |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
54 $samples |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
55 #end for |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
56 |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
57 ## If build cummerbund db |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
58 #if $build_cummerbund_db: |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
59 && echo 'library(cummeRbund)' > cuffData.r |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
60 #if $bias_correction.do_bias_correction == "Yes": |
3
f109453ecfa2
readCufflinks() parses feature file based on file extension: .gtf/.gff3 Also added gtfFile and genome input to cummerbund_wrapper.xml
Jim Johnson <jj@umn.edu>
parents:
2
diff
changeset
|
61 ## cummeRbund relies on the file extension to determine gff format |
f109453ecfa2
readCufflinks() parses feature file based on file extension: .gtf/.gff3 Also added gtfFile and genome input to cummerbund_wrapper.xml
Jim Johnson <jj@umn.edu>
parents:
2
diff
changeset
|
62 #set $gtf_link = '.'.join(['cuff',$gtf_input.extension]) |
f109453ecfa2
readCufflinks() parses feature file based on file extension: .gtf/.gff3 Also added gtfFile and genome input to cummerbund_wrapper.xml
Jim Johnson <jj@umn.edu>
parents:
2
diff
changeset
|
63 && ln -s $gtf_input $gtf_link |
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
64 #if $bias_correction.seq_source.index_source == "history": |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
65 ## Custom genome from history. |
3
f109453ecfa2
readCufflinks() parses feature file based on file extension: .gtf/.gff3 Also added gtfFile and genome input to cummerbund_wrapper.xml
Jim Johnson <jj@umn.edu>
parents:
2
diff
changeset
|
66 && echo 'cuff<-readCufflinks( dbFile = "cuffdata.db", gtfFile = "$gtf_link", genome = "$bias_correction.seq_source.ref_file", rebuild = T)' >> cuffData.r |
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
67 #else: |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
68 ## Built-in genome. |
3
f109453ecfa2
readCufflinks() parses feature file based on file extension: .gtf/.gff3 Also added gtfFile and genome input to cummerbund_wrapper.xml
Jim Johnson <jj@umn.edu>
parents:
2
diff
changeset
|
69 && echo 'cuff<-readCufflinks( dbFile = "cuffdata.db", gtfFile = "$gtf_link", genome = "${__get_data_table_entry__('sam_fa_indexes', 'value', $gtf_input.dbkey, 'path')}", rebuild = T)' >> cuffData.r |
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
70 #end if |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
71 #else |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
72 && echo 'cuff<-readCufflinks( dbFile = "cuffdata.db", rebuild = T)' >> cuffData.r |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
73 #end if |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
74 && Rscript --vanilla cuffData.r |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
75 && cp cuffdata.db $cummerbund_db |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
76 #end if |
0 | 77 </command> |
78 <inputs> | |
79 <param format="gtf,gff3" name="gtf_input" type="data" label="Transcripts" help="A transcript GFF3 or GTF file produced by cufflinks, cuffcompare, or other source."/> | |
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
80 |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
81 <repeat name="conditions" title="Condition" min="2"> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
82 <param name="name" title="Condition name" type="text" label="Name"/> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
83 <repeat name="samples" title="Replicate" min="1"> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
84 <param name="sample" label="Add replicate" type="data" format="sam,bam"/> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
85 </repeat> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
86 </repeat> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
87 <param name="time_series" type="boolean" checked="false" truevalue="--time-series" falsevalue="" optional="true" label="treat samples as a time-series"> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
88 <help> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
89 Instructs Cuffdiff to analyze the provided samples as a time series, rather than testing for differences between all pairs of samples. |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
90 Samples should be provided in increasing time order at the command line (e.g first time point SAM, second timepoint SAM, etc.) |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
91 </help> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
92 </param> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
93 |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
94 <param name="library_norm_method" type="select" label="Library normalization method"> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
95 <option value="geometric" selected="True">geometric</option> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
96 <option value="classic-fpkm">classic-fpkm</option> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
97 <option value="quartile">quartile</option> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
98 </param> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
99 |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
100 <param name="dispersion_method" type="select" label="Dispersion estimation method" help="If using only one sample per condition, you must use 'blind.'"> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
101 <option value="pooled" selected="True">pooled</option> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
102 <option value="per-condition">per-condition</option> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
103 <option value="blind">blind</option> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
104 </param> |
0 | 105 |
106 <param name="fdr" type="float" value="0.05" label="False Discovery Rate" help="The allowed false discovery rate."/> | |
107 | |
108 <param name="min_alignment_count" type="integer" value="10" label="Min Alignment Count" help="The minimum number of alignments in a locus for needed to conduct significance testing on changes in that locus observed between samples."/> | |
109 | |
110 <param name="multiread_correct" type="select" label="Use multi-read correct" help="Tells Cufflinks to do an initial estimation procedure to more accurately weight reads mapping to multiple locations in the genome."> | |
111 <option value="No" selected="true">No</option> | |
112 <option value="Yes">Yes</option> | |
113 </param> | |
114 | |
115 <conditional name="bias_correction"> | |
116 <param name="do_bias_correction" type="select" label="Perform Bias Correction" help="Bias detection and correction can significantly improve accuracy of transcript abundance estimates."> | |
117 <option value="No">No</option> | |
118 <option value="Yes">Yes</option> | |
119 </param> | |
120 <when value="Yes"> | |
121 <conditional name="seq_source"> | |
122 <param name="index_source" type="select" label="Reference sequence data"> | |
123 <option value="cached">Locally cached</option> | |
124 <option value="history">History</option> | |
125 </param> | |
126 <when value="cached"></when> | |
127 <when value="history"> | |
128 <param name="ref_file" type="data" format="fasta" label="Using reference file" /> | |
129 </when> | |
130 </conditional> | |
131 </when> | |
132 <when value="No"></when> | |
133 </conditional> | |
134 | |
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
135 <param name="include_read_group_files" type="select" label="Include Read Group Datasets" help="Read group datasets provide information on replicates."> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
136 <option value="No" selected="true">No</option> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
137 <option value="Yes">Yes</option> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
138 </param> |
3
f109453ecfa2
readCufflinks() parses feature file based on file extension: .gtf/.gff3 Also added gtfFile and genome input to cummerbund_wrapper.xml
Jim Johnson <jj@umn.edu>
parents:
2
diff
changeset
|
139 <param name="build_cummerbund_db" type="boolean" truevalue="yes" falsevalue="no" checked="false" label="Build cummeRbund database"> |
f109453ecfa2
readCufflinks() parses feature file based on file extension: .gtf/.gff3 Also added gtfFile and genome input to cummerbund_wrapper.xml
Jim Johnson <jj@umn.edu>
parents:
2
diff
changeset
|
140 <help>"Perform Bias Correction" must be selected to include gtf and genome reference in the cummeRbund database.</help> |
f109453ecfa2
readCufflinks() parses feature file based on file extension: .gtf/.gff3 Also added gtfFile and genome input to cummerbund_wrapper.xml
Jim Johnson <jj@umn.edu>
parents:
2
diff
changeset
|
141 </param> |
0 | 142 <conditional name="additional"> |
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
143 <param name="sAdditional" type="select" label="Set Additional Parameters? (not recommended for paired-end reads)"> |
0 | 144 <option value="No">No</option> |
145 <option value="Yes">Yes</option> | |
146 </param> | |
147 <when value="No"></when> | |
148 <when value="Yes"> | |
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
149 <param name="frag_mean_len" type="integer" value="" optional="true" label="Average Fragment Length Default: 200"> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
150 <help> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
151 Note: Cufflinks now learns the fragment length mean for each SAM file, |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
152 so using this option is no longer recommended with paired-end reads. |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
153 </help> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
154 </param> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
155 <param name="frag_len_std_dev" type="integer" value="" optional="true" label="Fragment Length Standard Deviation Default: 80"> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
156 <help> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
157 Note: Cufflinks now learns the fragment length mean for each SAM file, |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
158 so using this option is no longer recommended with paired-end reads. |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
159 </help> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
160 </param> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
161 <param name="max_bundle_frags" type="integer" value="" optional="true" label="--max-bundle-frags"> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
162 <help> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
163 Sets the maximum number of fragments a locus may have before being skipped. Skipped loci are listed in skipped.gtf. Default: 1000000 |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
164 </help> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
165 <validator type="in_range" message="Value greater than 0" min="1"/> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
166 </param> |
0 | 167 </when> |
168 </conditional> | |
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
169 </inputs> |
0 | 170 |
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
171 <stdio> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
172 <regex match="Error" source="both" level="fatal" description="Error"/> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
173 <regex match=".*" source="both" level="log" description="tool progress"/> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
174 </stdio> |
0 | 175 |
176 <outputs> | |
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
177 <!-- Optional read group datasets. --> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
178 <data format="cuffdatadb" name="cummerbund_db" label="${tool.name} on ${on_string}: cummeRbund sqlite Database" > |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
179 <filter>build_cummerbund_db</filter> |
0 | 180 </data> |
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
181 <data format="tabular" name="isoforms_read_group" label="${tool.name} on ${on_string}: isoforms read group tracking" from_work_dir="isoforms.read_group_tracking" > |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
182 <filter>(params['include_read_group_files'] == 'Yes'</filter> |
0 | 183 </data> |
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
184 <data format="tabular" name="genes_read_group" label="${tool.name} on ${on_string}: genes read group tracking" from_work_dir="genes.read_group_tracking" > |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
185 <filter>(params['include_read_group_files'] == 'Yes'</filter> |
0 | 186 </data> |
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
187 <data format="tabular" name="cds_read_group" label="${tool.name} on ${on_string}: CDs read group tracking" from_work_dir="cds.read_group_tracking" > |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
188 <filter>(params['include_read_group_files'] == 'Yes'</filter> |
0 | 189 </data> |
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
190 <data format="tabular" name="tss_groups_read_group" label="${tool.name} on ${on_string}: TSS groups read group tracking" from_work_dir="tss_groups.read_group_tracking" > |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
191 <filter>(params['include_read_group_files'] == 'Yes'</filter> |
0 | 192 </data> |
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
193 |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
194 <!-- Standard datasets. --> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
195 <data format="tabular" name="splicing_diff" label="${tool.name} on ${on_string}: splicing differential expression testing" from_work_dir="splicing.diff" /> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
196 <data format="tabular" name="promoters_diff" label="${tool.name} on ${on_string}: promoters differential expression testing" from_work_dir="promoters.diff" /> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
197 <data format="tabular" name="cds_diff" label="${tool.name} on ${on_string}: CDS overloading diffential expression testing" from_work_dir="cds.diff" /> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
198 <data format="tabular" name="cds_exp_fpkm_tracking" label="${tool.name} on ${on_string}: CDS FPKM differential expression testing" from_work_dir="cds_exp.diff" /> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
199 <data format="tabular" name="cds_fpkm_tracking" label="${tool.name} on ${on_string}: CDS FPKM tracking" from_work_dir="cds.fpkm_tracking" /> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
200 <data format="tabular" name="tss_groups_exp" label="${tool.name} on ${on_string}: TSS groups differential expression testing" from_work_dir="tss_group_exp.diff" /> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
201 <data format="tabular" name="tss_groups_fpkm_tracking" label="${tool.name} on ${on_string}: TSS groups FPKM tracking" from_work_dir="tss_groups.fpkm_tracking" /> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
202 <data format="tabular" name="genes_exp" label="${tool.name} on ${on_string}: gene differential expression testing" from_work_dir="gene_exp.diff" /> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
203 <data format="tabular" name="genes_fpkm_tracking" label="${tool.name} on ${on_string}: gene FPKM tracking" from_work_dir="genes.fpkm_tracking" /> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
204 <data format="tabular" name="isoforms_exp" label="${tool.name} on ${on_string}: transcript differential expression testing" from_work_dir="isoform_exp.diff" /> |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
205 <data format="tabular" name="isoforms_fpkm_tracking" label="${tool.name} on ${on_string}: transcript FPKM tracking" from_work_dir="isoforms.fpkm_tracking" /> |
0 | 206 </outputs> |
207 | |
208 <tests> | |
209 <test> | |
210 <!-- | |
211 cuffdiff cuffcompare_out5.gtf cuffdiff_in1.sam cuffdiff_in2.sam | |
212 --> | |
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
213 <!-- |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
214 NOTE: as of version 0.0.6 of the wrapper, tests cannot be run because multiple inputs to a repeat |
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
215 element are not supported. |
0 | 216 <param name="gtf_input" value="cuffcompare_out5.gtf" ftype="gtf" /> |
217 <param name="do_groups" value="No" /> | |
218 <param name="aligned_reads1" value="cuffdiff_in1.sam" ftype="sam" /> | |
219 <param name="aligned_reads2" value="cuffdiff_in2.sam" ftype="sam" /> | |
220 <param name="fdr" value="0.05" /> | |
221 <param name="min_alignment_count" value="0" /> | |
222 <param name="do_bias_correction" value="No" /> | |
223 <param name="do_normalization" value="No" /> | |
224 <param name="multiread_correct" value="No"/> | |
225 <param name="sAdditional" value="No"/> | |
226 <output name="splicing_diff" file="cuffdiff_out9.txt"/> | |
227 <output name="promoters_diff" file="cuffdiff_out10.txt"/> | |
228 <output name="cds_diff" file="cuffdiff_out11.txt"/> | |
229 <output name="cds_exp_fpkm_tracking" file="cuffdiff_out4.txt"/> | |
230 <output name="cds_fpkm_tracking" file="cuffdiff_out8.txt"/> | |
231 <output name="tss_groups_exp" file="cuffdiff_out3.txt" lines_diff="200"/> | |
232 <output name="tss_groups_fpkm_tracking" file="cuffdiff_out7.txt"/> | |
233 <output name="genes_exp" file="cuffdiff_out2.txt" lines_diff="200"/> | |
234 <output name="genes_fpkm_tracking" file="cuffdiff_out6.txt" lines_diff="200"/> | |
235 <output name="isoforms_exp" file="cuffdiff_out1.txt" lines_diff="200"/> | |
236 <output name="isoforms_fpkm_tracking" file="cuffdiff_out5.txt" lines_diff="200"/> | |
2
fdf01b3c1841
Update to new cuffdiff wrapper, add cuffdb_info.txt to cummerbund html output
Jim Johnson <jj@umn.edu>
parents:
0
diff
changeset
|
237 --> |
0 | 238 </test> |
239 </tests> | |
240 | |
241 <help> | |
242 **Cuffdiff Overview** | |
243 | |
244 Cuffdiff is part of Cufflinks_. Cuffdiff find significant changes in transcript expression, splicing, and promoter use. Please cite: Trapnell C, Williams BA, Pertea G, Mortazavi AM, Kwan G, van Baren MJ, Salzberg SL, Wold B, Pachter L. Transcript assembly and abundance estimation from RNA-Seq reveals thousands of new transcripts and switching among isoforms. Nature Biotechnology doi:10.1038/nbt.1621 | |
245 | |
246 .. _Cufflinks: http://cufflinks.cbcb.umd.edu/ | |
247 | |
248 ------ | |
249 | |
250 **Know what you are doing** | |
251 | |
252 .. class:: warningmark | |
253 | |
254 There is no such thing (yet) as an automated gearshift in expression analysis. It is all like stick-shift driving in San Francisco. In other words, running this tool with default parameters will probably not give you meaningful results. A way to deal with this is to **understand** the parameters by carefully reading the `documentation`__ and experimenting. Fortunately, Galaxy makes experimenting easy. | |
255 | |
256 .. __: http://cufflinks.cbcb.umd.edu/manual.html#cuffdiff | |
257 | |
258 ------ | |
259 | |
260 **Input format** | |
261 | |
262 Cuffdiff takes Cufflinks or Cuffcompare GTF files as input along with two SAM files containing the fragment alignments for two or more samples. | |
263 | |
264 ------ | |
265 | |
266 **Outputs** | |
267 | |
268 Cuffdiff produces many output files: | |
269 | |
270 1. Transcript FPKM expression tracking. | |
271 2. Gene FPKM expression tracking; tracks the summed FPKM of transcripts sharing each gene_id | |
272 3. Primary transcript FPKM tracking; tracks the summed FPKM of transcripts sharing each tss_id | |
273 4. Coding sequence FPKM tracking; tracks the summed FPKM of transcripts sharing each p_id, independent of tss_id | |
274 5. Transcript differential FPKM. | |
275 6. Gene differential FPKM. Tests difference sin the summed FPKM of transcripts sharing each gene_id | |
276 7. Primary transcript differential FPKM. Tests difference sin the summed FPKM of transcripts sharing each tss_id | |
277 8. Coding sequence differential FPKM. Tests difference sin the summed FPKM of transcripts sharing each p_id independent of tss_id | |
278 9. Differential splicing tests: this tab delimited file lists, for each primary transcript, the amount of overloading detected among its isoforms, i.e. how much differential splicing exists between isoforms processed from a single primary transcript. Only primary transcripts from which two or more isoforms are spliced are listed in this file. | |
279 10. Differential promoter tests: this tab delimited file lists, for each gene, the amount of overloading detected among its primary transcripts, i.e. how much differential promoter use exists between samples. Only genes producing two or more distinct primary transcripts (i.e. multi-promoter genes) are listed here. | |
280 11. Differential CDS tests: this tab delimited file lists, for each gene, the amount of overloading detected among its coding sequences, i.e. how much differential CDS output exists between samples. Only genes producing two or more distinct CDS (i.e. multi-protein genes) are listed here. | |
281 | |
282 ------- | |
283 | |
284 **Settings** | |
285 | |
286 All of the options have a default value. You can change any of them. Most of the options in Cuffdiff have been implemented here. | |
287 | |
288 ------ | |
289 | |
290 **Cuffdiff parameter list** | |
291 | |
292 This is a list of implemented Cuffdiff options:: | |
293 | |
294 -m INT Average fragement length; default 200 | |
295 -s INT Fragment legnth standard deviation; default 80 | |
296 -c INT The minimum number of alignments in a locus for needed to conduct significance testing on changes in that locus observed between samples. If no testing is performed, changes in the locus are deemed not significant, and the locus' observed changes don't contribute to correction for multiple testing. The default is 1,000 fragment alignments (up to 2,000 paired reads). | |
297 --FDR FLOAT The allowed false discovery rate. The default is 0.05. | |
298 --num-importance-samples INT Sets the number of importance samples generated for each locus during abundance estimation. Default: 1000 | |
299 --max-mle-iterations INT Sets the number of iterations allowed during maximum likelihood estimation of abundances. Default: 5000 | |
300 -N With this option, Cufflinks excludes the contribution of the top 25 percent most highly expressed genes from the number of mapped fragments used in the FPKM denominator. This can improve robustness of differential expression calls for less abundant genes and transcripts. | |
301 | |
302 </help> | |
303 </tool> |