code: def test_glm_ar():
      mulm, n, p, q = ar1_glm()
      assert len(mulm.labels_) == n
      assert len(mulm.results_) > 1
      tmp = sum(
          [mulm.results_[key].theta.shape[1] for key in mulm.results_]
      )
      assert tmp == n
  
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 13 for calling a global function with example usage: test_glm_ar() without return types