View source: R/getBMRatePrior.R
getBMRatePrior | R Documentation |
This function automatically calculates prior distributions for
the rate of trait evolution under the Brownian Motion (BM) model on a
discrete time-scale, at a given timeStep
, in the sense that that
variable is used with other TreEvo
functions like doRun_prc
.
getBMRatePrior(phy, traits, timeStep, verbose = TRUE)
phy |
A phylogenetic tree, in package |
traits |
Data matrix with rownames identical to |
timeStep |
time in a single iteration of the discrete-time simulation |
verbose |
If |
Returns a matrix of prior values that can be used in the doRun
functions.
Builds on functions in phylolm
to estimate distribution.
Returns a matrix of prior values
Brian O'Meara and Barb Banbury
data(simRunExample)
#timeStep = 0.1 -> effectively ~100 steps over the tree length
getBMRatePrior(phy = simPhyExample, traits = simCharExample,
timeStep = 0.01, verbose = TRUE)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.