def _pln_logpdf(x, alpha, nu, tau2):
    return np.log(alpha) + alpha * nu + alpha * tau2 / 2 - \
        (alpha + 1) * np.log(x) + \
        norm.logcdf((np.log(x) - nu - alpha * tau2) / np.sqrt(tau2))