bella-companion 0.0.47__py3-none-any.whl → 0.1.1__py3-none-any.whl

bella_companion/utils/explain/typeguards.py

@@ -0,0 +1,20 @@
+from typing import Any, TypeGuard
+
+import numpy as np
+from bella_lumiere.typing import Weights
+
+
+def is_list(obj: Any) -> TypeGuard[list[Any]]:
+    return isinstance(obj, list)
+
+
+def is_weights(obj: Any) -> TypeGuard[Weights]:
+    return is_list(obj) and all(isinstance(layer, np.ndarray) for layer in obj)
+
+
+def is_weights_list(weights: Any) -> TypeGuard[list[Weights]]:
+    return is_list(weights) and all(is_weights(w) for w in weights)
+
+
+def is_nested_weights_list(weights: Any) -> TypeGuard[list[list[Weights]]]:
+    return is_list(weights) and all(is_weights_list(w) for w in weights)

{bella_companion-0.0.47.dist-info → bella_companion-0.1.1.dist-info}/METADATA

@@ -1,12 +1,15 @@
-Metadata-Version: 2.4
+Metadata-Version: 2.3
 Name: bella-companion
-Version: 0.0.47
+Version: 0.1.1
 Summary: Bayesian Evolutionary Layered Learning Architectures (BELLA) companion
-Requires-Python: >=3.10
-Description-Content-Type: text/markdown
+Author: gabriele-marino
+Author-email: gabriele-marino <gabmarino.8601@gmail.com>
 Requires-Dist: arviz>=0.22.0
-Requires-Dist: bella-lumiere>=0.0.14
+Requires-Dist: bella-lumiere>=0.1.0
 Requires-Dist: dotenv>=0.9.9
 Requires-Dist: phylogenie>=3.1.4
 Requires-Dist: poetry>=2.2.1
 Requires-Dist: seaborn>=0.13.2
+Requires-Python: >=3.10
+Description-Content-Type: text/markdown
+

bella_companion-0.1.1.dist-info/RECORD

