Expected shortfall

From HandWiki
Short description: Measure of financial risk

Expected shortfall (ES) is a risk measure—a concept used in the field of financial risk measurement to evaluate the market risk or credit risk of a portfolio. The "expected shortfall at q% level" is the expected return on the portfolio in the worst [math]\displaystyle{ q\% }[/math] of cases. ES is an alternative to value at risk that is more sensitive to the shape of the tail of the loss distribution.

Expected shortfall is also called conditional value at risk (CVaR),[1] average value at risk (AVaR), expected tail loss (ETL), and superquantile.[2]

ES estimates the risk of an investment in a conservative way, focusing on the less profitable outcomes. For high values of [math]\displaystyle{ q }[/math] it ignores the most profitable but unlikely possibilities, while for small values of [math]\displaystyle{ q }[/math] it focuses on the worst losses. On the other hand, unlike the discounted maximum loss, even for lower values of [math]\displaystyle{ q }[/math] the expected shortfall does not consider only the single most catastrophic outcome. A value of [math]\displaystyle{ q }[/math] often used in practice is 5%.[citation needed]

Expected shortfall is considered a more useful risk measure than VaR because it is a coherent spectral measure of financial portfolio risk. It is calculated for a given quantile-level [math]\displaystyle{ q }[/math] and is defined to be the mean loss of portfolio value given that a loss is occurring at or below the [math]\displaystyle{ q }[/math]-quantile.

Formal definition

If [math]\displaystyle{ X \in L^p(\mathcal{F}) }[/math] (an Lp) is the payoff of a portfolio at some future time and [math]\displaystyle{ 0 \lt \alpha \lt 1 }[/math] then we define the expected shortfall as

[math]\displaystyle{ \operatorname{ES}_\alpha(X) = -\frac{1}{\alpha} \int_0^\alpha \operatorname{VaR}_\gamma(X) \, d\gamma }[/math]

where [math]\displaystyle{ \operatorname{VaR}_\gamma }[/math] is the value at risk. This can be equivalently written as

[math]\displaystyle{ \operatorname{ES}_\alpha(X) = -\frac{1}{\alpha} \left(\operatorname E[X \ 1_{\{X \leq x_{\alpha}\}}] + x_\alpha(\alpha - P[X \leq x_\alpha])\right) }[/math]

where [math]\displaystyle{ x_\alpha = \inf\{x \in \mathbb{R}: P(X \leq x) \geq \alpha\}=-\operatorname{VaR}_\alpha(X) }[/math] is the lower [math]\displaystyle{ \alpha }[/math]-quantile and [math]\displaystyle{ 1_A(x) = \begin{cases}1 &\text{if }x \in A\\ 0 &\text{else}\end{cases} }[/math] is the indicator function.[3] Note, that the second term vanishes for random variables with continuous distribution functions.

The dual representation is

[math]\displaystyle{ \operatorname {ES}_\alpha(X) = \inf_{Q \in \mathcal{Q}_\alpha} E^Q[X] }[/math]

where [math]\displaystyle{ \mathcal{Q}_\alpha }[/math] is the set of probability measures which are absolutely continuous to the physical measure [math]\displaystyle{ P }[/math] such that [math]\displaystyle{ \frac{dQ}{dP} \leq \alpha^{-1} }[/math] almost surely.[4] Note that [math]\displaystyle{ \frac{dQ}{dP} }[/math] is the Radon–Nikodym derivative of [math]\displaystyle{ Q }[/math] with respect to [math]\displaystyle{ P }[/math].

Expected shortfall can be generalized to a general class of coherent risk measures on [math]\displaystyle{ L^p }[/math] spaces (Lp space) with a corresponding dual characterization in the corresponding [math]\displaystyle{ L^q }[/math] dual space. The domain can be extended for more general Orlicz Hearts.[5]

If the underlying distribution for [math]\displaystyle{ X }[/math] is a continuous distribution then the expected shortfall is equivalent to the tail conditional expectation defined by [math]\displaystyle{ \operatorname{TCE}_{\alpha}(X) = E[-X\mid X \leq -\operatorname{VaR}_{\alpha}(X)] }[/math].[6]

Informally, and non-rigorously, this equation amounts to saying "in case of losses so severe that they occur only alpha percent of the time, what is our average loss".

Expected shortfall can also be written as a distortion risk measure given by the distortion function

[math]\displaystyle{ g(x) = \begin{cases}\frac{x}{1-\alpha} & \text{if }0 \leq x \lt 1-\alpha,\\ 1 & \text{if }1-\alpha \leq x \leq 1.\end{cases} \quad }[/math][7][8]

Examples

Example 1. If we believe our average loss on the worst 5% of the possible outcomes for our portfolio is EUR 1000, then we could say our expected shortfall is EUR 1000 for the 5% tail.

Example 2. Consider a portfolio that will have the following possible values at the end of the period:

probability
of event
ending value
of the portfolio
10% 0
30% 80
40% 100
20% 150

Now assume that we paid 100 at the beginning of the period for this portfolio. Then the profit in each case is (ending value−100) or:

probability
of event
profit
10% −100
30% −20
40% 0
20% 50

From this table let us calculate the expected shortfall [math]\displaystyle{ \operatorname{ES}_q }[/math] for a few values of [math]\displaystyle{ q }[/math]:

