R/vBellB12.R

In ActuarialM: Computation of Actuarial Measures Using Bell G Family

#' @export
#' @import stats
vBellB12<-function (p, a, b, k, lambda, log.p = FALSE, lower.tail = TRUE)
{
    if (log.p == TRUE)
        p = exp(p)
    if (lower.tail == FALSE)
        p = 1 - p
		t=-1/lambda*log(1-((log(1-p[p >= 0 & p <= 1]*(1-(exp(-exp(lambda)+1)))))/(-exp(lambda))))
    VaR=(a*(((1-t)^(-1/k) -1)^(1/b)))
    return(VaR)
}