bella-companion 0.0.25__py3-none-any.whl → 0.0.27__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 +7 -2
- {bella_companion-0.0.25.dist-info → bella_companion-0.0.27.dist-info}/METADATA +1 -1
- {bella_companion-0.0.25.dist-info → bella_companion-0.0.27.dist-info}/RECORD +9 -9
- {bella_companion-0.0.25.dist-info → bella_companion-0.0.27.dist-info}/WHEEL +0 -0
- {bella_companion-0.0.25.dist-info → bella_companion-0.0.27.dist-info}/entry_points.txt +0 -0
|
@@ -1,200 +1,449 @@
|
|
|
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="
|
|
36
|
-
<map name="
|
|
37
|
-
<map name="
|
|
38
|
-
<map name="
|
|
39
|
-
<map name="prior">beast.base.inference.distribution.Prior</map>
|
|
40
|
-
|
|
41
|
-
|
|
42
|
-
|
|
43
|
-
|
|
44
|
-
|
|
45
|
-
<alignment id="alignment" spec="AlignmentFromFasta"
|
|
46
|
-
|
|
47
|
-
|
|
48
|
-
<
|
|
49
|
-
|
|
50
|
-
|
|
51
|
-
|
|
52
|
-
|
|
53
|
-
|
|
54
|
-
|
|
55
|
-
|
|
2
|
+
<beast namespace="beast.base.inference
|
|
3
|
+
:beast.base.inference.parameter
|
|
4
|
+
:beast.base.evolution.branchratemodel
|
|
5
|
+
:beast.base.evolution.substitutionmodel
|
|
6
|
+
:beast.base.evolution.sitemodel
|
|
7
|
+
:beast.base.evolution.tree.coalescent
|
|
8
|
+
:beast.base.evolution.tree
|
|
9
|
+
:beast.base.evolution.likelihood
|
|
10
|
+
:beast.base.evolution.alignment
|
|
11
|
+
:beast.base.evolution.operator
|
|
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
|
+
:bdmmprime.util.operators
|
|
18
|
+
:bdmmprime.util.priors
|
|
19
|
+
:bdmmprime.distribution
|
|
20
|
+
:bdmmprime.parameterization
|
|
21
|
+
:bdmmprime.mapping
|
|
22
|
+
:bdmmprime.util
|
|
23
|
+
:bdmmprime.trajectories
|
|
24
|
+
:feast
|
|
25
|
+
:feast.parameter
|
|
26
|
+
:feast.fileio
|
|
27
|
+
:feast.function
|
|
28
|
+
:feast.expressions
|
|
29
|
+
:glmprior.util"
|
|
30
|
+
required="BEAST.base v2.7.7:BDMM-Prime v2.2.2" version="2.7">
|
|
31
|
+
|
|
32
|
+
|
|
33
|
+
<!-- Map elements -->
|
|
34
|
+
<map name="Uniform" >beast.base.inference.distribution.Uniform</map>
|
|
35
|
+
<map name="Exponential" >beast.base.inference.distribution.Exponential</map>
|
|
36
|
+
<map name="LogNormal" >beast.base.inference.distribution.LogNormalDistributionModel</map>
|
|
37
|
+
<map name="Normal" >beast.base.inference.distribution.Normal</map>
|
|
38
|
+
<map name="Poisson" >beast.base.inference.distribution.Poisson</map>
|
|
39
|
+
<map name="prior" >beast.base.inference.distribution.Prior</map>
|
|
40
|
+
|
|
41
|
+
<!-- Adapted from https://github.com/SarahNadeau/cov-europe-bdmm/blob/master/analyses/with_hubei_migration_decrease/with_hubei_migration_decrease_scrubbed.xml -->
|
|
42
|
+
<!-- Sequence data loaded from GISAID with the gisaid EPI ISLs -->
|
|
43
|
+
|
|
44
|
+
<!-- Alignment -->
|
|
45
|
+
<alignment id="alignment" spec= "AlignmentFromFasta"
|
|
46
|
+
fileName = "$(aligned_fasta)"/>
|
|
47
|
+
|
|
48
|
+
<typeSet id="typeSetBDMMPrime" spec="TypeSet">
|
|
49
|
+
<typeTraitSet id="typeTraitSet" spec="TraitSetFromTaxonSet"
|
|
50
|
+
traitname="type"
|
|
51
|
+
delimiter="|"
|
|
52
|
+
takeGroup="1">
|
|
53
|
+
<taxa spec="TaxonSet" alignment="@alignment"/>
|
|
54
|
+
</typeTraitSet>
|
|
55
|
+
</typeSet>
|
|
56
|
+
|
|
57
|
+
<trait id="dateTrait" spec="TraitSetFromTaxonSet"
|
|
58
|
+
traitname="date"
|
|
59
|
+
dateFormat="yyyy-M-dd"
|
|
60
|
+
delimiter= "|"
|
|
61
|
+
everythingAfterLast="true">
|
|
62
|
+
<taxa spec="TaxonSet" alignment="@alignment"/>
|
|
63
|
+
</trait>
|
|
64
|
+
|
|
65
|
+
<!-- BDMMPrime Population Model -->
|
|
66
|
+
|
|
67
|
+
<processLength id="originBDMMPrime" spec="RealParameter">10.0</processLength>
|
|
68
|
+
|
|
69
|
+
<Re id="ReSV" spec="SkylineVectorParameter"
|
|
70
|
+
timesAreAges="true" processLength="@originBDMMPrime" typeSet="@typeSetBDMMPrime">
|
|
71
|
+
<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>
|
|
72
|
+
<changeTimes id="ReRateChangeTimes" spec="RealParameter">0.123</changeTimes>
|
|
56
73
|
</Re>
|
|
57
74
|
|
|
58
|
-
<samplingProportion id="
|
|
59
|
-
|
|
75
|
+
<samplingProportion id="samplingProportionSV" spec="SkylineVectorParameter"
|
|
76
|
+
timesAreAges="true" processLength="@originBDMMPrime" typeSet="@typeSetBDMMPrime">
|
|
77
|
+
<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>
|
|
78
|
+
|
|
79
|
+
|
|
80
|
+
<changeTimes id="samplingProportionChangeTimes" spec="RealParameter">0.123 0.205</changeTimes>
|
|
60
81
|
</samplingProportion>
|
|
61
82
|
|
|
62
|
-
<becomeUninfectiousRate id="becomeUninfectiousRateSP" spec="SkylineVectorParameter" typeSet="@typeSet" skylineValues="36.5"/>
|
|
63
83
|
|
|
64
|
-
<
|
|
65
|
-
|
|
66
|
-
|
|
67
|
-
|
|
84
|
+
<becomeUninfectiousRate id="becomeUninfectiousRateSV" spec="SkylineVectorParameter"
|
|
85
|
+
typeSet="@typeSetBDMMPrime">
|
|
86
|
+
<skylineValues id="becomeUninfectiousValues" spec="RealParameter">36.5</skylineValues>
|
|
87
|
+
</becomeUninfectiousRate>
|
|
88
|
+
|
|
89
|
+
|
|
90
|
+
<migrationRate id="migrationRateSM" spec="SkylineMatrixParameter"
|
|
91
|
+
timesAreAges="true" processLength="@originBDMMPrime" typeSet="@typeSetBDMMPrime">
|
|
92
|
+
|
|
93
|
+
<skylineValues id="migrationRateValues" spec="glmprior.util.GLMLogLinear" transform="false">
|
|
94
|
+
<plate var="predictor" range="$(GLMpredictor_files)">
|
|
95
|
+
<predictor spec="RealParameterFromXSV" id="$(predictor)" fileName="$(predictor)"/>
|
|
68
96
|
</plate>
|
|
69
|
-
|
|
70
|
-
<
|
|
71
|
-
<
|
|
97
|
+
|
|
98
|
+
<parameter id="coefficientsGLM" spec="RealParameter" name="coefficients" value="0.1"/>
|
|
99
|
+
<parameter id="scaleFactorGLM" spec="RealParameter" name="scaleFactor" value="1"/>
|
|
100
|
+
<parameter id="indicatorGLM" spec="BooleanParameter" name="indicators" value="1"/>
|
|
101
|
+
<!-- <parameter id="errorTermGLM" spec="RealParameter" name="error" value="0.01" dimension="30"/> -->
|
|
102
|
+
|
|
72
103
|
</skylineValues>
|
|
73
|
-
|
|
104
|
+
|
|
105
|
+
<changeTimes id="migrationRateChangeTimes" spec="RealParameterFromXSV"
|
|
106
|
+
fileName="$(GLMpredictor_changetimes)"/>
|
|
74
107
|
</migrationRate>
|
|
75
108
|
|
|
76
|
-
<removalProb id="removalProbSP" spec="SkylineVectorParameter" skylineValues="1.0" processLength="@processLength" typeSet="@typeSet"/>
|
|
77
109
|
|
|
78
|
-
<startTypePriorProbs id="typeFrequencies" spec="RealParameter" value="1.0 0.0 0.0 0.0 0.0"/>
|
|
79
110
|
|
|
80
|
-
<
|
|
111
|
+
<removalProb id="removalProbSV" spec="SkylineVectorParameter" skylineValues="1.0"
|
|
112
|
+
processLength="@originBDMMPrime" typeSet="@typeSetBDMMPrime"/>
|
|
113
|
+
|
|
114
|
+
<startTypePriorProbs id="typeFrequencies" spec="RealParameter"
|
|
115
|
+
estimate="false" lower="0.0" upper="1.0">1.0 0.0 0.0 0.0 0.0</startTypePriorProbs>
|
|
116
|
+
|
|
117
|
+
|
|
118
|
+
<!-- Site Model -->
|
|
119
|
+
<siteModel id="siteModel" spec="SiteModel" gammaCategoryCount="4">
|
|
81
120
|
<parameter id="gammaShape" name="shape" value="1.0"/>
|
|
121
|
+
<parameter name="mutationRate" estimate="false" value="1.0"/>
|
|
122
|
+
<parameter name="proportionInvariant" estimate="false" value="0.0"/>
|
|
123
|
+
|
|
82
124
|
<substModel spec="HKY">
|
|
83
125
|
<kappa id="kappa" spec="RealParameter" value="2.0"/>
|
|
84
|
-
<frequencies spec="Frequencies" data="@alignment"/>
|
|
126
|
+
<frequencies id="empiricalFreqs" spec="Frequencies" data="@alignment"/>
|
|
85
127
|
</substModel>
|
|
86
128
|
</siteModel>
|
|
87
129
|
|
|
88
|
-
<branchRateModel id="branchRateModel" spec="StrictClockModel" clock.rate="8.0E-4"/>
|
|
89
130
|
|
|
90
|
-
|
|
131
|
+
<!-- Branch rate Model -->
|
|
132
|
+
<branchRateModel id="branchRateModel" spec="StrictClockModel">
|
|
133
|
+
<parameter id="clockRate" name="clock.rate" estimate="false">8.0E-4</parameter>
|
|
134
|
+
</branchRateModel>
|
|
135
|
+
|
|
91
136
|
|
|
92
|
-
|
|
137
|
+
<!-- Population Model for Tree initialiser -->
|
|
138
|
+
<popFunc id="popFunc" spec="ConstantPopulation">
|
|
139
|
+
<popSize spec="RealParameter" value="0.1"/>
|
|
140
|
+
</popFunc>
|
|
93
141
|
|
|
94
|
-
|
|
95
|
-
|
|
96
|
-
|
|
97
|
-
|
|
98
|
-
|
|
142
|
+
|
|
143
|
+
<!-- MCMC -->
|
|
144
|
+
<run id="mcmc" spec="MCMC" chainLength="10000000" numInitializationAttempts="100">
|
|
145
|
+
|
|
146
|
+
<!-- State -->
|
|
147
|
+
<state id="state" spec="State">
|
|
148
|
+
|
|
149
|
+
<stateNode id="treeBDMMPrime" spec="RandomTree"
|
|
150
|
+
taxa="@alignment"
|
|
151
|
+
populationModel="@popFunc"
|
|
152
|
+
trait="@dateTrait">
|
|
153
|
+
</stateNode>
|
|
154
|
+
|
|
155
|
+
<stateNode idref="originBDMMPrime"/>
|
|
156
|
+
<stateNode idref="ReValues"/>
|
|
157
|
+
<stateNode idref="samplingProportionValues"/>
|
|
99
158
|
<stateNode idref="coefficientsGLM"/>
|
|
100
159
|
<stateNode idref="scaleFactorGLM"/>
|
|
160
|
+
<stateNode idref="indicatorGLM"/>
|
|
161
|
+
|
|
101
162
|
<stateNode idref="kappa"/>
|
|
102
163
|
<stateNode idref="gammaShape"/>
|
|
164
|
+
|
|
103
165
|
</state>
|
|
104
166
|
|
|
167
|
+
|
|
168
|
+
<!-- Posterior Distribution -->
|
|
105
169
|
<distribution id="posterior" spec="CompoundDistribution">
|
|
170
|
+
|
|
171
|
+
<!-- Tree Likelihood -->
|
|
106
172
|
<distribution id="likelihood" spec="CompoundDistribution" useThreads="true">
|
|
107
|
-
|
|
173
|
+
|
|
174
|
+
<distribution id="treeLikelihood" spec="ThreadedTreeLikelihood"
|
|
175
|
+
data="@alignment"
|
|
176
|
+
tree="@treeBDMMPrime"
|
|
177
|
+
siteModel="@siteModel"
|
|
178
|
+
branchRateModel="@branchRateModel"/>
|
|
179
|
+
|
|
108
180
|
</distribution>
|
|
109
181
|
|
|
182
|
+
<!-- Prior distribution -->
|
|
110
183
|
<distribution id="prior" spec="CompoundDistribution">
|
|
111
|
-
|
|
112
|
-
|
|
184
|
+
|
|
185
|
+
<!-- BDMM-Prime tree prior distribution -->
|
|
186
|
+
<distribution id="BDMMPrime" spec="BirthDeathMigrationDistribution"
|
|
187
|
+
tree="@treeBDMMPrime"
|
|
188
|
+
conditionOnSurvival="true"
|
|
189
|
+
finalSampleOffset="0.0"
|
|
190
|
+
typeTraitSet="@typeTraitSet"
|
|
191
|
+
startTypePriorProbs="@typeFrequencies">
|
|
192
|
+
|
|
193
|
+
<!-- Parameterization BDMM-Prime -->
|
|
194
|
+
<parameterization id="EpiBDMMPrimeParameterization" spec="EpiParameterization"
|
|
195
|
+
processLength="@originBDMMPrime"
|
|
196
|
+
Re="@ReSV"
|
|
197
|
+
samplingProportion="@samplingProportionSV"
|
|
198
|
+
becomeUninfectiousRate="@becomeUninfectiousRateSV"
|
|
199
|
+
migrationRate="@migrationRateSM"
|
|
200
|
+
removalProb="@removalProbSV"
|
|
201
|
+
typeSet="@typeSetBDMMPrime"/>
|
|
202
|
+
|
|
113
203
|
</distribution>
|
|
114
204
|
|
|
115
|
-
|
|
205
|
+
|
|
206
|
+
<!-- Priors on parameters -->
|
|
207
|
+
<!-- BDMMPrime parameters priors -->
|
|
208
|
+
|
|
209
|
+
<prior id="originBDMMPrimePrior" name="distribution" x="@originBDMMPrime">
|
|
116
210
|
<LogNormal name="distr" M="-1.0" S="0.2"/>
|
|
117
211
|
</prior>
|
|
118
|
-
|
|
212
|
+
|
|
213
|
+
<!-- Reproductive number Prior -->
|
|
214
|
+
<distribution id="RePrior" spec="SmartZeroExcludingPrior" x="@ReValues">
|
|
119
215
|
<LogNormal name="distr" M="0.8" S="0.5"/>
|
|
120
216
|
</distribution>
|
|
121
|
-
|
|
217
|
+
|
|
218
|
+
<!-- Prior for sampling proportion for China -->
|
|
219
|
+
<distribution id="spChinaE2Prior" spec="SmartZeroExcludingPrior" x="@samplingProportionValues"
|
|
220
|
+
classesToExclude="1.2E-5 1.3E-5 1.4E-5 1.5E-5">
|
|
122
221
|
<Uniform name="distr" lower="0" upper="0.15"/>
|
|
123
222
|
</distribution>
|
|
124
|
-
|
|
223
|
+
<!-- Prior for sampling proportion for France -->
|
|
224
|
+
<distribution id="spFranceE1Prior" spec="SmartZeroExcludingPrior" x="@samplingProportionValues"
|
|
225
|
+
classesToExclude="1.1E-5 1.3E-5 1.4E-5 1.5E-5">
|
|
125
226
|
<Uniform name="distr" lower="0" upper="0.093"/>
|
|
126
227
|
</distribution>
|
|
127
|
-
|
|
228
|
+
<!-- Prior for sampling proportion for Germany Epoch -->
|
|
229
|
+
<distribution id="spGermanyE1Prior" spec="SmartZeroExcludingPrior" x="@samplingProportionValues"
|
|
230
|
+
classesToExclude="1.1E-5 1.2E-5 1.4E-5 1.5E-5">
|
|
128
231
|
<Uniform name="distr" lower="0" upper="0.10"/>
|
|
129
232
|
</distribution>
|
|
130
|
-
|
|
233
|
+
<!-- Prior for sampling proportion for Italy -->
|
|
234
|
+
<distribution id="spItalyE1Prior" spec="SmartZeroExcludingPrior" x="@samplingProportionValues"
|
|
235
|
+
classesToExclude="1.1E-5 1.2E-5 1.3E-5 1.5E-5">
|
|
131
236
|
<Uniform name="distr" lower="0" upper="0.005"/>
|
|
132
237
|
</distribution>
|
|
133
|
-
|
|
238
|
+
<!-- Prior for sampling proportion for OtherEuropean -->
|
|
239
|
+
<distribution id="spOEE1Prior" spec="SmartZeroExcludingPrior" x="@samplingProportionValues"
|
|
240
|
+
classesToExclude="1.1E-5 1.2E-5 1.3E-5 1.4E-5">
|
|
134
241
|
<Uniform name="distr" lower="0" upper="0.057"/>
|
|
135
242
|
</distribution>
|
|
136
|
-
|
|
137
|
-
|
|
243
|
+
|
|
244
|
+
|
|
245
|
+
<!-- Migration rate GLM Prior -->
|
|
246
|
+
<!-- <prior id="migRatePrior" name="distribution" x="@migrationRateValues">
|
|
247
|
+
<Uniform id="Uniform.migRate" name="distr" lower="0.0" upper="50.0"/>
|
|
248
|
+
</prior> -->
|
|
249
|
+
|
|
250
|
+
<!-- <prior id="migRatePrior" name="distribution" x="@migrationRateValues">
|
|
251
|
+
<LogNormal name="distr" M="0" S="1.5"/>
|
|
252
|
+
</prior> -->
|
|
253
|
+
|
|
254
|
+
<prior id="coefficientsPrior" name="distribution" x="@coefficientsGLM">
|
|
255
|
+
<Normal name="distr">
|
|
256
|
+
<parameter spec="RealParameter" estimate="false" name="mean">0.0</parameter>
|
|
257
|
+
<parameter spec="RealParameter" estimate="false" name="sigma">1.0</parameter>
|
|
258
|
+
</Normal>
|
|
138
259
|
</prior>
|
|
139
|
-
|
|
140
|
-
|
|
260
|
+
|
|
261
|
+
<prior id="scaleFactorPrior" name="distribution" x="@scaleFactorGLM">
|
|
262
|
+
<LogNormal name="distr">
|
|
263
|
+
<parameter spec="RealParameter" estimate="false" name="M">0.0</parameter>
|
|
264
|
+
<parameter spec="RealParameter" estimate="false" name="S">0.8</parameter>
|
|
265
|
+
</LogNormal>
|
|
141
266
|
</prior>
|
|
142
|
-
|
|
143
|
-
|
|
267
|
+
|
|
268
|
+
<!-- <prior id="nonZeroIndicatorGLMPrior" name="distribution">
|
|
269
|
+
<x id="indicatorSumGLM" spec="Sum" arg="@indicatorGLM"/>
|
|
270
|
+
<Poisson name="distr" lambda='0.693'/>
|
|
271
|
+
</prior> -->
|
|
272
|
+
|
|
273
|
+
<!-- Substitution model parameter priors -->
|
|
274
|
+
<prior id="GammaShapePrior" name="distribution" x="@gammaShape">
|
|
275
|
+
<Exponential name="distr">
|
|
276
|
+
<parameter spec="RealParameter" estimate="false" name="mean">0.5</parameter>
|
|
277
|
+
</Exponential>
|
|
144
278
|
</prior>
|
|
145
|
-
|
|
146
|
-
|
|
279
|
+
|
|
280
|
+
<prior id="KappaPrior" name="distribution" x="@kappa">
|
|
281
|
+
<LogNormal name="distr">
|
|
282
|
+
<parameter spec="RealParameter" estimate="false" name="M">1.0</parameter>
|
|
283
|
+
<parameter spec="RealParameter" estimate="false" name="S">1.25</parameter>
|
|
284
|
+
</LogNormal>
|
|
147
285
|
</prior>
|
|
286
|
+
|
|
148
287
|
</distribution>
|
|
149
288
|
</distribution>
|
|
150
289
|
|
|
151
|
-
|
|
152
|
-
|
|
153
|
-
|
|
154
|
-
<operator spec="
|
|
155
|
-
|
|
156
|
-
|
|
157
|
-
<operator spec="
|
|
158
|
-
|
|
159
|
-
|
|
160
|
-
<operator
|
|
161
|
-
|
|
162
|
-
|
|
163
|
-
<operator
|
|
164
|
-
|
|
165
|
-
|
|
290
|
+
|
|
291
|
+
<!-- Operators -->
|
|
292
|
+
<!-- Tree operators -->
|
|
293
|
+
<operator id="BDMMPrimeTreeRootScaler" spec="BactrianScaleOperator" tree="@treeBDMMPrime"
|
|
294
|
+
rootOnly="true" scaleFactor="0.25" weight="5.0"/>
|
|
295
|
+
|
|
296
|
+
<operator id="BDMMPrimeUniformOperator" spec="BactrianNodeOperator" tree="@treeBDMMPrime"
|
|
297
|
+
weight="30.0"/>
|
|
298
|
+
|
|
299
|
+
<operator id="BDMMPrimeSubtreeSlideScaler" spec="BactrianSubtreeSlide" tree="@treeBDMMPrime"
|
|
300
|
+
weight="15.0"/>
|
|
301
|
+
|
|
302
|
+
<operator id="BDMMPrimeNarrow" spec="Exchange" tree="@treeBDMMPrime"
|
|
303
|
+
weight="15.0"/>
|
|
304
|
+
|
|
305
|
+
<operator id="BDMMPrimeWide" spec="Exchange" tree="@treeBDMMPrime"
|
|
306
|
+
isNarrow="false" weight="3.0"/>
|
|
307
|
+
|
|
308
|
+
<operator id="BDMMPrimeWilsonBalding" spec="WilsonBalding" tree="@treeBDMMPrime"
|
|
309
|
+
weight="3.0"/>
|
|
310
|
+
|
|
311
|
+
<operator id="BDMMPrimeBICEPSEpochTop" spec="EpochFlexOperator" tree="@treeBDMMPrime"
|
|
312
|
+
scaleFactor="0.1" weight="2.0"/>
|
|
313
|
+
|
|
314
|
+
<operator id="BDMMPrimeBICEPSEpochAll" spec="EpochFlexOperator" tree="@treeBDMMPrime"
|
|
315
|
+
fromOldestTipOnly="false" scaleFactor="0.1" weight="2.0"/>
|
|
316
|
+
|
|
317
|
+
<operator id="BDMMPrimeBICEPSTreeFlex" spec="TreeStretchOperator" tree="@treeBDMMPrime"
|
|
318
|
+
scaleFactor="0.01" weight="2.0"/>
|
|
319
|
+
|
|
320
|
+
|
|
321
|
+
<!-- Population model operators -->
|
|
322
|
+
|
|
323
|
+
<operator id="originBDMMPrimeScaler" spec="BactrianScaleOperator" parameter="@originBDMMPrime"
|
|
324
|
+
scaleFactor="0.25" weight="3.0" />
|
|
325
|
+
|
|
326
|
+
<operator id="ReScaler" spec="SmartScaleOperator" parameter="@ReValues" scaleFactor="0.25"
|
|
327
|
+
weight="10.0" />
|
|
328
|
+
<!-- <operator id="ReScalerAll" spec="BactrianScaleOperator" parameter="@ReValues" scaleFactor="0.25"
|
|
329
|
+
weight="3.0" scaleAll="true" /> -->
|
|
330
|
+
|
|
331
|
+
<operator id="samplingProportionScaler" spec="SmartScaleOperator" parameter="@samplingProportionValues"
|
|
332
|
+
weight="3.0" />
|
|
333
|
+
|
|
334
|
+
<!-- <operator id="scaleFactorGLMScaler" spec="BactrianScaleOperator" parameter="@scaleFactorGLM" scaleFactor="0.25" weight="3.0" /> -->
|
|
335
|
+
|
|
336
|
+
<operator id="coefficientsGLMScaler" spec="BactrianRandomWalkOperator" parameter="@coefficientsGLM"
|
|
337
|
+
weight="5.0" scaleFactor="0.1"/>
|
|
338
|
+
|
|
339
|
+
<operator id="AVNMOperator" spec="beast.base.evolution.operator.kernel.AdaptableVarianceMultivariateNormalOperator" weight="5.0"
|
|
340
|
+
coefficient="1.0"
|
|
341
|
+
scaleFactor="0.5"
|
|
342
|
+
beta="0.05"
|
|
343
|
+
initial="800"
|
|
344
|
+
burnin="400"
|
|
345
|
+
every="1">
|
|
346
|
+
<transformations spec="beast.base.inference.operator.kernel.Transform$NoTransform">
|
|
347
|
+
<f idref="coefficientsGLM"/>
|
|
348
|
+
</transformations>
|
|
349
|
+
</operator>
|
|
350
|
+
|
|
351
|
+
<!-- <operator id="indicatorParamGLMFlip" spec="BitFlipOperator" parameter="@indicatorGLM"
|
|
352
|
+
weight="10.0" uniform="false" />
|
|
353
|
+
|
|
354
|
+
<operator id="indicatorSwap" spec="glmprior.operator.ExtendedSwapOperator" boolparameter="@indicatorGLM" weight="5.0" />
|
|
355
|
+
|
|
356
|
+
<operator id="indicatorParamGLMPick" spec="glmprior.operator.PickIndicatorOperator"
|
|
357
|
+
parameter="@indicatorGLM" weight="3.0" />
|
|
358
|
+
|
|
359
|
+
<operator id="BSSVSGLMoperator" spec="glmprior.operator.BitFlipBSSVSOperator" indicator="@indicatorGLM" mu="@scaleFactorGLM" weight="10"/>
|
|
360
|
+
|
|
361
|
+
<operator id="multiSwapOperator" spec="glmprior.operator.MultiSwapOperator" boolparameter="@indicatorGLM" parameter="@coefficientsGLM" weight="10.0" />
|
|
362
|
+
|
|
363
|
+
<operator id="JointSwapGLM1" spec="JointOperator" weight="3.0">
|
|
364
|
+
<operator idref="indicatorSwap"/>
|
|
365
|
+
<operator idref="AVNMOperator"/>
|
|
166
366
|
</operator>
|
|
167
|
-
<operator id="scaleFactorGLMScaler" spec="BactrianScaleOperator" parameter="@scaleFactorGLM" scaleFactor="0.25" weight="3.0" />
|
|
168
|
-
<operator spec="BactrianScaleOperator" parameter="@gammaShape" weight="1.0"/>
|
|
169
|
-
<operator spec="BactrianScaleOperator" parameter="@kappa" weight="1.0"/>
|
|
170
|
-
<operatorschedule spec="OperatorSchedule"/>
|
|
171
367
|
|
|
172
|
-
<
|
|
368
|
+
<operator id="JointSwapGLM2" spec="JointOperator" weight="3.0">
|
|
369
|
+
<operator idref="coefficientsGLMScaler"/>
|
|
370
|
+
<operator idref="scaleFactorGLMScaler"/>
|
|
371
|
+
</operator>
|
|
372
|
+
-->
|
|
373
|
+
<!-- <operator id="errorTermGLMScaler" spec="BactrianScaleOperator" parameter="@errorTermGLM"
|
|
374
|
+
weight="5.0" upper="20.0"/> -->
|
|
375
|
+
|
|
376
|
+
|
|
377
|
+
<!-- Substitution model operators -->
|
|
378
|
+
<operator id="gammaShapeScaler" spec="BactrianScaleOperator" parameter="@gammaShape"
|
|
379
|
+
weight="1.0"/>
|
|
380
|
+
|
|
381
|
+
<operator id="KappaScaler" spec="BactrianScaleOperator" parameter="@kappa"
|
|
382
|
+
weight="1.0"/>
|
|
383
|
+
|
|
384
|
+
|
|
385
|
+
<!-- Loggers -->
|
|
386
|
+
<logger id="tracelog" spec="Logger" fileName="$(filebase).log" logEvery="1000" model="@posterior" sanitiseHeaders="true" sort="smart">
|
|
173
387
|
<log idref="posterior"/>
|
|
174
388
|
<log idref="likelihood"/>
|
|
175
389
|
<log idref="prior"/>
|
|
176
|
-
<log idref="
|
|
177
|
-
<log idref="
|
|
178
|
-
<log idref="
|
|
179
|
-
<log
|
|
390
|
+
<log idref="treeLikelihood"/>
|
|
391
|
+
<log idref="gammaShape"/>
|
|
392
|
+
<log idref="kappa"/>
|
|
393
|
+
<log id="TreeHeight" spec="TreeStatLogger" tree="@treeBDMMPrime"/>
|
|
394
|
+
<log id="typedTreeStats" spec="TypedTreeStatsLogger" includeRootEdge="true"
|
|
395
|
+
typeLabel="type" typeSet="@typeSetBDMMPrime">
|
|
396
|
+
<typedTree id="typeMappedTree" spec="TypeMappedTree" bdmmDistrib="@BDMMPrime" startTypePriorProbs="@typeFrequencies" mapOnInit="false" remapOnLog="true" typeLabel="type" typeTraitSet="@typeTraitSet" untypedTree="@treeBDMMPrime"/>
|
|
397
|
+
</log>
|
|
398
|
+
<log idref="originBDMMPrime"/>
|
|
399
|
+
<log idref="ReSV"/>
|
|
400
|
+
<log idref="becomeUninfectiousRateSV"/>
|
|
401
|
+
<log idref="samplingProportionSV"/>
|
|
402
|
+
<log idref="migrationRateSM"/>
|
|
180
403
|
<log idref="coefficientsGLM"/>
|
|
181
404
|
<log idref="scaleFactorGLM"/>
|
|
182
|
-
<log
|
|
183
|
-
|
|
405
|
+
<log idref="indicatorGLM"/>
|
|
406
|
+
<log id="coefbyindGLM" spec="ExpCalculator" value="coefficientsGLM*indicatorGLM" useCaching="false">
|
|
407
|
+
<arg idref="coefficientsGLM"/>
|
|
408
|
+
<arg idref="indicatorGLM"/>
|
|
184
409
|
</log>
|
|
410
|
+
<plate var="predictor" range="$(GLMpredictor_files)">
|
|
411
|
+
<log idref="$(predictor)"/>
|
|
412
|
+
</plate>
|
|
413
|
+
<!-- <log idref="errorTermGLM"/> -->
|
|
414
|
+
<log idref="removalProbSV"/>
|
|
415
|
+
<log idref="typeFrequencies"/>
|
|
185
416
|
</logger>
|
|
186
417
|
|
|
187
|
-
<logger
|
|
188
|
-
<log idref="
|
|
418
|
+
<logger id="screenlog" spec="Logger" logEvery="1000">
|
|
419
|
+
<log idref="posterior"/>
|
|
420
|
+
<log arg="@posterior" id="ESS_posterior" spec="util.ESS"/>
|
|
421
|
+
<log idref="likelihood"/>
|
|
422
|
+
<log arg="@likelihood" id="ESS_likelihood" spec="util.ESS"/>
|
|
423
|
+
<log idref="prior"/>
|
|
424
|
+
<log arg="@prior" id="ESS_prior" spec="util.ESS"/>
|
|
425
|
+
<log idref="BDMMPrime"/>
|
|
426
|
+
<log arg="@BDMMPrime" id="ESS_BDMMPrime" spec="util.ESS"/>
|
|
189
427
|
</logger>
|
|
190
428
|
|
|
191
|
-
<logger
|
|
192
|
-
<log spec="
|
|
429
|
+
<logger id="treelog" spec="Logger" fileName="$(filebase).trees" logEvery="1000" mode="tree">
|
|
430
|
+
<log id="TreeWithMetaDataLogger" spec="TreeWithMetaDataLogger" tree="@treeBDMMPrime"/>
|
|
193
431
|
</logger>
|
|
194
|
-
|
|
195
|
-
<logger spec="OptionalLogger" enableLogger="true"
|
|
196
|
-
|
|
432
|
+
|
|
433
|
+
<logger id="typedTreeLogger" spec="OptionalLogger" enableLogger="true"
|
|
434
|
+
fileName="$(filebase).typed.trees" logEvery="10000" mode="tree">
|
|
435
|
+
<log idref="typeMappedTree"/>
|
|
436
|
+
</logger>
|
|
437
|
+
<logger id="nodeTypedTreeLogger.t:eubdmm_alignment" spec="OptionalLogger" enableLogger="true"
|
|
438
|
+
fileName="$(filebase).typed.node.trees" logEvery="10000" mode="tree">
|
|
439
|
+
<log id="nodeTypedTree" spec="TypedNodeTreeLogger" typedTree="@typeMappedTree"/>
|
|
197
440
|
</logger>
|
|
441
|
+
<logger id="trajLogger" spec="OptionalLogger" enableLogger="true" fileName="$(filebase).traj" logEvery="10000">
|
|
442
|
+
<log id="typedTraj" spec="SampledTrajectory" bdmmDistrib="@BDMMPrime"
|
|
443
|
+
startTypePriorProbs="@typeFrequencies" typeMappedTree="@typeMappedTree" useTauLeaping="true"/>
|
|
444
|
+
</logger>
|
|
445
|
+
|
|
446
|
+
<operatorschedule id="OperatorSchedule" spec="OperatorSchedule"/>
|
|
198
447
|
|
|
199
448
|
</run>
|
|
200
449
|
|