This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Extended Gini-type measures of risk and variability

\nameMohammed Berkhouch111Corresponding author. Phone number: (+212)6 10 88 46 86. LGII, ENSA, Ibn Zohr University, Agadir, Morocco. E-mail: mohammed.berkhouch@edu.uiz.ac.ma    \nameGhizlane Lakhnati LGII, ENSA, Ibn Zohr University, Agadir, Morocco. E-mail: g.lakhnati@uiz.ac.ma    \nameMarcelo Brutti Righi Federal University of Rio Grande do Sul, Porto Alegre, Brazil. E-mail: marcelo.righi@ufrgs.br
(August 14, 2025)
Abstract

The aim of this paper is to introduce a risk measure that extends the Gini-type measures of risk and variability, the Extended Gini Shortfall, by taking risk aversion into consideration. Our risk measure is coherent and catches variability, an important concept for risk management. The analysis is made under the Choquet integral representations framework. We expose results for analytic computation under well-known distribution functions. Furthermore, we provide a practical application.

JEL classification: C6, G10

keywords:
risk measures, variability measures, risk aversion, signed Choquet integral, Extended Gini Shortfall.

1 Introduction

In modern risk management, a large number of risk measures have been proposed in the literature. These measures are mappings from a set of random variables (financial losses) to real numbers. At first, the focus were on the variability over an expected return, as is the case for the well-known variance. After the collapses and crises in financial systems, a prominent trend associated with tail-based risk measures has emerged, especially with the most popular ones nowadays: the Value-at-Risk (VaR) and the Expected Shortfall (ES). However, this kind of risk measures does not capture the variability of a financial position, a primitive but relevant concept. In order to solution this issue, some authors have proposed and studied specific examples of risk measures.

In this sense, Fischer [8] considered combining the mean and semi-deviations. Regarding tail risk, Furman and Landsman [10] proposed a measure that weighs the mean and standard deviation in the truncated tail by VaR, while Righi and Ceretta [17] considered penalizing the ES by the dispersion of losses exceeding it. From a practical perspective, Righi and Borenstein [18] explored this concept, calling the approach as loss-deviation, for portfolio optimization. In a more general fashion, Righi [19] presents results and examples about compositions of risk and variability measures in order to ensure solid theoretical properties.

Recently, Furman et al. [11] introduced the Gini Shortfall (GS) risk measure which is coherent and satisfies co-monotonic additivity. GS is a composition between ES and tail based Gini coefficient. However, GS supposes that all individuals have the same attitude towards risk, while agents differ in the way they take personal decisions that involve risk because of discrepancies in their risk aversion. To incorporate such psychological behavior in tail risk analysis, we introduce a generalized version of the GS. This risk measure, called Extended Gini Shortfall (EGS), captures the notion of variability, satisfies the co-monotonic additivity property, and it is coherent under a necessary and sufficient condition for its loading parameter. The consideration of the decision-maker risk aversion, joined to these properties, is in consonance to what agents seek when searching for a suitable measure of risk. The approach followed in this article leads us to a new family of spectral risk measures, proposed by Acerbi [1], with an attractive weighting function.

In this sense, we discuss, in a separated manner, the properties from the variability term and our composed risk measure. Moreover, we discuss in details the role of each parameter in the mentioned weighting function. Furthermore, we expose results on analytic formulations for computation of EGS under known distribution functions. Our focus in this paper is on theoretical results, but this approach gives rise to further forthcoming investigations. In this sense, risk forecasting of our new family of risk measures is a subject that deserves a further and separate survey which will be made in a forthcoming work.

The rest of the paper is organized as follows. In Section 2, we present and discuss some preliminaries such as essential properties of measures of risk and variability, as well we elucidate the role of the signed Choquet integral. In section 3, we start with the Classical and Extended Gini functionals in order to introduce the concept and explore properties about what we call Tail Extended Gini functional and Extended Gini Shortfall. In Section 4, we give the closed-form of our risk measures class for elliptical distributions and then derive the uniform, Normal and Student-t cases. Section 5 illustrates an application of the introduced risk measures class in practice.

2 Preliminaries

We first introduce some basic notation. Let (Ω,𝒜,)(\Omega,\mathcal{A},\mathbb{P}) be an atomless probability space. All equations and inequalities are in the \mathbb{P} almost surely sense. Let LqL^{q}, q[0,)q\in[0,\infty), denote the set of all random variables (rv’s from now on) in (Ω,𝒜,)(\Omega,\mathcal{A},\mathbb{P}) with finite q-th moment and LL^{\infty} be the set of all essentially bounded rv’s. Throughout this paper, XL0X\in L^{0} is a rv modeling financial losses (profits) when it has positive (negative) values. For every XL0,FXX\in L^{0},\,F_{X} denotes the cdf of XX, and UXU_{X} denote any uniform [0,1][0,1] rv such that the equation FX1(UX)=XF^{-1}_{X}(U_{X})=X holds. The existence of such rv’s is assured in Rüschendorf ([21], Proposition 1.3). We denote xpx_{p} as the p-quantile of XX. Two rv’s XX and YY are co-monotonic when (X(ω)X(ω))(Y(ω)Y(ω))0 for (ω,ω)Ω×Ω(×)-almost surely(X(\omega)-X(\omega^{\prime}))(Y(\omega)-Y(\omega^{\prime}))\geq 0\text{ for }(\omega,\omega^{\prime})\in\Omega\times\Omega\quad(\mathbb{P}\times\mathbb{P})\text{-almost surely}. Throughout the present paper, we deal with several convex cones 𝒳\mathcal{X} of rv’s, of which 𝒳=L1\mathcal{X}=L^{1} is of particular importance and LL^{\infty} is always contained in 𝒳\mathcal{X}.

We begin by exposing definitions of both risk and variability measures. We assume throughout the paper that all functionals respect the following property, which is essential in order to obtain a functional directly from its distribution function.

Definition 2.1.

A functional ff is said to be law invariant if it fulfills the following property:

(A) Law Invariance: If X𝒳X\in\mathcal{X} and Y𝒳Y\in\mathcal{X} have the same distributions under \mathbb{P}, succinctly X=dYX\stackrel{{\scriptstyle d}}{{=}}Y, then f(X)=f(Y)f(X)=f(Y).

Definition 2.2.

A risk measure is a functional ρ:𝒳(,]\rho:\mathcal{X}\rightarrow(-\infty,\infty], which may fulfills the following properties:

(B1) Monotonicity: ρ(X)ρ(Y)\rho(X)\leq\rho(Y) when X,Y𝒳X,Y\in\mathcal{X} are such that XYX\leq Y \mathbb{P}-almost surely.

(B2) Translation invariance: ρ(X+m)=ρ(X)+m\rho(X+m)=\rho(X)+m for all mm\in\mathbb{R} and X𝒳X\in\mathcal{X}.

(A1) Positive homogeneity: ρ(λX)=λρ(X)\rho(\lambda X)=\lambda\rho(X) for all λ>0\lambda>0 and X𝒳X\in\mathcal{X}.

(A2) Subadditivity: ρ(X+Y)ρ(X)+ρ(Y)\rho(X+Y)\leq\rho(X)+\rho(Y) for all X,Y𝒳X,Y\in\mathcal{X}.

(A3) Co-monotonic additivity: ρ(X+Y)=ρ(X)+ρ(Y)\rho(X+Y)=\rho(X)+\rho(Y) for every co-monotonic pair X,Y𝒳X,Y\in\mathcal{X}.

A risk measure is monetary if it satisfies properties (B1) and (B2), and it is coherent if it satisfies furthermore (A1) and (A2).

Remark 2.3.

For interpretations of these properties, we refer the reader to Föllmer and Schied ([9], Chap 4), Delbaen [6], and McNeil et al. [16]. For example, both functionals VaR and ES are monetary and co-monotonically additive, whereas ES is even coherent.

Definition 2.4.

A functional ν:𝒳[0,]\nu:\mathcal{X}\rightarrow[0,\infty] is a measure of variability, which may fulfills the following properties222Inspired from the deviation measures notion of Rockafellar et al. [20].:

(C1) Standardization: ν(c)=0\nu(c)=0 for all cc\in\mathbb{R}.

(C2) Location invariance: ν(X+c)=ν(X)\nu(X+c)=\nu(X) for all cc\in\mathbb{R} and X𝒳X\in\mathcal{X}.

(A1) Positive homogeneity: ν(λX)=λν(X)\nu(\lambda X)=\lambda\nu(X) for all λ>0\lambda>0 and X𝒳X\in\mathcal{X}.

(A2) Subadditivity: ν(X+Y)ν(X)+ν(Y)\nu(X+Y)\leq\nu(X)+\nu(Y) for all X,Y𝒳X,Y\in\mathcal{X}

(A3) Co-monotonic additivity: ν(X+Y)=ν(X)+ν(Y)\nu(X+Y)=\nu(X)+\nu(Y) for every co-monotonic pair X,Y𝒳X,Y\in\mathcal{X}.

A measure of variability is coherent if it further satisfies (C1), (C2), (A1) and (A2).

Remark 2.5.

