beast-dev
Relx clock using Thorney beast
Dear BEAST developers,
Currently Im trying to delevop a large scale analysis of SARS-CoV-2 using the Thorney-BEAST aproach.
On december a was available, thanks for your help, to run a Thorney-BEAST applying a Relax clock model and using a 100 sample Dataset, however when a try to run similar data set using a large Dataset the next error message is appearing
Alternatively, it may be that the product of starting mutation rate and tree height is extremely small or extremely large.
Finally, it may be that the initial state is incompatible with one or more 'hard' constraints (on monophyly or bounds on parameter values. This will result in Priors with zero probability.
The individual components of the posterior are as follows: The initial posterior is zero: CompoundLikelihood(compoundModel)=( CTMCScalePrior(ctmcScalePrior)=5.6971, Exponential(ucld.stdev)=0.0986, Gamma(skygrid.precision)=-4.6139, GMRFSkygridLikelihood(gmrfSkyrideLikelihood[skygrid])=-71670.6436 Total = -71669.46187390713 ), CompoundLikelihood(compoundModel)=( ThorneyTreeLikelihood(thorneyTreeLikelihood[treeLikelihood])=-Inf Total = -Infinity ) Total = -Infinity For more information go to http://beast.community.
I should mention that when running the same dataset using a strict clock model that works without any problems.
Hope someone can help finding the error on my XML file to try to running the thorney-beast using a Relax clock model
Attached is my current XML file
Large-dataset-relax-clock-model.txt
