code: def test_mv_mixture_rvs_random(self):
      cov3 = np.array(
          [[1.0, 0.5, 0.75], [0.5, 1.5, 0.6], [0.75, 0.6, 2.0]]
      )
      mu = np.array([-1, 0.0, 2.0])
      mu2 = np.array([4, 2.0, 2.0])
      mvn3 = mvd.MVNormal(mu, cov3)
      mvn32 = mvd.MVNormal(mu2, cov3 / 2.0)
      np.random.seed(0)
      res = mv_mixture_rvs([0.4, 0.6], 5000, [mvn3, mvn32], 3)
      npt.assert_almost_equal(
          np.array([res.std(), res.mean(), res.var()]),
          np.array([1.874, 1.733, 3.512]),
          decimal=1,
      )
  
file_path: statsmodels/statsmodels/distributions/tests/test_mixture.py
imports: import numpy
  import numpy.testing
  from statsmodels.distributions.mixture_rvs import mv_mixture_rvs, MixtureDistribution
  import statsmodels.sandbox.distributions.mv_normal
  from scipy import stats
label_desc: Use image node_id 2 for calling the TestMixtureDistributions obj's underlying member method code with example usage: obj.test_mv_mixture_rvs_random() without return types