Issues with the Thorney Tree Likelihood Tutorial

143 views
Skip to first unread message

Clara Malekshahi

unread,
Jan 20, 2025, 2:54:03 PM1/20/25
to beast-users

Hello!

I have been trying to experiment with running BEAST starting with a tree to accommodate thousands of samples, as explained in this tutorial: Approaches for analyzing large phylogenetic datasets.

I first tried using the provided SG-thorney.template and beastgen to generate an XML file with my nwk tree. This give me the following error: 

Exception in thread "main" java.lang. NullPointerException

I was able to use beastgen and the exact code provided in the tutorial linked above to recreate the example XML file with a SARS-CoV-2 tree, but beastgen would not work with my XML or my tree. To see if the issue was with my file, I then attempted to run the example XML file on BEAST (I tried with both BEAST X v10.5.0-beta5 and the beast-mcmc release referenced in a previous post that seemed to have the same problem). Unfortunately, neither worked, and I received errors like the following:

Error in input BEAST XML file, 1K_SARS-CoV-2_SG-thorney.xml:

The '«thorneyTreeLikelihoods' element, with id, 'treeLikelihood', is incorrectly constructed.

The following was expected:

Exactly one ELEMENT of type MutationBranchMap REQUIRED

I did try to run an XML file (with my tree) I generated manually, and I got the same error when running BEAST. I was hoping someone would either know which BEAST release is most likely to work or whether I am using an incorrect version of the example SARS-CoV-2 XML file. I am happy to post more information or to try running the example XML file again to double check. Any suggestions would be wonderful! Thank you so much!

Clara M.



Clara Malekshahi

unread,
Feb 24, 2025, 3:54:21 PM2/24/25
to beast-users
Update to this post:

When I try to run the example XML file from the tutorial on BEAST 10.5.0, I get the following error: 

Random number seed: 1737735602371

Loading additional development parsers from development_parsers.properties, which is additional set of parsers only available for development version ...

Parsing XML file: 1K_SARS-CoV-2_SG-thorney.xml
  File encoding: UTF8
Looking for plugins in /Volumes/BEAST v10.5.0/plugins
Looking for plugins in /Users/claramal/Desktop/Moncla_Lab/plugins
Overriding checkpointing settings in the provided XML file

Taxon list 'taxa' created with 1000 taxa.
    most recent taxon date = 2020.5 years
Using node heights.

Creating the constrained tree model based on big fast tree model, 'treeModel'
  taxon count = 1000
  tree height = 0.4978826896494824
The gmrfSkyrideLikelihood has time aware smoothing
Creating a GMRF smoothed skyride model for multiple loci (SkyGrid)
Population sizes: 39

Feb 24, 2025 3:43:42 PM dr.app.beast.BeastMain <init>

SEVERE: 

******

Error in input BEAST XML file, 1K_SARS-CoV-2_SG-thorney.xml:

The '<thorneyTreeLikelihood>' element, with id, 'treeLikelihood', is incorrectly constructed.

The following was expected:

Exactly one ELEMENT of type MutationBranchMap REQUIRED

******

 
Does anyone know where the issue is coming from here? Am I using the wrong BEAST version?

Thank you very much!
Clara

JT McCrone

unread,
Mar 4, 2025, 5:24:09 PM3/4/25
to beast-users
Hi Clara,

The tutorial xml should run with the version linked to in the tutorial. (https://github.com/beast-dev/beast-mcmc/releases/tag/v1.10.5pre_thorney_v0.1.1).  (Please let me know if it doesn't!)

Overtime, the xml specification has been updated for better compatibility with standard analyses, but that work has not made it into a final release, yet. Once it does, the tutorial will be updated to reflect the changes. For now I would recommend using the tutorial specification with the prelease tagged in the tutorial. 

Please let me know if you run into any other difficulties. 

All the best,
JT
Reply all
Reply to author
Forward
0 new messages