"phylo"
object with a mapped binary or multistate trait (see read.simmap
) and data for a single continuously valued character. It then uses restricted maximum likelihood (REML) to fit the Brownian rate variation ("noncensored") model of O'Meara et al. (2006; Evolution). This function is similar to brownie.lite
but uses REML (which is faster and unbiased) instead of ML. REML optimization takes advantage of Felsenstein's (1985) contrasts algorithm.
brownieREML(tree, x, maxit=2000, ...)
"phylo"
format (see read.simmap
and make.simmap
).names(x)
should be the species names.optim
.O'Meara, B. C., C. Ane, M. J. Sanderson, and P. C. Wainwright. 2006. Testing for different rates of continuous trait evolution using likelihood. Evolution, 60, 922-933.
Revell, L. J. (2012) phytools: An R package for phylogenetic comparative biology (and other things). Methods Ecol. Evol., 3, 217-223.
brownie.lite
, evol.vcv
, evol.rate.mcmc