How do we predict on new unseen groups in a hierarchical model in PyMC3?

Good idea, done! Here’s the new issue