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)
"phylo"
format.minSplit
.$mcmc
from evol.rate.mcmc
(probably with burnin excluded).FALSE
).evol.rate.mcmc
, minSplit