code: def test_median_sigmaG_approx(axis, keepdims, atol=0.02):
      np.random.seed(0)
      a = np.random.normal(0, 1, size=(10, 10000))
  
      med, sigmaG = median_sigmaG(a, axis=axis, keepdims=keepdims)
      mu, sigma = mean_sigma(a, axis=axis, ddof=1, keepdims=keepdims)
  
      assert_allclose(med, mu, atol=atol)
      assert_allclose(sigmaG, sigma, atol=atol)
  
file_path: astroML/astroML/stats/tests/test_stats.py
imports: import pytest
  import numpy
  from numpy.testing import assert_array_almost_equal, assert_array_equal, assert_allclose
  from astroML.stats import mean_sigma, median_sigmaG, sigmaG, fit_bivariate_normal
  from astroML.stats.random import bivariate_normal, trunc_exp, linear
label_desc: Use image node_id 5 for calling a global function with example usage: test_median_sigmaG_approx(axis, keepdims, atol) without return types