code: def twosample_stat(
      Y1, V1, Y2, V2, stat_id, axis=0, Magics=None, niter=DEF_NITER
  ):
      """
      Wrapper for ts_stat and ts_stat_mfx
      """
      if stat_id.find("_mfx") < 0:
          return ts_stat(Y1, Y2, stat_id, axis, Magics)
      else:
          return ts_stat_mfx(
              Y1, V1, Y2, V2, stat_id, axis, Magics, niter
          )
  
imports: import numpy
  import warnings
  from nipy.algorithms.graph import wgraph_from_3d_grid
  from nipy.algorithms.graph.field import Field, field_from_graph_and_data
  from ..utils import zscore
  from .onesample import stat
  from .onesample import stat_mfx
  from .twosample import stat
  from .twosample import stat_mfx
label_desc: Use image node_id 8 for calling a global function with example usage: twosample_stat(Y1, V1, Y2, V2, stat_id, axis, Magics, niter) and returns: ts_stat, ts_stat_mfx