def _stats(self, p):
        mu = 1.0/p
        qr = 1.0-p
        var = p / p / qr
        g1 = (2.0-p) / sqrt(qr)
        g2 = np.polyval([1, -6, 6], p)/(1.0-p)
        return mu, var, g1, g2
dataset/ETHPy150Open scipy/scipy/scipy/stats/_discrete_distns.py/geom_gen._stats/swapped operands of `/`