For instance, the most classical measures of variability are the Variance and the Standard Deviation. The variance functional satisfies properties (A), (C1), (C2) but not (A1) or (A2), hence it is not coherent. On the other side, the standard deviation functional, since satisfying all aforementioned properties, is coherent. Neither the variance nor the standard deviation is co-monotonically additive.

The notion of signed Choquet integral plays a pivotal role thereafter. It originates from Choquet [4], in the framework of capacities, and is further characterized and studied in decision theory by Schmeidler ([22], [23]).

Definition 2.6.

A function h:[0,1]h:[0,1]\rightarrow\mathbb{R} is called a distortion function when it is non-decreasing and satisfies the boundary conditions h(0)=0h(0)=0 and h(1)=1h(1)=1. Let hh be a distortion function, the functional defined by the equation:

I(X)=0(1h(FX(x)))𝑑x0h(FX(x))𝑑xI(X)=\int_{0}^{\infty}(1-h(F_{X}(x)))dx-\int_{-\infty}^{0}h(F_{X}(x))dx (2.1)

for all X𝒳X\in\mathcal{X} is called the (increasing) Choquet integral. Whenever h:[0,1]h:[0,1]\rightarrow\mathbb{R} is of finite variation, II is called the signed Choquet integral.

Remark 2.7.

When hh is right-continuous, then equation(2.1) can be rewritten as (Wang et al. [25]):

I(X)=01FX1(t)𝑑h(t).I(X)=\int_{0}^{1}F^{-1}_{X}(t)dh(t). (2.2)

Furthermore, when hh is absolutely continuous, with ϕ\phi a function such that dh(t)=ϕ(t)dtdh(t)=\phi(t)dt, then equation(2.2) becomes:

I(X)=01FX1(t)ϕ(t)𝑑t.I(X)=\int_{0}^{1}F^{-1}_{X}(t)\phi(t)dt. (2.3)

In this case, ϕ\phi is called the weighting functional of the signed Choquet integral II.

Remark 2.8.

The signed Choquet integral is co-monotonically additive, as we can readily see from representation (2.2) (cf. Schmeidler [22]). Moreover we know from Yaari [26] and Föllmer and Schied ([9], Theorem 4.88), that any law-invariant risk measure is co-monotonically additive and monetary if and only if it can be represented as a Choquet integral. Furthermore, the functional II is sub-additive if and only if the function hh is convex (cf. Yaari [26] and Acerbi [1]). Moreover, as proved in Furman et al. [11], regarding the weighting functional ϕ\phi, the integral is monotone if and only if ϕ0\phi\geq 0 on [0,1][0,1] and it is sub-additive if and only if ϕ\phi is non-decreasing on [0,1][0,1]. The major difference between a (an increasing) Choquet integral and a signed one is that the latter, being more general, is not necessarily monotone.

One of the practical and theoretical reasons for what we are particularly interested in signed Choquet integral is that we know that a suitable risk measure should be monotone as argued by Artzner et al. [2], but this issue is irrelevant for a measure of variability. In other words, signed Choquet integral is relevant as long as a measure of variability is concerned. The following theorem is enunciated with a complete proof in Furman et al. [11], it gives the characterization for co-monotonically additive and coherent measures of variability.

Theorem 2.9.

Let ν:Lq\nu:L^{q}\rightarrow\mathbb{R} be any LqL^{q}-continuous functional. The following three statements are equivalent:

(i) ν\nu is a co-monotonically additive and coherent measure of variability.

(ii) There is a convex function h:[0,1],h(0)=h(1)=0h:[0,1]\rightarrow\mathbb{R},h(0)=h(1)=0, such that

ν(X)=01FX1(u)𝑑h(u),XLq.\nu(X)=\int_{0}^{1}F^{-1}_{X}(u)dh(u),\quad X\in L^{q}. (2.4)

(iii) There is a non-decreasing function g:[0,1]g:[0,1]\rightarrow\mathbb{R} such that

ν(X)=Cov[X,g(UX)],XLq.\nu(X)=Cov[X,g(U_{X})],\quad X\in L^{q}. (2.5)

Next, we recall a few partial orders of variability that have been popular in economics, insurance, finance and probability theory:

Definition 2.10.

For X,YL1X,Y\in L^{1}, XX is second-order stochastically dominated (SSD) by YY, succinctly XSSDYX\prec_{SSD}Y, if 𝔼[f(X)]𝔼[f(Y)]\mathbb{E}[f(X)]\leq\mathbb{E}[f(Y)] for all increasing convex functions ff.
If in addition, 𝔼[X]=𝔼[Y]\mathbb{E}[X]=\mathbb{E}[Y], then we say that XX is smaller than YY in convex order, succinctly XCXYX\prec_{CX}Y.333 We say equally YY is a Mean Preserving Spread of XX, succinctly YMPSXY\,MPS\,X.
Under this framework we have the following properties for risk and variability measures:

(B3) SSD-monotonicity: if XSSDYX\prec_{SSD}Y, then ρ(X)ρ(Y)\rho(X)\leq\rho(Y).

(C3) CX-monotonicity: if XCXYX\prec_{CX}Y, then ν(X)ν(Y)\nu(X)\leq\nu(Y).

Remark 2.11.

Let q[1,]q\in[1,\infty], on LqL^{q} all real-valued law-invariant coherent risk measures are SSD-monotone. We refer the reader to Dana [5], Grechuk et al. [14], and Föllmer and Schied [9] for proofs of the above assertions, and to Mao and Wang [15] for a characterization of SSD-monotone risk measures.

3 Extended Gini Shortfall

In this section, we expose our main contribution, which is based on the Gini coefficient, a free-center measure of variability that was introduced by Corrado Gini as an alternative to the variance measure (e.g., Giorgi ([12], [13]) and Ceriani and Verme [3]). The Gini coefficient has been remarkably influential in numerous research areas (e.g, Yitzhaki and Schechtman [30] and the references therein). Yitzhaki [28] lists more than a dozen alternative presentations of the Gini coefficient. We now present a formal definition.

Definition 3.1.

The Gini coefficient is a functional Gini:L1[0,]Gini:L^{1}\rightarrow[0,\infty] defined conform:

Gini(X)=𝔼[XX],Gini(X)=\mathbb{E}[\mid X^{*}-X^{**}\mid], (3.1)

where XX^{*} and XX^{**} are two independent copies of XX.

Remark 3.2.

The Gini coefficient can be written in terms of a signed Choquet integral:

Gini(X)=201FX1(u)(2u1)𝑑u.Gini(X)=2\int_{0}^{1}F^{-1}_{X}(u)(2u-1)du. (3.2)

From Theorem 2.9, it follows immediately that the Gini coefficient is a coherent measure of variability and it is CX-monotone. Moreover, equation (3.2) can be written in terms of covariance (which is the most common formula of the Gini coefficient):

Gini(X)=4Cov[FX1(U),U]=4Cov[X,UX].Gini(X)=4Cov[F^{-1}_{X}(U),U]=4Cov[X,U_{X}]. (3.3)

We recall that UU can be any uniformly on [0,1][0,1] distributed rv, and UXU_{X} is a uniform [0,1][0,1] rv such that the equation FX1(UX)=XF^{-1}_{X}(U_{X})=X holds.

The Gini functional supposes that all individuals have the same attitude towards risk. Nonetheless, the concept can be extended into a family of measures of variability differing from each other in the decision-maker’s degree of risk aversion, which is reflected in this paper by the parameter rr. The basic definition of the Extended Gini coefficient is based on the covariance term. We refer to Yitzhaki [27], Shalit and Yitzhaki [24], Yitzhaki and Schechtman [29], and Yitzhaki and Schechtman [30] for an overview of the Extended Gini properties. We now expose it in a formal sense.

Definition 3.3.

The Extended Gini coefficient is a functional EGinir:L1[0,]EGini_{r}:L^{1}\rightarrow[0,\infty] defined conform:

EGinir(X)=2rCov[X,(1FX(X))r1],r>1.EGini_{r}(X)=-2rCov[X,(1-F_{X}(X))^{r-1}],\>r>1. (3.4)
Remark 3.4.

There are special cases of interest for the Extended Gini:

\centerdot

For r=2r=2: the Extended Gini coefficient becomes the simple Gini.

\centerdot

For rr\rightarrow\infty: the Extended Gini reflects the attitude of a max-min investor who expresses risk only in terms of the worst outcome.

\centerdot

For r1r\rightarrow 1: the Extended Gini tends to zero and represents the attitude of a risk-neutral individual who does not care about variability.

We now explore the characterization of the Extended Gini coefficient as a signed Choquet integral. In this sense, note that from equation (2.5) in Theorem 2.9, if one sets gr(u)=r(1u)r1g_{r}(u)=-r(1-u)^{{}^{r-1}} for r>1r>1 and u[0,1]u\in[0,1], we run into ν(X)=rCov[X,(1FX(X))r1]\nu(X)=-r\,Cov[X,(1-F_{X}(X))^{r-1}]. With that in mind, we now state and prove the formal result.

Proposition 3.5.

The Extended Gini functional is a CX-monotone coherent measure of variability, represented by the signed Choquet integral

