code: def medcouple(y, axis=0):
      """
      Calculate the medcouple robust measure of skew.
  
      Parameters
      ----------
      y : array_like
          Data to compute use in the estimator.
      axis : {int, None}
          Axis along which the medcouple statistic is computed.  If `None`, the
          entire array is used.
  
      Returns
      -------
      mc : ndarray
          The medcouple statistic with the same shape as `y`, with the specified
          axis removed.
  
      Notes
      -----
      The current algorithm requires a O(N**2) memory allocations, and so may
      not work for very large arrays (N>10000).
  
      .. [*] M. Hubert and E. Vandervieren, "An adjusted boxplot for skewed
         distributions" Computational Statistics & Data Analysis, vol. 52, pp.
         5186-5201, August 2008.
      """
      y = np.asarray(y, dtype=np.double)  # GH 4243
      if axis is None:
          return _medcouple_1d(y.ravel())
  
      return np.apply_along_axis(_medcouple_1d, axis, y)
  
file_path: statsmodels/statsmodels/stats/stattools.py
imports: from scipy import stats
  import numpy
  from statsmodels.tools.sm_exceptions import ValueWarning
label_desc: Use image node_id 9 for calling a global function with example usage: medcouple(y, axis) and returns: np, _medcouple_1d