code: def test_stats(self):
      trkwds = {"moments": "mvsk"}
      trkwds.update(self.stkwds)
      stkwds = {"moments": "mvsk"}
      stkwds.update(self.stkwds)
      mvsktr = np.array(self.dist.stats(*self.trargs, **trkwds))
      mvskst = np.array(self.statsdist.stats(*self.stargs, **stkwds))
      assert_almost_equal(mvskst[:2], mvsktr[:2], 8)
      if np.any(np.abs(mvskst[2:]) < 1):
          assert_almost_equal(mvskst[2:], mvsktr[2:], 1)
      else:
          assert_almost_equal(mvskst[2:] / mvsktr[2:], np.ones(2), 0)
  
imports: import numpy
  from numpy.testing import assert_almost_equal
  from scipy import stats
  from statsmodels.sandbox.distributions.extras import ExpTransf_gen, LogTransf_gen, squarenormalg, absnormalg, negsquarenormalg, squaretg
label_desc: Use image node_id 5 for calling the CheckDistEquivalence obj's underlying member method code with example usage: obj.test_stats() without return types