code: def mean_squared_error(df, h0_intercept):
      return np.mean(
          (
              stats.t(df).cdf(np.sort(h0_intercept))
              - np.linspace(0, 1, h0_intercept.size + 1)[1:]
          )
          ** 2
      )
  
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 14 for calling a global function with example usage: mean_squared_error(df, h0_intercept) and returns: np