def METHOD_NAME(self):
    X = np.array([[1, 2, 3]])
    em_results = (0, 0), (4, 5), (0.5, 0.6)
    results = em.conditionals(X, em_results)
    np.testing.assert_allclose(stats.nbinom(n=4, p=0.5).pmf([[1, 2, 3]]), results[0])
    np.testing.assert_allclose(stats.nbinom(n=5, p=0.6).pmf([[1, 2, 3]]), results[1])