[math]\displaystyle{ q }[/math] expected shortfall [math]\displaystyle{ \operatorname{ES}_q }[/math]
5% 100
10% 100
20% 60
30% 46.6
40% 40
50% 32
60% 26.6
80% 20
90% 12.2
100% 6

To see how these values were calculated, consider the calculation of [math]\displaystyle{ \operatorname{ES}_{0.05} }[/math], the expectation in the worst 5% of cases. These cases belong to (are a subset of) row 1 in the profit table, which have a profit of −100 (total loss of the 100 invested). The expected profit for these cases is −100.

Now consider the calculation of [math]\displaystyle{ \operatorname{ES}_{0.20} }[/math], the expectation in the worst 20 out of 100 cases. These cases are as follows: 10 cases from row one, and 10 cases from row two (note that 10+10 equals the desired 20 cases). For row 1 there is a profit of −100, while for row 2 a profit of −20. Using the expected value formula we get

[math]\displaystyle{ \frac{ \frac{10}{100}(-100)+\frac{10}{100}(-20) }{ \frac{20}{100}} = -60. }[/math]

Similarly for any value of [math]\displaystyle{ q }[/math]. We select as many rows starting from the top as are necessary to give a cumulative probability of [math]\displaystyle{ q }[/math] and then calculate an expectation over those cases. In general, the last row selected may not be fully used (for example in calculating [math]\displaystyle{ -\operatorname{ES}_{0.20} }[/math] we used only 10 of the 30 cases per 100 provided by row 2).

As a final example, calculate [math]\displaystyle{ -\operatorname{ES}_1 }[/math]. This is the expectation over all cases, or

[math]\displaystyle{ 0.1(-100)+0.3(-20)+0.4\cdot 0+0.2\cdot 50 = -6. \, }[/math]

The value at risk (VaR) is given below for comparison.

[math]\displaystyle{ q }[/math] [math]\displaystyle{ \operatorname{VaR}_q }[/math]
[math]\displaystyle{ 0\% \le q \lt 10\% }[/math] −100
[math]\displaystyle{ 10\% \le q \lt 40\% }[/math] −20
[math]\displaystyle{ 40\% \le q \lt 80\% }[/math] 0
[math]\displaystyle{ 80\% \le q \le 100\% }[/math] 50

Properties

The expected shortfall [math]\displaystyle{ \operatorname{ES}_q }[/math] increases as [math]\displaystyle{ q }[/math] decreases.

The 100%-quantile expected shortfall [math]\displaystyle{ \operatorname{ES}_{1} }[/math] equals negative of the expected value of the portfolio.

For a given portfolio, the expected shortfall [math]\displaystyle{ \operatorname{ES}_q }[/math] is greater than or equal to the Value at Risk [math]\displaystyle{ \operatorname{VaR}_q }[/math] at the same [math]\displaystyle{ q }[/math] level.

Optimization of expected shortfall

Expected shortfall, in its standard form, is known to lead to a generally non-convex optimization problem. However, it is possible to transform the problem into a linear program and find the global solution.[9] This property makes expected shortfall a cornerstone of alternatives to mean-variance portfolio optimization, which account for the higher moments (e.g., skewness and kurtosis) of a return distribution.

Suppose that we want to minimize the expected shortfall of a portfolio. The key contribution of Rockafellar and Uryasev in their 2000 paper is to introduce the auxiliary function [math]\displaystyle{ F_{\alpha}(w,\gamma) }[/math] for the expected shortfall:[math]\displaystyle{ F_\alpha(w,\gamma) = \gamma + {1\over{1-\alpha}} \int_{\ell(w,x)\geq \gamma} \left[\ell(w,x)-\gamma\right] p(x) \, dx }[/math]Where [math]\displaystyle{ \gamma = \operatorname{VaR}_\alpha(X) }[/math] and [math]\displaystyle{ \ell(w,x) }[/math] is a loss function for a set of portfolio weights [math]\displaystyle{ w\in\mathbb{R}^p }[/math] to be applied to the returns. Rockafellar/Uryasev proved that [math]\displaystyle{ F_\alpha(w,\gamma) }[/math] is convex with respect to [math]\displaystyle{ \gamma }[/math] and is equivalent to the expected shortfall at the minimum point. To numerically compute the expected shortfall for a set of portfolio returns, it is necessary to generate [math]\displaystyle{ J }[/math] simulations of the portfolio constituents; this is often done using copulas. With these simulations in hand, the auxiliary function may be approximated by:[math]\displaystyle{ \widetilde{F}_\alpha(w,\gamma) = \gamma + {1\over{(1-\alpha)J}}\sum_{j=1}^J [\ell(w,x_j) - \gamma]_{+} }[/math]This is equivalent to the formulation:[math]\displaystyle{ \min_{\gamma,z,w} \; \gamma + {1\over{(1-\alpha)J}} \sum_{j=1}^J z_j, \quad \text{s.t. } z_j \geq \ell(w,x_j)-\gamma,\; z_j \geq 0 }[/math] Finally, choosing a linear loss function [math]\displaystyle{ \ell(w,x_{j}) = -w^T x_j }[/math] turns the optimization problem into a linear program. Using standard methods, it is then easy to find the portfolio that minimizes expected shortfall.

