add mcmc function