def _stats(self, mu1, mu2):
        mean = mu2 - mu1
        var = mu1 + mu2
        g1 = mean / sqrt((var)**3)
        g2 = 1 / var
        return mean, var, g1, g2
dataset/ETHPy150Open scipy/scipy/scipy/stats/_discrete_distns.py/skellam_gen._stats/swapped operands of `-`