Return lPosterior from gibbs sampler

Saves us from calculating it outside of the loop