Internal function to generate random numbers following over-dispersed Bernoulli distribution, based in the BetaBinomial distribution. It currently does not work and standard Bernoulli distribution is actually used instead.
Useful links