def _logpdf(self, x, dfn, dfd):
        n = 1.0 * dfn
        m = 1.0 * dfd
        lPx = m/2 * log(m) + n/2 * log(n) + (n/2 - 1) * log(x)
        lPx -= ((n+m)/2) * log(m + n*x) + special.betaln(n/2, m/2)
        return lPx
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/f_gen._logpdf/original