A non-negative whole number of the number of random samples to generate.
lambda
A non-negative numeric vector of means.
theta
A non-negative numeric vector of the dispersion parameters for the gamma mixing distribution. The variance for the resultant gamma-Poisson distribution is lambda + lambda^2 * theta.