Formulas for continuous probability distributions

Closed-form formulas exist for calculating the expected shortfall when the payoff of a portfolio [math]\displaystyle{ X }[/math] or a corresponding loss [math]\displaystyle{ L = -X }[/math] follows a specific continuous distribution. In the former case, the expected shortfall corresponds to the opposite number of the left-tail conditional expectation below [math]\displaystyle{ -\operatorname{VaR}_\alpha (X) }[/math]:

[math]\displaystyle{ \operatorname {ES}_\alpha(X) = E[-X\mid X \leq -\operatorname{VaR}_\alpha(X)] = -\frac{1}{\alpha}\int_0^\alpha \operatorname{VaR}_\gamma(X) \, d\gamma = -\frac{1}{\alpha} \int_{-\infty}^{-\operatorname{VaR}_\alpha(X)} xf(x) \, dx. }[/math]

Typical values of [math]\displaystyle{ \alpha }[/math] in this case are 5% and 1%.

For engineering or actuarial applications it is more common to consider the distribution of losses [math]\displaystyle{ L = -X }[/math], the expected shortfall in this case corresponds to the right-tail conditional expectation above [math]\displaystyle{ \operatorname{VaR}_\alpha (L) }[/math] and the typical values of [math]\displaystyle{ \alpha }[/math] are 95% and 99%:

[math]\displaystyle{ \operatorname {ES}_\alpha(L) = \operatorname E[L\mid L \geq \operatorname{VaR}_\alpha(L)] = \frac{1}{1-\alpha} \int^1_\alpha \operatorname{VaR}_\gamma(L) \, d\gamma = \frac{1}{1-\alpha} \int^{+\infty}_{\operatorname{VaR}_\alpha(L)} yf(y) \, dy. }[/math]

Since some formulas below were derived for the left-tail case and some for the right-tail case, the following reconciliations can be useful:

[math]\displaystyle{ \operatorname {ES}_\alpha(X) = -\frac{1}{\alpha} \operatorname E[X] + \frac{1-\alpha}{\alpha} \operatorname {ES}_\alpha(L) \text{ and } \operatorname{ES}_\alpha(L) = \frac{1}{1-\alpha} \operatorname E[L]+\frac{\alpha}{1-\alpha} \operatorname {ES}_\alpha(X). }[/math]

Normal distribution

If the payoff of a portfolio [math]\displaystyle{ X }[/math] follows the normal (Gaussian) distribution with p.d.f. [math]\displaystyle{ f(x) = \frac{1}{\sqrt{2\pi}\sigma}e^{-\frac{(x-\mu)^2}{2\sigma^2}} }[/math] then the expected shortfall is equal to [math]\displaystyle{ \operatorname{ES}_\alpha(X) = -\mu+\sigma\frac{\varphi(\Phi^{-1}(\alpha))}{\alpha} }[/math], where [math]\displaystyle{ \varphi(x)=\frac{1}{\sqrt{2\pi}}e^{-\frac{x^2}{2}} }[/math] is the standard normal p.d.f., [math]\displaystyle{ \Phi(x) }[/math] is the standard normal c.d.f., so [math]\displaystyle{ \Phi^{-1}(\alpha) }[/math] is the standard normal quantile.[10]

If the loss of a portfolio [math]\displaystyle{ L }[/math] follows the normal distribution, the expected shortfall is equal to [math]\displaystyle{ \operatorname{ES}_\alpha(L) = \mu+\sigma\frac{\varphi(\Phi^{-1}(\alpha))}{1-\alpha} }[/math].[11]

Generalized Student's t-distribution

If the payoff of a portfolio [math]\displaystyle{ X }[/math] follows the generalized Student's t-distribution with p.d.f. [math]\displaystyle{ f(x) = \frac{\Gamma\left(\frac{\nu+1}{2}\right)}{\Gamma\left(\frac{\nu}{2} \right) \sqrt{\pi\nu} \sigma} \left(1+\frac{1}{\nu}\left(\frac{x-\mu}{\sigma}\right)^2\right)^{-\frac{\nu+1}{2}} }[/math] then the expected shortfall is equal to [math]\displaystyle{ \operatorname{ES}_\alpha(X) = -\mu+\sigma\frac{\nu+(\Tau^{-1}(\alpha))^2}{\nu-1}\frac{\tau(\Tau^{-1}(\alpha))}{\alpha} }[/math], where [math]\displaystyle{ \tau(x)=\frac{\Gamma\bigl(\frac{\nu+1}{2}\bigr)}{\Gamma\bigl(\frac{\nu}{2}\bigr)\sqrt{\pi\nu}}\Bigl(1+\frac{x^2}{\nu}\Bigr)^{-\frac{\nu+1}{2}} }[/math] is the standard t-distribution p.d.f., [math]\displaystyle{ \Tau(x) }[/math] is the standard t-distribution c.d.f., so [math]\displaystyle{ \Tau^{-1}(\alpha) }[/math] is the standard t-distribution quantile.[10]

If the loss of a portfolio [math]\displaystyle{ L }[/math] follows generalized Student's t-distribution, the expected shortfall is equal to [math]\displaystyle{ \operatorname{ES}_\alpha(L) = \mu+\sigma\frac{\nu+(\Tau^{-1}(\alpha))^2}{\nu-1}\frac{\tau(\Tau^{-1}(\alpha))}{\alpha} }[/math].[11]

