R/ES.R

In GenHMM1d: Goodness-of-Fit for Zero-Inflated Univariate Hidden Markov Models

#'@title Expected shortfall function
#'
#'@description This function computes the expected shortfall of an univariate distribution, excluding zero-inflated.
#'
#'@param p  value (1  x 1) at which the expected shortfall needs to be computed; between 0 and 1; (e.g 0.01, 0.05)
#'@param param  parameters of the distribution; (1 x p)
#'@param family   distribution name; run the function distributions() for help
#'@param size additional parameter for some discrete distributions; run the command distributions() for help
#'@param Nsim number of simulations
#'
#'@return \item{es}{expected shortfall}
#'
#'@examples
#'family = "gaussian"
#'
#'theta = c(-1.5, 1.7) ;
#'es = ES( 0.01, theta, family)
#'print('Expected shortfall : ')
#'print(es$es)
#'
#'
#'@export



ES<-function(p, param, family,  size=0, Nsim=25000){

  var = QUANTILE(p, param, family,  size)

  param_sim = matrix(param, nrow = 1, ncol = length(param))

  sim = SimHMMGen(param_sim, size=0, 1, ZI=0, family,  Nsim)

  diffp = sim$SimData-var

  diffp[diffp>0] = 0

  es = var + (1/(1-p)) * mean(replace(diffp, diffp == 0, NA), na.rm = TRUE)

  out = list(es=es, var=var)
  return(out)


}