code: def test_robust_kurtosis_3d(self):
      x = np.tile(self.kurtosis_x, (10, 10, 1))
      kurtosis = np.array(robust_kurtosis(x, axis=2))
      for i, r in enumerate(self.expected_kurtosis):
          assert_almost_equal(r * np.ones((10, 10)), kurtosis[i])
  
imports: import numpy
  import pandas
  from numpy.testing import assert_almost_equal, assert_raises, assert_equal, assert_allclose
  from statsmodels.stats._adnorm import normal_ad
  from statsmodels.stats.stattools import omni_normtest, jarque_bera, durbin_watson, _medcouple_1d, medcouple, robust_kurtosis, robust_skewness
label_desc: Use image node_id 19 for calling the TestStattools obj's underlying member method code with example usage: obj.test_robust_kurtosis_3d() without return types