added mcmc sampling scripts for comparison with the message passing model