amkozlov / raxml-ng

RAxML Next Generation: faster, easier-to-use and more flexible
GNU Affero General Public License v3.0
374 stars 62 forks source link

Overcome the error: Constraint check failed after SPR round! #147

Closed padpadpadpad closed 1 year ago

padpadpadpad commented 1 year ago

Hi raxml-ng team

I am trying to set topological constraints on the tree based on the taxonomic assignment of some of the tips. The alignment is 246bp for 3541 taxa.

The command I am running is:

raxml-ng --msa alignment_myxo_97.7percent_new.fasta --model GTR+G --prefix myxo_97.7_constraint --threads 20 --seed 2 --outgroup otu_outgroup1 --tree-constraint contree.tre

However I am getting the error:

terminate called after throwing an instance of 'std::runtime_error'
  what():  ERROR in SPR round (LIBPLL-3584): Constraint check failed after SPR round!
Aborted (core dumped)

To us this hints that there is a problem with the constraint tree, but it is not obvious what this is. Any ideas what we could try to overcome this problem?

Many thanks Dan

Attached is the non-comprehensive constraint tree and the sequence alignment.

Archive.zip

amkozlov commented 1 year ago

Hi Dan,

apologies for the very late response.

I just tried to reproduce with raxml-ng v 1.1.0, but the analysis ran to completion without errors.

Which version of raxml-ng have you used? Could you please attach full log file?

Thanks, Alexey

padpadpadpad commented 1 year ago

Hi Alexey

I am very sorry for leaving this open, we solved the problem and I cannot remember how we did.

I have successfully ran raxml-ng using a constraint tree now.

Thank you for replying!

amkozlov commented 1 year ago

No problem, thanks for the confirmation!