def _stats(self, c):
        # See, for example, "A Statistical Study of Log-Gamma Distribution", by
        # Ping Shing Chan (thesis, McMaster University, 1993).
        mean = special.digamma(c)
        var = special.polygamma(1, c)
        skewness = np.power(var, 1.5) / special.polygamma(2, c)
        excess_kurtosis = special.polygamma(3, c) / (var*var)
        return mean, var, skewness, excess_kurtosis
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/loggamma_gen._stats/swapped operands of `/`