Laplace distribution

If the payoff of a portfolio [math]\displaystyle{ X }[/math] follows the Laplace distribution with the p.d.f.

[math]\displaystyle{ f(x) = \frac{1}{2b}e^{-|x-\mu|/b} }[/math]

and the c.d.f.

[math]\displaystyle{ F(x) = \begin{cases} 1 - \frac{1}{2} e^{-(x-\mu)/b} & \text{if }x \geq \mu,\\[4pt] \frac{1}{2} e^{(x-\mu)/b} & \text{if }x \lt \mu. \end{cases} }[/math]

then the expected shortfall is equal to [math]\displaystyle{ \operatorname{ES}_\alpha(X) = -\mu + b(1 - \ln 2\alpha) }[/math] for [math]\displaystyle{ \alpha \le 0.5 }[/math].[10]

If the loss of a portfolio [math]\displaystyle{ L }[/math] follows the Laplace distribution, the expected shortfall is equal to[11]

[math]\displaystyle{ \operatorname{ES}_\alpha(L) = \begin{cases} \mu + b \frac{\alpha}{1-\alpha} (1-\ln2\alpha) & \text{if }\alpha \lt 0.5,\\[4pt] \mu + b[1 - \ln(2(1-\alpha))] & \text{if }\alpha \ge 0.5. \end{cases} }[/math]

Logistic distribution

If the payoff of a portfolio [math]\displaystyle{ X }[/math] follows the logistic distribution with p.d.f. [math]\displaystyle{ f(x) = \frac{1}{s} e^{-\frac{x-\mu}{s}}\left(1+e^{-\frac{x-\mu}{s}}\right)^{-2} }[/math] and the c.d.f. [math]\displaystyle{ F(x) = \left(1+e^{-\frac{x-\mu}{s}}\right)^{-1} }[/math] then the expected shortfall is equal to [math]\displaystyle{ \operatorname{ES}_\alpha(X) = -\mu + s \ln\frac{(1-\alpha)^{1-\frac{1}{\alpha}}}{\alpha} }[/math].[10]

If the loss of a portfolio [math]\displaystyle{ L }[/math] follows the logistic distribution, the expected shortfall is equal to [math]\displaystyle{ \operatorname{ES}_\alpha(L) = \mu + s\frac{-\alpha\ln\alpha-(1-\alpha)\ln(1-\alpha)}{1-\alpha} }[/math].[11]

Exponential distribution

If the loss of a portfolio [math]\displaystyle{ L }[/math] follows the exponential distribution with p.d.f. [math]\displaystyle{ f(x) = \begin{cases}\lambda e^{-\lambda x} & \text{if }x \geq 0,\\ 0 & \text{if }x \lt 0.\end{cases} }[/math] and the c.d.f. [math]\displaystyle{ F(x) = \begin{cases}1 - e^{-\lambda x} & \text{if }x \geq 0,\\ 0 & \text{if }x \lt 0.\end{cases} }[/math] then the expected shortfall is equal to [math]\displaystyle{ \operatorname{ES}_\alpha(L) = \frac{-\ln(1-\alpha)+1}{\lambda} }[/math].[11]

Pareto distribution

If the loss of a portfolio [math]\displaystyle{ L }[/math] follows the Pareto distribution with p.d.f. [math]\displaystyle{ f(x) = \begin{cases} \frac{a x_m^a}{x^{a+1}} & \text{if }x \geq x_m,\\ 0 & \text{if }x \lt x_m. \end{cases} }[/math] and the c.d.f. [math]\displaystyle{ F(x) = \begin{cases} 1 - (x_m/x)^a & \text{if }x \geq x_m,\\ 0 & \text{if }x \lt x_m. \end{cases} }[/math] then the expected shortfall is equal to [math]\displaystyle{ \operatorname{ES}_\alpha(L) = \frac{x_m a}{(1-\alpha)^{1/a}(a-1)} }[/math].[11]

Generalized Pareto distribution (GPD)

If the loss of a portfolio [math]\displaystyle{ L }[/math] follows the GPD with p.d.f.

[math]\displaystyle{ f(x) = \frac{1}{s} \left( 1+\frac{\xi (x-\mu)}{s} \right)^{\left(-\frac{1}{\xi}-1\right)} }[/math]

and the c.d.f.

[math]\displaystyle{ F(x) = \begin{cases} 1 - \left(1+\frac{\xi(x-\mu)}{s}\right)^{-1 /\xi} & \text{if }\xi \ne 0,\\ 1-\exp \left( -\frac{x-\mu}{s} \right) & \text{if }\xi = 0. \end{cases} }[/math]

then the expected shortfall is equal to

[math]\displaystyle{ \operatorname{ES}_\alpha(L) = \begin{cases} \mu + s \left[ \frac{(1-\alpha)^{-\xi}}{1-\xi}+\frac{(1-\alpha)^{-\xi}-1}{\xi} \right] & \text{if }\xi \ne 0,\\ \mu + s \left[1 - \ln(1-\alpha) \right] & \text{if }\xi = 0, \end{cases} }[/math]

and the VaR is equal to[11]

