Greetings all,
Does anyone know a way to transform branching times into branch length of a
chronogram. Like the inverse function of ape::branching.time() ?
say you have a chronogram named chrono:
tr <- chrono
tr$edge.length <- NULL
bt <- branching.times(chrono)
chrono2 <- magical.function.here(bt, tr)
here I want :
all(chrono2$edge.length == chrono$edge.length)
This would be helpful in situations such as parsing the mcmctree output
(mcmc.txt) that are made of sampled branching times by the MCMC (a big
table with row=iteration and col=node age) and then, say, run the same
diversification analysis on the MCMC trees to get some credibility
intervals on the diversification parameters.
Thanks !
--
Tristan
[[alternative HTML version deleted]]
_______________________________________________
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/
