code: def mad(a, c=0.6745, axis=0):
      """
      Median Absolute Deviation:
  
      median(abs(a - median(a))) / c
      """
      _shape = a.shape
      a.shape = np.prod(a.shape, axis=0)
      m = np.median(np.fabs(a - np.median(a))) / c
      a.shape = _shape
      return m
  
file_path: nipy/nipy/algorithms/statistics/models/utils.py
imports: import numpy
  import scipy.interpolate
label_desc: Use image node_id 1 for calling a global function with example usage: mad(a, c, axis) and returns: m