[math]\displaystyle{ \operatorname{VaR}_\alpha(L) = \begin{cases} \mu + s \frac{(1-\alpha)^{-\xi}-1}{\xi} & \text{if }\xi \ne 0,\\ \mu - s \ln(1-\alpha) & \text{if }\xi = 0. \end{cases} }[/math]

Weibull distribution

If the loss of a portfolio [math]\displaystyle{ L }[/math] follows the Weibull distribution with p.d.f. [math]\displaystyle{ f(x) = \begin{cases} \frac{k}{\lambda} \left(\frac{x}{\lambda}\right)^{k-1} e^{-(x/\lambda)^k} & \text{if }x \geq 0,\\ 0 & \text{if }x \lt 0. \end{cases} }[/math] and the c.d.f. [math]\displaystyle{ F(x) = \begin{cases} 1 - e^{-(x/\lambda)^k} & \text{if }x \geq 0,\\ 0 & \text{if }x \lt 0. \end{cases} }[/math] then the expected shortfall is equal to [math]\displaystyle{ \operatorname{ES}_\alpha(L) = \frac{\lambda}{1-\alpha} \Gamma\left(1+\frac{1}{k},-\ln(1-\alpha)\right) }[/math], where [math]\displaystyle{ \Gamma(s,x) }[/math] is the upper incomplete gamma function.[11]

Generalized extreme value distribution (GEV)

If the payoff of a portfolio [math]\displaystyle{ X }[/math] follows the GEV with p.d.f. [math]\displaystyle{ f(x) = \begin{cases} \frac{1}{\sigma} \left( 1+\xi \frac{ x-\mu}{\sigma} \right)^{-\frac{1}{\xi}-1} \exp\left[-\left( 1 + \xi \frac{x-\mu}{\sigma} \right)^{-{1}/{\xi}}\right] & \text{if } \xi \ne 0,\\ \frac{1}{\sigma}e^{-\frac{x-\mu}{\sigma}}e^{-e^{-\frac{x-\mu}{\sigma}}} & \text{if } \xi = 0. \end{cases} }[/math] and c.d.f. [math]\displaystyle{ F(x) = \begin{cases} \exp\left(-\left(1+\xi\frac{x-\mu}{\sigma}\right)^{-{1}/{\xi}}\right) & \text{if }\xi \ne 0,\\ \exp\left(-e^{-\frac{x-\mu}{\sigma}}\right) & \text{if }\xi = 0. \end{cases} }[/math] then the expected shortfall is equal to [math]\displaystyle{ \operatorname{ES}_\alpha(X) = \begin{cases} -\mu - \frac{\sigma}{\alpha \xi} \big[ \Gamma(1-\xi,-\ln\alpha)-\alpha \big] & \text{if }\xi \ne 0,\\ -\mu - \frac{\sigma}{\alpha} \big[ \text{li}(\alpha) - \alpha \ln(-\ln \alpha) \big] & \text{if }\xi = 0. \end{cases} }[/math] and the VaR is equal to [math]\displaystyle{ \operatorname{VaR}_\alpha(X) = \begin{cases} -\mu - \frac{\sigma}{\xi} \left[(-\ln \alpha)^{-\xi}-1 \right] & \text{if }\xi \ne 0,\\ -\mu + \sigma \ln(-\ln\alpha) & \text{if }\xi = 0. \end{cases} }[/math], where [math]\displaystyle{ \Gamma(s,x) }[/math] is the upper incomplete gamma function, [math]\displaystyle{ \mathrm{li}(x) = \int \frac{dx}{\ln x} }[/math] is the logarithmic integral function.[12]

If the loss of a portfolio [math]\displaystyle{ L }[/math] follows the GEV, then the expected shortfall is equal to [math]\displaystyle{ \operatorname{ES}_\alpha(X) = \begin{cases} \mu + \frac{\sigma}{(1-\alpha) \xi} \bigl[ \gamma(1-\xi,-\ln\alpha)-(1-\alpha) \bigr] & \text{if }\xi \ne 0,\\ \mu + \frac{\sigma}{1-\alpha} \bigl[y - \text{li}(\alpha) + \alpha \ln(-\ln \alpha) \bigr] & \text{if }\xi = 0. \end{cases} }[/math], where [math]\displaystyle{ \gamma(s,x) }[/math] is the lower incomplete gamma function, [math]\displaystyle{ y }[/math] is the Euler-Mascheroni constant.[11]

Generalized hyperbolic secant (GHS) distribution

If the payoff of a portfolio [math]\displaystyle{ X }[/math] follows the GHS distribution with p.d.f. [math]\displaystyle{ f(x) = \frac{1}{2 \sigma} \operatorname{sech}\left(\frac{\pi}{2} \frac{x-\mu}{\sigma}\right) }[/math]and the c.d.f. [math]\displaystyle{ F(x) = \frac{2}{\pi}\arctan\left[\exp\left(\frac{\pi}{2}\frac{x-\mu}{\sigma}\right)\right] }[/math] then the expected shortfall is equal to [math]\displaystyle{ \operatorname{ES}_\alpha(X) = - \mu - \frac{2\sigma}{\pi} \ln\left( \tan \frac{\pi\alpha}{2} \right) - \frac{2\sigma}{\pi^2\alpha}i\left[\operatorname{Li}_2\left(-i\tan\frac{\pi\alpha}{2}\right)-\operatorname{Li}_2\left(i\tan\frac{\pi\alpha}{2}\right)\right] }[/math], where [math]\displaystyle{ \operatorname{Li}_2 }[/math] is the dilogarithm and [math]\displaystyle{ i=\sqrt{-1} }[/math] is the imaginary unit.[12]

