code: def test_roi_features():
      """ """
      mroi = make_subdomain()
      dshape = (8, 3)
      data = np.random.randn(*dshape)
      mroi.set_roi_feature("data_mean", data)
      assert mroi.roi_features["data_mean"].shape == dshape
  
imports: from os.path import dirname
  from os.path import join
  import numpy
  from nibabel import Nifti1Image, load
  from numpy.testing import assert_almost_equal, assert_array_equal
  from nipy.io.nibcompat import get_affine
  from ..discrete_domain import domain_from_binary_array, grid_domain_from_image
  from ..hroi import HROI_as_discrete_domain_blobs
  from ..mroi import subdomain_from_array, subdomain_from_balls
label_desc: Use image node_id 6 for calling a global function with example usage: test_roi_features() without return types