def _stats(self, a):
        mu2 = a*(a+1.0)
        return 0.0, mu2, 0.0, 3.0-(a+2.0)*(a+3.0)/mu2
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/dgamma_gen._stats/swapped operands of `-`