Implementation of Markov chain Monte Carlo sampling and the Metropolis-Hastings algorithm for multi-parameter Bayesian inference.