docs: def doane(data):
          """
          Modified Doane modified
          """
          from scipy.stats import skew
          n = len(data)
          sigma = np.sqrt(6. * (n - 2.) / (n + 1.) / (n + 3.))
          return 1 + np.log2(n) + \
              np.log2(1 + np.abs(skew(data)) / sigma)