This function fits a flexible Brownian multi-rate model using penalized likelihood.
multirateBM(tree, x, method=c("ML","REML"),
optim=c("L-BFGS-B","Nelder-Mead","BFGS","CG"),
maxit=NULL, n.iter=1, lambda=1, ...)an object of class "phylo".
a named numerical vector.
method of optimization. Currently only method="ML".
optimization routine to be used by optim. If more than one is specified and n.iter>1 then they will be alternated.
to be passed to optim. If set to maxit=NULL, the default value of maxit will be used, depending on the optimization method.
number of times to reiterate failed optimization.
lambda penalty term. High values of lambda correspond to high penalty for rate heterogeneity among edges. Low values of lambda correspond to low penalty.
optional arguments.
An object of class "multirateBM".
Revell, L. J. (2012) phytools: An R package for phylogenetic comparative biology (and other things). Methods Ecol. Evol., 3, 217-223.
Revell, L. J. (Accepted) A variable-rate quantitative trait evolution model using penalized-likelihood. PeerJ.