|
58
|
1 <tool id="iwtomics_loadandplot" name="IWTomics Load" version="@VERSION@.0">
|
|
18
|
2 <description>Smooth and Plot</description>
|
|
|
3 <macros>
|
|
|
4 <import>macros.xml</import>
|
|
|
5 </macros>
|
|
76
|
6 <!--<expand macro="requirements" />-->
|
|
18
|
7
|
|
74
|
8 <command detect_errors="exit_code">
|
|
18
|
9 <![CDATA[
|
|
70
|
10 Rscript '$__tool_directory__/loadandplot.R' outrdata='${outrdata}' outregions='${outregions}' outfeatures='${outfeatures}' outpdf='${outpdf}'
|
|
18
|
11
|
|
72
|
12 #set region_paths = '|'.join( [ str( $r ) for $r in $regions ] )
|
|
70
|
13 regionspaths='${region_paths}'
|
|
72
|
14 #set region_galaxyids = ','.join( [ str( $r.hid ) for $r in $regions ] )
|
|
70
|
15 regionsgalaxyids='c(${region_galaxyids})'
|
|
72
|
16 #set region_names = '|'.join( [ str( $r.name ) for $r in $regions ] )
|
|
70
|
17 regionsfilenames='${region_names}'
|
|
72
|
18 #set region_header = $regionsheadersection.regionsheader
|
|
18
|
19 #if $region_header != "":
|
|
70
|
20 regionsheaderfile='${region_header}'
|
|
18
|
21 #end if
|
|
70
|
22
|
|
72
|
23 #set feature_paths = '|'.join( [ str( $f ) for $f in $features ] )
|
|
70
|
24 featurespaths='${feature_paths}'
|
|
72
|
25 #set feature_galaxyids = ','.join( [ str( $f.hid ) for $f in $features ] )
|
|
70
|
26 featuresgalaxyids='c(${feature_galaxyids})'
|
|
72
|
27 #set feature_names = '|'.join( [ str( $f.name ) for $f in $features ] )
|
|
70
|
28 featuresfilenames='${feature_names}'
|
|
72
|
29 #set feature_header = $featuresheadersection.featuresheader
|
|
18
|
30 #if $feature_header != "":
|
|
70
|
31 featuresheaderfile='${feature_header}'
|
|
18
|
32 #end if
|
|
|
33
|
|
72
|
34 smoothing="'${conditionaltype.smoothing}'"
|
|
70
|
35 start.are.0based='${zerobased}'
|
|
18
|
36 #if $conditionaltype.smoothing == 'no':
|
|
72
|
37 alignment="'${conditionaltype.alignment}'"
|
|
18
|
38 #elif $conditionaltype.smoothing == 'locpoly':
|
|
70
|
39 bandwidth='${conditionaltype.locpolybandwidth}'
|
|
|
40 degree='${conditionaltype.locpolydegree}'
|
|
76
|
41 alignment="'${conditionaltype.conditionalscale.alignment}'"
|
|
18
|
42 #elif $conditionaltype.smoothing == 'kernel':
|
|
70
|
43 bandwidth='${conditionaltype.kernelbandwidth}'
|
|
76
|
44 alignment="'${conditionaltype.conditionalscale.alignment}'"
|
|
18
|
45 #elif $conditionaltype.smoothing == 'splines':
|
|
70
|
46 degree='${conditionaltype.splinesdegree}'
|
|
|
47 dist_knots='${conditionaltype.splinesdistknots}'
|
|
76
|
48 alignment="'${conditionaltype.conditionalscale.alignment}'"
|
|
|
49 #end if
|
|
|
50
|
|
|
51 #if $conditionaltype.smoothing != 'no':
|
|
|
52 #if $conditionaltype.conditionalscale.alignment == 'scale':
|
|
|
53 scale='${conditionaltype.conditionalscale.scalegrid}'
|
|
|
54 fill_gaps='TRUE'
|
|
|
55 #else
|
|
|
56 fill_gaps='${conditionaltype.conditionalscale.fillgaps}'
|
|
|
57 #end if
|
|
18
|
58 #end if
|
|
|
59
|
|
70
|
60 average='${plotres.average}'
|
|
|
61 size='${plotres.size}'
|
|
72
|
62 plottype="'${plotres.conditionalplottype.plottype}'"
|
|
18
|
63 #if $plotres.conditionalplottype.plottype == 'boxplot':
|
|
72
|
64 #set probs = ','.join( [ str( $p.prob ) for $p in $plotres.conditionalplottype.probabilitiessection.probabilities ] )
|
|
18
|
65 #if $probs != "":
|
|
70
|
66 probs='c(${plotres.conditionalplottype.probabilitiessection.prob0},${plotres.conditionalplottype.probabilitiessection.prob1},${plotres.conditionalplottype.probabilitiessection.prob2},${probs})'
|
|
18
|
67 #else:
|
|
70
|
68 probs='c(${plotres.conditionalplottype.probabilitiessection.prob0},${plotres.conditionalplottype.probabilitiessection.prob1},${plotres.conditionalplottype.probabilitiessection.prob2})'
|
|
18
|
69 #end if
|
|
74
|
70 #end if
|
|
18
|
71 ]]>
|
|
|
72 </command>
|
|
|
73
|
|
|
74 <inputs>
|
|
|
75 <!-- regions -->
|
|
|
76 <param format="tabular" name="regions" type="data" label="Select region files" multiple="true" help="" />
|
|
70
|
77 <section name="regionsheadersection" title="Insert header file for regions (optional)" expanded="false" help="Tabular file. Column 1: names of the region files; Column 2: unique identifiers of the region datasets (spaces and special characters are not allowed); Column 3: name of the region datasets to be used in the output.">
|
|
|
78 <param format="tabular" name="regionsheader" type="data" label="Select header file" optional="true" />
|
|
|
79 </section>
|
|
18
|
80
|
|
|
81 <!-- features -->
|
|
|
82 <param format="tabular" name="features" type="data" label="Select feature files" multiple="true" help="" />
|
|
70
|
83 <section name="featuresheadersection" title="Insert header file for features (optional)" expanded="false" help="Tabular file. Column 1: names of the feature files; Column 2: unique identifiers of the features (spaces and special characters are not allowed); Column 3: name of the features to be used in the output.">
|
|
|
84 <param format="tabular" name="featuresheader" type="data" label="Select header file" optional="true" />
|
|
|
85 </section>
|
|
18
|
86
|
|
|
87 <!-- zero-base one-base -->
|
|
|
88 <param name="zerobased" type="select" label="Start positions are:" multiple="false" help="">
|
|
|
89 <option value="TRUE">0-based</option>
|
|
|
90 <option value="FALSE">1-based</option>
|
|
|
91 </param>
|
|
|
92
|
|
|
93 <!-- conditional smoothing -->
|
|
|
94 <conditional name="conditionaltype">
|
|
|
95 <!-- smoothing -->
|
|
|
96 <param name="smoothing" type="select" label="Smoothing" help="Type of smoothing to be applied to the feature curves.">
|
|
|
97 <option value="no">No</option>
|
|
|
98 <option value="locpoly">Local Polynomials</option>
|
|
|
99 <option value="kernel">Gaussian Kernel</option>
|
|
|
100 <option value="splines">B-Splines</option>
|
|
|
101 </param>
|
|
|
102 <!-- conditional choice: smoothing=no -->
|
|
|
103 <when value="no">
|
|
76
|
104 <expand macro="alignment" />
|
|
18
|
105 </when>
|
|
|
106 <!-- conditional choice: smoothing=locpoly -->
|
|
|
107 <when value="locpoly">
|
|
|
108 <param name="locpolybandwidth" type="integer" value="5" min="2" label="Bandwidth for smoothing" help="" />
|
|
|
109 <param name="locpolydegree" type="integer" value="3" min="0" label="Degree of local polynomial" help="" />
|
|
76
|
110 <!-- conditional region alignment -->
|
|
|
111 <expand macro="alignment-s" />
|
|
18
|
112 </when>
|
|
|
113 <!-- conditional choice: smoothing=kernel -->
|
|
|
114 <when value="kernel">
|
|
76
|
115 <param name="kernelbandwidth" type="integer" value="5" min="2" label="Bandwidth for smoothing" help="" />
|
|
18
|
116 <!-- conditional region alignment -->
|
|
76
|
117 <expand macro="alignment-s" />
|
|
18
|
118 </when>
|
|
|
119 <!-- conditional choice: smoothing=splines -->
|
|
|
120 <when value="splines">
|
|
|
121 <param name="splinesdegree" type="integer" value="3" min="0" label="Degree of B-splines" help="" />
|
|
|
122 <param name="splinesdistknots" type="integer" value="10" min="2" label="Distance between nodes (approximate)" help="" />
|
|
76
|
123 <!-- conditional region alignment -->
|
|
|
124 <expand macro="alignment-s" />
|
|
18
|
125 </when>
|
|
|
126 </conditional>
|
|
|
127
|
|
|
128 <!-- plot IWTomics results -->
|
|
|
129 <section name="plotres" title="Plot data" expanded="True">
|
|
76
|
130 <expand macro="plot-params" />
|
|
18
|
131 </section>
|
|
|
132 </inputs>
|
|
|
133
|
|
|
134 <outputs>
|
|
|
135 <data format="rdata" name="outrdata" label="${tool.name} on ${on_string}: IWTomicsData Object" from_work_dir="iwtomics.loadandplot.RData" />
|
|
|
136 <data format="tabular" name="outregions" label="${tool.name} on ${on_string}: Region Dataset IDs" from_work_dir="iwtomics.loadandplot.regions.txt" />
|
|
|
137 <data format="tabular" name="outfeatures" label="${tool.name} on ${on_string}: Feature IDs" from_work_dir="iwtomics.loadandplot.features.txt" />
|
|
|
138 <data format="pdf" name="outpdf" label="${tool.name} on ${on_string}: Plotted Data" from_work_dir="iwtomics.loadandplot.pdf" />
|
|
|
139 </outputs>
|
|
|
140
|
|
56
|
141 <tests>
|
|
|
142 <test>
|
|
|
143 <param name="regions" value="input/Controls_regions.bed,input/Elements1_regions.bed,input/Elements2_regions.bed,input/Elements3_regions.bed" ftype="bed" />
|
|
|
144 <param name="regionsheader" value="input/regions.header.txt" ftype="tabular" />
|
|
|
145 <param name="features" value="input/Feature1.bed,input/Feature2.bed" ftype="bed" />
|
|
|
146 <param name="featuresheader" value="input/features.header.bed.txt" ftype="tabular" />
|
|
|
147 <param name="zerobased" value="TRUE" />
|
|
|
148 <param name="smoothing" value="kernel" />
|
|
|
149 <param name="alignment" value="center" />
|
|
|
150 <param name="kernelbandwidth" value="5" />
|
|
|
151 <param name="fillgaps" value="TRUE" />
|
|
|
152 <param name="average" value="TRUE" />
|
|
|
153 <param name="size" value="TRUE" />
|
|
|
154 <param name="plottype" value="boxplot" />
|
|
|
155 <param name="prob0" value="0.25" />
|
|
|
156 <param name="prob1" value="0.5" />
|
|
|
157 <param name="prob2" value="0.75" />
|
|
|
158 <output name="outrdata" file="output_loadandplot/iwtomics.loadandplot.RData" compare="sim_size" />
|
|
|
159 <output name="outregions" file="output_loadandplot/iwtomics.loadandplot.regions.txt" />
|
|
|
160 <output name="outfeatures" file="output_loadandplot/iwtomics.loadandplot.features.txt" />
|
|
|
161 <output name="outpdf" file="output_loadandplot/iwtomics.loadandplot.pdf" compare="sim_size" />
|
|
|
162 </test>
|
|
|
163 </tests>
|
|
|
164
|
|
68
|
165 <help><![CDATA[
|
|
18
|
166 This tool imports a collection of genomic region datasets, and associates to each region
|
|
|
167 multiple genomic feature measurements. It allows to align the regions in multiple ways
|
|
|
168 (center, left, right or scale alignment), to smooth the feature curves (possibly filling gaps in the
|
|
|
169 measurements) and to create a graphical representation of the feature measurements in each
|
|
|
170 region datasets (aligned curves or pointwise quantile curves).
|
|
|
171
|
|
|
172 -----
|
|
|
173
|
|
|
174 **Region datasets**
|
|
|
175
|
|
|
176 Each region dataset can be provided as a BED or Tabular file with tab delimited columns chr
|
|
|
177 start end (extra columns present in the input file are ignored). Regions can be of different
|
|
|
178 length::
|
|
|
179
|
|
|
180 chr2 49960150 50060150
|
|
|
181 chr2 55912445 56012445
|
|
|
182 ...
|
|
|
183
|
|
|
184 -----
|
|
|
185
|
|
|
186 **Feature measurements**
|
|
|
187
|
|
|
188 Feature measurements corresponding to all the regions can be provided as a BED or Tabular
|
|
|
189 file with tab delimited columns chr start end value::
|
|
|
190
|
|
|
191 chr2 49960150 49962150 0.9426
|
|
|
192 chr2 49962150 49964150 0.7816
|
|
|
193 ...
|
|
|
194
|
|
|
195 Each feature must be measured in windows of a fixed size inside all the regions (missing
|
|
|
196 values must be indicated as NA). Another way to import feature measurements is from a
|
|
|
197 Tabular file with the first three columns chr start end corresponding to the different genomic
|
|
|
198 regions, followed on the same row by all the measurements in fixed-size windows::
|
|
|
199
|
|
|
200 chr2 49960150 50060150 0.9426 0.7816 0.8921 ... ... 1.2063
|
|
|
201 chr2 55912445 56012445 0.8719 0.9975 1.1619 ... ... 0.9601
|
|
|
202 ...
|
|
|
203
|
|
|
204 -----
|
|
|
205
|
|
|
206 **Output**
|
|
|
207
|
|
|
208 The tool returns:
|
|
|
209
|
|
|
210 1. RData with the IWTomicsData object, that stores the aligned genomic region datasets, and their associated feature measurements;
|
|
|
211 2. Region dataset identifiers;
|
|
|
212 3. Feature identifiers;
|
|
|
213 4. PDF file with the plotted data.
|
|
|
214
|
|
|
215 1-3 can be used as input of the tool *IWTomics Test and Plot*
|
|
|
216
|
|
|
217 -----
|
|
|
218
|
|
|
219 .. class:: infomark
|
|
|
220
|
|
|
221 **Notes**
|
|
|
222
|
|
54
|
223 This Galaxy tool has been developed by Fabio Cumbo (Third University of Rome, Italy - fabio.cumbo@iasi.cnr.it) and Marzia A. Cremona (The Pennsylvania State University, USA - mac78@psu.edu).
|
|
18
|
224
|
|
|
225 It implements a simplified version of the methods *smooth* and *plot* for *IWTomicsData* objects.
|
|
|
226 The complete version can be found in the *R/Bioconductor* package *IWTomics* (see vignette_).
|
|
|
227
|
|
|
228 .. _vignette: https://bioconductor.org/packages/release/bioc/vignettes/IWTomics/inst/doc/IWTomics.pdf
|
|
54
|
229
|
|
|
230 Example data can be found at:
|
|
|
231
|
|
|
232 1. Simulated_data_
|
|
|
233 2. ETn_data_
|
|
|
234
|
|
|
235 .. _Simulated_data: https://usegalaxy.org/u/fabio-cumbo/h/iwtomics-example
|
|
|
236 .. _ETn_data: https://usegalaxy.org/u/fabio-cumbo/h/iwtomics-etn-example
|
|
68
|
237 ]]></help>
|
|
18
|
238
|
|
70
|
239 <expand macro="citations" />
|
|
18
|
240
|
|
|
241 </tool>
|
