code: def create_protocols():
      # block and event protocols
      block, bTcons, bFcons = protocol(
          descriptions["block"], "block", *delay.spectral
      )
      event, eTcons, eFcons = protocol(
          descriptions["event"], "event", *delay.spectral
      )
  
      # Now create the design matrices and contrasts
      # The 0 indicates that it will be these columns
      # convolved with the first HRF
      X = {}
      c = {}
      D = {}
      for f, cons, design_type in [
          (block, bTcons, "block"),
          (event, eTcons, "event"),
      ]:
          X[design_type], c[design_type] = f.design(t, contrasts=cons)
          D[design_type] = f.design(t, return_float=False)
      return X, c, D
  
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 4 for calling a global function with example usage: create_protocols() and returns: X,  c,  D