Add mcmc code