{"word_count": 247, "character_count": 2349}
Almost 6 years ago, I posted a brief comment on a sentence I found surprising, by that time, discovered in a report claiming thatthe expected shortfall [...] at the 99 % level corresponds quite closely to the [...] value-at-risk at a 99.6% levelwhich was inspired by a remark in Swiss Experience report,expected shortfall [...] on a 99% confidence level [...} corresponds to approximately 99.6% to 99.8% Value at RiskRecall thatwhileFor any (absolutely) continuous cumulative distribution function, strictly increasing, since both (the VaR and the TVaR) are continuous, and strictly increasing, it is possible to relate any TVaR to some VaR, with a different level.  I.e.Which is not the same asConsider for instance the lognormal distribution. Since there is no simple expression for the expected shortfall, use monte carlo simulation to approximate it. And then, use the cumulative distribution function to get the assocated level for the value at risk,&#62; n=1e7&#62; TVaR_VaR_LN=function(p){+     X=rlnorm(n)+     E=mean(X[X&#62;qlnorm(p)])+     return(plnorm(E))+ }E.g.&#62; TVaR_VaR_LN(.99)[1] 0.9967621In order to plot it, define&#62; prob=c(seq(.8,.99,by=.01),.995)&#62; P_ln=unlist(lapply(prob,TVaR_VaR_LN))Now, if we consider a distribution with lighter tails, like the exponential distribution&#62; TVaR_VaR_exp=function(p){+     X=rexp(n)+     E=mean(X[X&#62;qexp(p)])+     return(pexp(E))+ }&#62; P_exp=unlist(lapply(prob,TVaR_VaR_exp))or a distribution with heavier tails, like the Pareto one,&#62; qpareto=function(u,a=2){(1-u)^(-1/a)}&#62; rpareto=function(n,a=2){qpareto(runif(n),a)}&#62; ppareto=function(x,a=2){1-1/x^a}&#62; TVaR_VaR_par=function(p){+     X=rpareto(n)+     E=mean(X[X&#62;qpareto(p)])+     return(ppareto(E))+ }&#62; P_pareto=unlist(lapply(prob,TVaR_VaR_par))we have different probability levels.&#62; plot(prob,P_ln,type="l",xlab="TVaR probability level",ylab="VaR probability level")&#62; lines(prob,P_pareto,type="l",col="red")&#62; lines(prob,P_exp,type="l",col="blue")&#62; legend("topleft",+        c("Pareto","Log Normal","Exponential"),+        col=c("red","black","blue"),lty=1)Hence, the heavier the tail, the higher the probability level. So, always qppfoximating a 99% TVaR with 99.6% VaR might work in some cases, e.g.&#62; TVaR_VaR_exp(.99)[1] 0.9963071but I is usually an optimistic approximation.