First MCMC function + model