Hello colleagues,
I am calculating ace character mappings and model fit in a tree, and
noticed that rescaling the tree with root age in the range 0.5 - 120
produces identical results, but suddenly, for the root deeper than 130
or 140, all is wrong (e.g., all internal nodes have state probs 0.5/0.5
under ER).
I managed to reproduce this in some tutorial examples in phytools (e.g.,
http://www.phytools.org/Cordoba2017/ex/8/Anc-states-discrete.html),
simply by inserting a rescale line:
...
require(geiger)
...
eel.tree<-read.tree("elopomorph.tre")
eel.tree <- rescale(eel.tree, "depth", 130)
What am I doing wrong?
Any help will be much appreciated!
Martin
--
Martin J. Ramirez
Curador y Jefe, División Aracnología / Curator & Head, Division of Arachnology
Museo Argentino de Ciencias Naturales - CONICET
Av. Angel Gallardo 470
C1405DJR Buenos Aires
Argentina
rami...@macn.gov.ar
tel +54 11 4982-8370 ext. 169
_______________________________________________
R-sig-phylo mailing list - R-sig-phylo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
Searchable archive at http://www.mail-archive.com/r-sig-phylo@r-project.org/