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.
rbinomod(n, prob, d = 1)a vector with the simulated sample.
integer. Number of observations.
numeric. Probability of success on each trial.
numeric. Dispersion parameter.