mcmc_sample function

Sample from Stan model