def _stats(self, n, p):
        Q = 1.0 / p
        P = Q - 1.0
        mu = n*P
        var = n*P*Q
        g1 = sqrt(n*P*Q)/(Q+P)
        g2 = (1.0 + 6*P*Q) / (n*P*Q)
        return mu, var, g1, g2
dataset/ETHPy150Open scipy/scipy/scipy/stats/_discrete_distns.py/nbinom_gen._stats/swapped operands of `/`