This fucntion takes a phylogenetic tree, an average split position, and a raw MCMC output from evol.rate.mcmc
and returns a posterior sample of evolutionary rates rootward (\(\sigma_1^2\)) and tipward (\(\sigma_2^2\)) from the average split.
posterior.evolrate(tree, ave.shift, mcmc, tips, showTree=FALSE)
a phylogenetic tree in "phylo"
format.
mean or median shift-point from the posterior sample (see minSplit
.
matrix $mcmc
from evol.rate.mcmc
(probably with burnin excluded).
list of stips in state \(\sigma_1^2\) for each sampled generation of MCMC.
optional logical value indicating whether or not to plot the stretched and shrunken tree generated by the pre-processing algorithm implemented in this function (default is FALSE
).
A matrix containing the posterior sample of evolutionary rates and shift-points between rates.
Revell, L. J., D. L. Mahler, P. Peres-Neto, and B. D. Redelings (2012) A new method for identifying exceptional phenotypic diversification. Evolution, 66, 135-146.