evol.rate.mcmc(tree, x, ngen=10000, control=list())"phylo" format.names(x) should be the species names.sig1: starting value for $\sigma(1)^2$; sig2: starting value for $\sigma(2)^2$; a: starting value for a; sd1: standard deviation for the nocontrol are given in Revell et al. (2012).anc.Bayes, brownie.lite, evol.vcv, minSplit, posterior.evolrate