Hi Tristan, See compute.brtime in ape.
Cheers, Emmanuel Le 28/10/2015 15:31, Tristan Lefebure a écrit :
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/ Pour nous remonter une erreur de filtrage, veuillez vous rendre ici : http://f.security-mail.net/301eSIrvP0A
_______________________________________________ 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/