code: def test_model():
      # Check basics about the model fit
      # Check we fit the mean
      assert_array_almost_equal(RESULTS.theta[1], np.mean(Y))
      # Check we get the same as R
      assert_array_almost_equal(RESULTS.theta, [1.773, 2.5], 3)
      percentile = np.percentile
      pcts = percentile(RESULTS.residuals, [0, 25, 50, 75, 100])
      assert_array_almost_equal(
          pcts, [-1.6970, -0.6667, 0, 0.6667, 1.6970], 4
      )
  
file_path: nilearn/nilearn/glm/tests/test_model.py
imports: import numpy
  import pytest
  from numpy.testing import assert_array_almost_equal
  from nilearn.glm import OLSModel
label_desc: Use image node_id 1 for calling a global function with example usage: test_model() without return types