EGinir(X)=201FX1(u)(1+gr(u))𝑑u.EGini_{r}(X)=2\int_{0}^{1}F^{-1}_{X}(u)(1+g_{r}(u))du. (3.5)
Proof.

We recall that UU can be any uniformly distributed rv on [0,1][0,1] such that the equation FX1(U)=XF^{-1}_{X}(U)=X holds, where 𝔼[X]=m\mathbb{E}[X]=m. In order to obtain the proof, we need the claim that 𝔼[(1FX,p(X))r1|X>xp]=(1p)r1/r\mathbb{E}[(1-F_{X,p}(X))^{r-1}|X>x_{p}]=(1-p)^{r-1}/r for XL1X\in L^{1}, r(1,)r\in(1,\infty) and p[0,1)p\in[0,1). In order to prove it, we get that

𝔼[(1FX,p(X))r1|X>xp]\displaystyle\mathbb{E}[(1-F_{X,p}(X))^{r-1}|X>x_{p}] =𝔼[(1U)r1|U>p]\displaystyle=\mathbb{E}[(1-U)^{r-1}|U>p]
=11p(1u)r1𝟙[p,1](u)𝑑u\displaystyle=\frac{1}{1-p}\int_{\mathbb{R}}(1-u)^{r-1}\mathds{1}_{[p,1]}(u)du
=11pp1(1u)r1𝑑u\displaystyle=\frac{1}{1-p}\int_{p}^{1}(1-u)^{r-1}du
=11p[(1u)r/r]p1\displaystyle=\frac{1}{1-p}[-(1-u)^{r}/r]_{p}^{1}
=(1p)r1/r.\displaystyle=(1-p)^{r-1}/r.

Under this perspective, we can easily verify that: 𝔼[(1FX(X))r1]=1/r\mathbb{E}[(1-F_{X}(X))^{r-1}]=1/r. Thus, we get the following:

EGinir(X)\displaystyle EGini_{r}(X) =2rCov[X,(1FX(X))r1]\displaystyle=-2r\,Cov[X,(1-F_{X}(X))^{r-1}]
=2r𝔼[(X𝔼(X))((1FX(X))r1𝔼[(1FX(X))r1])]\displaystyle=-2r\,\mathbb{E}[(X-\mathbb{E}(X))((1-F_{X}(X))^{r-1}-\mathbb{E}[(1-F_{X}(X))^{r-1}])]
=2r𝔼[(FX1(U)m)((1U)r11r)]\displaystyle=-2r\,\mathbb{E}[(F^{-1}_{X}(U)-m)((1-U)^{r-1}-\frac{1}{r})]
=2r01(FX1(u)m)((1u)r11r)𝑑u\displaystyle=-2r\,\int_{0}^{1}(F^{-1}_{X}(u)-m)((1-u)^{r-1}-\frac{1}{r})du
=2r01FX1(u)((1u)r11r)𝑑u\displaystyle=-2r\,\int_{0}^{1}F^{-1}_{X}(u)((1-u)^{r-1}-\frac{1}{r})du
=201FX1(u)(1+gr(u))𝑑u.\displaystyle=2\,\int_{0}^{1}F^{-1}_{X}(u)(1+g_{r}(u))du.

From this signed Choquet integral representation, the other claims in the proposition follow immediately from Theorem 2.9 and the fact that all coherent measures of variability are CX-monotone. This concludes the proof. ∎

Remark 3.6.

In this case, hrh_{r} in equation (2.4) of Theorem 2.9 is given by:

hr(u)=u+(1u)r1,r>1h_{r}(u)=u+(1-u)^{r}-1,\>r>1 (3.6)

Based on the exposed content, we now turn the focus to an adaptation to the tails of the distribution function. In this sense, we now introduce the Tail Extended Gini (TEG) functional, as well formally prove its properties and Choquet integral representation.

Definition 3.7.

The Tail Extended Gini is a functional TEGinir,p:L1[0,]TEGini_{r,p}:L^{1}\rightarrow[0,\infty] defined conform:

TEGinir,p(X)=2r1pCov[X,(1FX(X))r1|X>xp],r>1, 0<p<1.TEGini_{r,p}(X)=\frac{-2r}{1-p}\,Cov[X,(1-F_{X}(X))^{r-1}|X>x_{p}],\>r>1,\>0<p<1. (3.7)
Proposition 3.8.

The Tail Extended Gini is standardized, location invariant, positively homogeneous and co-monotonic additive. Moreover, it is a signed Choquet integral conform:

TEGinir,p(X)=2(1p)2p1FX1(u)[gr(u)+(1p)r1]𝑑u.TEGini_{r,p}(X)=\frac{2}{(1-p)^{2}}\int_{p}^{1}F_{X}^{-1}(u)[g_{r}(u)+(1-p)^{r-1}]du. (3.8)
Proof.

The fact that the Tail Extended Gini is standardized, location invariant, positively homogeneous is easily realized from its definition. For the Choquet representation, we again recall that UU can be any uniformly distributed rv on [0,1][0,1] such that the equation FX1(U)=XF^{-1}_{X}(U)=X holds, where 𝔼[X]=m\mathbb{E}[X]=m. Thus, we obtain that

TEGinir,p(X)\displaystyle TEGini_{r,p}(X) =2r1pCov[X,(1FX(X))r1|X>xp]\displaystyle=\frac{-2r}{1-p}\,Cov[X,(1-F_{X}(X))^{r-1}|X>x_{p}]
=2r1p𝔼[(Xm)((1FX(X))r1(1p)r1/r)|X>xp]\displaystyle=\frac{-2r}{1-p}\,\mathbb{E}[(X-m)((1-F_{X}(X))^{r-1}-(1-p)^{r-1}/r)|X>x_{p}]
=2r1p𝔼[(Xm)((1FX(X))r1(1p)r1/r)|X>xp]\displaystyle=\frac{-2r}{1-p}\,\mathbb{E}[(X-m)((1-F_{X}(X))^{r-1}-(1-p)^{r-1}/r)|X>x_{p}]
=2r1p𝔼[(FX1(U)m)((1U)r1(1p)r1/r)|U>p]\displaystyle=\frac{-2r}{1-p}\,\mathbb{E}[(F_{X}^{-1}(U)-m)((1-U)^{r-1}-(1-p)^{r-1}/r)|U>p]
=2r(1p)2p1(FX1(u)m)((1u)r1(1p)r1/r)𝑑u\displaystyle=\frac{-2r}{(1-p)^{2}}\,\int_{p}^{1}(F^{-1}_{X}(u)-m)((1-u)^{r-1}-(1-p)^{r-1}/r)du
=2(1p)2[p1FX1(u)(r(1u)r1+(1p)r1)du\displaystyle=\frac{2}{(1-p)^{2}}[\int_{p}^{1}F^{-1}_{X}(u)(-r(1-u)^{r-1}+(1-p)^{r-1})du
+mrp1((1u)r1(1p)r1/r)du]\displaystyle\qquad\qquad+mr\,\int_{p}^{1}((1-u)^{r-1}-(1-p)^{r-1}/r)du]
=2(1p)2p1FX1(u)[gr(u)+(1p)r1]𝑑u.\displaystyle=\frac{2}{(1-p)^{2}}\,\int_{p}^{1}F^{-1}_{X}(u)[g_{r}(u)+(1-p)^{r-1}]du.

In fact, from the proof of the previous proposition we have that p1((1u)r1(1p)r1/r)𝑑u=0\int_{p}^{1}((1-u)^{r-1}-(1-p)^{r-1}/r)du=0. Finally, the Choquet representation implies co-monotonic additivity. This completes the proof. ∎

Remark 3.9.

However, as shown in a counter example (for r=2r=2) by Furman et al. [11], the Tail Extended Gini is not sub-additive. Therefore, unlike the Extended Gini functional, the tail counterpart is not a coherent measure of variability.

Despite the fact that TEG is not a coherent measure of variability, we will show now that a linear combination of the Expected Shortfall with the Tail Extended Gini gives rise to a coherent risk measure, the Extended Gini Shortfall, that quantifies both the magnitude and the variability of tail risks. We now define such a combination.

Definition 3.10.

Let p(0,1)p\in(0,1). We have that the Value at Risk and the Expected Shortfall are functionals VaRp:L0(,]VaR_{p}:L^{0}\rightarrow(-\infty,\infty] and ESp:L1(,]ES_{p}:L^{1}\rightarrow(-\infty,\infty] defined conform:

VaRp(X)=inf{x:FX(x)p},VaR_{p}(X)=\inf\{x\in\mathbb{R}:F_{X}(x)\geq p\}, (3.9)
ESp(X)=11pp1VaRq(X)𝑑q.ES_{p}(X)=\frac{1}{1-p}\int_{p}^{1}VaR_{q}(X)dq. (3.10)
Remark 3.11.

It is a well-established fact that ES is a SSD-monotone comonotonic additive coherent risk measure. When the cdf FXF_{X} is continuous, ES coincides with tail conditional expectation 𝔼[X|Xxp]\mathbb{E}[X|X\geq x_{p}].

Definition 3.12.

The Extended Gini Shortfall is a functional EGSr,pλ:L1(,]EGS_{r,p}^{\lambda}:L^{1}\rightarrow(-\infty,\infty] defined conform:

