R/eBellB12.R

In ActuarialM: Computation of Actuarial Measures Using Bell G Family

#' @export
#' @import stats
eBellB12<-function(p, a, b, k, lambda)
{
    fn = function(x) {
        vBellB12(x, a, b, k, lambda)
    }
    ES = p
    for (i in 1:length(p)) {
        ES[i] = (1/p[i]) * integrate(fn, lower = 0, upper = p[i],
            stop.on.error = FALSE)$value
    }
    return(ES)
}