|
0
|
1 <tool id="wtavg" name="Assign weighted-average" version="1.0.0">
|
|
|
2 <description> of the values of features overlapping an interval </description>
|
|
|
3 <command interpreter="python">WeightedAverage.py $genomic_interval $genomic_feature $out_file1 -1 ${genomic_interval.metadata.chromCol},${genomic_interval.metadata.startCol},${genomic_interval.metadata.endCol},${genomic_interval.metadata.strandCol} -2 ${genomic_feature.metadata.chromCol},${genomic_feature.metadata.startCol},${genomic_feature.metadata.endCol},${genomic_feature.metadata.strandCol},${genomic_feature.metadata.nameCol}</command>
|
|
|
4
|
|
|
5 <inputs>
|
|
|
6 <param format="interval" name="genomic_interval" type="data" label="Genomic intervals (first dataset)" help="Dataset missing? See Note below."/>
|
|
|
7 <param format="interval" name="genomic_feature" label="Genomic features (second dataset)" type="data" help="Make sure the value column is specified. See Note below." />
|
|
|
8 </inputs>
|
|
|
9 <outputs>
|
|
|
10 <data format="input" name="out_file1" metadata_source="genomic_interval" />
|
|
|
11 </outputs>
|
|
|
12 <tests>
|
|
|
13 <!-- Test data with valid values -->
|
|
|
14 <test>
|
|
|
15 <param name="genomic_interval" value="interval_interpolate.bed"/>
|
|
|
16 <param name="genomic_feature" value="value_interpolate.bed"/>
|
|
|
17 <output name="out_file1" file="interpolate_result.bed"/>
|
|
|
18 </test>
|
|
|
19
|
|
|
20 </tests>
|
|
|
21 <help>
|
|
|
22
|
|
|
23
|
|
|
24 .. class:: infomark
|
|
|
25
|
|
|
26 **What it does**
|
|
|
27
|
|
|
28 For each interval in your first dataset, this tool calculates the weighted average value of the overlapping features in your second dataset.
|
|
|
29
|
|
|
30 - When a genomic interval partially or totally overlaps a single genomic feature, the value of that genomic feature is assigned to the genomic interval.
|
|
|
31 - When a genomic interval partially or totally overlaps with more than one genomic features, the average of the values of the overlapping genomic features weighted by the corresponding number of overlapping bases is assigned to the genomic interval.
|
|
|
32 - When a genomic interval does not overlap with any genomic feature, 'NA' will be assigned as it's value.
|
|
|
33
|
|
|
34 -----
|
|
|
35
|
|
|
36 .. class:: warningmark
|
|
|
37
|
|
|
38 **Note**
|
|
|
39
|
|
|
40 The input datasets should be in **bed** or **interval** format. Please use "edit attributes"/pencil icon to specify the column containing the values for the features in the second dataset as **name/identifier** column.
|
|
|
41
|
|
|
42 The output will contain all the columns in the first input plus a new column containing the assigned value for each interval.
|
|
|
43
|
|
|
44 -----
|
|
|
45
|
|
|
46 **Example**
|
|
|
47
|
|
|
48 - Suppose our first dataset contains the following **genomic intervals**::
|
|
|
49
|
|
|
50 chr start stop
|
|
|
51 chr1 1000 2000
|
|
|
52 chr1 3000 5000
|
|
|
53 chr1 8000 9000
|
|
|
54
|
|
|
55 - and our second dataset contains the following **genomic features** each having an associated value (in fourth column) ::
|
|
|
56
|
|
|
57 chr start stop name
|
|
|
58 chr1 900 1200 0.5
|
|
|
59 chr1 2900 3100 0.2
|
|
|
60 chr1 4800 5100 0.8
|
|
|
61
|
|
|
62 - For each **genomic interval** in our first dataset, this tool calculates the weighted average value of the overlapping **genomic features** in our second dataset ::
|
|
|
63
|
|
|
64 chr1 1000 2000 0.5
|
|
|
65 chr1 3000 5000 0.6
|
|
|
66 chr1 8000 9000 NA
|
|
|
67
|
|
|
68
|
|
|
69
|
|
|
70 </help>
|
|
|
71 </tool> |