EGSr,pλ(X)=ESp(X)+λTEGinir,p(X),λ0.EGS_{r,p}^{\lambda}(X)=ES_{p}(X)+\lambda\,TEGini_{r,p}(X),\lambda\geq 0. (3.11)

In order to be a reasonable tool for risk management, the properties of coherent risk measures are desired. However, as mentioned in the previous section, TEG is not sub-additive, and as a measure of variability is not monotone. However, when λ\lambda is zero, then EGS obviously inherits all the properties of the ES which is coherent, but when λ\lambda is sufficiently large, then the TEG starts to dominate ES, and thus coherence of EGS cannot be expected. Intuitively, as suggested by Furman et al. [11], there might be a threshold that delineates the value of λ\lambda for which EGS is coherent. We now verify it in a formal way.

Proposition 3.13.

The Extended Gini shortfall is translation invariant, positively homogeneous, and co-monotonically additive, which can be represented as a signed Choquet integral conform:

EGSr,pλ(X)=01FX1(u)ϕr,pλ(u)𝑑uEGS_{r,p}^{\lambda}(X)=\int_{0}^{1}F^{-1}_{X}(u)\,\phi_{r,p}^{\lambda}(u)du (3.12)

where,

ϕr,pλ(u)=1(1p)2[1p+2λ(gr(u)+(1p)r1)] 1[p,1](u),u[0,1].\phi_{r,p}^{\lambda}(u)=\frac{1}{(1-p)^{2}}[1-p+2\lambda(g_{r}(u)+(1-p)^{r-1})]\,\mathds{1}_{[p,1]}(u),\quad u\in[0,1]. (3.13)

Moreover, the Extended Gini Shortfall is a SSD-monotone coherent risk measure for λ[0, 1/(2(r1)(1p)r2)]\lambda\in[0\,,\,1/(2(r-1)(1-p)^{r-2})].

Proof.

The translation invariance, positive homogeneity and co-monotonic additivity are easily verifiable from the properties of ES and TEG. Regarding the Choquet representation, we have that:

EGSr,pλ(X)\displaystyle EGS_{r,p}^{\lambda}(X) =ESp(X)+λTEGinir,p(X)\displaystyle=ES_{p}(X)+\lambda\,TEGini_{r,p}(X)
=11pp1FX1(u)𝑑u+2λ(1p)2p1FX1(u)[gr(u)+(1p)r1]𝑑u\displaystyle=\frac{1}{1-p}\int_{p}^{1}F^{-1}_{X}(u)du+\frac{2\lambda}{(1-p)^{2}}\,\int_{p}^{1}F^{-1}_{X}(u)[g_{r}(u)+(1-p)^{r-1}]du
=p1FX1(u)[11p+2λ(1p)2(gr(u)+(1p)r1)]𝑑u\displaystyle=\int_{p}^{1}F^{-1}_{X}(u)[\frac{1}{1-p}+\frac{2\lambda}{(1-p)^{2}}(g_{r}(u)+(1-p)^{r-1})]du
=01FX1(u)ϕr,pλ(u)𝑑u.\displaystyle=\int_{0}^{1}F^{-1}_{X}(u)\,\phi_{r,p}^{\lambda}(u)du.

For coherence, it remains to prove that EGS is monotone and sub-additive for λ[0, 1/(2(r1)(1p)r2)]\lambda\in[0\,,\,1/(2(r-1)(1-p)^{r-2})]. Note that ϕr,pλ\phi_{r,p}^{\lambda} is an increasing function on [0,1][0,1], therefore ϕr,pλ(u)\phi_{r,p}^{\lambda}(u) is non-negative if and only if ϕr,pλ(p)0\phi_{r,p}^{\lambda}(p)\geq 0. Thus, ϕr,pλ\phi_{r,p}^{\lambda} is non-negative if and only if
λ[0, 1/(2(r1)(1p)r2)]\lambda\in[0\,,\,1/(2(r-1)(1-p)^{r-2})]. This fact implies monotonicity and sub-additivity from the properties discussed on section 2. Thus, EGS is a coherent risk measure for this choice of λ\lambda. Finally, SSD-monotonicity for this choice of λ\lambda, it is directly implied by the fact that EGS is law invariant. This concludes the proof. ∎

Remark 3.14.

From the previous Proposition we can link the EGS with its acceptance set, which is defined as 𝒜EGSr,pλ={XL1:EGSr,pλ(X)0}\mathcal{A}_{EGS_{r,p}^{\lambda}}=\{X\in L^{1}:EGS_{r,p}^{\lambda}(X)\leq 0\}. It is well-known, see Föllmer and Schied [9] for instance, that this set is convex, law invariant, monotone, closed for multiplication with positive scalar and addition between co-monotonic variables. Moreover, we have that 𝒜EGSr,pλ\mathcal{A}_{EGS_{r,p}^{\lambda}} contains L+1L^{1}_{+} and has no intersection with {XL1:XL+1}\{X\in L^{1}:X\notin L^{1}_{+}\}. It is direct to verify, from Translation Invariance of EGS, that EGSr,pλ(X)=inf{m:X+m𝒜EGSr,pλ}EGS_{r,p}^{\lambda}(X)=\inf\{m:X+m\in\mathcal{A}_{EGS_{r,p}^{\lambda}}\}.

Remark 3.15.

We have that EGS can be represented as convex combination of ES at distinct levels of pp, a Kusuoka representation, conform EGSr,pλ(X)=01ESp(X)μ(dp)EGS_{r,p}^{\lambda}(X)=\int_{0}^{1}ES_{p}(X)\mu(dp), where ϕr,pλ(u)=[1u,1)1sμ(ds)\phi_{r,p}^{\lambda}(u)=\int_{[1-u,1)}\frac{1}{s}\mu(ds). Here, μ\mu is a probability measure over (0,1](0,1]. These representations are linked to the well-known dual representation, conform EGSr,pλ(X)=E[XQ]EGS_{r,p}^{\lambda}(X)=E[XQ], where FQ1(u)=ϕr,pλ(u)F^{-1}_{Q}(u)=\phi_{r,p}^{\lambda}(u). We can think about QQ as the relative density (Radon-Nikodym) of an alternative probability measure absolutely continuous in relation to \mathbb{P}.

Remark 3.16.

From the previous Proposition, we get that the Extended Gini Shortfall is part of the spectral risk measures class, introduced in Acerbi [1], characterized by the weighting function ϕr,pλ\phi_{r,p}^{\lambda}, which enables to reflect the individual’s subjective attitude toward risk. In Furman et al. [11] a specific case is introduced (r=2r=2), which assigns the same weighting function to all decision-makers. Thus, there is a connection to the individual’s risk aversion function. As a result, it is more legitimate for ϕr,pλ\phi_{r,p}^{\lambda} to be dependent on the parameters rr, pp and λ\lambda.

We now provide a result and interpretation about how this weighting function ϕ\phi behaves in relation to changes (partial derivatives) of each variable (parameter) among uu, rr, pp and λ\lambda. It is valid to point out that, by the spectral representation, results can be directly understood as the effect each parameter has over values for EGS.

Proposition 3.17.

Consider the weighting function

ϕr,pλ(u)=ϕ(u,r,p,λ)=1p+2λ[(1p)r1r(1u)r1](1p)2 1[p,1](u),\phi_{r,p}^{\lambda}(u)=\phi(u,r,p,\lambda)=\frac{1-p+2\lambda[(1-p)^{r-1}-r(1-u)^{r-1}]}{(1-p)^{2}}\,\mathds{1}_{[p,1]}(u),

where u[0,1]u\in[0,1], p(0,1)p\in(0,1), r>1r>1, and λ[0,1/(2(r1)(1p)r2)]\lambda\in[0,1/(2(r-1)(1-p)^{r-2})]. We have that:

(i) The interval for values of λ\lambda that make EGS subadditive has superior limit non-decreasing in pp, if r2r\geq 2 and non-increasing in pp otherwise. Moreover, the interval is non-decreasing in rr if, and only if, r11ln(1p)r\geq 1-\frac{1}{\ln(1-p)};

(ii) ϕ\phi is non-decreasing in uu;

(iii) ϕ\phi is non-decreasing in λ\lambda;

(iv) ϕ\phi is non-decreasing in pp if, and only if

u1(1p2λ(r3)(1p)r14λr)(1r1);u\geq 1-\left(\frac{1-p-2\lambda(r-3)(1-p)^{r-1}}{4\lambda r}\right)^{\left(\frac{1}{r-1}\right)};

(v) ϕ\phi is non-decreasing in rr if, and only if

(1p)(1p)(exp{(1u)r1[rln(1u)+1]})(1r1).(1-p)^{(1-p)}\geq\left(\exp\left\{(1-u)^{r-1}[r\ln(1-u)+1]\right\}\right)^{\left(\frac{1}{r-1}\right)}.
Proof.

For (i), we must remember that EGS is subadditive when λ[0,1/(2(r1)(1p)r2)]\lambda\in[0,1/(2(r-1)(1-p)^{r-2})]. Consider the functional

