def METHOD_NAME(self):
    rng = np.random.default_rng(2021)
    X = rng.negative_binomial(10, 0.5, (10, 10)) + rng.negative_binomial(100, 0.5, (10, 10))
    w, r, p = em.run_em(X, downsample=1e6, max_iter=100)
    # np.testing.assert_allclose([0.22986291449828194, 0.770137085501718], w)
    # np.testing.assert_allclose([22.612072509754206, 665.8608166342448], r)
    # np.testing.assert_allclose([0.1730657448451262, 0.8594447707151015], p)