Learn R Programming

phytools (version 0.5-38)

rateshift: Find the temporal position of one or more rate shifts

Description

Function finds the location of one or more rate shifts.

Usage

rateshift(tree, x, nrates=1, niter=10, method="ML", ...) "plot"(x, ...)

Arguments

tree
object of class "phylo".
x
vector of phenotypic trait values for species. names(x) should contain the species names and match tree$tip.label.
nrates
number of rates.
niter
number of iterations of optimization routine to ensure convergence.
method
optimization method. Can be "ML" (maximum likelihood) or "REML" (restricted maximum likelihood).
...
optional arguments. In the case of the plot method, these will be passed to plotSimmap. For rateshift, optional arguments include: tol, tolerance; plot & print, logical values indicating whether to plot or print the progress of the optimization (default to FALSE); quiet, logical argument indicating whether to suppress all notifications (defaults to FALSE); minL, numeric value; and fixed.shift, either a vector of fixed shift points, or a logical value.

Value

An object of class "rateshift".

References

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

See Also

brownie.lite