gibbs sampling ... (incomplete)