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
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