That’s a great idea, I’m glad that Nicole moved forward and started this project. I’m currently working on the Gaussian Mixture Models on my spare time and had a question: pyMC3 already has a GaussianMixture and Mixtureclass. I’m currently leaning towards reimplementing the marginalized model in the new package to have minimal coupling with pyMC3. Any argument against?