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

1 | ```
OU.vcv(phy, theta)
``` |

`phy` |
an ape format phylogeny |

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

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

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