code: def test_omni_normtest_axis(reset_randomstate):
      # test axis of omni_normtest
      x = np.random.randn(25, 3)
      nt1 = omni_normtest(x)
      nt2 = omni_normtest(x, axis=0)
      nt3 = omni_normtest(x.T, axis=1)
      assert_almost_equal(nt2, nt1, decimal=13)
      assert_almost_equal(nt3, nt1, decimal=13)
  
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 3 for calling a global function with example usage: test_omni_normtest_axis(reset_randomstate) without return types