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
.
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.