gibbs subsampler