code: def test_mean_sigma(a_shape, axis, ddof):
      np.random.seed(0)
  
      a = np.random.random(a_shape)
      mu1, sigma1 = mean_sigma(a, axis=axis, ddof=ddof)
  
      mu2 = np.mean(a, axis=axis)
      sigma2 = np.std(a, axis=axis, ddof=ddof)
  
      assert_array_almost_equal(mu1, mu2)
      assert_array_almost_equal(sigma1, sigma2)
  
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 1 for calling a global function with example usage: test_mean_sigma(a_shape, axis, ddof) without return types