How to avoid large number of unused random variables in hierarchical regression with 3 levels

numpyro
3

600 random variables flying around that serve absolutely no purpose at all

Could you elaborate on why this happens in your model?

For modeling, I would recommend reparameterizing your model (like this one).