|
48
|
1 <tool id="deeptools_heatmapper" name="heatmapper" version="@WRAPPER_VERSION@.0">
|
|
|
2 <description>creates a heatmap for a score associated to genomic regions</description>
|
|
|
3 <expand macro="requirements"/>
|
|
|
4 <expand macro="stdio" />
|
|
|
5 <macros>
|
|
|
6 <token name="@BINARY@">heatmapper</token>
|
|
|
7 <import>deepTools_macros.xml</import>
|
|
|
8 </macros>
|
|
|
9 <command>
|
|
52
|
10 <![CDATA[
|
|
48
|
11 heatmapper
|
|
|
12
|
|
|
13 --matrixFile $matrixFile
|
|
|
14 --outFileName $outFileName
|
|
|
15
|
|
|
16 #if $output.showOutputSettings == "yes"
|
|
|
17 --plotFileFormat $output.outFileFormat
|
|
|
18 #if $outFileNameData:
|
|
|
19 --outFileNameData '$outFileNameData'
|
|
|
20 #end if
|
|
|
21
|
|
|
22 #if $outFileNameMatrix:
|
|
|
23 --outFileNameMatrix '$outFileNameMatrix'
|
|
|
24 #end if
|
|
|
25
|
|
|
26 #if $outFileSortedRegions:
|
|
|
27 --outFileSortedRegions '$outFileSortedRegions'
|
|
|
28 #end if
|
|
|
29 #else
|
|
|
30 --plotFileFormat 'png'
|
|
|
31 #end if
|
|
|
32
|
|
|
33 #if $advancedOpt.showAdvancedOpt == "yes"
|
|
|
34 #if $advancedOpt.sortRegions:
|
|
|
35 --sortRegions '$advancedOpt.sortRegions'
|
|
|
36 #end if
|
|
|
37
|
|
|
38 #if $advancedOpt.sortUsing:
|
|
|
39 --sortUsing '$advancedOpt.sortUsing'
|
|
|
40 #end if
|
|
|
41
|
|
|
42 #if $advancedOpt.averageTypeSummaryPlot:
|
|
|
43 --averageTypeSummaryPlot '$advancedOpt.averageTypeSummaryPlot'
|
|
|
44 #end if
|
|
|
45
|
|
|
46 #if str($advancedOpt.missingDataColor.value) != "None":
|
|
|
47 --missingDataColor '$advancedOpt.missingDataColor'
|
|
|
48 #end if
|
|
|
49
|
|
|
50 --colorMap '$advancedOpt.colorMap'
|
|
|
51
|
|
|
52 #if str($advancedOpt.zMin).strip() != "":
|
|
|
53 --zMin $advancedOpt.zMin
|
|
|
54 #end if
|
|
|
55 #if $advancedOpt.zMax:
|
|
|
56 --zMax $advancedOpt.zMax
|
|
|
57 #end if
|
|
|
58
|
|
|
59 #if str($advancedOpt.yMin).strip() != "":
|
|
|
60 --yMin $advancedOpt.yMin
|
|
|
61 #end if
|
|
|
62 #if $advancedOpt.yMax:
|
|
|
63 --yMax $advancedOpt.yMax
|
|
|
64 #end if
|
|
|
65
|
|
|
66 --xAxisLabel '$advancedOpt.xAxisLabel'
|
|
|
67 --yAxisLabel '$advancedOpt.yAxisLabel'
|
|
|
68
|
|
|
69 --heatmapWidth $advancedOpt.heatmapWidth
|
|
|
70 --heatmapHeight $advancedOpt.heatmapHeight
|
|
|
71
|
|
|
72 --whatToShow '$advancedOpt.whatToShow'
|
|
|
73
|
|
|
74 --startLabel '$advancedOpt.startLabel'
|
|
|
75 --endLabel '$advancedOpt.endLabel'
|
|
|
76 --refPointLabel '$advancedOpt.referencePointLabel'
|
|
|
77 --regionsLabel '$advancedOpt.regionsLabel'
|
|
|
78
|
|
|
79 #if str($advancedOpt.plotTitle.value) != "None":
|
|
|
80 --plotTitle '$advancedOpt.plotTitle'
|
|
|
81 #end if
|
|
|
82
|
|
|
83 $advancedOpt.onePlotPerGroup
|
|
|
84
|
|
52
|
85 @KMEANS_CLUSTERING@
|
|
48
|
86
|
|
|
87 #end if
|
|
52
|
88 ]]>
|
|
48
|
89 </command>
|
|
|
90 <inputs>
|
|
|
91 <param name="matrixFile" format="bgzip" type="data" label="Matrix file from the computeMatrix tool"/>
|
|
|
92
|
|
|
93 <expand macro="input_graphic_output_settings">
|
|
|
94 <expand macro="input_image_file_format" />
|
|
|
95 <expand macro="input_save_matrix_values" />
|
|
|
96 </expand>
|
|
|
97
|
|
|
98 <conditional name="advancedOpt" >
|
|
|
99 <param name="showAdvancedOpt" type="select" label="Show advanced options" >
|
|
|
100 <option value="no" selected="true">no</option>
|
|
|
101 <option value="yes">yes</option>
|
|
|
102 </param>
|
|
|
103 <when value="no" />
|
|
|
104 <when value="yes">
|
|
|
105 <param name="sortRegions" type="select" label="Sort regions"
|
|
|
106 help="Whether the heatmap should present the regions sorted. The default is to sort in descending order based on the mean value per region.">
|
|
|
107 <option value="no">no ordering</option>
|
|
|
108 <option value="descend" selected="true">descending order</option>
|
|
|
109 <option value="ascend">ascending order</option>
|
|
|
110 </param>
|
|
|
111
|
|
|
112 <param name="sortUsing" type="select" label="Method used for sorting" help="For each row the method is computed." >
|
|
|
113 <option value="mean" selected="true">mean</option>
|
|
|
114 <option value="median">median</option>
|
|
|
115 <option value="min">min</option>
|
|
|
116 <option value="max">max</option>
|
|
|
117 <option value="sum">sum</option>
|
|
|
118 <option value="region_length">region length</option>
|
|
|
119 </param>
|
|
|
120
|
|
|
121 <param name="averageTypeSummaryPlot" type="select" label="Type of statistic that should be plotted in the summary image above the heatmap">
|
|
|
122 <option value="mean" selected="true">mean</option>
|
|
|
123 <option value="median">median</option>
|
|
|
124 <option value="min">min</option>
|
|
|
125 <option value="max">max</option>
|
|
|
126 <option value="sum">sum</option>
|
|
|
127 <option value="std">std</option>
|
|
|
128 </param>
|
|
|
129
|
|
52
|
130 <param name="missingDataColor" type="text" value="black" optional="true" label="Missing data color"
|
|
|
131 help="If 'Represent missing data as zero' is not set, such cases will be colored in black by default. By using this parameter a different color can be set. A value between 0 and 1 will be used for a gray scale (black is 0). Also color names can be used, see a list here: http://packages.python.org/ete2/reference/reference_svgcolors.html. Alternatively colors can be specified using the #rrggbb notation." />
|
|
48
|
132
|
|
|
133 <expand macro="colormap" />
|
|
|
134
|
|
52
|
135 <param name="zMin" type="float" value="" size="3"
|
|
|
136 label="Minimum value for the heatmap intensities. Leave empty for automatic values"/>
|
|
|
137 <param name="zMax" type="float" value="" size="3"
|
|
|
138 label="Maximum value for the heatmap intensities. Leave empty for automatic values"/>
|
|
|
139 <param name="yMin" type="float" value="" size="3"
|
|
|
140 label="Minimum value for the Y-axis of the summary plot. Leave empty for automatic values"/>
|
|
|
141 <param name="yMax" type="float" value="" size="3"
|
|
|
142 label="Maximum value for Y-axis of the summary plot. Leave empty for automatic values"/>
|
|
|
143 <param name="xAxisLabel" type="text" value="distance from TSS (bp)" size="200"
|
|
|
144 label="Description for the x-axis label" />
|
|
|
145 <param name="yAxisLabel" type="text" value="genes" size="30"
|
|
|
146 label="Description for the y-axis label for the top panel" />
|
|
48
|
147
|
|
52
|
148 <param name="heatmapWidth" type="float" value="7.5" min="1" max="100"
|
|
|
149 label="Heatmap width in cm" help="The minimum value is 1 and the maximum is 100."/>
|
|
|
150 <param name="heatmapHeight" type="float" value="25" min="3" max="100"
|
|
|
151 label="Heatmap height in cm" help="The minimum value is 3 and the maximum is 100."/>
|
|
48
|
152
|
|
52
|
153 <param name="whatToShow" type="select" label="What to show"
|
|
|
154 help ="The default is to include a summary or profile plot on top of the heatmap and a heatmap colorbar.">
|
|
48
|
155 <option value="plot, heatmap and colorbar" selected="true">summary plot, heatmap and colorbar</option>
|
|
|
156 <option value="plot and heatmap">summary plot and heatmap (no colorbar)</option>
|
|
|
157 <option value="heatmap only">heatmap only</option>
|
|
|
158 <option value="heatmap and colorbar">heatmap and colorbar</option>
|
|
|
159 <option value="colorbar only">colorbar only</option>
|
|
|
160 </param>
|
|
|
161
|
|
52
|
162 <param name="startLabel" type="text" value="TSS" size="10"
|
|
|
163 label="Label for the region start"
|
|
|
164 help ="[only for scale-regions mode] Label shown in the plot for the start of the region. Default is TSS (transcription start site), but could be changed to anything, e.g. "peak start"." />
|
|
|
165 <param name="endLabel" type="text" value="TES" size="10"
|
|
|
166 label="Label for the region end"
|
|
|
167 help="[only for scale-regions mode] Label shown in the plot for the region end. Default is TES (transcription end site)."/>
|
|
48
|
168
|
|
52
|
169 <param name="referencePointLabel" type="text" value="TSS" size="10"
|
|
|
170 label="Reference point label"
|
|
|
171 help ="[only for scale-regions mode] Label shown in the plot for the reference-point. Default is the same as the reference point selected (e.g. TSS), but could be anything, e.g. "peak start" etc." />
|
|
48
|
172 <param name="regionsLabel" type="text" value="genes" size="30"
|
|
|
173 label="Labels for the regions plotted in the heatmap"
|
|
|
174 help="If more than one region is being plotted a list of labels separated by comma and limited by quotes, is required. For example, label1, label2.">
|
|
|
175 <sanitizer>
|
|
|
176 <valid initial="string.printable">
|
|
|
177 </valid>
|
|
|
178 </sanitizer>
|
|
|
179 </param>
|
|
52
|
180 <param name="plotTitle" type="text" value="" size="30"
|
|
|
181 label="Title of the plot" help="Title of the plot, to be printed on top of the generated image. Leave blank for no title. (--plotTitle)" />
|
|
48
|
182 <param name="onePlotPerGroup" type="boolean" truevalue="--onePlotPerGroup" falsevalue=""
|
|
|
183 label="Do one plot per group"
|
|
|
184 help="When computeMatrix was used on more than one group of genes, the average plots for all the groups will be drawn in one panel by default. If this option is set, each group will get its own plot, stacked on top of each other."/>
|
|
|
185
|
|
|
186 <expand macro="kmeans_clustering" />
|
|
|
187 </when>
|
|
|
188 </conditional>
|
|
|
189 </inputs>
|
|
|
190 <outputs>
|
|
|
191 <expand macro="output_image_file_format" />
|
|
|
192 <expand macro="output_graphic_outputs" />
|
|
|
193 <expand macro="output_save_matrix_values" />
|
|
|
194 </outputs>
|
|
|
195 <tests>
|
|
|
196 <test>
|
|
52
|
197 <param name="matrixFile" value="computeMatrix_result1.gz" ftype="bgzip" />
|
|
|
198 <output name="outFileName" file="heatmapper_result1.png" ftype="png" compare="sim_size" delta="100" />
|
|
48
|
199 </test>
|
|
|
200 </tests>
|
|
|
201 <help>
|
|
52
|
202 <![CDATA[
|
|
48
|
203 **What it does**
|
|
|
204
|
|
|
205 The heatmapper visualizes scores associated with genomic regions, for example ChIP enrichment values around the TSS of genes.
|
|
|
206 Like profiler, it requires that computeMatrix was run first to calculate the values.
|
|
|
207
|
|
|
208 We implemented vast optional parameters to optimize the visual output and we encourage you to play around with the min/max values displayed in the heatmap as well as
|
|
|
209 with the different coloring options. The most powerful option is the k-means clustering where you simply need to indicate the number of
|
|
|
210 groups with similar read distributions that you expect and the algorithm will do the sorting for you.
|
|
|
211
|
|
|
212 Do check the examples on our help page with step-by-step protocols: https://github.com/fidelram/deepTools/wiki/Example-workflows
|
|
|
213
|
|
|
214
|
|
|
215 .. image:: $PATH_TO_IMAGES/visual_hm_DmelPolII.png
|
|
|
216 :alt: Heatmap of RNA Polymerase II ChIP-seq
|
|
|
217
|
|
|
218
|
|
|
219 You can find more details on the tool itself on the heatmapper wiki page: https://github.com/fidelram/deepTools/wiki/Visualizations#wiki-heatmapper
|
|
|
220
|
|
|
221
|
|
|
222 -----
|
|
|
223
|
|
|
224 @REFERENCES@
|
|
52
|
225 ]]>
|
|
48
|
226 </help>
|
|
|
227 <expand macro="citations" />
|
|
|
228 </tool>
|
