code: def test_fmri_example():
      # Run FMRI analysis using example data
      from nipy.utils import example_data
  
      fmri_files = [
          example_data.get_filename("fiac", "fiac0", run)
          for run in ["run1.nii.gz", "run2.nii.gz"]
      ]
      design_files = [
          example_data.get_filename("fiac", "fiac0", run)
          for run in ["run1_design.npz", "run2_design.npz"]
      ]
      mask = example_data.get_filename("fiac", "fiac0", "mask.nii.gz")
      multi_session_model = FMRILinearModel(
          fmri_files, design_files, mask
      )
      multi_session_model.fit()
      (z_image,) = multi_session_model.contrast([np.eye(13)[1]] * 2)
      # Check number of voxels with p < 0.001
      assert np.sum(z_image.get_fdata() > 3.09) == 671
  
imports: import numpy
  import pytest
  from nibabel import Nifti1Image, load, save
  from numpy.testing import assert_almost_equal, assert_array_almost_equal, assert_array_equal
  from nipy.io.nibcompat import get_affine
  from nipy.testing import funcfile
  from nipy.testing.decorators import if_example_data
  from ..glm import FMRILinearModel, GeneralLinearModel, data_scaling
label_desc: Use image node_id 28 for calling a global function with example usage: test_fmri_example() without return types