def _logpmf(self, x, n, p):
        coeff = gamln(n+x) - gamln(x+1) - gamln(n)
        return coeff + n*log(p) + special.xlog1py(x, -p)
dataset/ETHPy150Open scipy/scipy/scipy/stats/_discrete_distns.py/nbinom_gen._logpmf/original