Mercurial > repos > rnateam > mea
comparison mea.xml @ 8:e0f8bc01b84b draft
planemo upload commit 9099c0a9ff5f7f3d9407e2ec3598957e02850d88-dirty
author | rnateam |
---|---|
date | Thu, 10 Dec 2015 12:22:08 -0500 |
parents | 5f780b6d99d3 |
children | ee79a49755da |
comparison
equal
deleted
inserted
replaced
7:5f780b6d99d3 | 8:e0f8bc01b84b |
---|---|
1 <tool id="mea" name="MEA" version="0.6.4.0"> | 1 <tool id="mea" name="MEA" version="0.6.4.1"> |
2 <description>Predict MEA structures and compare structures of RNAs</description> | 2 <description>Predict MEA structures and compare structures of RNAs</description> |
3 | |
3 <requirements> | 4 <requirements> |
4 <requirement type="package" version="2.22.6">gengetopt</requirement> | 5 <requirement type="package" version="2.22.6" >gengetopt</requirement> |
5 <requirement type="package" version="0.6.4">mea</requirement> | 6 <requirement type="package" version="0.6.4" >mea</requirement> |
6 </requirements> | 7 </requirements> |
8 | |
7 <stdio> | 9 <stdio> |
8 <exit_code range=":-1"/> | 10 <exit_code range=":-1"/> |
9 <exit_code range="1:" /> | 11 <exit_code range="1:" /> |
10 </stdio> | 12 </stdio> |
13 | |
11 <version_command>mea --version</version_command> | 14 <version_command>mea --version</version_command> |
15 | |
12 <command><![CDATA[ | 16 <command><![CDATA[ |
13 mea | 17 mea |
14 #if str($predict_mode.predict_selector) == "predict": | 18 #if str($predict_mode.predict_selector) == "predict" |
15 $predict_mode.dotplot | 19 $predict_mode.dotplot |
16 #if str($predict_mode.mea_parameters.mea_parameters_selector) == "no_default" | 20 #if str($predict_mode.mea_parameters.mea_parameters_selector) == "no_default" |
17 --alpha $predict_mode.mea_parameters.alpha | 21 --alpha $predict_mode.mea_parameters.alpha |
18 --beta $predict_mode.mea_parameters.beta | 22 --beta $predict_mode.mea_parameters.beta |
19 --gamma $predict_mode.mea_parameters.gamma | 23 --gamma $predict_mode.mea_parameters.gamma |
20 --delta $predict_mode.mea_parameters.delta | 24 --delta $predict_mode.mea_parameters.delta |
21 #end if | 25 #end if |
22 #else | 26 #else |
23 --structure `cat $predict_mode.structure` | 27 --structure `cat $predict_mode.structure` |
24 #end if | 28 #end if |
25 #if str($compare_mode.compare_selector) == 'compare': | 29 #if $compare_mode.compare_selector |
26 --reference `cat $compare_mode.reference` | 30 --reference `cat $compare_mode.reference` |
27 $compare_mode.slide_rule | 31 $compare_mode.slide_rule |
28 $compare_mode.conflict_rule | 32 $compare_mode.conflict_rule |
29 #end if | 33 #end if |
30 > $stdout | 34 > $stdout |
31 ]]></command> | 35 ]]></command> |
36 | |
32 <inputs> | 37 <inputs> |
33 <conditional name="predict_mode"> | 38 <conditional name="predict_mode"> |
34 <param name="predict_selector" type="select" label="Predict MEA structure"> | 39 <param name="predict_selector" type="select" label="Predict MEA structure"> |
35 <option value="predict">Predict MEA structure</option> | 40 <option value="predict">Predict MEA structure</option> |
36 <option value="compare">Specify structure (for comparison)</option> | 41 <option value="compare">Specify structure (for comparison)</option> |
37 </param> | 42 </param> |
38 <when value="predict"> | 43 <when value="predict"> |
39 <param name="dotplot" type="data" format="rna_eps" label="Dotplot" optional="False" | 44 <param name="dotplot" type="data" format="rna_eps" label="Dotplot" |
40 help="Dotplot file (RNA base pair probabilities)"/> | 45 optional="false" help="Dotplot file (RNA base pair probabilities)"/> |
41 <conditional name="mea_parameters"> | 46 <conditional name="mea_parameters"> |
42 <param name="mea_parameters_selector" type="boolean" label="Use default parameters" | 47 <param name="mea_parameters_selector" |
43 truevalue="default" falsevalue="no_default" checked="True" | 48 type="boolean" |
44 help="The default parameters assign a base pair weight | 49 label="Use default parameters" |
45 of 0.5 and penalize long base pairs." /> | 50 truevalue="default" |
46 <when value="default" /> | 51 falsevalue="no_default" |
47 <when value="no_default"> | 52 checked="yes" |
48 <param name="alpha" label="Alpha value" type="float" | 53 help="The default parameters assign a base pair weight |
49 optional="False" value="0.012" | 54 of 0.5 and penalize long base pairs." |
50 help="Slope of base pair distance penalty. (--alpha)"/> | 55 /> |
51 <param name="beta" label="Beta value" type="float" | 56 <when value="default" /> |
52 optional="False" value="315" | 57 <when value="no_default"> |
53 help="Turning point of base pair distance penalty. (--beta)" /> | 58 <param name="alpha" label="Alpha" type="float" |
54 <param name="gamma" label="Gamma value" type="float" | 59 optional="false" value="0.012" |
55 optional="False" value="0.5" | 60 help="Slope of base pair distance penalty"/> |
56 help="Base pair weight factor. (--gamma)" /> | 61 <param name="beta" label="Beta" type="float" |
57 <param name="delta" label="Delta value" type="float" | 62 optional="false" value="315" |
58 optional="False" value="0.003" | 63 help="Turning point of base pair distance penalty" /> |
59 help="Minimum penalty factor for base pairs. (--delta)" /> | 64 <param name="gamma" label="Gamma" type="float" |
60 </when> | 65 optional="false" value="0.5" |
61 </conditional> | 66 help="Base pair weight factor" /> |
62 </when> | 67 <param name="delta" label="Delta" type="float" |
63 <when value="compare"> | 68 optional="false" value="0.003" |
64 <param name="structure" format="txt" type="data" label="Structure" optional="False" | 69 help="Minimum penalty factor for base pairs" /> |
65 help="(Predicted) RNA secondary structure for comparison to a reference structure. (--structure)" /> | 70 </when> |
66 </when> | 71 </conditional> |
67 </conditional> | 72 </when> |
68 <conditional name="compare_mode"> | 73 <when value="compare"> |
69 <param name="compare_selector" type="select" label="Compare to a reference structure" | 74 <param name="structure" format="txt" type="data" label="Structure" |
70 help="Whether to compare the predicted (or specified) structure to a reference structure."> | 75 optional="false" |
71 <option value="dont_compare">Do not compare</option> | 76 help="(Predicted) RNA secondary structure |
72 <option value="compare">Compare to reference structure</option> | 77 for comparison to a reference structure." /> |
73 </param> | 78 </when> |
74 <when value="dont_compare" /> | 79 </conditional> |
75 <when value="compare"> | 80 |
76 <param name="reference" | 81 |
77 optional="false" | 82 <conditional name="compare_mode"> |
78 format="txt" type="data" | 83 <param name="compare_selector" |
79 label="Reference" | 84 type="boolean" |
80 help="Reference structure in dot-bracket format"/> | 85 label="Compare to a reference structure" |
81 <param name="slide_rule" label="Slide Rule" type="boolean" | 86 optional="false" |
82 optional="False" checked="True" falsevalue="--no-slide-rule" truevalue="" help="Use slide rule. (--no-slide-rule)" /> | 87 checked="false" |
83 <param name="conflict_rule" label="Conflict Rule" type="boolean" | 88 truevalue="compare" |
84 optional="False" checked="True" falsevalue="--no-conflict-rule" truevalue="" help="Use onflict rule. (--no-conflict-rule)"/> | 89 falsevalue="dont_compare" |
85 </when> | 90 help="Whether to compare the predicted (or specified) structure |
86 </conditional> | 91 to a reference structure." |
92 /> | |
93 <when value="dont_compare" /> | |
94 <when value="compare"> | |
95 <param name="reference" | |
96 optional="false" | |
97 format="txt" type="data" | |
98 label="Reference" | |
99 help="Reference structure in dot-bracket format"/> | |
100 | |
101 <param name="slide_rule" label="Slide Rule" type="boolean" | |
102 optional="false" | |
103 checked="yes" | |
104 falsevalue="--no-slide-rule" truevalue="" | |
105 help="Use slide rule" /> | |
106 <param name="conflict_rule" label="Conflict Rule" type="boolean" | |
107 optional="false" | |
108 checked="yes" | |
109 falsevalue="--no-conflict-rule" truevalue="" | |
110 help="Use onflict rule"/> | |
111 </when> | |
112 </conditional> | |
87 </inputs> | 113 </inputs> |
114 | |
88 <outputs> | 115 <outputs> |
89 <data format="txt" name="stdout" label="${tool.name} on ${on_string}" /> | 116 <data format="txt" name="stdout" label="${tool.name} on ${on_string}" /> |
90 </outputs> | 117 </outputs> |
118 | |
91 <tests> | 119 <tests> |
92 <test> | 120 <test> |
93 <param name="mea_parameter_selector" value="default" /> | 121 <param name="mea_parameter_selector" value="default" /> |
94 <param name="predict_selector" value="predict" /> | 122 |
95 <param name="compare_selector" value="compare" /> | 123 <param name="predict_selector" value="predict" /> |
96 <param name="dotplot" value="test_dp.ps" /> | 124 <param name="compare_selector" value="true" /> |
97 <param name="reference" value="test_reference.txt" /> | 125 |
98 <output name="stdout" file="test_predict.out" /> | 126 <param name="dotplot" value="test_dp.ps" /> |
99 </test> | 127 <param name="reference" value="test_reference.txt" /> |
100 <test> | 128 |
101 <param name="mea_parameter_selector" value="default" /> | 129 <output name="stdout" file="test_predict.out" /> |
102 <param name="predict_selector" value="compare" /> | 130 </test> |
103 <param name="compare_selector" value="compare" /> | 131 |
104 <param name="structure" value="test_structure.txt" /> | 132 <test> |
105 <param name="reference" value="test_reference.txt" /> | 133 <param name="mea_parameter_selector" value="default" /> |
106 <output name="stdout" file="test_compare.out" /> | 134 <param name="predict_selector" value="compare" /> |
107 </test> | 135 <param name="compare_selector" value="true" /> |
136 | |
137 <param name="structure" value="test_structure.txt" /> | |
138 <param name="reference" value="test_reference.txt" /> | |
139 | |
140 <output name="stdout" file="test_compare.out" /> | |
141 </test> | |
108 </tests> | 142 </tests> |
143 | |
144 | |
109 <help><![CDATA[ | 145 <help><![CDATA[ |
110 ===== | 146 ===== |
111 MEA | 147 MEA |
112 ===== | 148 ===== |
113 | 149 |
114 `MEA <http://www.bioinf.uni-leipzig.de/Software/mea/>`_ predicts RNA maximum expected accuracy structures from RNA base pair probabilities and optionally compares them to a reference structure. In a special mode it skips the prediction and compares a | 150 MEA predicts RNA maximum expected accuracy structures from RNA base |
151 pair probabilities and optionally compares them to a reference | |
152 structure. In a special mode it skips the prediction and compares a | |
115 given structure to the reference. For the prediction, MEA allows to | 153 given structure to the reference. For the prediction, MEA allows to |
116 penalize long base pairs, using parameters alpha, beta, gamma, and | 154 penalize long base pairs, using parameters alpha, beta, gamma, and |
117 delta. For the comparison of secondary structures, several measures | 155 delta. For the comparison of secondary structures, several measures |
118 are computed from the confusion matrix of the RNA base pairs. | 156 are computed from the confusion matrix of the RNA base pairs. |
119 | 157 |
148 | 186 |
149 * FN = # false negatives | 187 * FN = # false negatives |
150 | 188 |
151 * TN = # true negatives | 189 * TN = # true negatives |
152 | 190 |
153 * SENS = TP/(TP+FN) 'Sensitivity'http://www.bioinf.uni-leipzig.de/Software/mea/ | 191 * SENS = TP/(TP+FN) 'Sensitivity' |
154 | 192 |
155 * PPV = TP/(TP+FP) 'Positive Predictive Value' | 193 * PPV = TP/(TP+FP) 'Positive Predictive Value' |
156 | 194 |
157 * F1 = PPV*SENS / (PPV+SENS), if PPV+SENS!=0; 0, otherwise 'F1-score' | 195 * F1 = PPV*SENS / (PPV+SENS), if PPV+SENS!=0; 0, otherwise 'F1-score' |
158 | 196 |
161 | 199 |
162 Special rules for prediction evaluation: | 200 Special rules for prediction evaluation: |
163 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ | 201 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
164 | 202 |
165 * Slide rule: tolerate shift of one base pair end by one base. This | 203 * Slide rule: tolerate shift of one base pair end by one base. This |
166 rule directly affects the number of true positives. | 204 rule directly affects the number of true positives. |
205 | |
167 * Conflict rule: predicted base pairs are false only if they | 206 * Conflict rule: predicted base pairs are false only if they |
168 conflict with the reference; two base pair conflict if and only if | 207 conflict with the reference; two base pair conflict if and only if |
169 they share one end This rule directly affects the number of false | 208 they share one end This rule directly affects the number of false |
170 positives. | 209 positives. |
171 | 210 |
211 -------- | |
212 Download | |
213 -------- | |
214 | |
215 The command line tool MEA is free software available for download and | |
216 local installation at | |
217 .. __: http://www.bioinf.uni-leipzig.de/Software/mea/ | |
172 ]]></help> | 218 ]]></help> |
173 <citations> | 219 <citations> |
174 <citation type="doi">10.1007/978-3-319-02624-4_1</citation> | 220 <citation type="doi">10.1007/978-3-319-02624-4_1</citation> |
175 </citations> | 221 </citations> |
222 | |
176 </tool> | 223 </tool> |