Learn R Programming

phytools (version 2.1-1)

sim.corrs: Multivariate Brownian simulation with multiple correlations and rates

Description

Simulates multivariate Brownian motion evolution on a tree with multiple evolutionary correlation/covariance matrices.

Usage

sim.corrs(tree, vcv, anc=NULL, internal=FALSE)

Value

A matrix containing the multivariate tip states for the n species in the tree (and nodes if internal=

TRUE).

Arguments

tree

an object of class "phylo" or "simmap".

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.

Author

Liam Revell liam.revell@umb.edu

Details

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.

References

Revell, L. J., and D. C. Collar (2009) Phylogenetic analysis of the evolutionary correlation using likelihood. Evolution, 63, 1090-1100.

Revell, L. J. (2024) phytools 2.0: an updated R ecosystem for phylogenetic comparative methods (and other things). PeerJ, 12, e16505.

See Also

fastBM, make.simmap, read.simmap, sim.history, sim.rates