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

Completely agree, do you want to open an issue here asking for that? GitHub - pymc-devs/pymc-examples: Examples of PyMC models, including a library of Jupyter notebooks.