bella-companion 0.0.25__py3-none-any.whl → 0.0.26__py3-none-any.whl
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Potentially problematic release.
This version of bella-companion might be problematic. Click here for more details.
- bella_companion/eucovid/beast_configs/GLM.xml +377 -128
- bella_companion/eucovid/beast_configs/Nonparametric.xml +331 -146
- bella_companion/eucovid/run_beast.py +15 -13
- bella_companion/platyrrhine/run_beast.py +1 -1
- bella_companion/utils/slurm.py +6 -1
- {bella_companion-0.0.25.dist-info → bella_companion-0.0.26.dist-info}/METADATA +1 -1
- {bella_companion-0.0.25.dist-info → bella_companion-0.0.26.dist-info}/RECORD +9 -9
- {bella_companion-0.0.25.dist-info → bella_companion-0.0.26.dist-info}/WHEEL +0 -0
- {bella_companion-0.0.25.dist-info → bella_companion-0.0.26.dist-info}/entry_points.txt +0 -0
|
@@ -1,213 +1,398 @@
|
|
|
1
1
|
<?xml version="1.0" encoding="UTF-8" standalone="no"?>
|
|
2
|
-
<beast
|
|
3
|
-
|
|
4
|
-
|
|
5
|
-
|
|
6
|
-
|
|
7
|
-
|
|
8
|
-
|
|
9
|
-
|
|
10
|
-
|
|
11
|
-
|
|
12
|
-
|
|
13
|
-
|
|
14
|
-
|
|
15
|
-
|
|
16
|
-
|
|
17
|
-
|
|
18
|
-
|
|
19
|
-
|
|
20
|
-
|
|
21
|
-
|
|
22
|
-
|
|
23
|
-
|
|
24
|
-
|
|
25
|
-
|
|
26
|
-
|
|
27
|
-
|
|
28
|
-
|
|
29
|
-
|
|
30
|
-
|
|
31
|
-
|
|
32
|
-
|
|
33
|
-
|
|
34
|
-
|
|
35
|
-
<map name="Uniform">beast.base.inference.distribution.Uniform</map>
|
|
36
|
-
<map name="Exponential">beast.base.inference.distribution.Exponential</map>
|
|
37
|
-
<map name="LogNormal">beast.base.inference.distribution.LogNormalDistributionModel</map>
|
|
38
|
-
<map name="prior">beast.base.inference.distribution.Prior</map>
|
|
39
|
-
|
|
40
|
-
|
|
41
|
-
|
|
42
|
-
|
|
43
|
-
|
|
44
|
-
<
|
|
45
|
-
|
|
46
|
-
|
|
47
|
-
<
|
|
48
|
-
|
|
49
|
-
|
|
50
|
-
|
|
51
|
-
|
|
52
|
-
|
|
53
|
-
|
|
2
|
+
<beast namespace="beast.base.inference
|
|
3
|
+
:beast.base.inference.parameter
|
|
4
|
+
:beast.base.inference.distribution
|
|
5
|
+
:beast.base.evolution.branchratemodel
|
|
6
|
+
:beast.base.evolution.substitutionmodel
|
|
7
|
+
:beast.base.evolution.sitemodel
|
|
8
|
+
:beast.base.evolution.tree.coalescent
|
|
9
|
+
:beast.base.evolution.tree
|
|
10
|
+
:beast.base.evolution.likelihood
|
|
11
|
+
:beast.base.evolution.alignment
|
|
12
|
+
:beast.base.evolution.operator.kernel
|
|
13
|
+
:beast.base.inference.operator.kernel
|
|
14
|
+
:beast.base.inference.operator
|
|
15
|
+
:beast.base.math.distributions
|
|
16
|
+
:beast.base.evolution
|
|
17
|
+
:beast.base.evolution.operator
|
|
18
|
+
:bdmmprime.util
|
|
19
|
+
:bdmmprime.util.operators
|
|
20
|
+
:bdmmprime.util.priors
|
|
21
|
+
:bdmmprime.distribution
|
|
22
|
+
:bdmmprime.parameterization
|
|
23
|
+
:bdmmprime.mapping
|
|
24
|
+
:bdmmprime.trajectories
|
|
25
|
+
:feast
|
|
26
|
+
:feast.parameter
|
|
27
|
+
:feast.fileio
|
|
28
|
+
:feast.function
|
|
29
|
+
:feast.expressions
|
|
30
|
+
:glmprior.util"
|
|
31
|
+
required="BEAST.base v2.7.7:BDMM-Prime v2.2.2" version="2.7">
|
|
32
|
+
|
|
33
|
+
|
|
34
|
+
<!-- Map elements -->
|
|
35
|
+
<map name="Uniform" >beast.base.inference.distribution.Uniform</map>
|
|
36
|
+
<map name="Exponential" >beast.base.inference.distribution.Exponential</map>
|
|
37
|
+
<map name="LogNormal" >beast.base.inference.distribution.LogNormalDistributionModel</map>
|
|
38
|
+
<map name="prior" >beast.base.inference.distribution.Prior</map>
|
|
39
|
+
|
|
40
|
+
<!-- Adapted from https://github.com/SarahNadeau/cov-europe-bdmm/blob/master/analyses/with_hubei_migration_decrease/with_hubei_migration_decrease_scrubbed.xml -->
|
|
41
|
+
<!-- Sequence data loaded from GISAID with the gisaid EPI ISLs -->
|
|
42
|
+
|
|
43
|
+
<!-- Alignment -->
|
|
44
|
+
<alignment id="alignment" spec= "AlignmentFromFasta"
|
|
45
|
+
fileName = "$(aligned_fasta)"/>
|
|
46
|
+
|
|
47
|
+
<typeSet id="typeSetBDMMPrime" spec="TypeSet">
|
|
48
|
+
<typeTraitSet id="typeTraitSet" spec="TraitSetFromTaxonSet"
|
|
49
|
+
traitname="type"
|
|
50
|
+
delimiter="|"
|
|
51
|
+
takeGroup="1">
|
|
52
|
+
<taxa spec="TaxonSet" alignment="@alignment"/>
|
|
53
|
+
</typeTraitSet>
|
|
54
|
+
</typeSet>
|
|
55
|
+
|
|
56
|
+
<trait id="dateTrait" spec="TraitSetFromTaxonSet"
|
|
57
|
+
traitname="date"
|
|
58
|
+
dateFormat="yyyy-M-dd"
|
|
59
|
+
delimiter= "|"
|
|
60
|
+
everythingAfterLast="true">
|
|
61
|
+
<taxa spec="TaxonSet" alignment="@alignment"/>
|
|
62
|
+
</trait>
|
|
63
|
+
|
|
64
|
+
<!-- BDMMPrime Population Model -->
|
|
65
|
+
|
|
66
|
+
<processLength id="originBDMMPrime" spec="RealParameter">1.0</processLength>
|
|
67
|
+
|
|
68
|
+
<Re id="ReSV" spec="SkylineVectorParameter"
|
|
69
|
+
timesAreAges="true" processLength="@originBDMMPrime" typeSet="@typeSetBDMMPrime">
|
|
70
|
+
<skylineValues id="ReValues" spec="RealParameter">1.0 1.1 1.2 1.3 1.4 1.01 1.1 1.2 1.3 1.4</skylineValues>
|
|
71
|
+
<changeTimes id="ReRateChangeTimes" spec="RealParameter">0.123</changeTimes>
|
|
54
72
|
</Re>
|
|
55
73
|
|
|
56
|
-
<samplingProportion id="
|
|
57
|
-
|
|
74
|
+
<samplingProportion id="samplingProportionSV" spec="SkylineVectorParameter"
|
|
75
|
+
timesAreAges="true" processLength="@originBDMMPrime" typeSet="@typeSetBDMMPrime">
|
|
76
|
+
<skylineValues id="samplingProportionValues" spec="RealParameter" dimension="15" lower="0.0" upper="1.0">1.11E-5 1.2E-5 1.3E-5 1.4E-5 1.5E-5 1.1E-5 1.2E-5 1.3E-5 1.4E-5 1.5E-5 0.0 0.0 0.0 0.0 0.0</skylineValues>
|
|
77
|
+
|
|
78
|
+
|
|
79
|
+
<changeTimes id="samplingProportionChangeTimes" spec="RealParameter">0.12 0.205</changeTimes>
|
|
58
80
|
</samplingProportion>
|
|
81
|
+
|
|
59
82
|
|
|
60
|
-
<becomeUninfectiousRate id="
|
|
61
|
-
|
|
62
|
-
|
|
63
|
-
|
|
64
|
-
|
|
65
|
-
|
|
66
|
-
|
|
67
|
-
|
|
68
|
-
|
|
69
|
-
|
|
70
|
-
|
|
71
|
-
|
|
72
|
-
|
|
73
|
-
|
|
74
|
-
|
|
75
|
-
|
|
76
|
-
|
|
77
|
-
|
|
78
|
-
|
|
79
|
-
|
|
80
|
-
|
|
81
|
-
|
|
82
|
-
|
|
83
|
-
|
|
84
|
-
|
|
85
|
-
|
|
86
|
-
|
|
87
|
-
|
|
88
|
-
|
|
89
|
-
|
|
90
|
-
|
|
91
|
-
|
|
83
|
+
<becomeUninfectiousRate id="becomeUninfectiousRateSV" spec="SkylineVectorParameter"
|
|
84
|
+
typeSet="@typeSetBDMMPrime">
|
|
85
|
+
<skylineValues id="becomeUninfectiousValues" spec="RealParameter">36.5</skylineValues>
|
|
86
|
+
</becomeUninfectiousRate>
|
|
87
|
+
|
|
88
|
+
|
|
89
|
+
<migrationRate id="migrationRateSM" spec="SkylineMatrixParameter"
|
|
90
|
+
timesAreAges="true" processLength="@originBDMMPrime" typeSet="@typeSetBDMMPrime">
|
|
91
|
+
<skylineValues id="migrationRateValues_scaled" name="skylineValues" spec="feast.expressions.ExpCalculator">
|
|
92
|
+
<![CDATA[
|
|
93
|
+
migrationRateValues *
|
|
94
|
+
{ f,f,f,f,
|
|
95
|
+
f, 1,1,1,
|
|
96
|
+
f,1, 1,1,
|
|
97
|
+
f,1,1, 1,
|
|
98
|
+
f,1,1,1,
|
|
99
|
+
|
|
100
|
+
1,1,1,1,
|
|
101
|
+
1, 1,1,1,
|
|
102
|
+
1,1, 1,1,
|
|
103
|
+
1,1,1, 1,
|
|
104
|
+
1,1,1,1 }
|
|
105
|
+
]]>
|
|
106
|
+
<arg id="migrationRateValues" spec="RealParameter" dimension="40">
|
|
107
|
+
0.11 0.15 0.2 0.24
|
|
108
|
+
0.0 0.14 0.18 0.22
|
|
109
|
+
0.0 0.1 0.19 0.23
|
|
110
|
+
0.0 0.12 0.16 0.25
|
|
111
|
+
0.0 0.13 0.17 0.21
|
|
112
|
+
|
|
113
|
+
0.11 0.15 0.2 0.24
|
|
114
|
+
0.0 0.14 0.18 0.22
|
|
115
|
+
0.0 0.1 0.19 0.23
|
|
116
|
+
0.0 0.12 0.16 0.25
|
|
117
|
+
0.0 0.13 0.17 0.21
|
|
118
|
+
</arg>
|
|
119
|
+
<arg id="f" spec="RealParameter">0.5</arg>
|
|
92
120
|
</skylineValues>
|
|
121
|
+
|
|
122
|
+
<changeTimes id="migrationRateChangeTimes" spec="RealParameter">0.123</changeTimes>
|
|
93
123
|
</migrationRate>
|
|
94
124
|
|
|
95
|
-
<removalProb id="removalProbSP" spec="SkylineVectorParameter" skylineValues="1.0" processLength="@processLength" typeSet="@typeSet"/>
|
|
96
125
|
|
|
97
|
-
<
|
|
98
|
-
|
|
99
|
-
|
|
100
|
-
|
|
126
|
+
<removalProb id="removalProbSV" spec="SkylineVectorParameter" skylineValues="1.0"
|
|
127
|
+
processLength="@originBDMMPrime" typeSet="@typeSetBDMMPrime"/>
|
|
128
|
+
|
|
129
|
+
<startTypePriorProbs id="typeFrequencies" spec="RealParameter"
|
|
130
|
+
estimate="false" lower="0.0" upper="1.0">1.0 0.0 0.0 0.0 0.0</startTypePriorProbs>
|
|
131
|
+
|
|
132
|
+
|
|
133
|
+
<!-- Site Model -->
|
|
134
|
+
<siteModel id="siteModel" spec="SiteModel" gammaCategoryCount="4">
|
|
135
|
+
<parameter id="gammaShape" name="shape" value="1.0"/>
|
|
136
|
+
<parameter name="mutationRate" estimate="false" value="1.0"/>
|
|
137
|
+
<parameter name="proportionInvariant" estimate="false" value="0.0"/>
|
|
138
|
+
|
|
101
139
|
<substModel spec="HKY">
|
|
102
140
|
<kappa id="kappa" spec="RealParameter" value="2.0"/>
|
|
103
|
-
<frequencies spec="Frequencies" data="@alignment"/>
|
|
141
|
+
<frequencies id="empiricalFreqs" spec="Frequencies" data="@alignment"/>
|
|
104
142
|
</substModel>
|
|
105
143
|
</siteModel>
|
|
106
144
|
|
|
107
|
-
<branchRateModel id="branchRateModel" spec="StrictClockModel" clock.rate="8.0E-4"/>
|
|
108
145
|
|
|
109
|
-
|
|
146
|
+
<!-- Branch rate Model -->
|
|
147
|
+
<branchRateModel id="branchRateModel" spec="StrictClockModel">
|
|
148
|
+
<parameter id="clockRate" name="clock.rate" estimate="false">8.0E-4</parameter>
|
|
149
|
+
</branchRateModel>
|
|
150
|
+
|
|
110
151
|
|
|
111
|
-
|
|
152
|
+
<!-- Population Model for Tree initialiser -->
|
|
153
|
+
<popFunc id="popFunc" spec="ConstantPopulation">
|
|
154
|
+
<popSize spec="RealParameter" value="0.1"/>
|
|
155
|
+
</popFunc>
|
|
112
156
|
|
|
157
|
+
|
|
158
|
+
<!-- MCMC -->
|
|
159
|
+
<run id="mcmc" spec="MCMC" chainLength="10000000" numInitializationAttempts="100">
|
|
160
|
+
|
|
161
|
+
<!-- State -->
|
|
113
162
|
<state id="state" spec="State" storeEvery="5000">
|
|
114
|
-
|
|
115
|
-
<stateNode
|
|
116
|
-
|
|
117
|
-
|
|
118
|
-
|
|
163
|
+
|
|
164
|
+
<stateNode id="treeBDMMPrime" spec="RandomTree"
|
|
165
|
+
taxa="@alignment"
|
|
166
|
+
populationModel="@popFunc"
|
|
167
|
+
trait="@dateTrait">
|
|
168
|
+
</stateNode>
|
|
169
|
+
|
|
170
|
+
<stateNode idref="originBDMMPrime"/>
|
|
171
|
+
<stateNode idref="ReValues"/>
|
|
172
|
+
<stateNode idref="samplingProportionValues"/>
|
|
173
|
+
<stateNode idref="migrationRateValues"/>
|
|
119
174
|
<stateNode idref="f"/>
|
|
175
|
+
|
|
120
176
|
<stateNode idref="kappa"/>
|
|
121
177
|
<stateNode idref="gammaShape"/>
|
|
178
|
+
|
|
122
179
|
</state>
|
|
180
|
+
|
|
123
181
|
|
|
182
|
+
<!-- Posterior Distribution -->
|
|
124
183
|
<distribution id="posterior" spec="CompoundDistribution">
|
|
184
|
+
|
|
185
|
+
<!-- Tree Likelihood -->
|
|
125
186
|
<distribution id="likelihood" spec="CompoundDistribution" useThreads="true">
|
|
126
|
-
|
|
187
|
+
|
|
188
|
+
<distribution id="treeLikelihood" spec="ThreadedTreeLikelihood"
|
|
189
|
+
data="@alignment"
|
|
190
|
+
tree="@treeBDMMPrime"
|
|
191
|
+
siteModel="@siteModel"
|
|
192
|
+
branchRateModel="@branchRateModel"/>
|
|
193
|
+
|
|
127
194
|
</distribution>
|
|
128
195
|
|
|
196
|
+
<!-- Prior distribution -->
|
|
129
197
|
<distribution id="prior" spec="CompoundDistribution">
|
|
130
|
-
|
|
131
|
-
|
|
198
|
+
|
|
199
|
+
<!-- BDMM-Prime tree prior distribution -->
|
|
200
|
+
<distribution id="BDMMPrime" spec="BirthDeathMigrationDistribution"
|
|
201
|
+
tree="@treeBDMMPrime"
|
|
202
|
+
conditionOnSurvival="true"
|
|
203
|
+
finalSampleOffset="0.0"
|
|
204
|
+
typeTraitSet="@typeTraitSet"
|
|
205
|
+
startTypePriorProbs="@typeFrequencies">
|
|
206
|
+
|
|
207
|
+
<!-- Parameterization BDMM-Prime -->
|
|
208
|
+
<parameterization id="EpiBDMMPrimeParameterization" spec="EpiParameterization"
|
|
209
|
+
processLength="@originBDMMPrime"
|
|
210
|
+
Re="@ReSV"
|
|
211
|
+
samplingProportion="@samplingProportionSV"
|
|
212
|
+
becomeUninfectiousRate="@becomeUninfectiousRateSV"
|
|
213
|
+
migrationRate="@migrationRateSM"
|
|
214
|
+
removalProb="@removalProbSV"
|
|
215
|
+
typeSet="@typeSetBDMMPrime"/>
|
|
216
|
+
|
|
132
217
|
</distribution>
|
|
133
218
|
|
|
134
|
-
|
|
219
|
+
|
|
220
|
+
<!-- Priors on parameters -->
|
|
221
|
+
<!-- BDMMPrime parameters priors -->
|
|
222
|
+
|
|
223
|
+
<prior id="originBDMMPrimePrior" name="distribution" x="@originBDMMPrime">
|
|
135
224
|
<LogNormal name="distr" M="-1.0" S="0.2"/>
|
|
136
225
|
</prior>
|
|
137
|
-
|
|
226
|
+
|
|
227
|
+
<!-- Reproductive number Prior -->
|
|
228
|
+
<distribution id="RePrior" spec="SmartZeroExcludingPrior" x="@ReValues">
|
|
138
229
|
<LogNormal name="distr" M="0.8" S="0.5"/>
|
|
139
230
|
</distribution>
|
|
140
|
-
|
|
231
|
+
|
|
232
|
+
<!-- Prior for sampling proportion for China -->
|
|
233
|
+
<distribution id="spChinaE2Prior" spec="SmartZeroExcludingPrior" x="@samplingProportionValues"
|
|
234
|
+
classesToExclude="1.2E-5 1.3E-5 1.4E-5 1.5E-5">
|
|
141
235
|
<Uniform name="distr" lower="0" upper="0.15"/>
|
|
142
236
|
</distribution>
|
|
143
|
-
|
|
237
|
+
<!-- Prior for sampling proportion for France -->
|
|
238
|
+
<distribution id="spFranceE1Prior" spec="SmartZeroExcludingPrior" x="@samplingProportionValues"
|
|
239
|
+
classesToExclude="1.1E-5 1.3E-5 1.4E-5 1.5E-5">
|
|
144
240
|
<Uniform name="distr" lower="0" upper="0.093"/>
|
|
145
241
|
</distribution>
|
|
146
|
-
|
|
242
|
+
<!-- Prior for sampling proportion for Germany Epoch -->
|
|
243
|
+
<distribution id="spGermanyE1Prior" spec="SmartZeroExcludingPrior" x="@samplingProportionValues"
|
|
244
|
+
classesToExclude="1.1E-5 1.2E-5 1.4E-5 1.5E-5">
|
|
147
245
|
<Uniform name="distr" lower="0" upper="0.10"/>
|
|
148
246
|
</distribution>
|
|
149
|
-
|
|
247
|
+
<!-- Prior for sampling proportion for Italy -->
|
|
248
|
+
<distribution id="spItalyE1Prior" spec="SmartZeroExcludingPrior" x="@samplingProportionValues"
|
|
249
|
+
classesToExclude="1.1E-5 1.2E-5 1.3E-5 1.5E-5">
|
|
150
250
|
<Uniform name="distr" lower="0" upper="0.005"/>
|
|
151
251
|
</distribution>
|
|
152
|
-
|
|
252
|
+
<!-- Prior for sampling proportion for OtherEuropean -->
|
|
253
|
+
<distribution id="spOEE1Prior" spec="SmartZeroExcludingPrior" x="@samplingProportionValues"
|
|
254
|
+
classesToExclude="1.1E-5 1.2E-5 1.3E-5 1.4E-5">
|
|
153
255
|
<Uniform name="distr" lower="0" upper="0.057"/>
|
|
154
256
|
</distribution>
|
|
155
|
-
|
|
257
|
+
|
|
258
|
+
<!-- Migration rate Prior -->
|
|
259
|
+
<distribution id="migRatePrior" name="distribution" spec="SmartZeroExcludingPrior" x="@migrationRateValues">
|
|
156
260
|
<LogNormal name="distr" M="0.0" S="1.0"/>
|
|
157
261
|
</distribution>
|
|
158
|
-
|
|
262
|
+
|
|
263
|
+
<distribution id="fPrior" name="distribution" spec="ZeroExcludingPrior" x="@f">
|
|
159
264
|
<Uniform name="distr" lower="0.0" upper="1.0"/>
|
|
160
265
|
</distribution>
|
|
161
|
-
|
|
162
|
-
|
|
266
|
+
|
|
267
|
+
<!-- Substitution model parameter priors -->
|
|
268
|
+
<prior id="GammaShapePrior" name="distribution" x="@gammaShape">
|
|
269
|
+
<Exponential name="distr">
|
|
270
|
+
<parameter spec="RealParameter" estimate="false" name="mean">0.5</parameter>
|
|
271
|
+
</Exponential>
|
|
163
272
|
</prior>
|
|
164
|
-
|
|
165
|
-
|
|
273
|
+
|
|
274
|
+
<prior id="KappaPrior" name="distribution" x="@kappa">
|
|
275
|
+
<LogNormal name="distr">
|
|
276
|
+
<parameter spec="RealParameter" estimate="false" name="M">1.0</parameter>
|
|
277
|
+
<parameter spec="RealParameter" estimate="false" name="S">1.25</parameter>
|
|
278
|
+
</LogNormal>
|
|
166
279
|
</prior>
|
|
280
|
+
|
|
167
281
|
</distribution>
|
|
168
282
|
</distribution>
|
|
169
283
|
|
|
170
|
-
|
|
171
|
-
|
|
172
|
-
|
|
173
|
-
<operator spec="
|
|
174
|
-
|
|
175
|
-
|
|
176
|
-
<operator spec="
|
|
177
|
-
|
|
178
|
-
|
|
179
|
-
<operator
|
|
180
|
-
|
|
181
|
-
|
|
182
|
-
<operator
|
|
183
|
-
|
|
184
|
-
|
|
185
|
-
<operator
|
|
186
|
-
|
|
187
|
-
|
|
284
|
+
|
|
285
|
+
<!-- Operators -->
|
|
286
|
+
<!-- Tree operators -->
|
|
287
|
+
<operator id="BDMMPrimeTreeRootScaler" spec="BactrianScaleOperator" tree="@treeBDMMPrime"
|
|
288
|
+
rootOnly="true" scaleFactor="0.25" weight="5.0"/>
|
|
289
|
+
|
|
290
|
+
<operator id="BDMMPrimeUniformOperator" spec="BactrianNodeOperator" tree="@treeBDMMPrime"
|
|
291
|
+
weight="30.0"/>
|
|
292
|
+
|
|
293
|
+
<operator id="BDMMPrimeSubtreeSlideScaler" spec="BactrianSubtreeSlide" tree="@treeBDMMPrime"
|
|
294
|
+
weight="15.0"/>
|
|
295
|
+
|
|
296
|
+
<operator id="BDMMPrimeNarrow" spec="Exchange" tree="@treeBDMMPrime"
|
|
297
|
+
weight="15.0"/>
|
|
298
|
+
|
|
299
|
+
<operator id="BDMMPrimeWide" spec="Exchange" tree="@treeBDMMPrime"
|
|
300
|
+
isNarrow="false" weight="3.0"/>
|
|
301
|
+
|
|
302
|
+
<operator id="BDMMPrimeWilsonBalding" spec="WilsonBalding" tree="@treeBDMMPrime"
|
|
303
|
+
weight="3.0"/>
|
|
304
|
+
|
|
305
|
+
<operator id="BDMMPrimeBICEPSEpochTop" spec="EpochFlexOperator" tree="@treeBDMMPrime"
|
|
306
|
+
scaleFactor="0.1" weight="2.0"/>
|
|
307
|
+
|
|
308
|
+
<operator id="BDMMPrimeBICEPSEpochAll" spec="EpochFlexOperator" tree="@treeBDMMPrime"
|
|
309
|
+
fromOldestTipOnly="false" scaleFactor="0.1" weight="2.0"/>
|
|
310
|
+
|
|
311
|
+
<operator id="BDMMPrimeBICEPSTreeFlex" spec="TreeStretchOperator" tree="@treeBDMMPrime"
|
|
312
|
+
scaleFactor="0.01" weight="2.0"/>
|
|
313
|
+
|
|
314
|
+
|
|
315
|
+
<!-- Population model operators -->
|
|
316
|
+
<operator id="originBDMMPrimeScaler" spec="BactrianScaleOperator" parameter="@originBDMMPrime"
|
|
317
|
+
scaleFactor="0.25" weight="5.0" />
|
|
318
|
+
|
|
319
|
+
<!-- <operator id="ReScaler" spec="BactrianScaleOperator" parameter="@ReValues"
|
|
320
|
+
weight="5.0" />
|
|
321
|
+
<operator id="ReScalerAll" spec="BactrianScaleOperator" parameter="@ReValues"
|
|
322
|
+
weight="5.0" scaleAll="true" /> -->
|
|
323
|
+
|
|
324
|
+
<operator id="ReScaler" spec="SmartScaleOperator" parameter="@ReValues"
|
|
325
|
+
weight="5.0" />
|
|
326
|
+
|
|
327
|
+
<operator id="samplingProportionScaler" spec="SmartScaleOperator" parameter="@samplingProportionValues"
|
|
328
|
+
weight="3.0" />
|
|
329
|
+
|
|
330
|
+
<operator id="migrationRateScaler" spec="SmartScaleOperator" parameter="@migrationRateValues"
|
|
331
|
+
weight="5.0" />
|
|
332
|
+
<operator id="fScaler" spec="BactrianScaleOperator" parameter="@f"
|
|
333
|
+
weight="5.0"/>
|
|
334
|
+
|
|
335
|
+
<!-- Substitution model operators -->
|
|
336
|
+
<operator id="gammaShapeScaler" spec="BactrianScaleOperator" parameter="@gammaShape"
|
|
337
|
+
weight="1.0"/>
|
|
338
|
+
|
|
339
|
+
<operator id="KappaScaler" spec="BactrianScaleOperator" parameter="@kappa"
|
|
340
|
+
weight="1.0"/>
|
|
341
|
+
|
|
342
|
+
|
|
343
|
+
<!-- Loggers -->
|
|
344
|
+
<logger id="tracelog" spec="Logger" fileName="$(filebase).log" logEvery="1000" model="@posterior" sanitiseHeaders="true" sort="smart">
|
|
188
345
|
<log idref="posterior"/>
|
|
189
346
|
<log idref="likelihood"/>
|
|
190
347
|
<log idref="prior"/>
|
|
191
|
-
<log idref="
|
|
192
|
-
<log idref="
|
|
193
|
-
<log idref="
|
|
194
|
-
<log
|
|
195
|
-
<log spec="TypedTreeStatsLogger" includeRootEdge="true"
|
|
348
|
+
<log idref="treeLikelihood"/>
|
|
349
|
+
<log idref="gammaShape"/>
|
|
350
|
+
<log idref="kappa"/>
|
|
351
|
+
<log id="TreeHeight" spec="TreeStatLogger" tree="@treeBDMMPrime"/>
|
|
352
|
+
<log id="typedTreeStats" spec="TypedTreeStatsLogger" includeRootEdge="true"
|
|
353
|
+
typeLabel="type" typeSet="@typeSetBDMMPrime">
|
|
196
354
|
<typedTree id="typeMappedTree" spec="TypeMappedTree" bdmmDistrib="@BDMMPrime" startTypePriorProbs="@typeFrequencies" mapOnInit="false" remapOnLog="true" typeLabel="type" typeTraitSet="@typeTraitSet" untypedTree="@treeBDMMPrime"/>
|
|
197
355
|
</log>
|
|
356
|
+
<log idref="originBDMMPrime"/>
|
|
357
|
+
<log idref="ReSV"/>
|
|
358
|
+
<log idref="becomeUninfectiousRateSV"/>
|
|
359
|
+
<log idref="samplingProportionSV"/>
|
|
360
|
+
<log idref="migrationRateSM"/>
|
|
361
|
+
<log idref="migrationRateValues"/>
|
|
362
|
+
<log idref="f"/>
|
|
363
|
+
<log idref="removalProbSV"/>
|
|
364
|
+
<log idref="typeFrequencies"/>
|
|
198
365
|
</logger>
|
|
199
366
|
|
|
200
|
-
<logger
|
|
201
|
-
<log idref="
|
|
367
|
+
<logger id="screenlog" spec="Logger" logEvery="1000">
|
|
368
|
+
<log idref="posterior"/>
|
|
369
|
+
<log arg="@posterior" id="ESS_posterior" spec="util.ESS"/>
|
|
370
|
+
<log idref="likelihood"/>
|
|
371
|
+
<log arg="@likelihood" id="ESS_likelihood" spec="util.ESS"/>
|
|
372
|
+
<log idref="prior"/>
|
|
373
|
+
<log arg="@prior" id="ESS_prior" spec="util.ESS"/>
|
|
374
|
+
<log idref="BDMMPrime"/>
|
|
375
|
+
<log arg="@BDMMPrime" id="ESS_BDMMPrime" spec="util.ESS"/>
|
|
202
376
|
</logger>
|
|
203
377
|
|
|
204
|
-
<logger
|
|
205
|
-
<log spec="
|
|
378
|
+
<logger id="treelog" spec="Logger" fileName="$(filebase).trees" logEvery="1000" mode="tree">
|
|
379
|
+
<log id="TreeWithMetaDataLogger" spec="TreeWithMetaDataLogger" tree="@treeBDMMPrime"/>
|
|
206
380
|
</logger>
|
|
207
|
-
|
|
208
|
-
<logger spec="OptionalLogger" enableLogger="true"
|
|
209
|
-
|
|
381
|
+
|
|
382
|
+
<logger id="typedTreeLogger" spec="OptionalLogger" enableLogger="true"
|
|
383
|
+
fileName="$(filebase).typed.trees" logEvery="10000" mode="tree">
|
|
384
|
+
<log idref="typeMappedTree"/>
|
|
385
|
+
</logger>
|
|
386
|
+
<logger id="nodeTypedTreeLogger.t:eubdmm_alignment" spec="OptionalLogger" enableLogger="true"
|
|
387
|
+
fileName="$(filebase).typed.node.trees" logEvery="10000" mode="tree">
|
|
388
|
+
<log id="nodeTypedTree" spec="TypedNodeTreeLogger" typedTree="@typeMappedTree"/>
|
|
210
389
|
</logger>
|
|
390
|
+
<logger id="trajLogger" spec="OptionalLogger" enableLogger="true" fileName="$(filebase).traj" logEvery="10000">
|
|
391
|
+
<log id="typedTraj" spec="SampledTrajectory" bdmmDistrib="@BDMMPrime"
|
|
392
|
+
startTypePriorProbs="@typeFrequencies" typeMappedTree="@typeMappedTree" useTauLeaping="true"/>
|
|
393
|
+
</logger>
|
|
394
|
+
|
|
395
|
+
<operatorschedule id="OperatorSchedule" spec="OperatorSchedule"/>
|
|
211
396
|
|
|
212
397
|
</run>
|
|
213
398
|
|