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

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.