How to replicate the hierarchical Poisson model in BDA

I am studying Bayesian data analysis written by Gelman et al. The second chapter discussed a hierarchical poisson model on kidney cancer rate. Do anyone know how to replicate the result in the book using Stan? Thank you in advance.

In the same book, a Poisson model with hierarchical parameters is discussed on p. 420-22.