powered by
Density function for a gamma-poisson distribution
dgpois(x, shape, rate, scale = 1/rate, log = FALSE)