code: def test_agreement():
      # The test: does Protocol manage to recreate the design of fMRIstat?
      X, c, D = create_protocols()
      for design_type in ["event", "block"]:
          dd = D[design_type]
          for i in range(X[design_type].shape[1]):
              _, cmax = matchcol(
                  X[design_type][:, i], fmristat[design_type]
              )
              if not dd.dtype.names[i].startswith("ns"):
                  assert np.greater(np.abs(cmax), 0.999)
  
file_path: nipy/nipy/modalities/fmri/fmristat/tests/test_FIAC.py
imports: import numpy
  from nipy.algorithms.statistics.formula import formulae
  from nipy.algorithms.statistics.models.regression import OLSModel
  from nipy.testing import assert_almost_equal
  from ...None import design, hrf, utils
  from ..None import hrf
  from .FIACdesigns import N_ROWS, altdescr, descriptions, fmristat, time_vector
label_desc: Use image node_id 7 for calling a global function with example usage: test_agreement() without return types