code: def test_glm_ols():
      mulm, n, p, q = ols_glm()
      assert_array_equal(mulm.labels_, np.zeros(n))
      assert list(mulm.results_) == [0.0]
      assert mulm.results_[0.0].theta.shape == (q, n)
      assert_almost_equal(mulm.results_[0.0].theta.mean(), 0, 1)
      assert_almost_equal(mulm.results_[0.0].theta.var(), 1.0 / p, 1)
  
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 9 for calling a global function with example usage: test_glm_ols() without return types