 Poisson-gamma conjugacy

For conditionally independent Poisson observations of mean $\theta$ and a gamma distribution prior of shape $r>0$ and rate $\lambda >0$, the likelihood function contributes ${\theta }^{{\sum }_i{x}_i}{e}^{-n\theta }$. Multiplying by the prior density gives a posterior distribution with shape $r+{\sum }_i{x}_i$ and rate $\lambda +n$. The single-observation update is shape $r+X$ and rate $\lambda +1$.

The posterior mean is the convex combinationThus the prior mean $r/\lambda$ and the sample maximum-likelihood estimate $\bar{x}$ bracket it. Integrating the next conditional Poisson distribution against this posterior distribution gives the Poisson-gamma mixture, a negative binomial distribution; its spread includes uncertainty in $\theta$ as well as the conditional count variation.