MCMC code seems to work -- decent likelihood match