Started work on gibbs sampling