code: def _create_design(rng, n_samples, n_descriptors, n_regressors):
      target_var = rng.standard_normal((n_samples, n_descriptors))
      tested_var = rng.standard_normal((n_samples, n_regressors))
  
      return target_var, tested_var, n_descriptors, n_regressors
  
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 4 for calling a global function with example usage: _create_design(rng, n_samples, n_descriptors, n_regressors) and returns: target_var,  tested_var,  n_descriptors,  n_regressors