code: def mvsk2mc(args):
      """convert mean, variance, skew, kurtosis to central moments"""
      X = _convert_to_multidim(args)
  
      def _local_counts(args):
          mu, sig2, sk, kur = args
          cnt = [None] * 4
          cnt[0] = mu
          cnt[1] = sig2
          cnt[2] = sk * sig2**1.5
          cnt[3] = (kur + 3.0) * sig2**2.0
          return tuple(cnt)
  
      res = np.apply_along_axis(_local_counts, 0, X)
      # for backward compatibility convert 1-dim output to list/tuple
      return _convert_from_multidim(res, tuple)
  
file_path: statsmodels/statsmodels/stats/moment_helpers.py
imports: import numpy
  from scipy.special import comb
label_desc: Use image node_id 8 for calling a global function with example usage: mvsk2mc(args) and returns: _convert_from_multidim, tuple