@@ -0,0 +1,64 @@
+bella_companion/.DS_Store,sha256=1lFlJ5EFymdzGAUAaI30vcaaLHt3F1LwpG7xILf9jsM,6148
+bella_companion/BELLA.jar,sha256=4mKDzPfz263rTqA9KS2f4Mc-UWGMuj8G48ZDonKSqVw,67408597
+bella_companion/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+bella_companion/cli.py,sha256=CrtV2CG6G1nKCCVu02gtNlpHNJHLnFtSaru-LFiHKKc,3014
+bella_companion/eucovid/__init__.py,sha256=09Ld_G7fZvCDCgEWsmaOyLVQ8pFc9QHAGnXKJ9Zm2AM,81
+bella_companion/eucovid/beast_configs/GLM.xml,sha256=-SZWxJgJdBthp132ckPpRGooKoYfKKREi6VV3pF22aY,11532
+bella_companion/eucovid/beast_configs/MLP.xml,sha256=C2d-Mcjw44CLFL-pCYgYdbZb1bNR-6OUBtaLWLUt5AU,10539
+bella_companion/eucovid/data/changetimes.csv,sha256=XEepIXKRwcKknJq4-ovKlorvCypGi67B8znuQ5f1LTg,60
+bella_companion/eucovid/data/flights.csv,sha256=tMR9-tuGtC3rVABXUtNdV9irUX-9aCCVb_fwksHTn8U,1068
+bella_companion/eucovid/data/msa.fasta,sha256=sWaVSINSJDDLjmwlpl31JWM5VnKF7yViSpERI69glQc,4839354
+bella_companion/eucovid/run_beast.py,sha256=O9jULAMYIJOtQ_o4EQlG0Xh8s1Z2i27wmnGjYbvNWeM,1331
+bella_companion/features.py,sha256=xj0iwGHhLV_4jxLmG0M-jXWejuKVFJqpB0dK8Skc9n8,283
+bella_companion/platyrrhine/__init__.py,sha256=gySgNPoKaEdhcZppKMsCmS-aAcYw_rDEw0viz72Udaw,259
+bella_companion/platyrrhine/beast_config.xml,sha256=bLzavt-Zd8cmSgQQnlCBfZ9KyM7d1inOGx12VNjPFb4,7056
+bella_companion/platyrrhine/data/change_times.csv,sha256=tedobHbaY7-6M3V6F6HSqGAU9dpeeajR5mu1XGJfp6w,68
+bella_companion/platyrrhine/data/traits.csv,sha256=N_mCnaM9Gy0Qot7-l_mK0vMvr7e_5BZ4bgUHpcoG9UM,3949
+bella_companion/platyrrhine/data/trees.nwk,sha256=lpSiGUG5s3zE--IKPOyWrJ5vpJb_xAPnfxychs3YNIs,488232
+bella_companion/platyrrhine/results.py,sha256=aRx9MkI43f5bhr-16UkShDH-_oMNFzhHxfU2WK1L3BY,7468
+bella_companion/platyrrhine/run_beast.py,sha256=t1gUmn_HEXKqxgGXYCdohGDWeL-QEpW7r3bKXlRfdog,2494
+bella_companion/platyrrhine/summarize.py,sha256=Qc3GE2vUspTqU_QjLdULvK0r7AZ3JTdqsBmd71izt7M,1356
+bella_companion/simulations/__init__.py,sha256=gjO_QB2wp335l7UO-zA9fFQ5vmZauGItBuO-QZFu93E,439
+bella_companion/simulations/beast_configs/epi-multitype/GLM.xml,sha256=VymAYsaZVuB3EPE0DCQSXpKepVq8pPjKaB5yhEpaVkE,4334
+bella_companion/simulations/beast_configs/epi-multitype/MLP.xml,sha256=AoJyjvVcwKtut64d7sB2gnnenx5sXyuTTNCivykUyqU,4077
+bella_companion/simulations/beast_configs/epi-multitype/Nonparametric.xml,sha256=DfCMcLM53b-XP9oIiQ9wJWOVTEP3GnoZaL5BMzLLfIQ,3208
+bella_companion/simulations/beast_configs/epi-skyline/GLM.xml,sha256=7LtJhU7o7aQsnWK7ZIKj-mzmHgCn_RjxvKKbLIlp1Go,3619
+bella_companion/simulations/beast_configs/epi-skyline/MLP.xml,sha256=ptoLM0FLDp54bUSxPm12qIvOYWyBkWR1e9VWVWjgU8E,3426
+bella_companion/simulations/beast_configs/epi-skyline/Nonparametric.xml,sha256=LTiatqeAMiFfdD1leQ5TrGaJhsYwlkbMDpxqhVFOgUQ,2561
+bella_companion/simulations/beast_configs/fbd-2traits/GLM.xml,sha256=Iff0eAShKLeOnlTmNu6044Wljtla3yckqo2LZ7mP6LY,5897
+bella_companion/simulations/beast_configs/fbd-2traits/MLP.xml,sha256=RvF70fcnvoBcJxjyT4nELlFm94atSZpKm9d1aS-cJZ4,5428
+bella_companion/simulations/beast_configs/fbd-2traits/Nonparametric.xml,sha256=ev6ZfMeCwMyK54G_4O8g-tbgoxoD3asxt71x91Cb1QY,3618
+bella_companion/simulations/beast_configs/fbd-no-traits/GLM.xml,sha256=Caq7Gxqx-wYAu3_YkMikVqni3DjBv1Wt00sc1Upo1X8,4743
+bella_companion/simulations/beast_configs/fbd-no-traits/MLP.xml,sha256=JE_cGh2DW5yPAXyXXr0D3DlovmV-aIbze7qrhENHK3w,4275
+bella_companion/simulations/beast_configs/fbd-no-traits/Nonparametric.xml,sha256=m2xaHZwskCFInvVY1QH1fwQqH6p_fgESlPX2BbBeUhk,2822
+bella_companion/simulations/figures/__init__.py,sha256=LeHPeBb4vxJo0pxLZi_9s6iRBKvWC32yhAZBByqK2_8,590
+bella_companion/simulations/figures/epi_multitype.py,sha256=ez-yCP_rs7J0CkfEjt4SFwaiI5T0r3GcLhUXjxc9F3M,2627
+bella_companion/simulations/figures/epi_skyline.py,sha256=SyFOsb8dwLuoN5IOCiQbG-ZCPDS8Wua_LL-cOuLtu_s,1673
+bella_companion/simulations/figures/fbd_2traits.py,sha256=MCzZLmVi-txEIk8Sdm3Tw6oUGMahnifzMuPzwaDw_sY,2726
+bella_companion/simulations/figures/fbd_no_traits.py,sha256=-_c-qnrFAYZ7vd4Ieeb-uYA3VUwURfrcP11BN9PG8JY,1907
+bella_companion/simulations/figures/scenarios.py,sha256=gbMz1TUxxT2RSIq2kQlFioNdgSHk-gQY2OQuf6-7Fww,2817
+bella_companion/simulations/figures/utils.py,sha256=0M5OrxaEuqcj9rR2uAc_O7utQvhEceZGH0sKrGRWaWs,3129
+bella_companion/simulations/generate_data.py,sha256=edfhXrs76hl30i_xizup4z6qgzXZyaDzX_PVC62Ytpw,849
+bella_companion/simulations/generate_figures.py,sha256=r3zGKnLGDpuVzzWVQIumkt1sTMl_4w1iKnVNy-fTJAg,371
+bella_companion/simulations/metrics.py,sha256=kSX7-B4a-LDm1V--S8CvEJDktlJ8qkpo_gY6V7XbOXw,1930
+bella_companion/simulations/run_beast.py,sha256=-L5PHkpPK-NGQoXxmeH7LMIWxtodvL4TOo-KkR2ISVQ,3299
+bella_companion/simulations/scenarios/__init__.py,sha256=3Kl1lKcFpfb3vLX64DmSW4XCF5kXU1ZoHtstFH-ZIzU,876
+bella_companion/simulations/scenarios/common.py,sha256=_ddaSuTvEVdttGkXB4HPc2B7IB1F_GBOCW3cVOPZ-ZM,807
+bella_companion/simulations/scenarios/epi_multitype.py,sha256=hgl2-m9gPm6OntnOZeRfquflbXFuu2Kt4RmKinSslYo,2427
+bella_companion/simulations/scenarios/epi_skyline.py,sha256=iLwsGw-DZK8QOIGhMRUyvFRDEwFrksSx46s2DB6sr0U,2353
+bella_companion/simulations/scenarios/fbd_2traits.py,sha256=-BNJMx-4GGW32IaTDxwVOHs2ZW7z781GJn-tMycygcI,3529
+bella_companion/simulations/scenarios/fbd_no_traits.py,sha256=OsJP2bJVswLgeerXSM0ULKnj5_lxRVHCCoENjs_fWPM,2508
+bella_companion/simulations/scenarios/scenario.py,sha256=oV9yxxD_De5YdCLbAwIFLmiRXqOaoflJjTIkY2yu5M8,608
+bella_companion/simulations/summarize.py,sha256=hfYpgmv2irD8wNQ-GYtJA324w94SKoYhx0bpDG2bb3w,1207
+bella_companion/utils/__init__.py,sha256=Ehu8rgO6T2ie_BHW25I0-huoBIHPtknkHYvoUr-37BY,446
+bella_companion/utils/beast.py,sha256=18LAAFkZeIvOGdA3VfTEjQxZKzDBW93OixGtMNZrRxU,2196
+bella_companion/utils/explain/__init__.py,sha256=ksk8HWIpdiMhT4lW7H_05p_s44ihMiHHNM2TJGO9vz4,234
+bella_companion/utils/explain/pdp.py,sha256=z7PhwnERDUnsM1XzYq971pjxadbgBhhaO7ZOwY0fKuM,5125
+bella_companion/utils/explain/shap.py,sha256=ueCymQiVOcLIBy_Pf0yq7aOlI8ATehzhoxnBPYk8r0Y,3453
+bella_companion/utils/explain/typeguards.py,sha256=U6N_lBosmXqFAiakB3zkjlHXzKd16e_cYkHLh76SSFY,591
+bella_companion/utils/slurm.py,sha256=xEyf-Jxk8fy71t3fko_Ic9WtUFSdLFE3w4tR17gaBPw,1868
+bella_companion/version.xml,sha256=Gr3nsYqFhVNP0Mszl_IfilxefsPN8ZpQlGlB8OJSwnE,522
+bella_companion-0.1.1.dist-info/WHEEL,sha256=M6du7VZflc4UPsGphmOXHANdgk8zessdJG0DBUuoA-U,78
+bella_companion-0.1.1.dist-info/entry_points.txt,sha256=Jmns5P9G0EXcV26VD5Na81SD2Fri_SiUISy9Rfz5ZYI,52
+bella_companion-0.1.1.dist-info/METADATA,sha256=f_dbV3ST1DaBgVxI2t-cMnTN9U5197hPA8yGf4HN38U,472
+bella_companion-0.1.1.dist-info/RECORD,,

