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))
  
imports: import numpy
  from numpy.testing import assert_array_almost_equal
  from nipy.algorithms.statistics.api import OLSModel
  from ..outputters import output_F, output_T
label_desc: Use image node_id 1 for calling a global function with example usage: test_model() without return types