B(r,p)=12(r1)(1p)r2B(r,p)=\frac{1}{2(r-1)(1-p)^{r-2}}

that represents the threshold that λ\lambda shall not exceed. Regarding pp, we get

Bp(r,p)=Bp(r,p)=(r2)(1p)1r2(r1).\frac{\partial B}{\partial p}(r,p)=B^{\prime}_{p}(r,p)=\frac{(r-2)(1-p)^{1-r}}{2(r-1)}.

It is direct that the sign of BpB^{\prime}_{p} depends on the sign of r2r-2. Thus, B(r,p)B(r,p) is non-decreasing in pp for r2r\geq 2 and non-increasing otherwise. Regarding rr, we thus have that

Br(r,p)=Br(r,p)=12(1p)r2[1+(r1)ln(1p)][(r1)(1p)r2]2.\frac{\partial B}{\partial r}(r,p)=B^{\prime}_{r}(r,p)=\frac{-1}{2}\frac{(1-p)^{r-2}[1+(r-1)\ln(1-p)]}{[(r-1)(1-p)^{r-2}]^{2}}.

The sign of BrB^{\prime}_{r} depends on the sign of Cp(r)=1+(r1)ln(1p)C_{p}(r)=1+(r-1)\ln(1-p). The decreasing function CpC_{p} maps (1,)(1,\infty) to (,1)(-\infty,1), then there exists a unique critical value r0=11ln(1p)r_{0}=1-\frac{1}{\ln(1-p)} such that BpB_{p} non-increases over (1,r0](1,r_{0}] and non-decreases on (r0,)(r_{0},\infty).

For (ii), the claim follows from the fact that EGS is a spectral risk measure. More specifically, we have that

ϕu(u,r,p,λ)=2λr(r1)(1u)r2(1p)2 1[p,1](u),\frac{\partial\phi}{\partial u}(u,r,p,\lambda)=\frac{2\lambda\,r(r-1)(1-u)^{r-2}}{(1-p)^{2}}\,\mathds{1}_{[p,1]}(u),

which is non-negative any case.

Regarding (iii), the claim is direct by definition of EGS. More specifically on the weighting function we obtain

ϕλ(u,r,p,λ)=2[(1p)r1r(1u)r1](1p)2 1[p,1](u).\frac{\partial\phi}{\partial\lambda}(u,r,p,\lambda)=\frac{2[(1-p)^{r-1}-r(1-u)^{r-1}]}{(1-p)^{2}}\,\mathds{1}_{[p,1]}(u).

We can note that this expression is non-decreasing in uu with critical value when u=1(1p)r1r1u=1-(1-p)r^{-\frac{1}{r-1}}. Thus, ϕλ0\frac{\partial\phi}{\partial\lambda}\geq 0 when u1(1p)r1r1pu\geq 1-(1-p)r^{-\frac{1}{r-1}}\geq p, since r1r11r^{\frac{1}{r-1}}\geq 1. But upu\geq p is the only case that matters because of the indicator function in ϕλ\frac{\partial\phi}{\partial\lambda}.

For item (iv) we begin by noticing that

ϕp(u,r,p,λ)=12λ(r3)(1p)r24λr(1u)r1(1p)1(1p)2 1[p,1](u),\frac{\partial\phi}{\partial p}(u,r,p,\lambda)=\frac{1-2\lambda(r-3)(1-p)^{r-2}-4\lambda r(1-u)^{r-1}(1-p)^{-1}}{(1-p)^{2}}\,\mathds{1}_{[p,1]}(u),

which is a non-decreasing expression in uu. Thus, we have that ϕp0\frac{\partial\phi}{\partial p}\geq 0 if, and only if,

u1(1p2λ(r3)(1p)r14λr)1r1.u\geq 1-\left(\frac{1-p-2\lambda(r-3)(1-p)^{r-1}}{4\lambda r}\right)^{\frac{1}{r-1}}.

Finally, regarding item (v), we follow the same reasoning to get

ϕr(u,r,p,λ)=2λ[ln(1p)(1p)r1(1u)r1rln(1u)(1u)r1](1p)2.\frac{\partial\phi}{\partial r}(u,r,p,\lambda)=\frac{2\lambda[\ln(1-p)(1-p)^{r-1}-(1-u)^{r-1}-r\ln(1-u)(1-u)^{r-1}]}{(1-p)^{2}}.

This is a complex expression which can assume both positive and negative values since some terms in numerator posses distinct signs without a domination. Moreover, isolating uu or rr is not trivial. Nonetheless, after some manipulation we obtain that ϕr0\frac{\partial\phi}{\partial r}\geq 0 if, and only if, (1p)(1p)(exp{(1u)r1[rln(1u)+1]})(1r1)(1-p)^{(1-p)}\geq\left(\exp\left\{(1-u)^{r-1}[r\ln(1-u)+1]\right\}\right)^{\left(\frac{1}{r-1}\right)}. This concludes the proof. ∎

Remark 3.18.

The non-linear behavior of the partial derivatives found by the previous proposition is related to the function gr(u)=r(1u)r1g_{r}(u)=-r(1-u)^{r-1}, which is central to the theory we develop, which is not trivial for r2r\neq 2, the case considered in Furman et al. [11]. As one can easily note, it recurrently appears in the expression for partial derivatives. Hence, our contribution is by extending the canonical case to a situation where more complexity regarding rr can be addressed.

Since pp reflects the prudence level, which is usually close to 11 in practice, then r0r_{0} in the proof of item (i) is in general close to 1. Thus, for practical matters, the superior limit for λ\lambda is non-decreasing in most relevant values of rr. The most complex parameter sensitivities are regarding the prudence level pp and the generalization term rr, respectively in items (iv) and (v), because both ϕp\frac{\partial\phi}{\partial p} and ϕr\frac{\partial\phi}{\partial r} are expressions that can assume both positive and negative values since some terms in numerator posses distinct signs without a domination. This can be explained due to the fact that ϕ\phi is a weighting function and changes in pp and rr alter how much ‘mass’ is put to any probability level uu. Because 01ϕr,pλ(u)𝑑u=1\int_{0}^{1}\phi_{r,p}^{\lambda}(u)du=1, it is necessary that the increase on ϕ(u)\phi(u) for some values of uu be compensated by some decrease in others.

Regarding prudence level, item (iv) in the previous proposition emphasizes that ϕ\phi non-decreases in pp for larger values of uu. This is in consonance with practical intuition because more weight is put to extreme probabilities. Moreover, this is corroborated when we verify the variations of ϕ\phi regarding uu and pp, where we obtain the non-negative expression

2ϕup=4λr(r1)(1u)r2(1p)3 1[p,1](u).\frac{\partial^{2}\phi}{\partial u\partial p}=\frac{4\lambda r(r-1)(1-u)^{r-2}}{(1-p)^{3}}\,\mathds{1}_{[p,1]}(u).

When we consider the sensibility of variations regarding uu and rr, in a case we get

2ϕur=2λ(1u)r2[(2r1)+(r2r)ln(1u)](1p)2 1[p,1](u),\frac{\partial^{2}\phi}{\partial u\partial r}=\frac{2\lambda(1-u)^{r-2}[(2r-1)+(r^{2}-r)\ln(1-u)]}{(1-p)^{2}}\,\mathds{1}_{[p,1]}(u),

there is divergence about sign – first term inside brackets is non-negative while the second is non-positive, corroborating with the previous argument. Nonetheless, ϕ\phi is non-decreasing in rr, a situation where rr behaves more like a risk aversion coefficient, when a non-decreasing function of pp, (1p)(1p)(1-p)^{(1-p)} is greater than a threshold that is an expression depending on both uu and rr. Repeating the argument of practical choices for the prudence level, pp will be close to 11 and ϕ\phi will be increasing in rr for more values of uu. In this sense, rr must be understood as a generalization parameter for a family of EGS risk measures rather than a single linear risk aversion coefficient.

4 Extended Gini Shortfall for usual distributions

In this section we provide analytical formulations to compute the proposed Extended Gini Shortfall to known and very used distribution functions. In that sense, a location-scale family is a family of probability distributions parameterized by a location parameter and a non-negative scale parameter. Suppose that ZZ is a fixed rv taking values in \mathbb{R}. For α\alpha\in\mathbb{R} and β(0,)\beta\in(0,\infty), let X=α+βZX=\alpha+\beta Z. The two-parameter family of distributions associated with XX is called the location-scale family associated with the given distribution of ZZ; α\alpha is called the location parameter and β\beta the scale parameter. The standard form of any distribution is the form whose location and scale parameters are 0 and 11, respectively. In this section, we restrain our attention into standardized rv’s. In general, when X=α+βZX=\alpha+\beta Z for α\alpha\in\mathbb{R} and β(0,)\beta\in(0,\infty), we have, directly from their properties, both ESp(X)=α+βESp(Z)ES_{p}(X)=\alpha+\beta ES_{p}(Z) and TEGinir,p(X)=βTEGinir,p(Z)TEGini_{r,p}(X)=\beta TEGini_{r,p}(Z).

