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)
A matrix containing the multivariate tip states for the n
species in the tree (and nodes if internal=
TRUE
).
is a phylogenetic tree in 'phylo' format; or a modified 'phylo' tree with a mapped discrete character.
is a square covariance matrix or named list of matrices (one for each mapped state on the tree).
optional vector of values for the root state.
logical value indicating whether to return states at internal nodes.
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.
fastBM
, make.simmap
, read.simmap
, sim.history
, sim.rates