gibbs_betabin function

Gibbs sampling of the beta-binomial distribution