|
3
|
1 <?xml version="1.0"?>
|
|
|
2 <tool id="mrbayes" name="MrBayes">
|
|
|
3 <description>with options and commands</description>
|
|
4
|
4 <command> mb $command_file > $outfile
|
|
|
5 #silent sys.stderr.write("\n \n \n !!!! Cheetah Template Variables !!!!\n")
|
|
|
6 #for k,v in $searchList[2].items()
|
|
|
7 #silent sys.stderr.write(" %s = %s\n" % (str(k), str(v) ))
|
|
|
8 #end for
|
|
|
9 #silent sys.stderr.write("!!!! end-of-list !!!!\n \n \n")
|
|
|
10 </command>
|
|
3
|
11 <inputs>
|
|
|
12 <param name="data" type="data" format="nex" label="Nexus input file"/>
|
|
|
13
|
|
|
14 <param name="outgroup" type="text" label="Outgroup" size="10" value="1"
|
|
|
15 help="The name or number identifying one taxon from the data." />
|
|
|
16 <param name="model" type="select" label="Choose model">
|
|
|
17 <sanitizer sanitize="False" />
|
|
|
18 <option value="lset nst=6" >GRT</option>
|
|
|
19 <option value="lset nst=6; prset statefreqpr=fixed(equal)">SYM</option>
|
|
|
20 <option value="lset nst=2" >HKY</option>
|
|
|
21 <option value="lset nst=2; prset statefreqpr=fixed(equal)">K2P</option>
|
|
|
22 <option value="lset nst=1" >F81</option>
|
|
|
23 <option value="lset nst=1; prset statefreqpr=fixed(equal)">JC </option>
|
|
|
24 <option value="lset nst=mixed" >Mixed </option>
|
|
|
25 <help>Models are defined below and are only valid for nucleotide data.</help>
|
|
|
26 </param>
|
|
|
27 <param name="rates" type="select" label="Choose rates">
|
|
|
28 <option value="equal"/>
|
|
|
29 <option value="gamma"/>
|
|
|
30 <option value="propinv"/>
|
|
|
31 <option value="invgamma"/>
|
|
|
32 <option value="adgamma"/>
|
|
|
33 </param>
|
|
|
34 <param name="ngen" type="integer" label="Number of generations" value="1000" >
|
|
|
35 <validator type="in_range" min="1" max="inf" message="Must be greater than or equal to 1"/>
|
|
|
36 </param>
|
|
|
37 <param name="nchain" type="integer" label="Number of chains" value="2" >
|
|
|
38 <validator type="in_range" min="1" max="inf" message="Must be greater than or equal to 1"/>
|
|
|
39 </param>
|
|
|
40 <param name="nrun" type="integer" label="Number of runs" value="2" >
|
|
|
41 <validator type="in_range" min="1" max="inf" message="Must be greater than or equal to 1"/>
|
|
|
42 </param>
|
|
|
43 <param name="checkfreq" type="integer" label="Checkpoint frequency" value="100000" >
|
|
|
44 <validator type="in_range" min="100" max="inf" message="Must be greater than or equal to 100"/>
|
|
4
|
45 <help>Frequency with which checkpoints are written. Used for continuing analysis with append.</help>
|
|
3
|
46 </param>
|
|
|
47 <param name="samplefreq" type="integer" label="Sample frequency" value="500" >
|
|
|
48 <validator type="in_range" min="1" max="inf" message="Must be greater than or equal to 1"/>
|
|
|
49 <help>Frequency with which output it written to files." </help>
|
|
|
50 </param>
|
|
|
51 <param name="printfreq" type="integer" label="Print frequency" value="500" >
|
|
|
52 <validator type="in_range" min="1" max="inf" message="Must be greater than or equal to 1"/>
|
|
|
53 <help>Frequency with which output is printed to the screen.</help>
|
|
|
54 </param>
|
|
|
55 <param name="burninfrac" type="float" label="Burn-in fraction" value="0.25">
|
|
|
56 <validator type="in_range" min="0.01" max=".50" message="Must be between 0.01 and .50"/>
|
|
|
57 <help>Fraction of samples to discard when summarizing. </help>
|
|
|
58 </param>
|
|
|
59 <param name="stoprule" type="select" label="Stop rule?" >
|
|
|
60 <option value="no">No</option>
|
|
|
61 <option value="yes">Yes</option>
|
|
|
62 </param>
|
|
|
63
|
|
|
64 <param name="sump" type="boolean" truevalue="sump;" falsevalue="" label="Include sump?"
|
|
|
65 help="Summarizes sampled parameter value."/>
|
|
|
66 <param name="sumt" type="boolean" truevalue="sumt;" falsevalue="" label="Include sumt?"
|
|
|
67 help="Summarizes the trees."/>
|
|
|
68 <param name="plot" type="boolean" truevalue="plot;" falsevalue="" label="Include plot?"
|
|
|
69 help="Creates an x-y graph of the parameter over the course of the chain."/>
|
|
|
70
|
|
4
|
71 <param name="append" type="boolean" truevalue="mcmcp append=yes" falsevalue="mcmcp append=no" label="Append to last checkpoint?"
|
|
|
72 help="Requires .ckp file (i.e. that ngen > checkfreq)." />
|
|
|
73
|
|
3
|
74 </inputs>
|
|
|
75
|
|
|
76 <configfiles>
|
|
|
77 <configfile name="command_file">
|
|
|
78 begin mrbayes;
|
|
|
79 set quitonerror=no;
|
|
|
80 execute $data;
|
|
|
81 outgroup $outgroup;
|
|
|
82 $model;
|
|
|
83 lset rates=$rates;
|
|
|
84 mcmcp ngen=$ngen nrun=$nrun nchain=$nchain checkfreq=$checkfreq samplefreq=$samplefreq printfreq=$printfreq stoprule=$stoprule burninfrac=$burninfrac;
|
|
4
|
85 #if int($ngen) >= int($checkfreq)
|
|
|
86 $append;
|
|
|
87 #end if#
|
|
3
|
88 mcmc;
|
|
|
89 $sump
|
|
|
90 $sumt
|
|
|
91 $plot
|
|
|
92 end;
|
|
|
93 </configfile>
|
|
|
94 </configfiles>
|
|
|
95
|
|
|
96 <outputs>
|
|
|
97 <data name="outfile" format="txt"/>
|
|
|
98 </outputs>
|
|
|
99
|
|
|
100 <requirements>
|
|
|
101 <requirement version="3.2.1" type="binary">mb</requirement>
|
|
|
102 </requirements>
|
|
|
103 <help>
|
|
|
104
|
|
|
105 .. class:: warningmark
|
|
|
106
|
|
4
|
107 This tool takes as input a Nexus file (.nex), described here_, that should only contain a data block, without commands. If a command block exists in the file, it will be executed first, followed by the ones specified above.
|
|
3
|
108
|
|
|
109 .. _here: http://en.wikipedia.org/wiki/Nexus_file.
|
|
|
110
|
|
|
111 -----
|
|
|
112
|
|
|
113 **What it does**
|
|
|
114
|
|
|
115 MrBayes is a program for Bayesian inference and model choice across a wide range of phylogenetic and evolutionary models. MrBayes uses Markov chain Monte Carlo (MCMC) methods to estimate the posterior distribution of model parameters.
|
|
|
116
|
|
4
|
117 See the manual_ or the `command reference`_ for more information.
|
|
|
118
|
|
|
119 .. _manual: http://mrbayes.sourceforge.net/manual.php
|
|
|
120 .. _command reference: http://mrbayes.sourceforge.net/commref_mb3.2.pdf
|
|
|
121
|
|
3
|
122 -----
|
|
|
123
|
|
|
124 **Model Definitions**
|
|
|
125
|
|
|
126 Models are defined by the following MrBayes options::
|
|
|
127
|
|
|
128 GRT: lset nst=6;
|
|
|
129 SYM: lset nst=6; prset statefreqpr=fixed(equal);
|
|
|
130 HKY: lset nst=2;
|
|
|
131 K2P: lset nst=2; prset statefreqpr=fixed(equal);
|
|
|
132 F81: lset nst=1;
|
|
|
133 JC: lset nst=1; prset statefreqpr=fixed(equal);
|
|
|
134 Mixed: lset nst=mixed;
|
|
|
135
|
|
|
136 Each model makes different assumptions about nucelotide substitution rates and state frequencies.
|
|
|
137
|
|
|
138 The mixed model performs Markov chain sampling over the space of all possible reversible substitution models, by grouping the six rates in various combinations.
|
|
|
139
|
|
|
140 -----
|
|
|
141
|
|
|
142 .. class:: infomark
|
|
|
143
|
|
4
|
144 **Tip for using append**
|
|
3
|
145
|
|
4
|
146 If you want to use the append functionality, all parameters will have to be identical to your initial analysis. The easiest way to ensure this is by using the blue "Run this job again" arrow on the dataset, and checking the additional append box. Note also that ngen will be the total number of generations, not the number of generations to add, so you should increase this value except if a run was interrupted. Selecting append while checkfreq > ngen will have no effect.
|
|
3
|
147 </help>
|
|
|
148 </tool>
|
|
|
149
|
