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.