def METHOD_NAME(self, p):
    lambda_ = self.lambda_
    sigma_ = exp(self.ln_sigma_)
    if lambda_ > 0:
        return exp(sigma_ * log(gammainccinv(1 / lambda_ ** 2, p) * lambda_ ** 2) / lambda_) * exp(self.mu_)
    return exp(sigma_ * log(gammaincinv(1 / lambda_ ** 2, p) * lambda_ ** 2) / lambda_) * exp(self.mu_)