Markov chain Monte Carlo (MCMC) in julia