Johnson's SU-distribution

If the payoff of a portfolio [math]\displaystyle{ X }[/math] follows Johnson's SU-distribution with the c.d.f. [math]\displaystyle{ F(x) = \Phi\left[\gamma+\delta\sinh^{-1}\left(\frac{x-\xi}{\lambda}\right)\right] }[/math] then the expected shortfall is equal to [math]\displaystyle{ \operatorname{ES}_\alpha(X) = -\xi - \frac{\lambda}{2\alpha} \left[ \exp\left(\frac{1-2\gamma\delta}{2\delta^2}\right) \; \Phi\left(\Phi^{-1}(\alpha)-\frac{1}{\delta}\right) - \exp\left(\frac{1+2\gamma\delta}{2\delta^2}\right) \; \Phi\left(\Phi^{-1}(\alpha)+\frac{1}{\delta}\right) \right] }[/math], where [math]\displaystyle{ \Phi }[/math] is the c.d.f. of the standard normal distribution.[13]

Burr type XII distribution

If the payoff of a portfolio [math]\displaystyle{ X }[/math] follows the Burr type XII distribution the p.d.f. [math]\displaystyle{ f(x) = \frac{ck}{\beta} \left(\frac{x-\gamma}{\beta}\right)^{c-1} \left[1+\left(\frac{x-\gamma}{\beta} \right)^c\right]^{-k-1} }[/math] and the c.d.f. [math]\displaystyle{ F(x) = 1-\left[1+\left(\frac{x-\gamma}{\beta} \right)^c \right]^{-k} }[/math], the expected shortfall is equal to [math]\displaystyle{ \operatorname{ES}_\alpha(X) = - \gamma - \frac{\beta}{\alpha} \left( (1-\alpha)^{-1/k}-1 \right)^{1/c} \left[ \alpha -1+{_2F_1}\left(\frac{1}{c},k;1+\frac{1}{c};1-(1-\alpha)^{-1/k}\right) \right] }[/math], where [math]\displaystyle{ _2F_1 }[/math] is the hypergeometric function. Alternatively, [math]\displaystyle{ \operatorname{ES}_\alpha(X) = - \gamma - \frac{\beta}{\alpha} \frac{ck}{c+1} \left( (1-\alpha)^{-1/k}-1 \right)^{1+\frac{1}{c}} {_2F_1}\left(1+\frac{1}{c}, k+1;2+\frac{1}{c};1-(1-\alpha)^{-1/k}\right) }[/math].[12]

Dagum distribution

If the payoff of a portfolio [math]\displaystyle{ X }[/math] follows the Dagum distribution with p.d.f. [math]\displaystyle{ f(x) = \frac{ck}{\beta} \left(\frac{x-\gamma}{\beta}\right)^{ck-1} \left[1+\left(\frac{x-\gamma}{\beta}\right)^c\right]^{-k-1} }[/math] and the c.d.f. [math]\displaystyle{ F(x) = \left[1+\left(\frac{x-\gamma}{\beta}\right)^{-c}\right]^{-k} }[/math], the expected shortfall is equal to [math]\displaystyle{ \operatorname{ES}_\alpha(X) = - \gamma - \frac{\beta}{\alpha} \frac{ck}{ck+1} \left( \alpha^{-1/k}-1 \right)^{-k-\frac{1}{c}} {_2F_1}\left(k+1,k+\frac{1}{c};k+1+\frac{1}{c};-\frac{1}{\alpha^{-1/k}-1}\right) }[/math], where [math]\displaystyle{ _2F_1 }[/math] is the hypergeometric function.[12]

Lognormal distribution

If the payoff of a portfolio [math]\displaystyle{ X }[/math] follows lognormal distribution, i.e. the random variable [math]\displaystyle{ \ln(1+X) }[/math] follows the normal distribution with p.d.f. [math]\displaystyle{ f(x) = \frac{1}{\sqrt{2\pi}\sigma}e^{-\frac{(x-\mu)^2}{2\sigma^2}} }[/math], then the expected shortfall is equal to [math]\displaystyle{ \operatorname{ES}_\alpha(X) = 1 - \exp\left(\mu+\frac{\sigma^2}{2}\right) \frac{\Phi\left(\Phi^{-1}(\alpha)-\sigma\right)}{\alpha} }[/math], where [math]\displaystyle{ \Phi(x) }[/math] is the standard normal c.d.f., so [math]\displaystyle{ \Phi^{-1}(\alpha) }[/math] is the standard normal quantile.[14]

Log-logistic distribution

