rod function

Internal function to generate random numbers following over-dispersed Bernoulli or Poisson distribution. It wraps rpoisod and rbinomod functions.