Dear all,

I am trying to compute the diversification rates as Jetz et al. 2012 (
https://www.nature.com/articles/nature11631), where the authors have said
that the formula is the inverse of the equal splits distinctiveness value
of Redding and Moors 2006 (
https://onlinelibrary.wiley.com/doi/full/10.1111/j.1523-1739.2006.00555.x).
But, in another paper, Pyron & Jetz 2018 (
https://www.nature.com/articles/s41559-018-0515-5) said “is closely
inversely related to the equal splits metric”. Looking in both formulas
they seem to be the “same” but only the inverse of equal splits, as Jetz et
al. said in 2012.

I would like to know if some of you have already calculated this
diversification rate and if it is in fact just the inverse, or it is a kind
of "closely" related formula. I wrote a small example of what I’ve done,
and would like to ask if you think that this is the right way to do that.


library(phytools)

tr <- rtree(20)

plot(tr)

tr <- compute.brlen(tr)

plot(tr)

tr$edge.length

tr$edge

plot(tr, show.tip.label = TRUE, show.node.label = TRUE)



ES = evol.distinct(tr, type = "equal.splits")$w

DR = evol.distinct(tr, type = "equal.splits")$w^-1


Thank you in advance,


-- 
Lilian Sayuri Ouchi de Melo
(PhD Student in Animal Biology at UNESP/IBILCE)

        [[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/

Reply via email to