def _fitstart(self, data):  # pab
        loc = data.min() - 1e-4
        scale = (data - loc).max()
        m = np.mean((data - loc) / scale)
        fun = lambda c: (c - sc.log1p(c)) / (c * sc.log1p(c)) - m