Learn R Programming

phytools (version 0.5-64)

evol.vcv: Likelihood test for variation in the evolutionary VCV matrix

Description

This function takes a modified "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).

Usage

evol.vcv(tree, X, maxit=2000, vars=FALSE, ...)

Arguments

tree
a phylogenetic tree in modified "phylo" format (see read.simmap).
X
an n x m matrix of tip values for m continuously valued traits in n species - row names should be species names.
maxit
an optional integer value indicating the maximum number of iterations for optimization - may need to be increased for large trees.
vars
an optional logical value indicating whether or not to estimate the variances of the parameter estimates from the Hessian matrix.
...
optional arguments.

Value

A list with the following components:

Details

This function performs optimization by first optimizing the likelihood with respect to the Cholesky matrices using 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.

References

Revell, L. J., and D. C. Collar (2009) Phylogenetic analysis of the evolutionary correlation using likelihood. Evolution, 63, 1090-1100.

Revell, L. J. (2012) phytools: An R package for phylogenetic comparative biology (and other things). Methods Ecol. Evol., 3, 217-223.

See Also

evol.rate.mcmc, brownie.lite