run gibbs sampling