In what follows, we start with the general elliptical family and then specialize the obtained result to uniform, normal and Student-t distribution cases. We recall that XX is an elliptical distribution if X=dα+βZX\stackrel{{\scriptstyle d}}{{=}}\alpha+\beta Z, where ZZ is a spherical distribution. Let ZZ be a spherical rv with characteristic generator ψ:[0,)\psi:[0,\infty)\rightarrow\mathbb{R}; succinctly ZS(ψ)Z\sim S(\psi). When ZZ has a probability density function (pdf), then there is a density generator g:[0,)[0,)g:[0,\infty)\rightarrow[0,\infty) such that 0z1/2g(z)𝑑z<\int_{0}^{\infty}z^{-1/2}g(z)dz<\infty, we succinctly write ZS(g)Z\sim S(g). We can express the pdf f:[0,]f:\mathbb{R}\rightarrow[0,\infty] of ZZ by f(z)=cg(z2/2)f(z)=c\,g(z^{2}/2), where c>0c>0 is the normalizing constant. The mean 𝔼[Z]\mathbb{E}[Z] is finite when 0g(z)𝑑z<\int_{0}^{\infty}g(z)dz<\infty, in which case we have 𝔼[Z]=0\mathbb{E}[Z]=0 because the pdf ff is symmetric around 0. Under this condition, we define the function G¯:[0,)[0,)\overline{G}:[0,\infty)\rightarrow[0,\infty) by G¯(y)=cyg(x)𝑑x\overline{G}(y)=c\int_{y}^{\infty}g(x)dx, which is called the tail generator of ZZ. We now state and prove the main result in this section.

Proposition 4.1.

Let ZS(g)Z\sim S(g), 𝔼[Z]\mathbb{E}[Z] finite, and p(0,1)p\in(0,1). Then, we have:

ESp(Z)=G¯(zp2/2)1p,ES_{p}(Z)=\frac{\overline{G}(z_{p}^{2}/2)}{1-p}, (4.1)
TEGinir,p(Z)=2r(r1)1p𝔼[(1FZ(Z))r2G¯(Z2/2)|Z>zp]+2[1r(1p)r2]ESp(Z).TEGini_{r,p}(Z)=\frac{2r(r-1)}{1-p}\mathbb{E}\left[(1-F_{Z}(Z))^{r-2}\overline{G}(Z^{2}/2)|Z>z_{p}\right]+2[1-r(1-p)^{r-2}]ES_{p}(Z). (4.2)
Proof.

For ES. we have that:

ESp(Z)\displaystyle ES_{p}(Z) =𝔼[Z|Z>zp]\displaystyle=\mathbb{E}[Z|Z>z_{p}]
=11pzpzf(z)𝑑z\displaystyle=\frac{1}{1-p}\int_{z_{p}}^{\infty}zf(z)dz
=c1pzpzg(z2/2)𝑑z\displaystyle=\frac{c}{1-p}\int_{z_{p}}^{\infty}zg(z^{2}/2)dz
=c1pzp2/2g(x)𝑑x\displaystyle=\frac{c}{1-p}\int_{z_{p}^{2}/2}^{\infty}g(x)dx
=G¯(zp2/2)1p.\displaystyle=\frac{\overline{G}(z_{p}^{2}/2)}{1-p}.

Concerning to TEG, we get:

TEGinir,p(Z)=2r1p[𝔼[Z(1FZ(Z))r1|Z>zp]𝔼[Z|Z>zp]𝔼[(1FZ(Z))r1|Z>zp]]TEGini_{r,p}(Z)=\frac{-2r}{1-p}\left[\mathbb{E}[Z(1-F_{Z}(Z))^{r-1}|Z>z_{p}]-\mathbb{E}[Z|Z>z_{p}]\mathbb{E}[(1-F_{Z}(Z))^{r-1}|Z>z_{p}]\right]
𝔼[Z(1FZ(Z))r1|Z>zp]\displaystyle\mathbb{E}[Z(1-F_{Z}(Z))^{r-1}|Z>z_{p}] =11p𝔼[Z(1FZ(Z))r1 1{Z>zp}]\displaystyle=\frac{1}{1-p}\mathbb{E}[Z(1-F_{Z}(Z))^{r-1}\,\mathds{1}_{\{Z>z_{p}\}}]
=11pzpz(1FZ(z))r1f(z)𝑑z\displaystyle=\frac{1}{1-p}\int_{z_{p}}^{\infty}z(1-F_{Z}(z))^{r-1}f(z)dz

Note that zf(z)dz=dG¯(z2/2)zf(z)dz=-d\overline{G}(z^{2}/2) and G¯(z2/2)=(1p)ESp(Z)\overline{G}(z^{2}/2)=(1-p)ES_{p}(Z). Integration by parts leads to:

𝔼[Z(1FZ(Z))r1|Z>zp]\displaystyle\mathbb{E}[Z(1-F_{Z}(Z))^{r-1}|Z>z_{p}] =11p(1p)r1G¯(z2/2)(r1)𝔼[(1FZ(Z))r2G¯(Z2/2)|Z>zp]\displaystyle=\frac{1}{1-p}(1-p)^{r-1}\overline{G}(z^{2}/2)-(r-1)\mathbb{E}[(1-F_{Z}(Z))^{r-2}\overline{G}(Z^{2}/2)|Z>z_{p}]
=(1p)r1ESp(Z)(r1)𝔼[(1FZ(Z))r2G¯(Z2/2)|Z>zp].\displaystyle=(1-p)^{r-1}ES_{p}(Z)-(r-1)\mathbb{E}[(1-F_{Z}(Z))^{r-2}\overline{G}(Z^{2}/2)|Z>z_{p}].
𝔼[Z|Z>zp]=ESp(Z).\mathbb{E}[Z|Z>z_{p}]=ES_{p}(Z).

Finally, similarly to the proof of Proposition 3.5, we have:

𝔼[(1FZ(Z))r1|Z>zp]=(1p)r1r.\mathbb{E}[(1-F_{Z}(Z))^{r-1}|Z>z_{p}]=\frac{(1-p)^{r-1}}{r}.

This completes the proof. ∎

Remark 4.2.

Note that Var(Z)Var(Z) is finite whenever 0z1/2g(z)𝑑z<\int_{0}^{\infty}z^{1/2}g(z)dz<\infty, in which case Var(Z)Var(Z) is equal to G¯(z2/2)𝑑z\int_{-\infty}^{\infty}\overline{G}(z^{2}/2)dz (by an integration by parts). Hence, we have that f(z)=G¯(zp2/2)/Var(Z)f^{*}(z)=\overline{G}(z_{p}^{2}/2)/Var(Z) is a pdf. In this case we get the expression TEGinir,p(Z)=2r(r1)1pVar(Z)𝔼[(1FZ(Z))r2f(Z)|Z>zp]+2[1r(1p)r2]ESp(Z)TEGini_{r,p}(Z)=\frac{2r(r-1)}{1-p}Var(Z)\mathbb{E}\left[(1-F_{Z}(Z))^{r-2}f^{*}(Z)|Z>z_{p}\right]+2[1-r(1-p)^{r-2}]ES_{p}(Z).

We now focus on the case of a standard uniform rv ZU[1,1]Z\sim U[-1,1]. In this sense, we have the following Corollary.

Corollary 4.3.

Let ZU[1,1]Z\sim U[-1,1] and p(0,1)p\in(0,1). Then, we have:

ESp(Z)=1zp24(1p),ES_{p}(Z)=\frac{1-z_{p}^{2}}{4(1-p)}, (4.3)
TEGinir,p(Z)=r3(1p)2(1zp2)r1+2[1r(1p)r2]ESp(Z).TEGini_{r,p}(Z)=\frac{r}{3(1-p)^{2}}\left(\frac{1-z_{p}}{2}\right)^{r-1}+2[1-r(1-p)^{r-2}]ES_{p}(Z). (4.4)
Proof.

The standard uniform is a spherical distribution with density fZ(z)=12𝟙[1,1](z)=12𝟙[0,1/2](z2/2),z[1,1]f_{Z}(z)=\frac{1}{2}\mathds{1}_{[-1,1]}(z)=\frac{1}{2}\mathds{1}_{[0,1/2]}(z^{2}/2),\quad z\in[-1,1]. Thus, we obtain g(z)=𝟙[0,1/2](z) and c=12g(z)=\mathds{1}_{[0,1/2]}(z)\quad\text{ and }\quad c=\frac{1}{2}. In this case, the tail generator is given by:

G¯(y)=y12g(t)𝑑t=14y2.\overline{G}(y)=\int_{y}^{\infty}\frac{1}{2}g(t)dt=\frac{1}{4}-\frac{y}{2}.

Hence, we obtain

G¯(zp2/2)=1zp24,ESp(Z)=1zp24(1p)\overline{G}(z_{p}^{2}/2)=\frac{1-z_{p}^{2}}{4},\>ES_{p}(Z)=\frac{1-z_{p}^{2}}{4(1-p)}

In order to prove the result for TEGini, we use Remark 4.2. Thus, we obtain:

TEGinir,p(Z)=2r(r1)3(1p)𝔼[(1Z2)r2f(Z)|Z>zp]+2[1r(1p)r2]ESp(Z),TEGini_{r,p}(Z)=\frac{2r(r-1)}{3(1-p)}\mathbb{E}\left[\left(\frac{1-Z}{2}\right)^{r-2}f^{*}(Z)|Z>z_{p}\right]+2[1-r(1-p)^{r-2}]ES_{p}(Z),

