Parameterized as mu (the conditional mean) and shape (the same quantity
as size in rnbinom()). Thus Var(y) = mu + mu^2 / shape, which approaches
the Poisson variance as shape approaches infinity.
Parameterized as mu (the conditional mean) and shape (the same quantity
as size in rnbinom()). Thus Var(y) = mu + mu^2 / shape, which approaches
the Poisson variance as shape approaches infinity.