code: def compare_to_ref_score(
      own_score, tested_var, target_var, covars=None
  ):
      reference = ref_score(tested_var, target_var, covars)
      assert_array_almost_equal(own_score, reference, decimal=6)
      return reference
  
imports: import nibabel
  import numpy
  import pytest
  from numpy.testing import assert_array_almost_equal, assert_array_less, assert_equal
  from scipy import stats
  from nilearn.conftest import _rng
  from nilearn.maskers import NiftiMasker
  from nilearn.mass_univariate import permuted_ols
label_desc: Use image node_id 2 for calling a global function with example usage: compare_to_ref_score(own_score, tested_var, target_var, covars) and returns: reference