with

𝔼[(1Z2)r2f(Z)|Z>zp]\displaystyle\mathbb{E}\left[\left(\frac{1-Z}{2}\right)^{r-2}f^{*}(Z)|Z>z_{p}\right] =14(1p)zp1(1z2)r2𝑑z\displaystyle=\frac{1}{4(1-p)}\int_{z_{p}}^{1}\left(\frac{1-z}{2}\right)^{r-2}dz
=12(1p)(r1)(1zp2)r1.\displaystyle=\frac{1}{2(1-p)(r-1)}\left(\frac{1-z_{p}}{2}\right)^{r-1}.

This completes the proof. ∎

In this next corollary, we deal with the standard Normal rv Z𝒩(0,1)Z\sim\mathcal{N}(0,1) whose cdf we denote by Φ\Phi.

Corollary 4.4.

Let Z𝒩(0,1)Z\sim\mathcal{N}(0,1) and p(0,1)p\in(0,1). Then, we have:

ESp(Z)=Φ(zp)1p,ES_{p}(Z)=\frac{\Phi^{\prime}(z_{p})}{1-p}, (4.5)
TEGinir,p(Z)=2r(r1)1p𝔼[(1Φ(Z))r2Φ(Z)|Z>zp]+2[1r(1p)r2]ESp(Z)TEGini_{r,p}(Z)=\frac{2r(r-1)}{1-p}\mathbb{E}\left[(1-\Phi(Z))^{r-2}\Phi^{\prime}(Z)|Z>z_{p}\right]+2[1-r(1-p)^{r-2}]ES_{p}(Z) (4.6)
Proof.

The standard normal is a spherical distribution with density generator g(z)=exp(z),c=1/2π and G¯(z2/2)=Φ(z)g(z)=exp(-z),\,c=1/\sqrt{2\pi}\,\text{ and }\,\overline{G}(z^{2}/2)=\Phi^{\prime}(z). Thus, equation (4.5) follows immediately from equation (4.1). To establish the second part, we use Remark 4.2 with Var(Z)=1Var(Z)=1. ∎

Finally, we expose a corollary concerning to the Student-t distribution. Let n+n\in\mathbb{R}_{+}^{*}, we say that ZZ has a standard Student-t distribution with nn degree of freedom if its pdf is:

fn(z)=1nBeta(n/2,1/2)(1+z2n)(n+12),zf_{n}(z)=\frac{1}{\sqrt{n}\,Beta(n/2,1/2)}\left(1+\frac{z^{2}}{n}\right)^{-(\frac{n+1}{2})},\>z\in\mathbb{R}

Where Beta(a,b)Beta(a,b) stands for the Beta distribution. We set θ=n+12\theta=\frac{n+1}{2} and kθ=n2k_{\theta}=\frac{n}{2}. Then, we denote Zt(θ)Z\sim t(\theta) with parameter θ>12\theta>\frac{1}{2} and the pdf can be rewritten as:

fθ(z)=cθ(1+z22kθ)θ,zf_{\theta}(z)=c_{\theta}\left(1+\frac{z^{2}}{2k_{\theta}}\right)^{-\theta},\>z\in\mathbb{R}

where cθ=(2kθBeta(θ1/2,1/2))1c_{\theta}=(\sqrt{2k_{\theta}}Beta(\theta-1/2,1/2))^{-1}. The expected value of ZZ is well-defined only for θ>1\theta>1 and the variance of ZZ is finite only if θ>32\theta>\frac{3}{2}. Let FθF_{\theta} be the cdf of ZZ.

Corollary 4.5.

Let Zt(θ)Z\sim t(\theta), θ>1\theta>1, and p(0,1)p\in(0,1). Then, we have:

ESp(Z)=cθkθ(1p)(θ1)(1+zp22kθ)(θ1),ES_{p}(Z)=\frac{c_{\theta}k_{\theta}}{(1-p)(\theta-1)}\left(1+\frac{z^{2}_{p}}{2k_{\theta}}\right)^{-(\theta-1)}, (4.7)
TEGinir,p(Z)\displaystyle TEGini_{r,p}(Z) =2r(r1)(1p)(θ1)cθkθcθ1𝔼[(1Fθ(Z))r2fθ1(kθ1kθZ)|Z>zp]\displaystyle=\frac{2r(r-1)}{(1-p)(\theta-1)}\frac{c_{\theta}k_{\theta}}{c_{\theta-1}}\mathbb{E}\left[(1-F_{\theta}(Z))^{r-2}f_{\theta-1}\left(\sqrt{\frac{k_{\theta-1}}{k_{\theta}}}Z\right)|Z>z_{p}\right]\newline
+2[1r(1p)r2]ESp(Z).\displaystyle+2[1-r(1-p)^{r-2}]ES_{p}(Z). (4.8)
Proof.

The tail generator of the standard Student-t is given by:

g(z)=(1+zkθ)θ,zg(z)=\left(1+\frac{z}{k_{\theta}}\right)^{-\theta},\qquad z\in\mathbb{R}

Hence, we get that:

G¯(z)\displaystyle\overline{G}(z) =cθz(1+tkθ)θ𝑑t\displaystyle=c_{\theta}\int_{z}^{\infty}\left(1+\frac{t}{k_{\theta}}\right)^{-\theta}dt
=cθkθθ1(1+zkθ)(θ1),\displaystyle=\frac{c_{\theta}k_{\theta}}{\theta-1}\left(1+\frac{z}{k_{\theta}}\right)^{-(\theta-1)},

which leads to:

G¯(Z2/2)=cθkθcθ1(θ1)fθ1(kθ1kθZ)\overline{G}(Z^{2}/2)=\frac{c_{\theta}k_{\theta}}{c_{\theta-1}(\theta-1)}f_{\theta-1}\left(\sqrt{\frac{k_{\theta-1}}{k_{\theta}}}Z\right)

Then, equation (4.6) follows immediately from equation (4.2). ∎

5 From theory to practice

In this section we are going to provide an illustration for the practical usage of EGSEGS. In this sense, consider the EGSr,pλEGS_{r,p}^{\lambda} definition:

EGSr,pλ(X)=01FX1(u)ϕr,pλ(u)𝑑u,EGS_{r,p}^{\lambda}(X)=\int_{0}^{1}F^{-1}_{X}(u)\,\phi_{r,p}^{\lambda}(u)du,

where,

ϕr,pλ(u)=1p+2λ[(1p)r1r(1u)r1](1p)2 1[p,1](u),\phi_{r,p}^{\lambda}(u)=\frac{1-p+2\lambda[(1-p)^{r-1}-r(1-u)^{r-1}]}{(1-p)^{2}}\,\mathds{1}_{[p,1]}(u),

with u[0,1]u\in[0,1], p(0,1)p\in(0,1), r>1r>1, and λ[0,1/(2(r1)(1p)r2)].\lambda\in[0,1/(2(r-1)(1-p)^{r-2})].

In practice, the assessment of EGSr,pλEGS_{r,p}^{\lambda} can be reduced to evaluating its discrete version proposed by Acerbi[1] as a consistent estimator for spectral risk measures:

EGS^r,pλ(X)=i=1NX(i)ϕi,\widehat{EGS}_{r,p}^{\lambda}(X)=\sum_{i=1}^{N}X_{(i)}\,\phi_{i}, (5.1)

where, {X(i);i=1,,N}\{X_{(i)};i=1,...,N\} are the ordered statistics given by the N-tuple {X1,,XN}\{X_{1},...,X_{N}\} of observations, and ϕi\phi_{i} is the natural choice for a suitable weighting function given by:

ϕi=ϕr,pλ(i/N)k=1Nϕr,pλ(k/N)i=1,,N\phi_{i}=\frac{\phi_{r,p}^{\lambda}(i/N)}{\sum_{k=1}^{N}\phi_{r,p}^{\lambda}(k/N)}\qquad i=1,...,N (5.2)

and satisfying iϕi=1.\sum_{i}\phi_{i}=1.

According to the expression of ϕr,pλ\phi_{r,p}^{\lambda}, the EGSr,pλEGS_{r,p}^{\lambda} is concerned only with losses beyond the VaRpVaR_{p}, thus we are interested in tail risks.

Remark 5.1.

When rr is sufficiently high, (1p)r1r(1u)r10(1-p)^{r-1}-r(1-u)^{r-1}\rightarrow 0 thus ϕr,pλ11p\phi_{r,p}^{\lambda}\rightarrow\displaystyle\frac{1}{1-p}.
That is to say, under this condition, EGSr,pλEGS_{r,p}^{\lambda} is confounded with ESpES_{p}.

This remark ensures that, for a highly risk averse investor, we can simply use ESpES_{p}. In other words, ESpES_{p} can be considered as a limit of EGSr,pλEGS_{r,p}^{\lambda} for a quite high risk aversion degree. Furthermore, since more risk averse the investor is less risk he takes then we can claim that EGSr,pλEGS_{r,p}^{\lambda} is at least equal to ESpES_{p} i.e., EGSr,pλESpEGS_{r,p}^{\lambda}\geq ES_{p}.

