Hi Jingfei
I don’t think you can mix sklearn BayesianGaussianMixture with PyMC3. The SciPy library usually calculates the single point which maximises the likelihood whereas PyMC (and other bayesian PPL languages) give a full posterior distribution.
Maybe if you could describe what you want to achieve then we can see if PyMC can help.