Sampler occasionally gets static in pymc3.1, but not in pymc3.0

FYI, if you are doing hierarchical linear model or linear mixed model, bambi is a really good library build on top of PyMC3 (and PyStan). It’s similar to brm in R. In my experience they reparameterize the model really well.