In the following, we illustrate the above approach with the use of a numerical example. Our dataset consists of daily returns from the MASI index covering the period of 15th November 2016 to the 15th November 2017, which includes a total of N=250N=250 observations. The MASI (Moroccan All Shares Index) is a stock index that tracks the performance of all securities listed in the Casablanca Stock Exchange located at Casablanca in Morocco.444http://www.casablanca-bourse.com

The proposed methodology does not make any assumptions about the distribution that describes the data, except that an Augmented Unit Root test is performed to make sure that our data series is stationary. Moreover, as shown in Figure 1, the return graph validates this verification since the series fluctuates around 0 and has no trend.

Refer to caption
Figure 1: Graph of the daily observed MASI return.
Remark 5.2.

When defining the EGSr,pλEGS_{r,p}^{\lambda} we were making the convention that our rv XX represents financial losses (profits) when it has positive (negative) values. However, to be in compliance with the real world data the latter convention is updated.

By using the sorted returns series we calculate the VaRpVaR_{p} value to identify the concerned losses. Then, we affect to each value its own weight in accordance with the parameters p,r,p,r, and λ\lambda as shown below in Table 1.
To fulfill the condition λ[0,1/(2(r1)(1p)r2)]\lambda\in[0,1/(2(r-1)(1-p)^{r-2})] we take λ\lambda arbitrarily as the midpoint of the interval, thus λ=14(r1)(1p)r2\displaystyle\lambda=\frac{1}{4(r-1)(1-p)^{r-2}}.

Sorted Returns(%) Weights
-1.146 0.04
-1.151 0.05
-1.160 0.05
-1.181 0.06
-1.270 0.06
-1.273 0.07
-1.304 0.08
-1.621 0.08
-1.685 0.09
-1.690 0.10
-1.880 0.10
-1.924 0.11
-2.090 0.11
Table 1: Weighted losses beyond VaR for p=95% and r=2.

In Table 2 below, we report the calculation results of EGS^r,pλ\widehat{EGS}_{r,p}^{\lambda} for different values of pp and rr:

EGS^r,pλ\widehat{EGS}_{r,p}^{\lambda} r=2(GSpλ)r=2\,(GS^{\lambda}_{p}) r=3r=3 r=6r=6 r=20r=20 r=30r=30
p=90%p=90\%
VaR=0.73%VaR=0.73\% 1.32% 1.28% 1.25% 1.22% 1.22%
ES=1.21%ES=1.21\%
p=95%p=95\%
VaR=1.11%VaR=1.11\% 1.58% 1.55% 1.52% 1.50% 1.49%
ES=1.49%ES=1.49\%
p=99%p=99\%
VaR=1.79%VaR=1.79\% 1.99% 1.98% 1.97% 1.96% 1.96%
ES=1.96%ES=1.96\%
Table 2: Outcomes of the EGSEGS estimator according to pp and rr.

This empirical exercise using the daily returns for the MASI index between 15th November 2016 and 15th November 2017 is a historical approach that illustrates the practical use of EGSr,pλEGS_{r,p}^{\lambda} in the real world by considering the psychological attitude of the investor. The obtained results confirm earlier remarks in the previous subsection: first, more the investor is risk averse less risks he takes and then smaller is the amount of capital required to hedge his position; furthermore, we have EGSr,pλESpEGS_{r,p}^{\lambda}\geq ES_{p} and in a highly risk averse context (r20r\geq 20 for this survey) both risk measures can be confounded.

References

  • [1] Acerbi, C. (2002). Spectral measures of risk: A coherent representation of subjective risk aversion. Journal of Banking and Finance, 26, 1505-1518. doi: 10.1016/S0378-4266(02)00281-9
  • [2] Artzner, P., Delbaen, F., Eber, J.-M. and Heath, D. (1999). Coherent measures of risk. Mathematical Finance, 9, 203-228. doi: 10.1111/1467-9965.00068
  • [3] Ceriani, L. and Verme, P. (2012). The origins of the Gini index: extracts from Variabilità e Mutabilità (1912) by Corrado Gini. Journal of Economic Inequality, 10, 421-443. doi: 10.1007/s10888-011-9188-x
  • [4] Choquet, G. (1954). Theory of capacities. Annales de l’institut Fourier, 5, 131-295. doi: 10.5802/aif.53
  • [5] Dana, R.-A. (2005). A representation result for concave Schur-concave functions. Mathematical Finance, 15, 613-634. doi: 10.1111/j.1467-9965.2005.00253.x
  • [6] Delbaen, F. (2012): Monetary Utility Functions. Osaka University Press, Osaka.
  • [7] Embrechts, P., Puccetti, G., Rüschendorf, L., Wang, R. and Beleraj, A. (2014). An academic response to Basel 3.5. Risks, 2, 25-48. doi:10.3390/risks2010025
  • [8] Fischer, T. (2003). Risk capital allocation by coherent risk measures based on one-sided moments. Insurance: Mathematics and Economics 32, 135-146. doi: 10.1016/S0167-6687(02)00209-3
  • [9] Föllmer, H. and Schied, A. (2011). Stochastic Finance: An Introduction in Discrete Time. (Third Edition.) Walter de Gruyter, Berlin.
  • [10] Furman, E. and Landsman, Z. (2006). Tail variance premium with applications for elliptical portfolio of risks. ASTIN Bulletin: Journal of the International Actuarial Association, 36, 433-462. doi: 10.2143/AST.36.2.2017929
  • [11] Furman, E., Wang, R. and Zitikis, R. (2017). Gini-type measures of risk and variability: Gini shortfall, capital allocation and heavy-tailed risks. Journal of Banking and Finance. In Press, Accepted Manuscript. doi: 10.1016/j.jbankfin.2017.06.013
  • [12] Giorgi, G.M. (1990) Bibliographic portrait of the Gini concentration ratio. Metron, 48, 183-221.
  • [13] Giorgi, G.M. (1993) A fresh look at the topical interest of the Gini concentration ratio. Metron, 51, 83-98.
  • [14] Grechuk, B., Molyboha, A. and Zabarankin, M. (2009). Maximum entropy principle with general deviation measures. Mathematics of Operations Research, 34, 445-467. doi: 10.1287/moor.1090.0377
  • [15] Mao, T. and Wang, R. (2016). Risk aversion in regulatory capital principles. Working paper. doi: 10.2139/ssrn.2658669
  • [16] McNeil, A.J., Frey, R. and Embrechts, P. (2005). Quantitative Risk Management: Concepts, Techniques and Tools. (Revised Edition.) Princeton University Press, Princeton, NJ.
  • [17] Righi, M. and Ceretta, P. (2016). Shortfall Deviation Risk: an alternative to risk measurement. Journal of Risk 19, 81-116. doi: 10.21314/JOR.2016.349
  • [18] Righi, M. and Borenstein, D. (2017). A simulation comparison of risk measures for portfolio optimization. Finance Research Letters, in press. doi: 10.1016/j.frl.2017.07.013
  • [19] Righi, M. (2017). A composition between risk and deviation measures. Working Paper.
  • [20] Rockafellar, R. T., Uryasev, S. and Zabarankin, M. (2006). Generalized deviation in risk analysis. Finance and Stochastics, 10, 51-74. doi: 10.1007/s00780-005-0165-8
  • [21] Rüschendorf, L. (2013). Mathematical Risk Analysis. Dependence, Risk Bounds, Optimal Allocations and Portfolios. Springer, Heidelberg.
  • [22] Schmeidler, D. (1986). Integral representation without additivity. Proceedings of the American Mathematical Society, 97, 255-261. doi: 10.2307/2046508
  • [23] Schmeidler, D. (1989). Subject probability and expected utility without additivity. Econometrica, 57(3), 571-587. doi: 10.2307/1911053
  • [24] Shalit, H., and Yitzhaki, S. (1984). Mean-Gini, portfolio theory and the pricing of risky assets. Journal of Finance, 39(5 (December)), 1449-1468. doi: 10.1111/j.1540-6261.1984.tb04917.x
  • [25] Wang, R., Wei, Y. and Willmot., G.E. (2017). Characterization, Robustness and Aggregation of Signed Choquet Integrals. Working paper.
  • [26] Yaari, M.E. (1987). The dual theory of choice under risk. Econometrica, 55, 95-115. doi: 10.2307/1911158
  • [27] Yitzhaki, S. (1983) On an extension of Gini inequality index. International Economic Review, 24, no. 3, 617-628. doi: 10.2307/2648789
  • [28] Yitzhaki, S. (1998) More than a dozen alternative ways of spelling Gini. Research on Economic Inequality, 8, 13-30. doi: 10.1007/978-1-4614-4720-7
  • [29] Yitzhaki, S. and Schechtman, E. (2005). The properties of the extended Gini measures of variability and inequality. Metron, LXIII(3), 401-443. doi: 10.2139/ssrn.815564
  • [30] Yitzhaki, S. and Schechtman, E. (2013). The Gini Methodology: A Primer on a Statistical Methodology. Springer, New York, NY.