powered by
Simple Gibbs sampler demonstration on conditional exponentials from Chapter 1 (pages 25-27).
expo.gibbs(B,k,m)
an upper bound
length of the subchains
number of iterations
Jeff Gill