OU.vcv: The variance covariance matrix for an Ornstein-Uhlenbeck process

Description

This function returns the variance-covariance matrix corresponding to an Ornstein-Uhlenbeck process run along a phylogeny

Usage

1
OU.vcv(phy, theta)

Arguments

phy

an ape format phylogeny

theta

the constraint parameter of the Ornstein-Uhlenbeck process

Value

a matrix with nrow = ncol = length(phy$tip.label), where the i,jth entry corresponds to the covariance between species i and j. NB: this is computed assuming that the rate of evolution is equal to one, and can be rescaled simply by multiplying


Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.