I’m new to both Bayesian statistics and PyMC3, so if you can assume I don’t know anything, it would be much appreciated.

I’ve seen the example in the PyMC3 docs of implementing Gelman’s radon problem as a centered hierarchical linear model. I’ve also seen Tom Wiecki’s example of reparameterizing the same problem as a non-centered model, which seems like a clever way of specifying this model.

What I’m struggling with is, what if there were more levels to the hierarchical structure in the non-centered example? What would be an effective way of specifying that model? For example, suppose we had data not only for Minnesota, but now for all 50 states. We could then specify a hierarchy of state > county > town. I think seeing an example of how to make the hierarchical structure extensible would be extremely helpful to me. Thanks for any help you can provide.