comparison iqtree.xml @ 3:a396ccf6079c draft default tip

"planemo upload for repository https://github.com/quadram-institute-bioscience/galaxy-tools/tree/master/tools/iqtree"
author leomrtns
date Wed, 18 Dec 2019 13:46:22 +0000
parents 87daf702e477
children
comparison
equal deleted inserted replaced
2:87daf702e477 3:a396ccf6079c
14 #if $general_options.s 14 #if $general_options.s
15 -s '$general_options.s' 15 -s '$general_options.s'
16 #end if 16 #end if
17 17
18 ## file 18 ## file
19 #if $general_options.t 19 #if str($initial_tree.mode.tree) == "file"
20 -t '$general_options.t' 20 -t '$initial_tree.mode.t'
21 $tree_parameters.constructing_consensus.con 21 $tree_parameters.constructing_consensus.con
22 $tree_parameters.constructing_consensus.net 22 $tree_parameters.constructing_consensus.net
23 #if str($tree_parameters.constructing_consensus.bi) != '' 23 #if str($tree_parameters.constructing_consensus.bi) != ''
24 -bi '$tree_parameters.constructing_consensus.bi' 24 -bi '$tree_parameters.constructing_consensus.bi'
25 #end if 25 #end if
37 #if $tree_parameters.computing_robinson_foulds.rf 37 #if $tree_parameters.computing_robinson_foulds.rf
38 -rf '$tree_parameters.computing_robinson_foulds.rf' 38 -rf '$tree_parameters.computing_robinson_foulds.rf'
39 $tree_parameters.computing_robinson_foulds.rf_all 39 $tree_parameters.computing_robinson_foulds.rf_all
40 $tree_parameters.computing_robinson_foulds.rf_adj 40 $tree_parameters.computing_robinson_foulds.rf_adj
41 #end if 41 #end if
42 #end if 42 #else if str($initial_tree.mode.tree) == "fixed"
43 43 -te '$initial_tree.mode.te'
44 ## file 44 #else if str($initial_tree.mode.tree) == "default"
45 #if $general_options.te 45 #if str($initial_tree.mode.ninit) != ''
46 -te '$general_options.te' 46 -ninit '$initial_tree.mode.ninit'
47 #end if
48 #if str($initial_tree.mode.sprrad) != ''
49 -sprrad '$initial_tree.mode.sprrad'
50 #end if
51 #else if str($initial_tree.mode.tree) == "fast"
52 -fast -nbest 1
53 #else if str($initial_tree.mode.tree) == "pars"
54 -t PARS -ninit 2 -nbest 1
55 #else
56 -t BIONJ -ninit 2 -nbest 1
57 #end if
58
59 #if str($initial_tree.mode.tree) in ["fixed", "default", "file"]:
60 #if str($tree_parameters.tree_search.nbest) != ''
61 -nbest '$tree_parameters.tree_search.nbest'
62 #end if
47 #end if 63 #end if
48 64
49 #if str($general_options.st) != '' 65 #if str($general_options.st) != ''
50 -st '$general_options.st' 66 -st '$general_options.st'
51 #end if 67 #end if
152 #end if 168 #end if
153 169
154 $modelling_parameters.site_specific_frequency.fmax 170 $modelling_parameters.site_specific_frequency.fmax
155 171
156 172
157
158 #if str($tree_parameters.tree_search.ninit) != ''
159 -ninit '$tree_parameters.tree_search.ninit'
160 #end if
161
162 #if str($tree_parameters.tree_search.ntop) != '' 173 #if str($tree_parameters.tree_search.ntop) != ''
163 -ntop '$tree_parameters.tree_search.ntop' 174 -ntop '$tree_parameters.tree_search.ntop'
164 #end if 175 #end if
165 176
166 #if str($tree_parameters.tree_search.nbest) != '' 177
167 -nbest '$tree_parameters.tree_search.nbest'
168 #end if
169 178
170 #if str($tree_parameters.tree_search.nstop) != '' 179 #if str($tree_parameters.tree_search.nstop) != ''
171 -nstop '$tree_parameters.tree_search.nstop' 180 -nstop '$tree_parameters.tree_search.nstop'
172 #end if 181 #end if
173 182
174 #if str($tree_parameters.tree_search.n) != '' 183 #if str($tree_parameters.tree_search.n) != ''
175 -n '$tree_parameters.tree_search.n' 184 -n '$tree_parameters.tree_search.n'
176 #end if 185 #end if
177 186
178 #if str($tree_parameters.tree_search.sprrad) != ''
179 -sprrad '$tree_parameters.tree_search.sprrad'
180 #end if
181
182 #if str($tree_parameters.tree_search.pers) != '' 187 #if str($tree_parameters.tree_search.pers) != ''
183 -pers '$tree_parameters.tree_search.pers' 188 -pers '$tree_parameters.tree_search.pers'
184 #end if 189 #end if
185 190
186 $tree_parameters.tree_search.allnni 191 $tree_parameters.tree_search.allnni
187 $tree_parameters.tree_search.djc 192 $tree_parameters.tree_search.djc
188 $tree_parameters.tree_search.fast
189 193
190 ## file 194 ## file
191 #if $tree_parameters.tree_search.g 195 #if $tree_parameters.tree_search.g
192 -g '$tree_parameters.tree_search.g' 196 -g '$tree_parameters.tree_search.g'
193 #end if 197 #end if
276 </command> 280 </command>
277 <inputs> 281 <inputs>
278 <section name="general_options" expanded="True" title="General options"> 282 <section name="general_options" expanded="True" title="General options">
279 <param argument="-s" type="data" format="txt" optional="true" label="Specify input alignment file in PHYLIP, FASTA, NEXUS, CLUSTAL or MSF format."/> 283 <param argument="-s" type="data" format="txt" optional="true" label="Specify input alignment file in PHYLIP, FASTA, NEXUS, CLUSTAL or MSF format."/>
280 <param argument="-st" type="select" label="Specify sequence type as either of DNA, AA, BIN, MORPH, CODON or NT2AA for DNA, amino-acid, binary, morphological, codon or DNA-to-AA-translated sequences"> 284 <param argument="-st" type="select" label="Specify sequence type as either of DNA, AA, BIN, MORPH, CODON or NT2AA for DNA, amino-acid, binary, morphological, codon or DNA-to-AA-translated sequences">
281 <help><![CDATA[
282 Note that -st CODON is always necessary when using codon models and you also need to specify a genetic code like this if differed from the standard genetic code.
283 <br/><i>-st</i> NT2AA tells IQ-TREE to translate protein-coding DNA into AA sequences and then subsequent analysis will work on the AA sequences. You can also use a genetic code like -st NT2AA5 for the Invertebrate Mitochondrial Code (see genetic code table).]]>
284 </help>
285 <option value="DNA">DNA</option> 285 <option value="DNA">DNA</option>
286 <option value="AA">AA</option> 286 <option value="AA">AA</option>
287 <option value="BIN">BIN</option> 287 <option value="BIN">BIN</option>
288 <option value="MORPH">MORPH</option> 288 <option value="MORPH">MORPH</option>
289 <option value="CODON">CODON</option> 289 <option value="CODON">CODON</option>
290 <option value="NT2AA">NT2AA</option> 290 <option value="NT2AA">NT2AA</option>
291 <help><![CDATA[
292 Note that -st CODON is always necessary when using codon models and you also need to specify a genetic code like this if differed from the standard genetic code.
293 <br/><i>-st</i> NT2AA tells IQ-TREE to translate protein-coding DNA into AA sequences and then subsequent analysis will work on the AA sequences.
294 You can also use a genetic code like -st NT2AA5 for the Invertebrate Mitochondrial Code (see genetic code table).
295 ]]>
296 </help>
291 </param> 297 </param>
292 <param argument="-t" type="data" format="nhx" optional="true" label="Specify a file containing starting tree for tree search"/>
293 <param argument="-te" type="data" format="nhx" optional="true" label="Like -t but fixing user tree" help="That means, no tree search is performed and IQ-TREE computes the log-likelihood of the fixed user tree."/>
294 <param argument="-seed" type="integer" optional="true" label="Specify a random number seed to reproduce a previous run (leave blank to randomize)"/> 298 <param argument="-seed" type="integer" optional="true" label="Specify a random number seed to reproduce a previous run (leave blank to randomize)"/>
295 <param argument="-keep_ident" type="boolean" truevalue="-keep-ident" falsevalue="" checked="false" label="Keep identical sequences in the alignment" help="By default: IQ-TREE will remove them during the analysis and add them in the end."/> 299 <param argument="-keep_ident" type="boolean" truevalue="-keep-ident" falsevalue="" checked="false" label="Keep identical sequences in the alignment" help="By default: IQ-TREE will remove them during the analysis and add them in the end."/>
296 <param argument="-safe" type="boolean" truevalue="-safe" falsevalue="" checked="false" label="Turn on safe numerical mode to avoid numerical underflow for large data sets with many sequences (typically in the order of thousands)" help="This mode is automatically turned on when having more than 2000 sequences."/> 300 <param argument="-safe" type="boolean" truevalue="-safe" falsevalue="" checked="false" label="Turn on safe numerical mode to avoid numerical underflow for large data sets with many sequences (typically in the order of thousands)" help="This mode is automatically turned on when having more than 2000 sequences."/>
297 </section> 301 </section>
302
303 <section name="initial_tree" expanded="True" title="Input and initial trees">
304 <conditional name="mode">
305 <param name="tree" type="select" label="Initial tree search">
306 <option value="default" selected="true">parsimony using PLL library (default) </option>
307 <option value="fast">fasttree emulation (option "-fast")</option>
308 <option value="pars">fast parsimony, using IQTREE internal algorithm (option "-t PARS")</option>
309 <option value="bionj">bioNJ tree (option "-t BIONJ")</option>
310 <option value="fixed">fixed tree, no optimisation (option "-te"). Used in e.g. ancestral reconstruction</option>
311 <option value="file">initial trees from file (option "-t"). Also used for consensus and for tree distance calculations</option>
312 <help><![CDATA[
313 The options for the initial tree estimation are parsimony (default), fast parsimony, fasttree emulation, bioNJ, fixed (no optimisation), and from file. <br/>
314 The default parsimony estimation is through the PLL library, which is linear on the number of <i>sites</i> and not <i>patterns</i> (so slower than other algorithms).
315 Notice that this library is also used to populate the "initial parsimony trees" set ("<i>-ninit</i>" below), so keep that in mind. <br/>
316
317 <br/>According to the documentation (see http://www.iqtree.org/doc/Command-Reference#tree-search-parameters), the option "-fast"
318 "Turn on the fast tree search mode, where IQ-TREE will just construct two starting trees: maximum parsimony and BIONJ, which are then optimized
319 by nearest neighbor interchange (NNI)."
320 If you chose fast parsimony, fasttree, or bioNJ, it is assumed that the initial set of trees is two ("-ninit 2") since doing otherwise would defeat the purpose
321 of avoiding the default PLL.<br/>
322 <br/>You can also select a tree file as input to IQ_TREE. In most cases this tree file is used as
323 initial state for optimisation (or as a fixed topology). But in some cases it is used for consensus tree calculation or Robinson-Foulds distances.
324 There are a few reasons why you should provide trees to IQ-TREE:
325 <ul>
326 <li>Through a <b>constrained tree search</b>, where you give a (potentially small, multifurcating) tree which should be respected when searching for optimal trees.
327 This is given through option <i>"-g"</i>. You use this option when you "know" or want to test a specific grouping of taxa. (Section "Tree search parameters")</li>
328 <li>Input trees for <b>tree fit comparison</b>, as in section "Tree topology tests" below. IQ-TREE may still do model optimisation, initial tree search, etc. since
329 it compares the user-provided trees assuming a common model. </li>
330 <li>Suggest the <b>initial trees</b> (option "<i>-t</i>"). From these, IQ-TREE will start exploring trees (to speed up the tree search in case you already have
331 good candidate trees). This same command (i.e. same tree file) can also be used for calculating consensus trees or distances between trees. The other options will
332 inform IQ-TREE how it should use this file (here, the default is to use it as starting tree). </li>
333 <li>By <b>fixing the tree</b> (option "<i>-te</i>"). This means that IQ-TREE will not try any optimisation with other topologies, and can be used e.g. in ancestral reconstruction. </li>
334 </ul>
335 ]]></help>
336 </param>
337 <when value="default">
338 <param argument="-ninit" type="integer" value="100" optional="true" label="Specify number of initial parsimony trees (default=100)"
339 help="This value must be strictly higher than the number of best trees to maintain (i.e. higher than 5 under default settings)."/>
340 <param argument="-sprrad" type="integer" value="6" optional="true" label="Specify SPR radius for the initial parsimony tree search"/>
341 </when>
342 <when value="file">
343 <param argument="-t" type="data" format="nhx" optional="false" label="Specify a file containing the starting tree for tree search"
344 help="also used by consensus or distance calculations, if you set the appropriate variables in 'Tree Parameters' below"/>
345 </when>
346 <when value="fixed">
347 <param argument="-te" type="data" format="nhx" optional="false" label="file with fixed user tree" help="No tree search is performed and IQ-TREE
348 computes the log-likelihood of the fixed user tree."/>
349 </when>
350 <when value="fast"/>
351 <when value="pars"/>
352 <when value="bionj"/>
353 </conditional>
354 </section>
355
298 <section name="likelihood_mapping" expanded="False" title="Likelihood mapping analysis"> 356 <section name="likelihood_mapping" expanded="False" title="Likelihood mapping analysis">
299 <param argument="-lmap" type="integer" optional="true" label="Specify the number of quartets to be randomly drawn" help="If you specify -lmap ALL, all unique quartets will be drawn, instead."/> 357 <param argument="-lmap" type="integer" optional="true" label="Specify the number of quartets to be randomly drawn" help="If you specify -lmap ALL, all unique quartets will be drawn, instead."/>
300 <param argument="-lmclust" type="data" format="txt" optional="true" label="Specify a NEXUS file containing taxon clusters (see below for example) for quartet mapping analysis."/> 358 <param argument="-lmclust" type="data" format="txt" optional="true" label="Specify a NEXUS file containing taxon clusters (see below for example) for quartet mapping analysis."/>
301 <param argument="-wql" type="boolean" truevalue="-wql" falsevalue="" checked="false" label="Write quartet log-likelihoods into .lmap.quartetlh file (typically not needed)."/> 359 <param argument="-wql" type="boolean" truevalue="-wql" falsevalue="" checked="false" label="Write quartet log-likelihoods into .lmap.quartetlh file (typically not needed)."/>
302 </section> 360 </section>
303 <section name="modelling_parameters" title="Modelling Parameters"> 361 <section name="modelling_parameters" title="Modelling Parameters">
304 <section name="automatic_model" expanded="False" title="Automatic model selection"> 362 <section name="automatic_model" expanded="False" title="Automatic model selection">
305 <conditional name="cond_model" > 363 <conditional name="cond_model" >
306 <param name="opt_custommodel" type="boolean" checked="false" label="Use Custom Model" help="See http://www.iqtree.org/doc/Substitution-Models"/> 364 <param name="opt_custommodel" type="boolean" checked="false" label="Use Custom Model">
365 <help><![CDATA[
366 With this option, you can restrict the analysis to a fixed model instead of letting iq-tree select the best one (slow model selection).
367 Suggestions include "HKY+G" (good compromise speed-performance) or "GTR+G+I" (parameter-rich model).
368 <br>See http://www.iqtree.org/doc/Substitution-Models for a full list]]>
369 </help></param>
307 <when value="true"> 370 <when value="true">
308 <param argument="-m" type="text" label="Model"> 371 <param argument="-m" type="text" label="Model">
309 <sanitizer> 372 <sanitizer>
310 <valid initial="string.ascii_uppercase"> 373 <valid initial="string.ascii_uppercase">
311 <add value="&#43;" /> 374 <add value="&#43;" />
424 <param argument="-fmax" type="boolean" truevalue="-fmax" falsevalue="" checked="false" label="Switch to posterior maximum mode for obtaining site-specific profiles" help="Default: posterior mean."/> 487 <param argument="-fmax" type="boolean" truevalue="-fmax" falsevalue="" checked="false" label="Switch to posterior maximum mode for obtaining site-specific profiles" help="Default: posterior mean."/>
425 </section> 488 </section>
426 </section> 489 </section>
427 <section name="tree_parameters" title="Tree Parameters"> 490 <section name="tree_parameters" title="Tree Parameters">
428 <section name="tree_search" expanded="False" title="Tree search parameters"> 491 <section name="tree_search" expanded="False" title="Tree search parameters">
429 <param argument="-fast" type="boolean" truevalue="-fast" falsevalue="" checked="false" label="fast tree search mode"> 492 <param argument="-nstop" type="integer" value="100" optional="true" label="Specify number of unsuccessful iterations to stop
430 <help><![CDATA[ 493 (-nstop and -pers are the most relevant when increasing exploration space, according to IQ_TREE reference manual)"/>
431 <br/>Emulates fasttree. According to the documentation (see http://www.iqtree.org/doc/Command-Reference#tree-search-parameters): 494 <param argument="-pers" type="float" value="0.5" optional="true" label="Specify perturbation strength (between 0 and 1) for randomized NNI
432 <br/>"Turn on the fast tree search mode, where IQ-TREE will just construct two starting trees: maximum parsimony and BIONJ, which are then optimized 495 (-nstop and -pers are the most relevant when increasing exploration space, according to IQ_TREE reference manual)"/>
433 by nearest neighbor interchange (NNI). Introduced in version 1.6."
434 <br/> Some options below may therefore be unused (e.g. the initial number of parsimony trees...)
435 ]]>
436 </help>
437 </param>
438 <param argument="-n" type="integer" optional="true" label="Skip subsequent tree search (n = 0), useful when you only want to assess the phylogenetic information of the alignment. Otherwise specify number of iterations to stop (this option overrides -nstop criterion)." /> 496 <param argument="-n" type="integer" optional="true" label="Skip subsequent tree search (n = 0), useful when you only want to assess the phylogenetic information of the alignment. Otherwise specify number of iterations to stop (this option overrides -nstop criterion)." />
439 <param argument="-ninit" type="integer" value="100" optional="true" label="Specify number of initial parsimony trees"/> 497 <param argument="-nbest" type="integer" value="5" optional="true" label="Specify number of trees in the candidate set to maintain during ML tree search (default=5)"
498 help="This value must be strictly smaller than the number of initial pasimony trees (i.e. less than 100 for default settings)"/>
440 <param argument="-ntop" type="integer" value="20" optional="true" label="Specify number of top initial parsimony trees to optimize with ML nearest neighbor interchange (NNI) search to initialize the candidate set"/> 499 <param argument="-ntop" type="integer" value="20" optional="true" label="Specify number of top initial parsimony trees to optimize with ML nearest neighbor interchange (NNI) search to initialize the candidate set"/>
441 <param argument="-nbest" type="integer" value="5" optional="true" label="Specify number of trees in the candidate set to maintain during ML tree search"/>
442 <param argument="-nstop" type="integer" value="100" optional="true" label="Specify number of unsuccessful iterations to stop"/>
443 <param argument="-sprrad" type="integer" value="6" optional="true" label="Specify SPR radius for the initial parsimony tree search"/>
444 <param argument="-pers" type="float" value="0.5" optional="true" label="Specify perturbation strength (between 0 and 1) for randomized NNI"/>
445 <param argument="-allnni" type="boolean" truevalue="-allnni" falsevalue="" checked="false" label="Turn on more thorough and slower NNI search"/> 500 <param argument="-allnni" type="boolean" truevalue="-allnni" falsevalue="" checked="false" label="Turn on more thorough and slower NNI search"/>
446 <param argument="-djc" type="boolean" truevalue="-djc" falsevalue="" checked="false" label="Avoid computing ML pairwise distances and BIONJ tree."/> 501 <param argument="-djc" type="boolean" truevalue="-djc" falsevalue="" checked="false" label="Avoid computing ML pairwise distances and BIONJ tree (when generating tree set)"
447 <param argument="-g" type="data" format="txt" optional="true" label="Specify a topological constraint tree file in NEWICK format" help="The constraint tree can be a multifurcating tree and need not to include all taxa."/> 502 help="This is independent on how the initial tree/trees were found, and is added to set of initial trees (since it was estimated using likelihood parameters) "/>
503 <param argument="-g" type="data" format="txt" optional="true" label="Specify a topological constraint tree file in NEWICK format" help="The constraint tree can be a multifurcating tree and need not to include all taxa. (see https://github.com/Cibiv/IQ-TREE/wiki/Advanced-Tutorial#constrained-tree-search)"/>
448 </section> 504 </section>
449 <section name="single_branch" expanded="False" title="Single branch tests"> 505 <section name="single_branch" expanded="False" title="Single branch tests">
450 <param argument="-alrt" type="integer" optional="true" label="Specify number of replicates (&gt;=1000) to perform SH-like approximate likelihood ratio test (SH-aLRT) (Guindon et al., 2010)" help="If number of replicates is set to 0 (-alrt 0), then the parametric aLRT test (Anisimova and Gascuel 2006) is performed, instead of SH-aLRT."/> 506 <param argument="-alrt" type="integer" optional="true" label="Specify number of replicates (&gt;=1000) to perform SH-like approximate likelihood ratio test (SH-aLRT) (Guindon et al., 2010)" help="If number of replicates is set to 0 (-alrt 0), then the parametric aLRT test (Anisimova and Gascuel 2006) is performed, instead of SH-aLRT."/>
451 <param argument="-abayes" type="boolean" truevalue="-abayes" falsevalue="" checked="false" label="Perform approximate Bayes test (Anisimova et al., 2011)."/> 507 <param argument="-abayes" type="boolean" truevalue="-abayes" falsevalue="" checked="false" label="Perform approximate Bayes test (Anisimova et al., 2011)."/>
452 <param argument="-lbp" type="integer" optional="true" label="Specify number of replicates (&gt;=1000) to perform fast local bootstrap probability method (Adachi and Hasegawa, 1996)."/> 508 <param argument="-lbp" type="integer" optional="true" label="Specify number of replicates (&gt;=1000) to perform fast local bootstrap probability method (Adachi and Hasegawa, 1996)."/>
453 </section> 509 </section>
454 <section name="tree_topology" expanded="False" title="Tree topology tests"> 510 <section name="tree_topology" expanded="False" title="Tree topology tests">
455 <param argument="-z" type="data" format="txt" optional="true" label="Specify a file containing a set of trees" help="IQ-TREE will compute the log-likelihoods of all trees."/> 511 <param argument="-z" type="data" format="txt" optional="true" label="Specify a file containing a set of trees">
512 <help><![CDATA[
513 IQ-TREE will compute and compare the log-likelihoods of all trees in this file. It may still search for an optimal tree, in order to find best model and parameters
514 (unless you fixed those as well). It is strongly suggested to also perform some tests below (like RELL, AU, etc.). See more info at
515 https://github.com/Cibiv/IQ-TREE/wiki/Advanced-Tutorial#tree-topology-tests
516 ]]>
517 </help>
518 </param>
456 <param argument="-zb" type="integer" optional="true" label="Specify the number of RELL (Kishino et al., 1990) replicates (&gt;=1000) to perform several tree topology tests for all trees passed via -z" help="The tests include bootstrap proportion (BP), KH test (Kishino and Hasegawa, 1989), SH test (Shimodaira and Hasegawa, 1999) and expected likelihood weights (ELW) (Strimmer and Rambaut, 2002)."/> 519 <param argument="-zb" type="integer" optional="true" label="Specify the number of RELL (Kishino et al., 1990) replicates (&gt;=1000) to perform several tree topology tests for all trees passed via -z" help="The tests include bootstrap proportion (BP), KH test (Kishino and Hasegawa, 1989), SH test (Shimodaira and Hasegawa, 1999) and expected likelihood weights (ELW) (Strimmer and Rambaut, 2002)."/>
457 <param argument="-zw" type="boolean" truevalue="-zw" falsevalue="" checked="false" label="Used together with -zb to additionally perform the weighted-KH and weighted-SH tests."/> 520 <param argument="-zw" type="boolean" truevalue="-zw" falsevalue="" checked="false" label="Used together with -zb to additionally perform the weighted-KH and weighted-SH tests."/>
458 <param argument="-au" type="boolean" truevalue="-au" falsevalue="" checked="false" label="Used together with -zb to additionally perform the approximately unbiased (AU) test (Shimodaira, 2002)" help="Note that you have to specify the number of replicates for the AU test via -zb."/> 521 <param argument="-au" type="boolean" truevalue="-au" falsevalue="" checked="false" label="Used together with -zb to additionally perform the approximately unbiased (AU) test (Shimodaira, 2002)" help="Note that you have to specify the number of replicates for the AU test via -zb."/>
459 </section> 522 </section>
460 <section name="constructing_consensus" expanded="False" title="Constructing consensus tree"> 523 <section name="constructing_consensus" expanded="False" title="Constructing consensus tree">