Completed gibbs sampling ... need to test