Learn R Programming

phytools (version 0.7-80)

multirateBM: Function to fit a multi-rate Brownian evolution model

Description

This function fits a flexible Brownian multi-rate model using penalized likelihood.

Usage

multirateBM(tree, x, method=c("ML","REML"), 
	optim=c("L-BFGS-B","Nelder-Mead","BFGS","CG"),
	maxit=NULL, n.iter=1, lambda=1, ...)

Arguments

tree

an object of class "phylo".

x

a named numerical vector.

method

method of optimization. Currently only method="ML".

optim

optimization routine to be used by optim. If more than one is specified and n.iter>1 then they will be alternated.

maxit

to be passed to optim. If set to maxit=NULL, the default value of maxit will be used, depending on the optimization method.

n.iter

number of times to reiterate failed optimization.

lambda

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.

Value

An object of class "multirateBM".

References

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.

See Also

brownie.lite