code: def test_EC2(shape):
      fs = 0
      es = 0
      vs = 0
      ec = 0
  
      for f in decompose2d(shape, dim=3):
          ec += 1
          fs += 1
      for e in decompose2d(shape, dim=2):
          ec -= 1
          es += 1
      for v in decompose2d(shape, dim=1):
          ec += 1
          vs += 1
      return fs, es, vs, ec
  
file_path: nipy/nipy/algorithms/statistics/utils.py
imports: from itertools import combinations
  import numpy
  from scipy.stats import norm
label_desc: Use image node_id 10 for calling a global function with example usage: test_EC2(shape) and returns: fs,  es,  vs,  ec