fitBayes(tree, x, ngen=10000, model="BM", method="reduced", control=list())
"phylo"
.names(x)
should contain the species names (not individual IDs)."BM"
or "lambda"
."reduced"
or "full"
.sig2
: starting value for $\sigma^2$ (BM rate); lambda
: starting value for the $\lambda$ parameter; a
: starting for the state at the root node; xbar
: starting values for the states at the tips; intV
: starting value for the intraspecific variance (reduced method); or v
: starting value for the vector of intraspecific variances for all species (full method); pr.mean
: means for the prior distributions in the following order - sig2
, lambda
(if applicable), a
, xbar
, intV
or v
(if applicable), note that the prior probability distribution is exponential for sig2
and normal for a
and y
; pr.var
: variances on the prior distributions, same order as pr.mean
.ngen/control$sample+1
containing the posterior sample and likelihoods. Matrix columns are labeled by species (for species means and variances), or by the corresponding evolutionary parameter.
Revell, L. J. and R. G. Reynolds. (2012) A new Bayesian method for fitting evolutionary models to comparative data with intraspecific variation. Evolution, 66, 2697-2707.
anc.Bayes
, brownie.lite
, evol.rate.mcmc