Gamma–Poisson hierarchical model

A Gamma–Poisson hierarchical model uses $Y_i\mid {\theta }_i\sim \mathrm{Poisson}(t_i{\theta }_i)$, ${\theta }_i\mid b\sim \mathrm{Gamma}(a,b)$, and $b\sim \mathrm{Gamma}(c,r)$, with gamma distributions in shape-rate convention. Conditional independence gives full conditional distributions ${\theta }_i\mid b,y\sim \mathrm{Gamma}(a+y_i,b+t_i)$ and $b\mid \theta ,y\sim \mathrm{Gamma}(c+na,r+{\sum }_i{\theta }_i)$.