No reason at all to go the other way! Thanks for the advice. I’m a bit self taught and new to PyMC3/Bayesian Inference, do you know of any good articles/resources that provide a guide to creating a multinomial hierarchical model that actually divvys up the first level predictions into second level predictions?