"phylo"
object with a mapped binary or multistate trait and data for an arbitrary number of continuously valued character. It then fits the multiple evolutionary variance-covariance matrix (rate matrix) model of Revell & Collar (2009; Evolution).
evol.vcv(tree, X, maxit=2000, vars=FALSE, ...)
"phylo"
format (see read.simmap
).n
x m
matrix of tip values for m
continuously valued traits in n
species - row names should be species names.optim
. Optimization is by method="Nelder-Mead"
. Using box constraints does not make sense here as they would be applied to the Cholesky matrix rather than the target parameters. May have to increase maxit
for large trees and more than 2 traits.
Revell, L. J. (2012) phytools: An R package for phylogenetic comparative biology (and other things). Methods Ecol. Evol., 3, 217-223.
evol.rate.mcmc
, brownie.lite