def _stats(self, skew):
        ans, x, transx, skew, mask, invmask, beta, alpha, zeta = (
            self._preprocess([1], skew))
        m = zeta + alpha / beta
        v = alpha / (beta**2)
        s = 2.0 / (alpha**0.5) * np.sign(beta)
        k = 6.0 / alpha
        return m, v, s, k
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/pearson3_gen._stats/original