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