comparison ngs.plot.xml @ 0:b848d0eb3841 draft default tip

Uploaded
author kosrou
date Mon, 02 Mar 2015 08:16:45 -0500
parents
children
comparison
equal deleted inserted replaced
-1:000000000000 0:b848d0eb3841
1 <tool id="ngs_plot_1" name="ngs.plot">
2 <command interpreter="bash">
3 #if $gene_list_source_type.gene_list == "plot_whole_genome":
4 #if $genomic_region_source_type.genomic_region == "tss":
5 ngs.plot.tss_tes.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $genomic_region_source_type.flanking_region_size1 $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area
6 #else if $genomic_region_source_type.genomic_region == "tes":
7 ngs.plot.tss_tes.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $genomic_region_source_type.flanking_region_size2 $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area
8 #else if $genomic_region_source_type.genomic_region == "genebody":
9 #if $genomic_region_source_type.flanking_region_option_source_type1.flanking_region_option3 == "flanking_region_size":
10 ngs.plot.genebody_exon_cgi.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type1.flanking_region_option3 $genomic_region_source_type.flanking_region_option_source_type1.flanking_region_size3 $genomic_region_source_type.further_information1 $genomic_region_source_type.interval_size1
11 #else if $genomic_region_source_type.flanking_region_option_source_type1.flanking_region_option3 == "flanking_floating_size":
12 ngs.plot.genebody_exon_cgi.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type1.flanking_region_option3 $genomic_region_source_type.flanking_region_option_source_type1.specify_flanking_region3 $genomic_region_source_type.further_information1 $genomic_region_source_type.interval_size1
13 #end if
14 #else if $genomic_region_source_type.genomic_region == "exon":
15 #if $genomic_region_source_type.flanking_region_option_source_type2.flanking_region_option4 == "flanking_region_size":
16 ngs.plot.genebody_exon_cgi.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type2.flanking_region_option4 $genomic_region_source_type.flanking_region_option_source_type2.flanking_region_size4 $genomic_region_source_type.further_information2 $genomic_region_source_type.interval_size2
17 #else if $genomic_region_source_type.flanking_region_option_source_type2.flanking_region_option4 == "flanking_floating_size":
18 ngs.plot.genebody_exon_cgi.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type2.flanking_region_option4 $genomic_region_source_type.flanking_region_option_source_type2.specify_flanking_region4 $genomic_region_source_type.further_information2 $genomic_region_source_type.interval_size2
19 #end if
20 #else if $genomic_region_source_type.genomic_region == "cgi":
21 #if $genomic_region_source_type.flanking_region_option_source_type3.flanking_region_option5 == "flanking_region_size":
22 ngs.plot.genebody_exon_cgi.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type3.flanking_region_option5 $genomic_region_source_type.flanking_region_option_source_type3.flanking_region_size5 $genomic_region_source_type.further_information3 $genomic_region_source_type.interval_size3
23 #else if $genomic_region_source_type.flanking_region_option_source_type3.flanking_region_option5 == "flanking_floating_size":
24 ngs.plot.genebody_exon_cgi.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type3.flanking_region_option5 $genomic_region_source_type.flanking_region_option_source_type3.specify_flanking_region5 $genomic_region_source_type.further_information3 $genomic_region_source_type.interval_size3
25 #end if
26 #else if $genomic_region_source_type.genomic_region == "bed":
27 #if $genomic_region_source_type.flanking_region_option_source_type4.flanking_region_option6 == "flanking_region_size":
28 ngs.plot.bed.sh $genome_name $bed_file $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type4.flanking_region_option6 $genomic_region_source_type.flanking_region_option_source_type4.flanking_region_size6 $genomic_region_source_type.interval_size4
29 #else if $genomic_region_source_type.flanking_region_option_source_type4.flanking_region_option6 == "flanking_floating_size":
30 ngs.plot.bed.sh $genome_name $bed_file $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type4.flanking_region_option6 $genomic_region_source_type.flanking_region_option_source_type4.specify_flanking_region6 $genomic_region_source_type.interval_size4
31 #end if
32 #end if
33 #else if $gene_list_source_type.gene_list == "offer_gene_list":
34 #if $genomic_region_source_type.genomic_region == "tss":
35 ngs.plot.tss_tes_genelist.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $genomic_region_source_type.flanking_region_size1 $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $gene_list_source_type.gene_list_detail
36 #else if $genomic_region_source_type.genomic_region == "tes":
37 ngs.plot.tss_tes_genelist.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $genomic_region_source_type.flanking_region_size2 $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $gene_list_source_type.gene_list_detail
38 #else if $genomic_region_source_type.genomic_region == "genebody":
39 #if $genomic_region_source_type.flanking_region_option_source_type1.flanking_region_option3 == "flanking_region_size":
40 ngs.plot.genebody_exon_cgi_genelist.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type1.flanking_region_option3 $genomic_region_source_type.flanking_region_option_source_type1.flanking_region_size3 $genomic_region_source_type.further_information1 $genomic_region_source_type.interval_size1 $gene_list_source_type.gene_list_detail
41 #else if $genomic_region_source_type.flanking_region_option_source_type1.flanking_region_option3 == "flanking_floating_size":
42 ngs.plot.genebody_exon_cgi_genelist.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type1.flanking_region_option3 $genomic_region_source_type.flanking_region_option_source_type1.specify_flanking_region3 $genomic_region_source_type.further_information1 $genomic_region_source_type.interval_size1 $gene_list_source_type.gene_list_detail
43 #end if
44 #else if $genomic_region_source_type.genomic_region == "exon":
45 #if $genomic_region_source_type.flanking_region_option_source_type2.flanking_region_option4 == "flanking_region_size":
46 ngs.plot.genebody_exon_cgi_genelist.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type2.flanking_region_option4 $genomic_region_source_type.flanking_region_option_source_type2.flanking_region_size4 $genomic_region_source_type.further_information2 $genomic_region_source_type.interval_size2 $gene_list_source_type.gene_list_detail
47 #else if $genomic_region_source_type.flanking_region_option_source_type2.flanking_region_option4 == "flanking_floating_size":
48 ngs.plot.genebody_exon_cgi_genelist.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type2.flanking_region_option4 $genomic_region_source_type.flanking_region_option_source_type2.specify_flanking_region4 $genomic_region_source_type.further_information2 $genomic_region_source_type.interval_size2 $gene_list_source_type.gene_list_detail
49 #end if
50 #else if $genomic_region_source_type.genomic_region == "cgi":
51 #if $genomic_region_source_type.flanking_region_option_source_type3.flanking_region_option5 == "flanking_region_size":
52 ngs.plot.genebody_exon_cgi_genelist.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type3.flanking_region_option5 $genomic_region_source_type.flanking_region_option_source_type3.flanking_region_size5 $genomic_region_source_type.further_information3 $genomic_region_source_type.interval_size3 $gene_list_source_type.gene_list_detail
53 #else if $genomic_region_source_type.flanking_region_option_source_type3.flanking_region_option5 == "flanking_floating_size":
54 ngs.plot.genebody_exon_cgi_genelist.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type3.flanking_region_option5 $genomic_region_source_type.flanking_region_option_source_type3.specify_flanking_region5 $genomic_region_source_type.further_information3 $genomic_region_source_type.interval_size3 $gene_list_source_type.gene_list_detail
55 #end if
56 #else if $genomic_region_source_type.genomic_region == "bed":
57 #if $genomic_region_source_type.flanking_region_option_source_type4.flanking_region_option6 == "flanking_region_size":
58 ngs.plot.bed_genelist.sh $genome_name $bed_file $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type4.flanking_region_option6 $genomic_region_source_type.flanking_region_option_source_type4.flanking_region_size6 $genomic_region_source_type.interval_size4 $gene_list_source_type.gene_list_detail
59 #else if $genomic_region_source_type.flanking_region_option_source_type4.flanking_region_option6 == "flanking_floating_size":
60 ngs.plot.bed_genelist.sh $genome_name $bed_file $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type4.flanking_region_option6 $genomic_region_source_type.flanking_region_option_source_type4.specify_flanking_region6 $genomic_region_source_type.interval_size4 $gene_list_source_type.gene_list_detail
61 #end if
62 #end if
63 #end if
64 </command>
65 <inputs>
66 <param type="data" name="bam_file" label="BAM File for Plotting"/>
67 <param type="select" name="genome_name" label="Choose Genome">
68 <option value="hg19">hg19</option>
69 <option value="mm9">mm9</option>
70 <option value="rn4">rn4</option>
71 </param>
72 <conditional name="genomic_region_source_type">
73 <param type="select" name="genomic_region" label="Genomic Region interested in">
74 <option value="tss">tss</option>
75 <option value="tes">tes</option>
76 <option value="genebody">genebody</option>
77 <option value="exon">exon</option>
78 <option value="cgi">cgi</option>
79 <option value="bed">customized bed file</option>
80 </param>
81 <when value="tss">
82 <!-- <param type="select" name="flanking_region_option1" label="Choose one method of below two options for plotting flanking region. If choosing one option, You need give this option's value, and ignore the other option and value">
83 <option value="flanking_region_size">Using flanking region size</option>
84 <option value="flanking_floating_size">Using floating size with interval size</option>
85 </param>
86 -->
87 <param name="flanking_region_size1" size="30" type="text" value="2000" label="Flanking Region Size"/>
88 <!-- <param name="specify_flanking_region1" size="30" type="text" value="0.33" label="2) Specify flanking region as a multiple(or fraction) of interval size, can be less than 1. In this case, flanking region has floating size with interval size. default: 33%"/>
89 -->
90 </when>
91 <when value="tes">
92 <!-- <param type="select" name="flanking_region_option2" label="Choose one method of below two options for plotting flanking region. If choosing one option, You need give this option's value, and ignore the other option and value">
93 <option value="flanking_region_size">Using flanking region size</option>
94 <option value="flanking_floating_size">Using floating size with interval size</option>
95 </param>
96 -->
97 <param name="flanking_region_size2" size="30" type="text" value="2000" label="Flanking Region Size"/>
98 <!-- <param name="specify_flanking_region2" size="30" type="text" value="0.33" label="2) Specify flanking region as a multiple(or fraction) of interval size, can be less than 1. In this case, flanking region has floating size with interval size. default: 33%"/>
99 -->
100 </when>
101 <when value="genebody">
102 <param type="select" name="further_information1" label="If you select genebody, further information can be provided for specific regions to plot: ChIP-seq(default), RNA-seq">
103 <option value="chipseq">ChIP-seq</option>
104 <option value="rnaseq">RNA-seq</option>
105 </param>
106 <param name="interval_size1" type="select" label="Internal Region Size">
107 <option value="automatic">Automatic</option>
108 <option value="0">0</option>
109 <option value="1">1</option>
110 </param>
111
112 <conditional name="flanking_region_option_source_type1">
113 <param type="select" name="flanking_region_option3" label="Choose one method below for plotting flanking region">
114 <option value="flanking_region_size">Using flanking region size</option>
115 <option value="flanking_floating_size">Using floating size with interval size</option>
116 </param>
117 <when value="flanking_region_size">
118 <param name="flanking_region_size3" size="30" type="text" value="2000" label="Flanking Region Size"/>
119 </when>
120 <when value="flanking_floating_size">
121 <param name="specify_flanking_region3" size="30" type="text" value="0.33" label="Flanking Region Factor"/>
122 </when>
123 </conditional>
124 </when>
125 <when value="exon">
126 <param type="select" name="further_information2" label="If you select exon, further information can be provided for specific regions to plot: canonical(default), variant, promoter, polyA, altAcceptor, altDonor, altBoth">
127 <option value="canonical">canonical</option>
128 <option value="variant">variant</option>
129 <option value="promoter">promoter</option>
130 <option value="polyA">polyA</option>
131 <option value="altAcceptor">altAcceptor</option>
132 <option value="altDonor">altDonor</option>
133 <option value="altBoth">altBoth</option>
134 </param>
135 <param name="interval_size2" type="select" label="Internal Region Size">
136 <option value="automatic">Automatic</option>
137 <option value="0">0</option>
138 <option value="1">1</option>
139 </param>
140
141 <conditional name="flanking_region_option_source_type2">
142 <param type="select" name="flanking_region_option4" label="Choose one method below for plotting flanking region">
143 <option value="flanking_region_size">Using flanking region size</option>
144 <option value="flanking_floating_size">Using floating size with interval size</option>
145 </param>
146 <when value="flanking_region_size">
147 <param name="flanking_region_size4" size="30" type="text" value="500" label="Flanking Region Size"/>
148 </when>
149 <when value="flanking_floating_size">
150 <param name="specify_flanking_region4" size="30" type="text" value="0.33" label="Flanking Region Factor"/>
151 </when>
152 </conditional>
153 </when>
154 <when value="cgi">
155 <param type="select" name="further_information3" label="If you select cgi, further information can be provided for specific regions to plot: proximalpromoter(default), genebody, genedesert, otherintergenic, pericentromere, promoter1k, promoter3k">
156 <option value="ProximalPromoter">proximalpromoter</option>
157 <option value="Genebody">genebody</option>
158 <option value="Genedesert">genedesert</option>
159 <option value="OtherIntergenic">otherintergenic</option>
160 <option value="Pericentromere">pericentromere</option>
161 <option value="Promoter1k">promoter1k</option>
162 <option value="Promoter3k">promoter3k</option>
163 </param>
164 <param name="interval_size3" type="select" label="Internal Region Size">
165 <option value="automatic">Automatic</option>
166 <option value="0">0</option>
167 <option value="1">1</option>
168 </param>
169
170 <conditional name="flanking_region_option_source_type3">
171 <param type="select" name="flanking_region_option5" label="Choose one method below for plotting flanking region">
172 <option value="flanking_region_size">Using flanking region size</option>
173 <option value="flanking_floating_size">Using floating size with interval size</option>
174 </param>
175 <when value="flanking_region_size">
176 <param name="flanking_region_size5" size="30" type="text" value="500" label="Flanking Region Size"/>
177 </when>
178 <when value="flanking_floating_size">
179 <param name="specify_flanking_region5" size="30" type="text" value="0.33" label="Flanking Region Factor"/>
180 </when>
181 </conditional>
182 </when>
183
184
185 <when value="bed">
186 <param type="data" name="bed_file" label="customized *.bed file for plotting"/>
187 <param name="interval_size4" type="select" label="Internal Region Size">
188 <option value="automatic">Automatic</option>
189 <option value="0">0</option>
190 <option value="1">1</option>
191 </param>
192
193
194 <conditional name="flanking_region_option_source_type4">
195 <param type="select" name="flanking_region_option6" label="Choose one method below for plotting flanking region">
196 <option value="flanking_region_size">Using flanking region size</option>
197 <option value="flanking_floating_size">Using floating size with interval size</option>
198 </param>
199 <when value="flanking_region_size">
200 <param name="flanking_region_size6" size="30" type="text" value="1000" label="Flanking Region Size"/>
201 </when>
202 <when value="flanking_floating_size">
203 <param name="specify_flanking_region6" size="30" type="text" value="0.33" label="Flanking Region Factor"/>
204 </when>
205 </conditional>
206 </when>
207
208 </conditional>
209 <param type="select" name="gene_database" label="Gene Database">
210 <option value="refseq">refseq</option>
211 <option value="ensembl">ensembl</option>
212 </param>
213 <param name="image_title" size="30" type="text" value="Noname" label="Image title"/>
214
215 <conditional name="gene_list_source_type">
216 <param type="select" name="gene_list" label="Choose scale for ploting">
217 <option value="plot_whole_genome">Plot whole genome</option>
218 <option value="offer_gene_list">Offer gene list</option>
219 </param>
220 <when value="plot_whole_genome">
221 </when>
222 <when value="offer_gene_list">
223 <param name="gene_list_detail" type="data" label="Choose the uploaded gene list for plotting"/>
224 </when>
225 </conditional>
226
227 <!-- <param name="interval_size" size="30" type="text" value="3000" label="Interval size(optional, default varies by type: genebody=3000;exon=250;cgi=500)"/>
228 <param name="flanking_region_size" size="30" type="text" value="1000" label="Flanking region size(optional, default varies by type: tss,tes,genebody=1000;exon=500;cgi=500)"/>
229 <param name="specify_flanking_region" size="30" type="text" value="specify_flanking_region" label="Specify flanking region as a multiple(or fraction) of interval size, can be less than 1. In this case, flanking region has floating size with interval size."/> -->
230 <param name="randomly_sample" size="10" type="text" value="1" label="Randomly sample the regions for plotting"/>
231 <param type="select" name="gene_order" label="Gene order algorithm">
232 <option value="total">total</option>
233 <option value="hc">hc</option>
234 <option value="max">max</option>
235 <option value="prod">prod</option>
236 <option value="diff">diff</option>
237 <option value="pca">pca</option>
238 <option value="none">none</option>
239 </param>
240 <param name="chunk_size" size="10" type="text" value="100" label="Chunk size for loading genes in batch"/>
241 <param name="fragement_length" size="10" type="text" value="150" label="Expected Fragement Length"/>
242 <param name="quality_requirement" size="10" type="text" value="20" label="Mapping Quality Requirement"/>
243 <param type="select" name="standard_error" label="Boolean tag for plotting standard errors">
244 <option value="1">On</option>
245 <option value="0">Off</option>
246 </param>
247
248 <param name="radius_size" size="10" type="text" value="0" label="The fraction of extreme values to be trimmed on both ends"/>
249 <param name="flooding_fraction" size="10" type="text" value="0.02" label="Flooding fraction"/>
250 <param type="select" name="smooth_method" label="Moving window width to smooth avg. profiles">
251 <option value="1">no</option>
252 <option value="3">slightly</option>
253 <option value="5">somewhat</option>
254 <option value="9">quite</option>
255 <option value="13">super</option>
256 </param>
257 <param name="shaded_area" size="10" type="text" value="0" label="Opacity of shaded area, suggested value:[0, 0.5]"/>
258 <!-- <param name="weighted_coverage" size="30" type="text" value="weighted_coverage" label="By default, do NOT calculate weighted coverage for splined curves using gene length. However, this can be useful, e.g. when comparing enrichment of a histone mark under two conditions. Can be turned on by setting -E to 1."/>
259 -->
260 </inputs>
261 <outputs>
262 <data format="pdf" name="out_avg_png" metadata_source="bam_file" />
263 <data format="pdf" name="out_hm_png" />
264 <data format="zip" name="out_zip"/>
265
266 </outputs>
267 <help>
268 **NGS PLOT Overview**
269
270 ngs.plot_ needs an indexed bam file or a configuration file as an input to plot the peak trends and patterns cross the genomic region or users' interested region. ngs.plot.r will generate multiple files including avgprof file, heatmap file, and a zip file for replotting.
271
272 .. _ngs.plot: http://code.google.com/p/ngsplot/
273
274 ------
275
276 **Input formats**
277
278 ngs.plot needs an indexed bam file or a configuration file as an input to plot the peak trends and patterns cross the genomic region or users' interested region.
279
280 ------
281
282 **Outputs**
283
284 ngs.plot generates multiple files including avgprof file, heatmap file, and a zip file for replotting.
285
286 ------
287
288 **ngs.plot parameter list**
289
290 This is a list of implemented ngs.plot options::
291
292 -G Genome name, currently supported: mm9, hg19, rn4, sacCer3(ensembl only)
293 -R Genomic regions to plot: tss, tes, genebody, exon, cgi or *.bed
294 -C Indexed bam file or a configuration file for multiplot
295 -E Gene list to subset regions
296 -T Image title(default=Noname), will be used in figure legend.
297 -F If you select genebody, exon or cgi, further information can be provided: for genebody: ChIP-seq(default), RNA-seq. RNA-seq is handled differently from ChIP-seq to cope with intronic regions. for exon: canonical(default), variant, promoter, polyA, altAcceptor, altDonor, altBoth. the alt's mean exons with alternative boundaries. for cgi, choose one based on location: ProximalPromoter(default), Promoter1k, Promoter3k, Genebody, Genedesert, OtherIntergenic, Pericentromere
298 -D Gene database: refseq(default), ensembl
299 -I Shall interval be larger than flanking in plot?(0 or 1)
300 -L Flanking region size (Size in bps. By default, when -R=tss, tes, genebody, -L=2000; when -R=exon, cgi, -L=500; when -R=*.bed, -L=1000)
301 -N Flanking region factor(will override flanking size) -S Randomly sample the regions for plot, must be:(0, 1]
302 -S Randomly sample the regions for plot, must be:(0, 1] (default=1, i.e. all regions)
303 -GO Gene order algorithm used in heatmaps: total(default), hc, max, prod, diff, pca and none(according to gene list supplied)
304 -CS Chunk size for loading genes in batch(default=100)
305 -FL Fragment length used to calculate physical coverage(default=150)
306 -MQ Mapping quality cutoff to filter reads(default=20)
307 -SE Shall standard errors be plotted?(0 or 1)
308 -RB The fraction of extreme values to be trimmed on both ends (default=0, 0.05 means 5% of extreme values will be trimmed)
309 -FC Flooding fraction:[0, 1), default=0.02
310 -MW Moving window width to smooth avg. profiles, must be integer:1=no(default); 3=slightly; 5=somewhat; 9=quite; 13=super.
311 -H Opacity of shaded area, suggested value:[0, 0.5] (default=0, i.e. no shading, just curves)
312 </help>
313 </tool>
314