R/eexp.R

In AntoineUC/Expectrem: Bias-reduced estimation of extreme value index and extreme quantiles, expectiles and Expected Shortfalls

eexp=function(probs,lambda=1){
  
  if(any(probs>=1) || any(probs<=0)){
    stop("only asymmetries between 0 and 1 allowed.")
  }
  
  if(length(lambda)>1){
    stop("lambda must be of size 1.")
  }
  
  if(lambda<=0){
    stop("lambda must be strictly positive.")
  }
  
  return(1/lambda*(1+lambertWp((2*probs-1)/(1-probs)*exp(-1))))
}