Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Cannot compute beta density at 1 when beta is negative #19

Closed
richelbilderbeek opened this issue Mar 2, 2018 · 2 comments
Closed

Cannot compute beta density at 1 when beta is negative #19

richelbilderbeek opened this issue Mar 2, 2018 · 2 comments

Comments

@richelbilderbeek
Copy link
Member

[1] "ERROR"
 [1] ""                                                                                  
 [2] "                        BEAST v2.4.7, 2002-2017"                                   
 [3] "             Bayesian Evolutionary Analysis Sampling Trees"                        
 [4] "                       Designed and developed by"                                  
 [5] " Remco Bouckaert, Alexei J. Drummond, Andrew Rambaut & Marc A. Suchard"            
 [6] "                                    "                                              
 [7] "                     Department of Computer Science"                               
 [8] "                         University of Auckland"                                   
 [9] "                        remco@cs.auckland.ac.nz"                                   
[10] "                        alexei@cs.auckland.ac.nz"                                  
[11] "                                    "                                              
[12] "                   Institute of Evolutionary Biology"                              
[13] "                        University of Edinburgh"                                   
[14] "                           a.rambaut@ed.ac.uk"                                     
[15] "                                    "                                              
[16] "                    David Geffen School of Medicine"                               
[17] "                 University of California, Los Angeles"                            
[18] "                           msuchard@ucla.edu"                                      
[19] "                                    "                                              
[20] "                      Downloads, Help & Resources:"                                
[21] "                           http://beast2.org/"                                     
[22] "                                    "                                              
[23] "  Source code distributed under the GNU Lesser General Public License:"            
[24] "                   http://github.com/CompEvol/beast2"                              
[25] "                                    "                                              
[26] "                           BEAST developers:"                                      
[27] "   Alex Alekseyenko, Trevor Bedford, Erik Bloomquist, Joseph Heled, "              
[28] " Sebastian Hoehna, Denise Kuehnert, Philippe Lemey, Wai Lok Sibon Li, "            
[29] "Gerton Lunter, Sidney Markowitz, Vladimir Minin, Michael Defoin Platel, "          
[30] "                 Oliver Pybus, Chieh-Hsi Wu, Walter Xie"                           
[31] "                                    "                                              
[32] "                               Thanks to:"                                         
[33] "          Roald Forsberg, Beth Shapiro and Korbinian Strimmer"                     
[34] ""                                                                                  
[35] "Random number seed: 1519987266495"                                                 
[36] ""                                                                                  
[37] "61430_aco: 960 4"                                                                  
[38] "626029_aco: 960 4"                                                                 
[39] "630116_aco: 960 4"                                                                 
[40] "630210_aco: 960 4"                                                                 
[41] "AMNH13589_aco: 960 4"                                                              
[42] "B25702_aco: 960 4"                                                                 
[43] "B41613_aco: 960 4"                                                                 
[44] "B431_aco: 960 4"                                                                   
[45] "B48218_aco: 960 4"                                                                 
[46] "B87109_aco: 960 4"                                                                 
[47] "FALK1_aco: 960 4"                                                                  
[48] "KU21673_aco: 960 4"                                                                
[49] "KU25127_aco: 960 4"                                                                
[50] "KU3604_aco: 960 4"                                                                 
[51] "KU9813_aco: 960 4"                                                                 
[52] "UWBM54394_aco: 960 4"                                                              
[53] "UWBM54511_aco: 960 4"                                                              
[54] "UWBM54556_aco: 960 4"                                                              
[55] "bas3_aco: 960 4"                                                                   
[56] "chacoensis_aco: 960 4"                                                             
[57] "dabbenei_aco: 960 4"                                                               
[58] "meridae_aco: 960 4"                                                                
[59] "Alignment(anthus_aco)"                                                             
[60] "  22 taxa"                                                                         
[61] "  960 sites"                                                                       
[62] "  70 patterns"                                                                     
[63] ""                                                                                  
[64] ""                                                                                  
[65] "Error 110 parsing the xml input file"                                              
[66] ""                                                                                  
[67] "validate and intialize error: Cannot compute beta density at 1 when beta = %.3g"   
[68] ""                                                                                  
[69] "Error detected about here:"                                                        
[70] "  <beast>"                                                                         
[71] "      <run id='mcmc' spec='MCMC'>"                                                 
[72] "          <distribution id='posterior' spec='util.CompoundDistribution'>"          
[73] "              <distribution id='prior' spec='util.CompoundDistribution'>"          
[74] "                  <prior id='YuleBirthRatePrior.t:anthus_aco' name='distribution'>"
[75] ""                                                                                  
attr(,"status")
[1] 1
[1] "site model:"
$name
[1] "JC69"

$id
[1] NA

$gamma_site_model
$gamma_site_model$gamma_cat_count
[1] 0

$gamma_site_model$gamma_shape
[1] 0.8105238

$gamma_site_model$prop_invariant
[1] 0.7663484

$gamma_site_model$gamma_shape_prior_distr
$gamma_site_model$gamma_shape_prior_distr$name
[1] "inv_gamma"

$gamma_site_model$gamma_shape_prior_distr$id
[1] NA

$gamma_site_model$gamma_shape_prior_distr$alpha
$gamma_site_model$gamma_shape_prior_distr$alpha$name
[1] "alpha"

$gamma_site_model$gamma_shape_prior_distr$alpha$id
[1] NA

$gamma_site_model$gamma_shape_prior_distr$alpha$estimate
[1] FALSE

$gamma_site_model$gamma_shape_prior_distr$alpha$value
[1] 92.08046


