PyMC Discourse
ExpectationMax