If the payoff of a portfolio [math]\displaystyle{ X }[/math] follows log-logistic distribution, i.e. the random variable [math]\displaystyle{ \ln(1+X) }[/math] follows the logistic distribution with p.d.f. [math]\displaystyle{ f(x) = \frac{1}{s} e^{-\frac{x-\mu}{s}} \left(1+e^{-\frac{x-\mu}{s}}\right)^{-2} }[/math], then the expected shortfall is equal to [math]\displaystyle{ \operatorname{ES}_\alpha(X) = 1-\frac{e^\mu}{\alpha}I_\alpha(1+s,1-s)\frac{\pi s}{\sin\pi s} }[/math], where [math]\displaystyle{ I_\alpha }[/math] is the regularized incomplete beta function, [math]\displaystyle{ I_\alpha(a,b)=\frac{\Beta_\alpha(a,b)}{\Beta(a,b)} }[/math].

As the incomplete beta function is defined only for positive arguments, for a more generic case the expected shortfall can be expressed with the hypergeometric function: [math]\displaystyle{ \operatorname{ES}_\alpha(X) = 1-\frac{e^\mu \alpha^s}{s+1} {_2F_1}(s,s+1;s+2;\alpha) }[/math].[14]

If the loss of a portfolio [math]\displaystyle{ L }[/math] follows log-logistic distribution with p.d.f. [math]\displaystyle{ f(x) = \frac{\frac{b}{a}(x/a)^{b-1}}{(1+(x/a)^b)^2} }[/math] and c.d.f. [math]\displaystyle{ F(x) = \frac{1}{1+(x/a)^{-b}} }[/math], then the expected shortfall is equal to [math]\displaystyle{ \operatorname{ES}_\alpha(L) = \frac{a}{1-\alpha} \left[ \frac{\pi}{b} \csc\left(\frac{\pi}{b}\right) - \Beta_\alpha \left(\frac{1}{b}+1,1-\frac{1}{b}\right) \right] }[/math], where [math]\displaystyle{ B_\alpha }[/math] is the incomplete beta function.[11]

Log-Laplace distribution

If the payoff of a portfolio [math]\displaystyle{ X }[/math] follows log-Laplace distribution, i.e. the random variable [math]\displaystyle{ \ln(1+X) }[/math] follows the Laplace distribution the p.d.f. [math]\displaystyle{ f(x) = \frac{1}{2b}e^{-\frac{|x-\mu|}{b}} }[/math], then the expected shortfall is equal to

[math]\displaystyle{ \operatorname{ES}_\alpha(X) = \begin{cases} 1 - \frac{e^\mu (2\alpha)^b}{b+1} & \text{if }\alpha \le 0.5,\\ 1 - \frac{e^\mu 2^{-b}}{\alpha(b-1)} \left[(1-\alpha)^{(1-b)}-1\right] & \text{if } \alpha \gt 0.5. \end{cases} }[/math][14]

Log-generalized hyperbolic secant (log-GHS) distribution

If the payoff of a portfolio [math]\displaystyle{ X }[/math] follows log-GHS distribution, i.e. the random variable [math]\displaystyle{ \ln(1+X) }[/math] follows the GHS distribution with p.d.f. [math]\displaystyle{ f(x) = \frac{1}{2 \sigma} \operatorname{sech} \left(\frac{\pi}{2}\frac{x-\mu}{\sigma}\right) }[/math], then the expected shortfall is equal to

[math]\displaystyle{ \operatorname{ES}_\alpha(X) = 1 - \frac{1}{\alpha(\sigma+{\pi/2})} \left(\tan\frac{\pi \alpha}{2}\exp\frac{\pi \mu}{2\sigma}\right)^{2\sigma/\pi} \tan\frac{\pi \alpha}{2} {_2F_1}\left(1,\frac{1}{2}+\frac{\sigma}{\pi};\frac{3}{2}+\frac{\sigma}{\pi};-\tan\left(\frac{\pi \alpha}{2}\right)^2\right), }[/math]

where [math]\displaystyle{ _2F_1 }[/math] is the hypergeometric function.[14]

Dynamic expected shortfall

The conditional version of the expected shortfall at the time t is defined by

[math]\displaystyle{ \operatorname{ES}_\alpha^t(X) = \operatorname{ess\sup}_{Q \in \mathcal{Q}_{\alpha}^t} E^Q[-X \mid \mathcal{F}_t] }[/math]

where [math]\displaystyle{ \mathcal{Q}_{\alpha}^t = \left\{Q = P\,\vert_{\mathcal{F}_t}: \frac{dQ}{dP} \leq \alpha_t^{-1} \text{ a.s.}\right\} }[/math].[15][16]

This is not a time-consistent risk measure. The time-consistent version is given by

[math]\displaystyle{ \rho_{\alpha}^t(X) = \operatorname{ess\sup}_{Q \in \tilde{\mathcal{Q}}_{\alpha}^t} E^Q[-X\mid\mathcal{F}_t] }[/math]

such that[17]

[math]\displaystyle{ \tilde{\mathcal{Q}}_{\alpha}^t = \left\{Q \ll P: \operatorname{E}\left[\frac{dQ}{dP} \mid \mathcal{F}_{\tau+1} \right] \leq \alpha_t^{-1} \operatorname{E}\left[\frac{dQ}{dP} \mid \mathcal{F}_{\tau}\right] \; \forall \tau \geq t \text{ a.s.}\right\}. }[/math]

See also

Methods of statistical estimation of VaR and ES can be found in Embrechts et al.[18] and Novak.[19] When forecasting VaR and ES, or optimizing portfolios to minimize tail risk, it is important to account for asymmetric dependence and non-normalities in the distribution of stock returns such as auto-regression, asymmetric volatility, skewness, and kurtosis.[20]

