code: def test_jarque_bera():
      # tests against R fBasics
      st_pv_R = np.array([1.9662677226861689, 0.3741367669648314])
      jb = jarque_bera(x)[:2]
      assert_almost_equal(jb, st_pv_R, 14)
  
      st_pv_R = np.array([78.329987305556, 0.000000000000])
      jb = jarque_bera(x**2)[:2]
      assert_almost_equal(jb, st_pv_R, 13)
  
      st_pv_R = np.array([5.7135750796706670, 0.0574530296971343])
      jb = jarque_bera(np.log(x**2))[:2]
      assert_almost_equal(jb, st_pv_R, 14)
  
      st_pv_R = np.array([2.6489315748495761, 0.2659449923067881])
      jb = jarque_bera(np.exp(-(x**2)))[:2]
      assert_almost_equal(jb, st_pv_R, 14)
  
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 4 for calling a global function with example usage: test_jarque_bera() without return types