Multivariate Brownian simulation with multiple correlations and rates
This function conducts BM simulation on a tree with multiple rates and/or multiple evolutionary correlations between characters. If vcv
is a single matrix, instead of a list of matrices, sim.corrs
will simulate multivariate BM with a single rate matrix.
sim.corrs(tree, vcv, anc=NULL, internal=FALSE)
tree |
is a phylogenetic tree in 'phylo' format; or a modified 'phylo' tree with a mapped discrete character. |
vcv |
is a square covariance matrix or named list of matrices (one for each mapped state on the tree). |
anc |
optional vector of values for the root state. |
internal |
logical value indicating whether to return states at internal nodes. |
A matrix containing the multivariate tip states for the n
species in the tree (and nodes if internal=
TRUE
).
Liam Revell liam.revell@umb.edu
Revell, L. J., and D. C. Collar (2009) Phylogenetic analysis of the evolutionary correlation using likelihood. Evolution, 63, 1090-1100.
Revell, L. J. (2012) phytools: An R package for phylogenetic comparative biology (and other things). Methods Ecol. Evol., 3, 217-223.
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.