$gamma_site_model$gamma_shape_prior_distr$beta
$gamma_site_model$gamma_shape_prior_distr$beta$name
[1] "beta"

$gamma_site_model$gamma_shape_prior_distr$beta$id
[1] NA

$gamma_site_model$gamma_shape_prior_distr$beta$estimate
[1] FALSE

$gamma_site_model$gamma_shape_prior_distr$beta$value
[1] -41.83111




[1] "clock model:"
$name
[1] "strict"

$id
[1] NA

$clock_rate_param
$clock_rate_param$name
[1] "clock_rate"

$clock_rate_param$id
[1] NA

$clock_rate_param$estimate
[1] FALSE

$clock_rate_param$value
[1] 66.6412


$clock_rate_distr
$clock_rate_distr$name
[1] "laplace"

$clock_rate_distr$id
[1] NA

$clock_rate_distr$mu
$clock_rate_distr$mu$name
[1] "mu"

$clock_rate_distr$mu$id
[1] NA

$clock_rate_distr$mu$estimate
[1] TRUE

$clock_rate_distr$mu$value
[1] 9.882322


$clock_rate_distr$scale
$clock_rate_distr$scale$name
[1] "scale"

$clock_rate_distr$scale$id
[1] NA

$clock_rate_distr$scale$estimate
[1] TRUE

$clock_rate_distr$scale$value
[1] -98.63774



[1] "tree prior:"
$name
[1] "yule"

$id
[1] NA

$birth_rate_distr
$birth_rate_distr$name
[1] "beta"

$birth_rate_distr$id
[1] NA

$birth_rate_distr$alpha
$birth_rate_distr$alpha$name
[1] "alpha"

$birth_rate_distr$alpha$id
[1] NA

$birth_rate_distr$alpha$estimate
[1] TRUE

$birth_rate_distr$alpha$value
[1] -4.153554


$birth_rate_distr$beta
$birth_rate_distr$beta$name
[1] "beta"

$birth_rate_distr$beta$id
[1] NA

$birth_rate_distr$beta$estimate
[1] FALSE

$birth_rate_distr$beta$value
[1] -93.9874
@richelbilderbeek
Copy link
Member Author

In the XML,

<?xml version="1.0" encoding="UTF-8" standalone="no"?><beast beautitemplate='Standard' beautistatus='' namespace="beast.core:beast.evolution.alignment:beast.evolution.tree.coalescent:beast.core.util:beast.evolution.nuc:beast.evolution.operators:beast.evolution.sitemodel:beast.evolution.substitutionmodel:beast.evolution.likelihood" required="" version="2.4">


    <data
