R/risk_gomp.R

In shannon: Computation of Entropy Measures and Relative Loss

#' @export
#' @import VaRES
#' @import stats
rlse_gomp <- function(p, alpha, beta) {
    SE <- function(p, alpha, beta) {
        fn = function(x) {
            (extraDistr::dgompertz(x,alpha, beta)) * ((extraDistr::dgompertz(x,alpha, beta, log = TRUE)) - log(extraDistr::pgompertz(p[i],
                alpha, beta)))
        }
        ES = p
        for (i in 1:length(p)) {
            ES[i] = (stats::integrate(fn, lower = 0, upper = p[i], stop.on.error = FALSE)$value)
        }
        return(ES)
    }
    sinc <- -1 * (1/extraDistr::pgompertz(p, alpha, beta)) * SE(p, alpha, beta)
    ((se_gomp(alpha, beta) - sinc))/(se_gomp(alpha, beta))
}

#' @export
rlre_gomp <- function(p, alpha, beta, delta) {
    RE <- function(p, alpha, beta, delta) {
        fn = function(x) {
            (extraDistr::dgompertz(x, alpha, beta))^delta
        }
        ES = p
        for (i in 1:length(p)) {
            ES[i] = (stats::integrate(fn, lower = 0, upper = p[i], stop.on.error = FALSE)$value)
        }
        return(ES)
    }

    rinc <- (1/(1 - delta)) * log((1/extraDistr::pgompertz(p, alpha, beta))^delta * RE(p,
        alpha, beta, delta))

    (re_gomp(alpha, beta, delta) - rinc)/re_gomp(alpha, beta, delta)
}

#' @export
rlhce_gomp <- function(p, alpha, beta, delta) {
    HC <- function(p, alpha, beta, delta) {
        fn = function(x) {
            (extraDistr::dgompertz(x, alpha, beta))^delta
        }
        ES = p
        for (i in 1:length(p)) {
            ES[i] = (stats::integrate(fn, lower = 0, upper = p[i], stop.on.error = FALSE)$value)
        }
        return(ES)
    }
    hcinc <- (1/(2^(1 - delta) - 1)) * (((1/extraDistr::pgompertz(p, alpha, beta))^delta *
        HC(p, alpha, beta, delta)) - 1)
    (hce_gomp(alpha, beta, delta) - hcinc)/hce_gomp(alpha, beta, delta)
}

#' @export
rlae_gomp <- function(p, alpha, beta, delta) {
    AE <- function(p, alpha, beta, delta) {
        fn = function(x) {
            (extraDistr::dgompertz(x, alpha, beta))^delta
        }
        ES = p
        for (i in 1:length(p)) {
            ES[i] = (stats::integrate(fn, lower = 0, upper = p[i], stop.on.error = FALSE)$value)
        }
        return(ES)
    }
    ainc <- (delta/(1 - delta)) * (((1/extraDistr::pgompertz(p, alpha, beta))^delta * AE(p,
        alpha, beta, delta))^(1/delta) - 1)

    (ae_gomp(alpha, beta, delta) - ainc)/ae_gomp(alpha, beta, delta)
}