The function sim.rates
conducts BM simulation on a tree with multiple rates.
The function multiOU
conducts multi-regime OU simulations on the tree under a range of conditions.
sim.rates(tree, sig2, anc=0, nsim=1, internal=FALSE, plot=FALSE)
multiOU(tree, alpha, sig2, theta=NULL, a0=NULL, nsim=1, internal=FALSE, ...)
A vector (for nsim=1
) or matrix containing the tip states for the n
species in the tree.
is a stochastic map format phylogenetic tree in modified "phylo"
format (e.g., see make.simmap
).
a named vector containing the rates for each state; names should be states in mtree
.
optional value for the root state.
number of simulations.
logical value indicating whether to return states at internal nodes.
logical value indicating whether or not to visual the rate heterogeneity (default value is FALSE
.
single value or vector of values of the OU \(\alpha\) parameter.
single value or vector of values of the OU \(\theta\) parameter.
optional value of the root state. Defaults to zero.
optional arguments.
Liam Revell liam.revell@umb.edu
multiOU
uses a difference equation approximation of the OU process.
Revell, L. J. (2012) phytools: An R package for phylogenetic comparative biology (and other things). Methods Ecol. Evol., 3, 217-223.
fastBM
, make.simmap
, read.simmap
, sim.history