id="anthus_aco"
name="alignment">
                    <sequence id="seq_61430_aco" taxon="61430_aco" totalcount="4" value="acaggttagaaactactctgttttctggctccttgtttaatgccctgtcctattttattgcgaaaattgtctgtttttcacagaaaactgagagtagtcaagggattccttgtcctttgctttggtctgcacagctgtcttgttttaaggccaagtggaatgagacagctgactcttcaggtgtgaaaacttggatgtaggtgtagattggctctcagttgacctccagctggtgcagattcttcagtttgtttgatggagctttgagaagtcctttcagactgaaggatactctgaattttagctatggagataatgtggatgttgtgtgtttaagtcctgtttgcagtttttttctgttcagtcagttattttactgtgtgagtcaggacccttagaagccctcagtggcaaccacagagcgcacagttaattttctgtgcaagaaaattaagatcatactctgtgtccaggaaagtcaagaatattcctggttttctctactgtaaaattttatcttgtaacttgtgtttgggtctgcatgattattcaaaaatcttagtagatttggaaggatgttgcatattatggaaacaaagttggaaaaagtttgtatcagttgcagtatttcttcacatcatttnttaacnncntnnnnnnnnngcttctgccacttgaaaagacaaattaaaaacnaatttataatgcttatatgctttagttacattngggtctttcagtaactttagtgcttttgatagccatacctgtgagnttgacagtgtctaaaattagaagtgttccttttcttctgctcttcccattctcgtgtgtcttcaatagtttctgcaaataatgatgtgcagacttagcattgattcaacagcagaggtaagcatacctgtggcttacttggcttcagcttatccagcagtgccaaccactctctgtttgtcttac"/>
                    <sequence id="seq_626029_aco" taxon="626029_aco" totalcount="4" value="acaggttagaaactactctgttttctggctgcttgtttaatgccctctcctattttattgtgacgattgtctgtttttcacagaaaactgagagtagtcaagggattccttgtcctttgctttggtgtgcacagctgtcttgttttaaggcccagtggaatgagacagctgactcttcaggtgtgaaaacttggatgtacgtgtagattggctctcagttgacctccagctggtgcagattcttcagtttgtttgatggagctttgagaagtcctttcagactgaaggatactctgaattttagctatggaaataatgtggatgttgtgtgtttaagcactgtttgcagtttttttctgttcagtcagttattttactgtgtgagtcaggactcttagaagccctctgtggcaaccacagagcgcatagttaattttctgtacaagaaaattaagatcctactcagtgttcaggaaagtcaagaatattcctggttttctctactgtaaaattttatcttgtaacttgtgtttgggtctgcatgattattcaaaaatcttagtagatttggaaggatgttgcattttatggaaacaaagttgggaaaagtttgtatcagttccagtatttcttcacatcatttnttaacnncntnnnnnnnnngcttctaccacttgaaaagacaaattaaaaacnaatttataatgcttatatgctttagttacattngggtctttcagtaactttagtgctttacttagccatacctgtgagcntggcagtgtctaaaattagaagtgttccttttcttctgctcttcccattctcgtgtgtcttcaatagtttctgcaaataatgatgtgcagacttagcattggttcaacagcagaggtaaaaatacctgtggcttacccggcttcagcttatccagcagtgccaaccactctctgtttgtcttac"/>
                    <sequence id="seq_630116_aco" taxon="630116_aco" totalcount="4" value="acaggttagaaactactctgttttctggctccttgtttaatgccctgtcctattttattgcgaaaattgtctgtttttcacagaaaactgagagtagtcaagggattccttgtcctttgctttggtctgcacagctgtcttgttttaaggccaagtggaatgagacagctgactcttcaggtgtgaaaacttggatgtaggtgtagattggctctcagttgacctccagctggtgcagattcttcagtttgtttgatggagctttgagaagtcctttcagactgaaggatactctgaattttagctatggagataatgtggatgttgtgtgtttaagtcctgtttgcagtttttttctgttcagtcagttattttactgtgtgagtcaggacccttagaagccctcagtggcaaccacagagcgcacagttaattttctgtgcaagaaaattaagatcatactctgtgtccaggaaagtcaagaatattcctggttttctctactgtaaaattttatcttgtaacttgtgtttgggtctgcatgattattcaaaaatcttagtagatttggaaggatgttgcatattatggaaacaaagttggaaaaagtttgtatcagttgcagtatttcttcacatcatttnttaacnncntnnnnnnnnngcttctgccacttgaaaagacaaattaaaaacnaatttataatgcttatatgctttagttacattngggtctttcagtaactttagtgcttttgatagccatacctgtgagnttgacagtgtctaaaattagaagtgttccttttcttctgctcttcccattctcgtgtgtcttcaatagtttctgcaaataatgatgtgcagacttagcattgattcaacagcagaggtaagcatacctgtggcttacttggcttcagcttatccagcagtgccaaccactctctgtttgtcttac"/>
                    <sequence id="seq_630210_aco" taxon="630210_aco" totalcount="4" value="acaggttagaaactactctgttttctggctccttgtttaatgccctgtcctattttattgtgacaattgtctgtttttcacagaaaactgagagtagtcaagggattccttgtcctttgctttggtctgcacagctgtcttgttttaaggccaagtggaatgagacagctgactcttcaggtgtgaaaacttggatgtaggtgtagattggctcttagttgacctccagctggtgcagattcttcagtttgtttgatggagctttgagaagtcctttcagactgaaggatactctgaattttagctatggagataatgtggatgttgtgtgtttaagccctgtttgcagtttttttctgttcagtcagttattttactgtgtgagtcaggacccttagaagccctcagtggcaaccacagagcgcacagttaattttctgtgcaagaaaattaacaacatactctgtgtccaggaaagtcaagaatattcctggttttctctactgtaaaattttatcttgtaacttgtgtttgggtctgcatgattattcaaaaatcttagtagatttggaaggatgttgcatattatggaaacaaagttggaaaaagtttgtatcagttccagtagttcnncangtnatttcttcacatcatttttaacctgcttctaccacttgaaaagacaaattaaaaacnaatttataatgcttatatgctttagttacattngggtctttcagtaactttagtgcttttgatagccatacctgtgagnttgacagtgtctaaaattagaagtgttccttttcttctgctcttcccattctygtgtgtcttcaatagtttctgcaaataatgatgtgcagacttagcattgattcaacagcagaggtaascatacctgtggcttacttggcttcagcttatccagcagtgccaaccactctctgtttgtcttac"/>
                    <sequence id="seq_AMNH13589_aco" taxon="AMNH13589_aco" totalcount="4" value="acaggttagaaactactctgttttctggctccttgtttaatgccctgtcctattttattgcgaaaattgtctgtttttcacagaaaactgagagtagtcaagggattccttgtcctttgctttggtctgcacagctgtcttgttttaaggccaagtggaatgagacagctgactcttcaggtgtgaaaacttggatgtaggtgtagattggctctcagttgacctccagctggtgcagattcttcagtttgtttgatggagctttgagaagtcctttcagactgaaggatactctgaattttagctatggagataatgtggatgttgtgtgtttaagtcctgtttgcagtttttttctgttcagtcagttattttactgtgtgagtcaggacccttagaagccctcagtggcaaccacagagcgcacagttaattttctgtgcaagaaaattaagatcatactctgtgtccaggaaagtcaagaatattcctggttttctctactgtaaaattttatcttgtaacttgtgtttgggtctgcatgattattcaaaaatcttagtagatttggaaggatgttgcatattatggaaacaaagttggaaaaagtttgtatcagttgcagtatttcttcacatcatttnttaacnncntnnnnnnnnngcttctgccacttgaaaagacaaattaaaaacnaatttataatgcttatatgctttagttacattngggtctttcagtaactttagtgcttttgatagccatacctgtgagnttgacagtgtctaaaattagaagtgttccttttcttctgctcttcccattctcgtgtgtcttcaatagtttctgcaaataatgatgtgcagacttagcattgattcaacagcagaggtaagcatacctgtggcttacttggcttcagcttatccagcagtgccaaccactctctgtttgtcttac"/>
                    <sequence id="seq_B25702_aco" taxon="B25702_aco" totalcount="4" value="acaggttagaaactactctgttttctggctccttgtttaatgccctgtcctattttattgtgacaattgtctgtttttcacagaaaactgagagtagtcaagggattccttgtcctttgctttggtctgcacagctgtcttgttttaaggcccagtggaatgagacagctgactcttcaggtgtgaaaacttggatgtaggtgtagattggctctcagttgacctccagctggtgcagattcttcactttgtttggtggagctttgagaagtcctttcagactgaaggatactctgaattttagctatggaaataatgtggatgttgtgtgtttaagccctgtttgcagtttttttctgttcagtcagttattttactgtgtgagtcaggactcttagaagccctctgtggcaaccacagagcgcatagttaattttctgtgcaagaaaattaagatcctactcaatgttcaggaaagtcaagaatattcctgtttttttctactgtaaaattttatcttgtaacttgtgtttgggtctgcatgattattcaaaaatcttagtagatttggaaggatgttgcatattatggaaacaaagttggaaaaagtttgtatcagttccagtatttcttcacatcatttnttaacnncntnnnnnnnnngcttctaccacttgaaaagacaaattaaaaacnaatttataatgcttatatgctttagttacattngggtctttcagtaactttagtgcttttgatagccatacctgtgagcttggcagtgtctaaaattagaagtgttccttttcttctgctcttcccattcttgtgtgtcttcaatagtttctgcaaataatgatgtgcagacttagcattggttcaacagcagaggtaaaaatacctgtggcttacttggcttcagcttatccagcagtgccaaccactctctgtttgtcttac"/>
                    <sequence id="seq_B41613_aco" taxon="B41613_aco" totalcount="4" value="acaggttagaaactactctgttttctggctgcttgtttaatgccctctcctattttattgtgacaattgtctgtttttcacagaaaactgagagtagtcaagggattccttgtcctttgctttggtctgcacagctgtcttgttttaaggcccagtggaatgagacagctgactcttcaggtgtgaaaacttggatgtaggtgtagattggctctcagttgacctccagctggtgcagattcttcagtttgtttgatggagctttgagaagtcctttcagactgaaggatactctgaattttagctatggaaataatgtggatgttgtgtgtttaagcactgtttgcagtttttttctgttcagtcagttattttactgtgtgagtcaggactcttagaagccctctgtggcaaccacagagcgcatagttaattttctgtacaagaaaattaagatcctactcagtgttcaggaaagtcaagaatattcctggttttctctactgtaaaattttatcttgtaacttgtgtttgggtctgcatgattattcaaaaatcttagtagatttggaaggatgttgcattttatggaaacaaagttgggaaaagtttgtatcagttccagtatttcttcacatcatttnttaacnncntnnnnnnnnngcttctaccacttgaaaagacaaattaaaaacnaatttataatgcttatatgctttagttacattngggtctttcagtaactttagtgcttttcatagccatacctgtgagcntggcagtgtctaaaattagaagtgttccttttcttctgctcttcccattctcgtgtgtcttcaatagtttctgcaaataatgatgtgcagacttagcattggttcaacagcagaggtaaaaatacctgtggcttacttggcttcagcttatccagcagtgccaaccactctctgtttgtcttac"/>
                    <sequence id="seq_B431_aco" taxon="B431_aco" totalcount="4" value="acaggttagaaactactctgttttctggctccttgtttaatgccctgtcctattttattgcgacaattgtctgtttttcacagaaaactgagagtagtcaagggattccttgtcctttgctttggtctgcacagctgtcttgttttaaggcgaagtggaatgagacagctgactcttcaggtgtgaaaacttggatgtaggtgtagattggctctcagttgacctccagctggtgcagattcttcagtttgtttgatggagctttgagaagtcctttcagactgaaggatactctgaattttagctatggagataatgtggatgttgtgtgtttaagccctgtttgcagtttttttctgttcagtcagttattttactgtgtgagtcaggacccttagaagccctcagtggcaaccacagagcgcacatttaattttctgtgcaagaaaattaacaacatactctgtgtccaggaaagtcaagaatattcctggttttctctactgtaaaattttatcttgtaacttgtgtttgggtctgcatgattattcaaaaatcttagtagatttggaaggatgttgcatattatggaaacaaagttggaaaaagtttgtatcagttccagtatttcttcacatcatttnttaacnncntnnnnnnnnngcttctaccacttcaaaagacaaattaaaaacnaatttataatgcttatatgctttagttacattagggtctttcagtaactttagtgcttttgatagccatacctgtgagnttgacagtgtctaaaattagaagtattccttttcttctgctcttcccattctcgtgtgtcttcaatagtttctgcaaataatgatgtgcagacttaacattgattcaccagcagaggtaagcatacctgtggcttacttggcttcagcttatccagcagtgccagccactctctgtttgtcttac"/>
                    <sequence id="seq_B48218_aco" taxon="B48218_aco" totalcount="4" value="acaggttagaaactactctgttttctggctccttgtttattgccctgtcctattttattgcgacaattgtctgtttttcacagaaaactgagagtagtcaagggattccttgtcctttgctttggtctgcacagctgtcttgttttaaggccaagtggaatgagacagctgactcttcaggtgtgaaaacttggatgtaggtgtagattggctctcagttgacttccagctggtgcagattcttcagtttgtttgatggagctttgagaagtcctttcagactgaaggatactctgaattttagctatggagataatgtggatgttgtgtgtttaagccctgtttgcagtttttttctgttcagtcagttattttactgtgtgagtcaggacccttagaagccctcagtggcaaccacagagcgcacagttaattttctgtgcaagaaaattaagatcatactctgtgtccaggaaagtcaagaatattcctggttttctctactgtaaaattttatcttataacttgtgtttgggtctgcatgattattcaaaaatcttagtagatttggaaggatgtngcatattatggaaacaaagttggaaaaagtttgtatcagttccagtatttcttcacatcatttnttaacnncntnnnnnnnnngcttctaccacttgaaaagacaaattaaaaacnaatttataatgcttatatgctttagttacattngggtctttcagtaactttagtgcttttgatagccatacctgtgagnttgacagtgtctaaaattagaagtgttccttttcttctgctcttcccattcttgtgtgtcttcaatagtttctgcaaataatgatgtgcagacttagcattgattcancagcagaggtaagcatacctgtggcttacttggcttcagcttatccagcagtgccaaccactctctgtttgtcttac"/>
                    <sequence id="seq_B87109_aco" taxon="B87109_aco" totalcount="4" value="acaggttagaaactactctgttttctggctgcttgtttaatgccctctcctattttattgtgacaattgtctgtttttcacagaaaactgagagtagtcaagggattccttgtcctttgctttggtctgcacagctgtcttgttttaaggcccagtggaatgagacagctgactcttcaggtgtgaaaacttggatgtaggtgtagattggctctcagttgacctccagctggtgcagatttttcagtttgtttgatggagctttgagaagtcctttcagactgaaggatactctgaattttagctatggaaataatgtggatgttgtgtgtttaagcactgtttgcagtttttttctgttcagtcagttattttactgtgtgagtcaggactcttagaagccctctgtggcaaccacagagcgcatagttaattttctgtgcaagaaaattaagatcctactcagtgttcaggaaagtcaagaatattcctggttttctctactgtaaaattttatcttgtaacttgtgtttgggtctgcatgattattcaaaaatcttagtagatttggaaggatgttgcatattatggaaacaaagttgggaaaagtttgtatcagttccagtatttcttcacatcatttnttaacnncntnnnnnnnnngcttctaccacttgaaaagacaaattaaaaacnaatttataatgcttatatgctttagttacattngggtctttcagtaactttagtgcttttcatagccatacctgtgagcntggcagtgtctaaaattagaagtgttccttttcttctgctcttcccattctcgtgtgtcttcaatagtttctgcaaataatgatgtgcagacttagcattggttcaacagcagaggtaaaaatacctgtggcttacttggcttcagcttatccagcagtgccaaccactctctgtttgtcttac"/>
                    <sequence id="seq_FALK1_aco" taxon="FALK1_aco" totalcount="4" value="acaggttagaaactactctgttttctggctccttgtttaatgccctgtcctattttattgcgacaattgtctgtttttcacagaaaactgagagtagtcaagggattccttgtcctttgctttggtctgcacagctgtcttgttttaaggccaagtggaatgagacagctgactcttcaggtgtgaaaacttggatgtaggtgtagattggctctcagttgacctccagctggtgcagattcttcagtttgtttgatggagctttgagaagtcctttcagactgaaggatactctgaattttagctatggagataatgtggatgttgtgtgtttaagtcctgtttgcagtttttttctgttcagtcagttattttactgtgtgagtcaggacccttagaagccctcagtggcaaccacagagcgcacagttaattttctgtgcaagaaaattaagatcatactctgtgtccaggaaagtcaagaatattcctggttttctctactgtaaaattttatcttgtaacttgtgtttgggtctgcatgattattcaaaaatcttagtagatttggaaggatgttgcatattatggaaacaaagttggaaaaagtttgtatcagttgcagtatttcttcacatcatttnttaacnncntnnnnnnnnngcttctgccacttgaaaagacaaattaaaaacnaatttataatgcttatatgctttagttacattngggtctttcagtaactttagtgcttttgatagccatacctgtgagnttgacagtgtctaaaattagaagtgttccttttcttctgctcttcccattctcgtgtgtcttcaatagtttctgcaaataatgatgtgcagacttagcattgattcaacagcagaggtaagcatacctgtggcttacttggcttcagcttatccagcagtgccaaccactctctgtttgtcttac"/>
                    <sequence id="seq_KU21673_aco" taxon="KU21673_aco" totalcount="4" value="acaggttagaaactactctgttttctggctccttgtttaatgccctgtcctattttattgtgacaattgtctgtttttcacagaaaactgagagtagtcaagggattccttgtcctttgctttggtctgcacagctgtcttgttttaaggcccagtggaatgagacagctgactcttcaggtgtgaaaacttggatgtaggtgtagattggctctcagttgacctccagctggtgcagattcttcagtttgtttgatggagctttgggaagtcctttcagactgaaggatactctgatttttagctatggaaataatgtggatgttgtgtgtttaagccctgtttgcagtttttttctgttcagtcagttattttactgtgtgagtcaggactcttagaagccctctgtggcaaccacagagcgcatagttaattttctgtgcaagaaaattaagatcctactcagtgttcaggaaagtcaagaatattcctggttttctctactgtaaaattttatcttgtaacttgtgtttgggtctgcatgattattcaaaaatcttagtagatttggaaggatgttgcatattatggaaacaaagttggaaaaagtttgtatcagttccagtatttcttcacatcatttnttaacnncntnnnnnnnnngcttctaccacttgaaaagacaaattaaaaaccaatttataatgcttatatgctttagttacattngggtctttcagtaactttagtgcttttgatagccatacctgtgagcttggcagtgtctaaaattagaagtgttccttttcttctgctcttcccattctcgtgtgtcttcaatagtttctgcaaataatgatgtgcagacttagcattggttcaacagcagaggtaaaaatacctgtggcttacttggcttcagcttatccagcagtgccaaccactctctgtttgtcttac"/>
                    <sequence id="seq_KU25127_aco" taxon="KU25127_aco" totalcount="4" value="acaggttagaaactactctgttttctggctccttgtttaatgccctgtcctattttattgcgacaattgtctgtttttcacagaaaactgagagtagtcaagggattccttgtcctttgctttggtctgcacagctgtcttgttttaaggccaagtggaatgagacagctgactcttcaggtgtgaaaacttggatgtaggtgtagattggctctcagttgacctccagctggtgcagattcttcagtttgtttgatggagctttgagaagtcctttcagactgaaggatactctgaattttagctatggagataatgtggatgttgtgtgtttaagccctgtttgcagtttttttctgttcagtcagttattttactgtgtgagtcaggacccttagaagccctcagtggcaaccacagagcgcacatttaattttctgtgcaagaaaattaacaacatactctgtgtccaggaaagtcaagaatattcctggttttctctactgtaaaattttatcttgtaacttgtgtttgggtctgcatgattattcaaaaatcttagtagatttggaaggatgttgcatattatggaaacaaagttggaaaaagtttgtatcagttccagtatttcttcacatcatttnttaacnncntnnnnnnnnngcttctaccacttcaaaagacaaattaaaaacnaatttataatgcttatatgctttagttacattagggtctttcagtaactttagtgcttttgatagccatacctgtgagnttgacagtgtctaaaattagaagtattccttttcttctgctcttcccattctngtgtgtcttcaatagtttctgcaaataatgatgtgcagacttaacattgattcaccagcagaggtaagcatacctgtggcttacttggcttcagcttatccagcagtgccagccactctctgtttgtcttac"/>
                    <sequence id="seq_KU3604_aco" taxon="KU3604_aco" totalcount="4" value="acaggttagaaactactctgttttctggctccttgtttaatgccctgtcctattttattgcgacaattgtctgtttttcacagaaaactgagagtagtcaagggattccttgtcctttgctttggtctgcacagctgtcttgttttaaggccnagtggaatgagacagctgactcttcaggtgtgaaaacttggatgtaggtgtagattggctctcagttgacctccagctggtgcagattcttcattttgtttgatggagctttgagaagtcctttcagactgaaggatactctgaattttagctatggagataatgtgaatgttgtgtgtttaagccctgtttgcagtttttttctgttcagtcagttattttactgtgtcagtcaggacccttagaagccctcagtggcaaccacagagcgcacagttaattttctgtgcaagaaaattaagatcatactctgtgtccaggaaagtcaagaatattcctggttttctctactgtaaaattttatcttgtaacttgtgtttgggtctgcatgattattcaaaaatcttagtagatttggaaggatgttgcatattatggaaacaaagttggaaaaagtttgtatcagttctagtatttcttcacatcatttnttaacnncntnnnnnnnnngcttctaccacttgaaaagacaaattaaaaacnaatttataatgcttatatgctttagttacattngggtctttcagtaactttagtgcttttgatagccatacctgtgagcttgacagtgtctaaaattagaagtattccttttcttctgctcttcccattctagtgtgtcttcaatagtttctgcaaataatgatgtgcagacttaccattgattcaacagcagaggtaannatacctgtggcttacttggcttcagcttatccagcagtgccaaccactctctgtttgtcttac"/>
                    <sequence id="seq_KU9813_aco" taxon="KU9813_aco" totalcount="4" value="acaggttagaaactactctgttttctggctccttgtttaatgccctgtcctattttattgtgacaattgtctgtttttcacagaaaactgagagtagtcaagggattccttgtcctttgctttggtctgcacagctgtcttgttttaaggccaagtggaatgagacagctgactcttcaggtgtgaaaacttggatgtaggtgtagattggctcttagttgacctccagctggtgcagattcttcagtttgtttgatggagctttgagaagtcctttcagactgaaggatactctgaattttagctatggagataatgtggatgttgtgtgtttaagccctgtttgcagtttttttctgttcagtcagttattttactgtgtgagtcaggacccttagaagccctcagtggcaaccacagagcgcacagttaattttctgtgcaagaaaattaacaacatactctgtgtccaggaaagtcaagaatattcctggttttctctactgtaaaattttatcttgtaacttgtgtttgggtctgcatgattattcaaaaatcttagtagatttggaaggatgttgcatattatggaaacaaagttggaaaaagtttgtatcagttccagtagttcnncangtnatttcttcacatcatttttaacctgcttctaccacttgaaaagacaaattaaaaacnaatttataatgcttatatgctttagttacattngggtctttcagtaactttagtgcttttgatagccatacctgtgagnttgacagtgtctaaaattagaagtgttccttttcttctgctcttcccattctcgtgtgtcttcaatagtttctgcaaataatgatgtgcagacttagcattgattcaacagcagaggtaagcatacctgtggcttacttggcttcagcttatccagcagtgccaaccactctctgtttgtcttac"/>
                    <sequence id="seq_UWBM54394_aco" taxon="UWBM54394_aco" totalcount="4" value="acaggttagaaactactctgttttctggctccttgtttaatgccctgtcctattttattgcgacaattgtctgtttttcacagaaaactgagagtagtcaagggattccttgtcctttgctttggtctgcacagctgtcttgttttaaggcgaagtggaatgagacagctgactcttcaggtgtgaaaacttggatgtaggtgtagattggctctcagttgacctccagctggtgcagattcttcagtttgtttgatggagctttgagaagtcctttcagactgaaggatactctgaattttagctatggagataatgtggatgttgtgtgtttaagccctgtttgcagtttttttctgttcagtcagttattttactgtgtgagtcaggacccttagaagccctcagtggcaaccacagagcgcacatttaattttctgtgcaagaaaattaacaacatactctgtgtccaggaaagtcaagaatattcctggttttctctactgtaaaattttatcttgtaacttgtgtttgggtctgcatgattattcaaaaatcttagtagatttggaaggatgttgcatattatggaaacaaagttggaaaaagtttgtatcagttccagtatttcttcacatcatttnttaacnncntnnnnnnnnngcttctaccacttgaaaagacaaattaaaaacnaatttataatgcttatatgctttagttacattagggtctttcagtaactttagtgcttttgatagccatacctgtgagnttgacagtgtctaaaattagaagtattccttttcttctgctcttcccattctcgtgtgtcttcaatagtttctgcaaataatgatgtgcagacttagcattgattcaccagcagaggtaagcatacctgtggcttacttggcttcagcttatccagcagtgccagccactctctgtttgtcttac"/>
                    <sequence id="seq_UWBM54511_aco" taxon="UWBM54511_aco" totalcount="4" value="acaggttagaaactactctgttttctggctccttgtttaatgccctgtcctattttattgcgaaaattgtctgtttttcacagaaaactgagagtagtcaagggattccttgtcctttgctttggtctgcacagctgtcttgttttaaggccaagtggaatgagacagctgactcttcaggtgtgaaaacttggatgtaggtgtagattggctctcagttgacctccagctggtgcagattcttcagtttgtttgatggagctttgagaagtcctttcagactgaaggatactctgaattttagctatggagataatgtggatgttgtgtgtttaagtcctgtttgcagtttttttctgttcagtcagttattttactgtgtgagtcaggacccttagaagccctcagtggcaaccacagagcgcacagttaattttctgtgcaagaaaattaagatcatactctgtgtccaggaaagtcaagaatattcctggttttctctactgtaaaattttatcttgtaacttgtgtttgggtctgcatgattattcaaaaatcttagtagatttggaaggatgttgcatattatggaaacaaagttggaaaaagtttgtatcagttgcagtatttcttcacatcatttnttaacnncntnnnnnnnnngcttctgccacttgaaaagacaaattaaaaacnaatttataatgcttatatgctttagttacattngggtctttcagtaactttagtgcttttgatagccatacctgtgagnttgacagtgtctaaaattagaagtgttccttttcttctgctcttcccattctcgtgtgtcttcaatagtttctgcaaataatgatgtgcagacttagcattgattcaacagcagaggtaagcatacctgtggcttacttggcttcagcttatccagcagtgccaaccactctctgtttgtcttac"/>
                    <sequence id="seq_UWBM54556_aco" taxon="UWBM54556_aco" totalcount="4" value="acaggttagaaactactctgttttctggctgcttgtttaatgccctgtcctattttattgtgacaattgtctgtttttcacagaaaactgagagtagtcaagggattccttgtcctttgctttggtctgcacagctgtcttgttttaaggcccagtggaatgagacagctgactcttcaggtgtgaaaacttggatgtaggtgtagattggctctcagttgacctccagctggtgcagattcttcagtttgtttgatggagctttgggaagtcctttcagactgaaggatactctgatttttagctatggaaataatgtggatgttgtgtgtttaagccctgtttgcagtttttttctgttcagtcagttattttactgtgtgagtcaggactcttagaagccctctgtggcaaccacagagcgcatagtttattttctgtgcaagaaaattaagatcctactcagtgttcaggaaagtcaagaatattcctggttttctctactgtaaaattttatcttgtaacttgtgtttgggtctgcatgattattcaaaaatcttagtagatttggaaggatgttgcatattatggaaacaaagttggaaaaagtttgtatcagttccagtatttcttcacatcatttnttaacnncntnnnnnnnnngcttctaccacttgaaaagacaaattaaaaaccaatttataatgcttatatgctttagttacattngagtctttcagtaactttagtgcttttgatagccatacctgtgagcntggcagtgtctaaaattagaagtgttccttttcttctgctcttcccattctcgtgtgtcttcaatagtttctgcaaataatgatgtgcagacttagcattggttcaacagcagaggtaaaaatacctgtggcttacttggcttcagcttatccagcagtgccaaccactctctgtttgtcttac"/>
                    <sequence id="seq_bas3_aco" taxon="bas3_aco" totalcount="4" value="nnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnn"/>
                    <sequence id="seq_chacoensis_aco" taxon="chacoensis_aco" totalcount="4" value="nnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnn"/>
                    <sequence id="seq_dabbenei_aco" taxon="dabbenei_aco" totalcount="4" value="nnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnn"/>
                    <sequence id="seq_meridae_aco" taxon="meridae_aco" totalcount="4" value="nnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnn"/>
                </data>


    


    


    
