powered by
posterior.evolrate
ave.rates(tree, shift, tips, sig1, sig2, ave.shift, showTree=TRUE)
shift
evol.rate.mcmc
minSplit