bella_companion-0.1.1.dist-info/WHEEL

@@ -0,0 +1,4 @@
+Wheel-Version: 1.0
+Generator: uv 0.9.5
+Root-Is-Purelib: true
+Tag: py3-none-any

{bella_companion-0.0.47.dist-info → bella_companion-0.1.1.dist-info}/entry_points.txt

@@ -1,2 +1,3 @@
 [console_scripts]
 bella = bella_companion.cli:main
+

bella_companion/eucovid/beast_configs/Nonparametric.xml

@@ -1,399 +0,0 @@
-<?xml version="1.0" encoding="UTF-8" standalone="no"?>
-<beast namespace="beast.base.inference
-    :beast.base.inference.parameter
-    :beast.base.inference.distribution
-    :beast.base.evolution.branchratemodel
-    :beast.base.evolution.substitutionmodel
-    :beast.base.evolution.sitemodel
-    :beast.base.evolution.tree.coalescent
-    :beast.base.evolution.tree
-    :beast.base.evolution.likelihood
-    :beast.base.evolution.alignment
-    :beast.base.evolution.operator.kernel
-    :beast.base.inference.operator.kernel
-    :beast.base.inference.operator
-    :beast.base.math.distributions
-    :beast.base.evolution
-    :beast.base.evolution.operator
-    :bdmmprime.util
-    :bdmmprime.util.operators
-    :bdmmprime.util.priors
-    :bdmmprime.distribution
-    :bdmmprime.parameterization
-    :bdmmprime.mapping
-    :bdmmprime.trajectories
-    :feast
-    :feast.parameter
-    :feast.fileio
-    :feast.function
-    :feast.expressions
-    :glmprior.util"
-    required="BEAST.base v2.7.7:BDMM-Prime v2.2.2" version="2.7">
-
-
-<!-- Map elements -->
-    <map name="Uniform" >beast.base.inference.distribution.Uniform</map>
-    <map name="Exponential" >beast.base.inference.distribution.Exponential</map>
-    <map name="LogNormal" >beast.base.inference.distribution.LogNormalDistributionModel</map>
-    <map name="prior" >beast.base.inference.distribution.Prior</map>
-
-<!-- Adapted from https://github.com/SarahNadeau/cov-europe-bdmm/blob/master/analyses/with_hubei_migration_decrease/with_hubei_migration_decrease_scrubbed.xml -->
-<!-- Sequence data loaded from GISAID with the gisaid EPI ISLs -->
-
-<!-- Alignment -->
-    <alignment id="alignment" spec= "AlignmentFromFasta"
-        fileName = "$(aligned_fasta)"/>
-
-    <typeSet id="typeSetBDMMPrime" spec="TypeSet">
-        <typeTraitSet id="typeTraitSet" spec="TraitSetFromTaxonSet" 
-            traitname="type" 
-            delimiter="|" 
-            takeGroup="1"> 
-        <taxa spec="TaxonSet" alignment="@alignment"/>
-        </typeTraitSet>
-    </typeSet>
-
-    <trait id="dateTrait" spec="TraitSetFromTaxonSet" 
-           traitname="date" 
-           dateFormat="yyyy-M-dd" 
-           delimiter= "|" 
-           everythingAfterLast="true">
-        <taxa spec="TaxonSet" alignment="@alignment"/>
-    </trait>
-
- <!-- BDMMPrime Population Model -->
-
-    <processLength id="originBDMMPrime" spec="RealParameter">1.0</processLength>
-
-    <Re id="ReSV" spec="SkylineVectorParameter" 
-        timesAreAges="true" processLength="@originBDMMPrime" typeSet="@typeSetBDMMPrime">
-        <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>
-         <changeTimes id="ReRateChangeTimes" spec="RealParameter">0.123</changeTimes>
-    </Re>
-
-    <samplingProportion id="samplingProportionSV" spec="SkylineVectorParameter" 
-                        timesAreAges="true" processLength="@originBDMMPrime" typeSet="@typeSetBDMMPrime">
-        <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>
-
-
-        <changeTimes id="samplingProportionChangeTimes" spec="RealParameter">0.12 0.205</changeTimes>
-    </samplingProportion>
-    
-
-    <becomeUninfectiousRate id="becomeUninfectiousRateSV" spec="SkylineVectorParameter"
-                            typeSet="@typeSetBDMMPrime">
-        <skylineValues id="becomeUninfectiousValues" spec="RealParameter">36.5</skylineValues>
-    </becomeUninfectiousRate>
-
-
-    <migrationRate id="migrationRateSM" spec="SkylineMatrixParameter" 
-                   timesAreAges="true"  processLength="@originBDMMPrime" typeSet="@typeSetBDMMPrime">
-        <skylineValues id="migrationRateValues_scaled" name="skylineValues" spec="feast.expressions.ExpCalculator">
-                            <![CDATA[
-                            migrationRateValues *
-                            {  f,f,f,f,
-                             f,  1,1,1,
-                             f,1,  1,1,
-                             f,1,1,  1,
-                             f,1,1,1,
-
-                               1,1,1,1,
-                             1,  1,1,1,
-                             1,1,  1,1,
-                             1,1,1,  1,
-                             1,1,1,1  }
-                            ]]>
-            <arg id="migrationRateValues" spec="RealParameter" dimension="40">
-                0.11 0.15 0.2 0.24 
-            0.0      0.14 0.18 0.22 
-            0.0 0.1       0.19 0.23 
-            0.0 0.12 0.16      0.25 
-            0.0 0.13 0.17 0.21 
-            
-                0.11 0.15 0.2 0.24 
-            0.0      0.14 0.18 0.22 
-            0.0 0.1       0.19 0.23 
-            0.0 0.12 0.16      0.25 
-            0.0 0.13 0.17 0.21
-        </arg>
-            <arg id="f" spec="RealParameter">0.5</arg>
-        </skylineValues>
-
-        <changeTimes id="migrationRateChangeTimes" spec="RealParameter">0.123</changeTimes>
-    </migrationRate> 
-
-
-    <removalProb id="removalProbSV" spec="SkylineVectorParameter" skylineValues="1.0"
-        processLength="@originBDMMPrime" typeSet="@typeSetBDMMPrime"/>
-
-    <startTypePriorProbs id="typeFrequencies" spec="RealParameter" 
-        estimate="false" lower="0.0" upper="1.0">1.0 0.0 0.0 0.0 0.0</startTypePriorProbs>
-
-
-<!-- Site Model -->
-    <siteModel id="siteModel" spec="SiteModel" gammaCategoryCount="4">
-        <parameter id="gammaShape" name="shape" value="1.0"/>
-        <parameter name="mutationRate"  estimate="false" value="1.0"/>
-        <parameter name="proportionInvariant" estimate="false" value="0.0"/>
-
-        <substModel spec="HKY">
-            <kappa id="kappa" spec="RealParameter" value="2.0"/>
-            <frequencies id="empiricalFreqs" spec="Frequencies" data="@alignment"/>
-        </substModel>
-    </siteModel>
-
-
-<!-- Branch rate Model -->
-    <branchRateModel id="branchRateModel" spec="StrictClockModel">
-        <parameter id="clockRate" name="clock.rate"  estimate="false">8.0E-4</parameter>
-    </branchRateModel>
-
-
-<!-- Population Model for Tree initialiser -->
-    <popFunc id="popFunc" spec="ConstantPopulation">
-        <popSize spec="RealParameter" value="0.1"/>
-    </popFunc>
-
-
-<!-- MCMC -->
-    <run id="mcmc" spec="MCMC" chainLength="10000000" numInitializationAttempts="100">
-
-    <!-- State -->
-        <state id="state" spec="State" storeEvery="5000">
-            
-            <stateNode id="treeBDMMPrime" spec="RandomTree" 
-                       taxa="@alignment" 
-                       populationModel="@popFunc"
-                       trait="@dateTrait">
-            </stateNode>
-            
-            <stateNode idref="originBDMMPrime"/>
-            <stateNode idref="ReValues"/>
-            <stateNode idref="samplingProportionValues"/>
-            <stateNode idref="migrationRateValues"/>
-            <stateNode idref="f"/>
-
-            <stateNode idref="kappa"/>
-            <stateNode idref="gammaShape"/>
-
-        </state>
-    
-
-    <!-- Posterior Distribution -->
-        <distribution id="posterior" spec="CompoundDistribution">
-
-             <!-- Tree Likelihood -->
-            <distribution id="likelihood" spec="CompoundDistribution" useThreads="true">
-                
-                <distribution id="treeLikelihood" spec="ThreadedTreeLikelihood" 
-                              data="@alignment" 
-                              tree="@treeBDMMPrime"
-                              siteModel="@siteModel"
-                              branchRateModel="@branchRateModel"/>
-                    
-            </distribution>
-
-            <!-- Prior distribution -->
-            <distribution id="prior" spec="CompoundDistribution">
-
-                <!-- BDMM-Prime tree prior distribution -->
-                <distribution id="BDMMPrime" spec="BirthDeathMigrationDistribution" 
-                              tree="@treeBDMMPrime"
-                              conditionOnSurvival="true"
-                              finalSampleOffset="0.0"
-                              typeTraitSet="@typeTraitSet"
-                              startTypePriorProbs="@typeFrequencies">
-
-                    <!-- Parameterization BDMM-Prime -->
-                    <parameterization id="EpiBDMMPrimeParameterization" spec="EpiParameterization" 
-                                      processLength="@originBDMMPrime"
-                                      Re="@ReSV"
-                                      samplingProportion="@samplingProportionSV"
-                                      becomeUninfectiousRate="@becomeUninfectiousRateSV"
-                                      migrationRate="@migrationRateSM"
-                                      removalProb="@removalProbSV"
-                                      typeSet="@typeSetBDMMPrime"/>
-
-                </distribution> 
-
-
-            <!-- Priors on parameters -->
-            <!-- BDMMPrime parameters priors -->
-
-                <prior id="originBDMMPrimePrior" name="distribution" x="@originBDMMPrime">
-                    <LogNormal name="distr" M="-1.0" S="0.2"/>
-                </prior>
-
-                <!-- Reproductive number Prior -->
-                <distribution id="RePrior" spec="SmartZeroExcludingPrior" x="@ReValues">
-                    <LogNormal name="distr" M="0.8" S="0.5"/>
-                </distribution>
-
-                <!-- Prior for sampling proportion for China -->
-                <distribution id="spChinaE2Prior" spec="SmartZeroExcludingPrior" x="@samplingProportionValues" 
-                    classesToExclude="1.2E-5 1.3E-5 1.4E-5 1.5E-5">
-                    <Uniform name="distr" lower="0" upper="0.15"/>
-                </distribution>
-                 <!-- Prior for sampling proportion for France -->
-                <distribution id="spFranceE1Prior" spec="SmartZeroExcludingPrior" x="@samplingProportionValues" 
-                    classesToExclude="1.1E-5 1.3E-5 1.4E-5 1.5E-5">
-                    <Uniform name="distr" lower="0" upper="0.093"/>
-                </distribution>
-                <!-- Prior for sampling proportion for Germany Epoch -->
-                <distribution id="spGermanyE1Prior" spec="SmartZeroExcludingPrior" x="@samplingProportionValues" 
-                    classesToExclude="1.1E-5 1.2E-5 1.4E-5 1.5E-5">
-                    <Uniform name="distr" lower="0" upper="0.10"/>
-                </distribution>
-                <!-- Prior for sampling proportion for Italy -->
-                <distribution id="spItalyE1Prior" spec="SmartZeroExcludingPrior" x="@samplingProportionValues" 
-                    classesToExclude="1.1E-5 1.2E-5 1.3E-5 1.5E-5">
-                    <Uniform name="distr" lower="0" upper="0.005"/>
-                </distribution>
-                <!-- Prior for sampling proportion for OtherEuropean -->
-                <distribution id="spOEE1Prior" spec="SmartZeroExcludingPrior" x="@samplingProportionValues" 
-                    classesToExclude="1.1E-5 1.2E-5 1.3E-5 1.4E-5">
-                    <Uniform name="distr" lower="0" upper="0.057"/>
-                </distribution>
-
-                <!-- Migration rate  Prior -->
-                <distribution id="migRatePrior" name="distribution" spec="SmartZeroExcludingPrior" x="@migrationRateValues">
-                    <LogNormal name="distr" M="0.0" S="1.0"/>
-                </distribution> 
-
-                <distribution id="fPrior" name="distribution" spec="ZeroExcludingPrior" x="@f">
-                    <Uniform name="distr" lower="0.0" upper="1.0"/>
-                </distribution>
-
-            <!-- Substitution model parameter priors -->
-                <prior id="GammaShapePrior" name="distribution" x="@gammaShape">
-                    <Exponential name="distr">
-                        <parameter spec="RealParameter" estimate="false" name="mean">0.5</parameter>
-                    </Exponential>
-                </prior>
-
-                <prior id="KappaPrior" name="distribution" x="@kappa">
-                    <LogNormal name="distr">
-                        <parameter spec="RealParameter" estimate="false" name="M">1.0</parameter>
-                        <parameter spec="RealParameter" estimate="false" name="S">1.25</parameter>
-                    </LogNormal>
-                </prior>
-
-            </distribution>
-        </distribution>
-
-
-    <!-- Operators -->
-        <!-- Tree operators -->
-        <operator id="BDMMPrimeTreeRootScaler" spec="BactrianScaleOperator" tree="@treeBDMMPrime" 
-            rootOnly="true" scaleFactor="0.25"  weight="5.0"/>
-
-        <operator id="BDMMPrimeUniformOperator" spec="BactrianNodeOperator" tree="@treeBDMMPrime" 
-            weight="30.0"/>
-
-        <operator id="BDMMPrimeSubtreeSlideScaler" spec="BactrianSubtreeSlide" tree="@treeBDMMPrime" 
-            weight="15.0"/>
-
-        <operator id="BDMMPrimeNarrow" spec="Exchange" tree="@treeBDMMPrime" 
-            weight="15.0"/>
-
-        <operator id="BDMMPrimeWide" spec="Exchange" tree="@treeBDMMPrime" 
-            isNarrow="false" weight="3.0"/>
-
-        <operator id="BDMMPrimeWilsonBalding" spec="WilsonBalding" tree="@treeBDMMPrime" 
-            weight="3.0"/>
-
-        <operator id="BDMMPrimeBICEPSEpochTop" spec="EpochFlexOperator" tree="@treeBDMMPrime" 
-            scaleFactor="0.1" weight="2.0"/>
-
-        <operator id="BDMMPrimeBICEPSEpochAll" spec="EpochFlexOperator" tree="@treeBDMMPrime" 
-            fromOldestTipOnly="false" scaleFactor="0.1" weight="2.0"/>
-
-        <operator id="BDMMPrimeBICEPSTreeFlex" spec="TreeStretchOperator" tree="@treeBDMMPrime" 
-            scaleFactor="0.01" weight="2.0"/>
-
-
-        <!-- Population model operators -->
-        <operator id="originBDMMPrimeScaler" spec="BactrianScaleOperator" parameter="@originBDMMPrime" 
-            scaleFactor="0.25" weight="5.0" />
-
-<!--         <operator id="ReScaler" spec="BactrianScaleOperator" parameter="@ReValues" 
-                weight="5.0" />
-        <operator id="ReScalerAll" spec="BactrianScaleOperator" parameter="@ReValues" 
-                weight="5.0" scaleAll="true" /> -->
-
-        <operator id="ReScaler" spec="SmartScaleOperator" parameter="@ReValues" 
-        weight="5.0" />
-
-        <operator id="samplingProportionScaler" spec="SmartScaleOperator" parameter="@samplingProportionValues" 
-                weight="3.0" />
-
-        <operator id="migrationRateScaler" spec="SmartScaleOperator" parameter="@migrationRateValues" 
-                weight="5.0" />
-        <operator id="fScaler" spec="BactrianScaleOperator" parameter="@f" 
-             weight="5.0"/>
-
-        <!-- Substitution model operators -->
-        <operator id="gammaShapeScaler" spec="BactrianScaleOperator" parameter="@gammaShape" 
-             weight="1.0"/>
-
-        <operator id="KappaScaler" spec="BactrianScaleOperator" parameter="@kappa" 
-             weight="1.0"/>
-
-
-    <!-- Loggers -->
-        <logger id="tracelog" spec="Logger" fileName="$(filebase).log" logEvery="1000" model="@posterior" sanitiseHeaders="true" sort="smart">
-            <log idref="posterior"/>
-            <log idref="likelihood"/>
-            <log idref="prior"/>
-            <log idref="treeLikelihood"/>
-            <log idref="gammaShape"/>   
-            <log idref="kappa"/>
-            <log id="TreeHeight" spec="TreeStatLogger" tree="@treeBDMMPrime"/>
-            <log id="typedTreeStats" spec="TypedTreeStatsLogger" includeRootEdge="true" 
-                typeLabel="type" typeSet="@typeSetBDMMPrime">
-                <typedTree id="typeMappedTree" spec="TypeMappedTree" bdmmDistrib="@BDMMPrime" startTypePriorProbs="@typeFrequencies" mapOnInit="false" remapOnLog="true" typeLabel="type" typeTraitSet="@typeTraitSet" untypedTree="@treeBDMMPrime"/>
-            </log>
-            <log idref="originBDMMPrime"/>
-            <log idref="ReSV"/>
-            <log idref="becomeUninfectiousRateSV"/>
-            <log idref="samplingProportionSV"/>
-            <log idref="migrationRateSM"/>
-            <log idref="migrationRateValues"/>
-            <log idref="f"/>
-            <log idref="removalProbSV"/>
-            <log idref="typeFrequencies"/>
-        </logger>
-
-        <logger id="screenlog" spec="Logger" logEvery="1000">
-            <log idref="posterior"/>
-            <log arg="@posterior" id="ESS_posterior" spec="util.ESS"/>
-            <log idref="likelihood"/>
-            <log arg="@likelihood" id="ESS_likelihood" spec="util.ESS"/>
-            <log idref="prior"/>
-            <log arg="@prior" id="ESS_prior" spec="util.ESS"/>
-            <log idref="BDMMPrime"/>
-            <log arg="@BDMMPrime" id="ESS_BDMMPrime" spec="util.ESS"/>
-        </logger>
-
-        <logger id="treelog" spec="Logger" fileName="$(filebase).trees" logEvery="1000" mode="tree">
-            <log id="TreeWithMetaDataLogger" spec="TreeWithMetaDataLogger" tree="@treeBDMMPrime"/>
-        </logger>
- 
-        <logger id="typedTreeLogger" spec="OptionalLogger" enableLogger="true" 
-            fileName="$(filebase).typed.trees" logEvery="10000" mode="tree">
-            <log idref="typeMappedTree"/>
-        </logger>
-        <logger id="nodeTypedTreeLogger.t:eubdmm_alignment" spec="OptionalLogger" enableLogger="true" 
-            fileName="$(filebase).typed.node.trees" logEvery="10000" mode="tree">
-            <log id="nodeTypedTree" spec="TypedNodeTreeLogger" typedTree="@typeMappedTree"/>
-        </logger>
-        <logger id="trajLogger" spec="OptionalLogger" enableLogger="true" fileName="$(filebase).traj" logEvery="10000">
-            <log id="typedTraj" spec="SampledTrajectory" bdmmDistrib="@BDMMPrime" 
-                startTypePriorProbs="@typeFrequencies" typeMappedTree="@typeMappedTree" useTauLeaping="true"/>
-        </logger>
-
-        <operatorschedule id="OperatorSchedule" spec="OperatorSchedule"/>
-
-    </run>
-
-</beast>