<map name="Uniform" >beast.math.distributions.Uniform</map>
<map name="Exponential" >beast.math.distributions.Exponential</map>
<map name="LogNormal" >beast.math.distributions.LogNormalDistributionModel</map>
<map name="Normal" >beast.math.distributions.Normal</map>
<map name="Beta" >beast.math.distributions.Beta</map>
<map name="Gamma" >beast.math.distributions.Gamma</map>
<map name="LaplaceDistribution" >beast.math.distributions.LaplaceDistribution</map>
<map name="prior" >beast.math.distributions.Prior</map>
<map name="InverseGamma" >beast.math.distributions.InverseGamma</map>
<map name="OneOnX" >beast.math.distributions.OneOnX</map>


<run id="mcmc" spec="MCMC" chainLength="10000000">
    <state id="state" storeEvery="5000">
        <tree id="Tree.t:anthus_aco" name="stateNode">
            <taxonset id="TaxonSet.anthus_aco" spec="TaxonSet">
                <alignment idref="anthus_aco"/>
            </taxonset>
        </tree>
        <parameter id="birthRate.t:anthus_aco" name="stateNode">1.0</parameter>
    </state>

    <init id="RandomTree.t:anthus_aco" spec="beast.evolution.tree.RandomTree" estimate="false" initial="@Tree.t:anthus_aco" taxa="@anthus_aco">
        <populationModel id="ConstantPopulation0.t:anthus_aco" spec="ConstantPopulation">
            <parameter id="randomPopSize.t:anthus_aco" name="popSize">1.0</parameter>
        </populationModel>
    </init>

    <distribution id="posterior" spec="util.CompoundDistribution">
        <distribution id="prior" spec="util.CompoundDistribution">
            <distribution id="YuleModel.t:anthus_aco" spec="beast.evolution.speciation.YuleModel" birthDiffRate="@birthRate.t:anthus_aco" tree="@Tree.t:anthus_aco"/>
            <prior id="YuleBirthRatePrior.t:anthus_aco" name="distribution" x="@birthRate.t:anthus_aco">
                <Beta id="Beta.100" name="distr">
                    <parameter id="RealParameter.200" estimate="true" name="alpha">-4.1535540483892</parameter>
                    <parameter id="RealParameter.201" estimate="false" name="beta">93.9874036703259</parameter>
                </Beta>
            </prior>
        </distribution>
        <distribution id="likelihood" spec="util.CompoundDistribution" useThreads="true">
            <distribution id="treeLikelihood.anthus_aco" spec="ThreadedTreeLikelihood" data="@anthus_aco" tree="@Tree.t:anthus_aco">
                <siteModel id="SiteModel.s:anthus_aco" spec="SiteModel">
                    <parameter id="mutationRate.s:anthus_aco" estimate="false" name="mutationRate">1.0</parameter>
                    <parameter id="gammaShape.s:anthus_aco" estimate="false" name="shape">1.0</parameter>
                    <parameter id="proportionInvariant.s:anthus_aco" estimate="false" lower="0.0" name="proportionInvariant" upper="1.0">0.766348364297301</parameter>
                    <substModel id="JC69.s:anthus_aco" spec="JukesCantor"/>
                </siteModel>
                <branchRateModel id="StrictClock.c:anthus_aco" spec="beast.evolution.branchratemodel.StrictClockModel">
                    <parameter id="clockRate.c:anthus_aco" estimate="false" name="clock.rate">66.6411960497499</parameter>
                </branchRateModel>
            </distribution>
        </distribution>
    </distribution>

    <operator id="YuleBirthRateScaler.t:anthus_aco" spec="ScaleOperator" parameter="@birthRate.t:anthus_aco" scaleFactor="0.75" weight="3.0"/>

    <operator id="YuleModelTreeScaler.t:anthus_aco" spec="ScaleOperator" scaleFactor="0.5" tree="@Tree.t:anthus_aco" weight="3.0"/>

    <operator id="YuleModelTreeRootScaler.t:anthus_aco" spec="ScaleOperator" rootOnly="true" scaleFactor="0.5" tree="@Tree.t:anthus_aco" weight="3.0"/>

    <operator id="YuleModelUniformOperator.t:anthus_aco" spec="Uniform" tree="@Tree.t:anthus_aco" weight="30.0"/>

    <operator id="YuleModelSubtreeSlide.t:anthus_aco" spec="SubtreeSlide" tree="@Tree.t:anthus_aco" weight="15.0"/>

    <operator id="YuleModelNarrow.t:anthus_aco" spec="Exchange" tree="@Tree.t:anthus_aco" weight="15.0"/>

    <operator id="YuleModelWide.t:anthus_aco" spec="Exchange" isNarrow="false" tree="@Tree.t:anthus_aco" weight="3.0"/>

    <operator id="YuleModelWilsonBalding.t:anthus_aco" spec="WilsonBalding" tree="@Tree.t:anthus_aco" weight="3.0"/>

    <logger id="tracelog" fileName="anthus_aco.log" logEvery="1000" model="@posterior" sanitiseHeaders="true" sort="smart">
        <log idref="posterior"/>
        <log idref="likelihood"/>
        <log idref="prior"/>
        <log idref="treeLikelihood.anthus_aco"/>
        <log id="TreeHeight.t:anthus_aco" spec="beast.evolution.tree.TreeHeightLogger" tree="@Tree.t:anthus_aco"/>
        <log idref="YuleModel.t:anthus_aco"/>
        <log idref="birthRate.t:anthus_aco"/>
    </logger>

    <logger id="screenlog" logEvery="1000">
        <log idref="posterior"/>
        <log id="ESS.0" spec="util.ESS" arg="@posterior"/>
        <log idref="likelihood"/>
        <log idref="prior"/>
    </logger>

    <logger id="treelog.t:anthus_aco" fileName="$(tree).trees" logEvery="1000" mode="tree">
        <log id="TreeWithMetaDataLogger.t:anthus_aco" spec="beast.evolution.tree.TreeWithMetaDataLogger" tree="@Tree.t:anthus_aco"/>
    </logger>

</run>

</beast>

change

                    <parameter id="RealParameter.201" estimate="false" name="beta">-93.9874036703259</parameter>

to

                    <parameter id="RealParameter.201" estimate="false" name="beta">93.9874036703259</parameter>

and everything is OK again.

@richelbilderbeek
Copy link
Member Author

Fixed

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant