ratematrixMCMC function

Estimate the evolutionary rate matrix using Markov-chain Monte Carlo