Added mcmc script and model, needs diagnostics