R/eBellW.R

In ActuarialM: Computation of Actuarial Measures Using Bell G Family

#' @export
#' @import stats
eBellW<-function(p, alpha , beta, lambda)
{
    fn = function(x) {
        vBellW(x, alpha , beta, 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)
}