This function returns the variance-covariance matrix corresponding to an Ornstein-Uhlenbeck process run along a phylogeny
an ape format phylogeny
the constraint parameter of the Ornstein-Uhlenbeck process
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
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.