skewt.py 文件源码

def logpdf_internal_prior(x, df, loc=0.0, scale=1.0, gamma = 1.0):
        if x-loc < 0.0:
            return np.log(2.0) - np.log(gamma + 1.0/gamma) + ss.t.logpdf(x=gamma*x, loc=loc*gamma,df=df, scale=scale)
        else:
            return np.log(2.0) - np.log(gamma + 1.0/gamma) + ss.t.logpdf(x=x/gamma, loc=loc/gamma,df=df, scale=scale)