ratematrixJointMCMC function

Estimate the evolutionary rate matrix together with the regimes using Markov-chain Monte Carlo