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)