code: def test_cochransq2():
      # from an example found on web, verifies 13.286
      data = np.array(
          """
          0 0 0 1
          0 0 0 1
          0 0 0 1
          1 1 1 1
          1 0 0 1
          0 1 0 1
          1 0 0 1
          0 0 0 1
          0 1 0 0
          0 0 0 0
          1 0 0 1
          0 0 1 1""".split(),
          int,
      ).reshape(-1, 4)
  
      res = cochrans_q(data)
      assert_allclose(
          [res.statistic, res.pvalue],
          [13.2857143, 0.00405776],
          rtol=1e-6,
      )
  
imports: from statsmodels.compat.python import lzip
  import numpy
  from numpy.testing import assert_allclose, assert_almost_equal, assert_approx_equal, assert_
  from scipy import stats
  import pytest
  from statsmodels.stats.contingency_tables import mcnemar, cochrans_q, SquareTable
  from statsmodels.sandbox.stats.runs import Runs, runstest_1samp, runstest_2samp
  from statsmodels.sandbox.stats.runs import mcnemar
  from statsmodels.stats.nonparametric import rank_compare_2indep, rank_compare_2ordinal, prob_larger_continuous, cohensd2problarger
  from statsmodels.tools.testing import Holder
label_desc: Use image node_id 7 for calling a global function with example usage: test_cochransq2() without return types