conflicting requirements for --tree when using comprehensive constraint tree

17 views
Skip to first unread message

Steven Tavis

unread,
Jun 10, 2021, 5:31:28 PMJun 10
to raxml

Hello,

I am attempting to estimate branch lengths based on a comprehensive constraint tree however I am running into what seems to be mutually incompatible requirements. If I use a random starting tree I get an error telling me that I need a newick tree as an input:

/home/cades/SLT02/programs/ParGenes/pargenes/pargenes_src/../pargenes_binaries/raxml-ng --msa $msa --model JTT+G4 --prefix $prefix --tree rand{1} --seed 1 --threads 1 --evaluate --tree-constraint $constrainttree
ERROR: Please provide a valid Newick file as an argument of --tree option.

However If I start with a newick tree it tells me that this option does not work with a constraint tree:

/home/cades/SLT02/programs/ParGenes/pargenes/pargenes_src/../pargenes_binaries/raxml-ng --msa $msa --model JTT+G4 --prefix $prefix --tree $besttree --seed 1 --threads 1 --evaluate --tree-constraint $constrainttree

RAxML-NG v. 0.9.0 released on 20.05.2019 by The Exelixis Lab.
Developed by: Alexey M. Kozlov and Alexandros Stamatakis.
Contributors: Diego Darriba, Tomas Flouri, Benoit Morel, Sarah Lutteropp, Ben Bettisworth.
Latest version: https://github.com/amkozlov/raxml-ng
Questions/problems/suggestions? Please visit: https://groups.google.com/forum/#!forum/raxml

RAxML-NG was called at 10-Jun-2021 16:52:39 as follows:

/home/cades/SLT02/programs/ParGenes/pargenes/pargenes_src/../pargenes_binaries/raxml-ng --msa /home/cades/SLT02/data/corrEvo/SLT01_corrEvoPgenesFirstRun/parse_run/results/PP_5460_faa/PP_5460_faa.raxml.rba --model JTT+G4 --prefix /home/cades/SLT02/data/corrEvo/SLT01_corrEvoPgenesFixedTopology/PP_5460_faa/PP_5460_faa --tree /home/cades/SLT02/data/corrEvo/SLT01_corrEvoPgenesFirstRun/mlsearch_run/results/PP_5460_faa/PP_5460_faa.raxml.bestTree --seed 1 --threads 1 --evaluate --tree-constraint /home/cades/SLT02/data/corrEvo/SLT01_corrEvoPgenesFixedTopology/PP_5460_faa/PP_5460_faa_constriant.newick

Analysis options:
  run mode: Evaluate tree likelihood
  start tree(s): user
  topological constraint: /home/cades/SLT02/data/corrEvo/SLT01_corrEvoPgenesFixedTopology/PP_5460_faa/PP_5460_faa_constriant.newick
  random seed: 1
  tip-inner: OFF
  pattern compression: ON
  per-rate scalers: OFF
  site repeats: ON
  branch lengths: proportional (ML estimate, algorithm: NR-FAST)
  SIMD kernels: AVX2
  parallelization: NONE/sequential


ERROR:  User and parsimony starting trees are not supported in combination with constrained tree inference.
       Please use random starting trees instead.

Both commands work if I remove the --evaluate and --tree-constraint $constrainttree options. Is there a way I can fix this?

Thanks for your help,
SLT

Alexey Kozlov

unread,
Jun 10, 2021, 7:44:11 PMJun 10
to ra...@googlegroups.com
Hi Steven,

raxml-ng can estimate branch lengths on binary trees only. So if your constraint tree contains
multifurcations, you will have to run constrained tree search instead of calling --evaluate.

Best,
Alexey
> --
> You received this message because you are subscribed to the Google Groups "raxml" group.
> To unsubscribe from this group and stop receiving emails from it, send an email to
> raxml+un...@googlegroups.com <mailto:raxml+un...@googlegroups.com>.
> To view this discussion on the web visit
> https://groups.google.com/d/msgid/raxml/9d761c3c-267d-4f93-af7b-7e3a521d9b13n%40googlegroups.com
> <https://groups.google.com/d/msgid/raxml/9d761c3c-267d-4f93-af7b-7e3a521d9b13n%40googlegroups.com?utm_medium=email&utm_source=footer>.

Steven Tavis

unread,
Jun 11, 2021, 9:52:39 AMJun 11
to raxml
Hi Alexi,

Thank you.

Steven
Reply all
Reply to author
Forward
0 new messages