estDiversity | R Documentation |
This function estimates the lineage density at each node in the tree based on a biogeographic model (based on Mahler et al. 2010).
estDiversity(tree, x, method=c("asr","simulation"), model="ER", ...)
tree |
is a phylogenetic tree in |
x |
a vector containing the biogeographic area for each of the tip taxa. |
method |
method for reconstructing ancestral biogeography. |
model |
model for ancestral character estimation. In theory, any model from |
... |
optional arguments. So far, this includes only |
Two different methods are implemented in the current version. For method="asr"
the state at the current node, and at each position along each co-extant internal edge, is computed as the marginal (empirical Bayesian) ancestral state reconstruction using the re-rooting method of Yang (2006). The lineage density is then computed as the sum of the marginal reconstructions (posterior probabilities) times the summed marginal ancestral reconstructions across co-extant edges. In method="simulation"
, stochastic character mapping is used to generate optional argument nsim
stochastic maps of ancestral biogeography. Then the lineage density at each node is computed as the number of co-existing lineages with the same biogeography as the focal node, averaged acrossed stochastic maps. The importance of this distinction may depend on the degree to which reconstructions at internal nodes are independent, which relates to the distinction between marginal and joint reconstruction (e.g., see Yang 2006).
A vector containing the estimated lineage density at each node
Liam Revell liam.revell@umb.edu
Mahler, D. L, L. J. Revell, R. E. Glor, and J. B. Losos. (2010) Ecological opportunity and the rate of morphological evolution in the diversification of Greater Antillean anoles. Evolution, 64, 2731-2745.
Revell, L. J. (2012) phytools: An R package for phylogenetic comparative biology (and other things). Methods Ecol. Evol., 3, 217-223.
Yang, Z. (2006) Computational Molecular Evolution. Oxford University Press.
fitDiversityModel
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.