References

  1. Rockafellar, R. Tyrrell; Uryasev, Stanislav (2000). "Optimization of conditional value-at-risk". Journal of Risk 2 (3): 21–42. doi:10.21314/JOR.2000.038. https://sites.math.washington.edu/~rtr/papers/rtr179-CVaR1.pdf. 
  2. Rockafellar, R. Tyrrell; Royset, Johannes (2010). "On Buffered Failure Probability in Design and Optimization of Structures". Reliability Engineering and System Safety 95 (5): 499–510. doi:10.1016/j.ress.2010.01.001. https://sites.math.washington.edu/~rtr/papers/rtr211-BufferedProb.pdf. 
  3. Carlo Acerbi; Dirk Tasche (2002). "Expected Shortfall: a natural coherent alternative to Value at Risk". Economic Notes 31 (2): 379–388. doi:10.1111/1468-0300.00091. http://www.bis.org/bcbs/ca/acertasc.pdf. Retrieved April 25, 2012. 
  4. Föllmer, H.; Schied, A. (2008). Convex and coherent risk measures. http://wws.mathematik.hu-berlin.de/~foellmer/papers/CCRM.pdf. Retrieved October 4, 2011. 
  5. Patrick Cheridito; Tianhui Li (2008). "Dual characterization of properties of risk measures on Orlicz hearts". Mathematics and Financial Economics 2: 2–29. doi:10.1007/s11579-008-0013-7. 
  6. "Average Value at Risk". https://statistik.ets.kit.edu/download/doc_secure1/7_StochModels.pdf. 
  7. Julia L. Wirch; Mary R. Hardy. "Distortion Risk Measures: Coherence and Stochastic Dominance". http://pascal.iseg.utl.pt/~cemapre/ime2002/main_page/papers/JuliaWirch.pdf. 
  8. Balbás, A.; Garrido, J.; Mayoral, S. (2008). "Properties of Distortion Risk Measures". Methodology and Computing in Applied Probability 11 (3): 385. doi:10.1007/s11009-008-9089-z. http://e-archivo.uc3m.es/bitstream/10016/14071/1/properties_balbas_MCAP_2009_ps.pdf. 
  9. Rockafellar, R. Tyrrell; Uryasev, Stanislav (2000). "Optimization of conditional value-at-risk". Journal of Risk 2 (3): 21–42. doi:10.21314/JOR.2000.038. https://sites.math.washington.edu/~rtr/papers/rtr179-CVaR1.pdf. 
  10. 10.0 10.1 10.2 10.3 Khokhlov, Valentyn (2016). "Conditional Value-at-Risk for Elliptical Distributions". Evropský časopis Ekonomiky a Managementu 2 (6): 70–79. 
  11. 11.0 11.1 11.2 11.3 11.4 11.5 11.6 11.7 11.8 11.9 Norton, Matthew; Khokhlov, Valentyn; Uryasev, Stan (2018-11-27). "Calculating CVaR and bPOE for Common Probability Distributions With Application to Portfolio Optimization and Density Estimation". arXiv:1811.11301 [q-fin.RM].
  12. 12.0 12.1 12.2 12.3 Khokhlov, Valentyn (2018-06-21). Conditional Value-at-Risk for Uncommon Distributions. doi:10.2139/ssrn.3200629. 
  13. Stucchi, Patrizia (2011-05-31). Moment-Based CVaR Estimation: Quasi-Closed Formulas. doi:10.2139/ssrn.1855986. 
  14. 14.0 14.1 14.2 14.3 Template:Cite SSRN
  15. Detlefsen, Kai; Scandolo, Giacomo (2005). "Conditional and dynamic convex risk measures". Finance Stoch. 9 (4): 539–561. doi:10.1007/s00780-005-0159-6. http://www.dmd.unifi.it/scandolo/pdf/Scandolo-Detlefsen-05.pdf. Retrieved October 11, 2011. [|permanent dead link|dead link}}]
  16. Acciaio, Beatrice; Penner, Irina (2011). Dynamic convex risk measures. http://wws.mathematik.hu-berlin.de/~penner/Acciaio_Penner.pdf. Retrieved October 11, 2011. 
  17. Cheridito, Patrick; Kupper, Michael (May 2010). "Composition of time-consistent dynamic monetary risk measures in discrete time". International Journal of Theoretical and Applied Finance. http://wws.mathematik.hu-berlin.de/~kupper/papers/comp2010.pdf. Retrieved February 4, 2011. 
  18. Embrechts P., Kluppelberg C. and Mikosch T., Modelling Extremal Events for Insurance and Finance. Springer (1997).
  19. Novak S.Y., Extreme value methods with applications to finance. Chapman & Hall/CRC Press (2011). ISBN:978-1-4398-3574-6.
  20. Low, R.K.Y.; Alcock, J.; Faff, R.; Brailsford, T. (2013). "Canonical vine copulas in the context of modern portfolio management: Are they worth it?". Journal of Banking & Finance 37 (8): 3085–3099. doi:10.1016/j.jbankfin.2013.02.036. https://espace.library.uq.edu.au/view/UQ:297895/EC15UQ297895.pdf. 

External links