Thanks for the reply @ricardoV94. The pymc.Mixture example you linked appears to be multidimensional mixture in the sense that each component of the vector is a mixture of Gaussians. I’m looking to fit data that looks like this:
Where each 2-dimensional vector is a mixture of a multivariate Gaussian.