def _pdf(self, x, a, b):
        trm = _norm_pdf(a + b*log(x/(1.0-x)))
        return b*1.0/(x*(x-1))*trm
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/johnsonsb_gen._pdf/swapped operands of `-`