Next Article in Journal
The Determinants of CDS Spreads in Multiple Industry Sectors: A Comparison between the US and Europe
Next Article in Special Issue
Robust Estimations for the Tail Index of Weibull-Type Distribution
Previous Article in Journal
On the Evaluation of the Distribution of a General Multivariate Collective Model: Recursions versus Fast Fourier Transform
Previous Article in Special Issue
Bayesian Adjustment for Insurance Misrepresentation in Heavy-Tailed Loss Regression
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Some Results on Measures of Interaction between Paired Risks

1
Department of Mathematics, Shantou University, Shantou 515063, China
2
Department of Mathematical Sciences, Stevens Institute of Technology, Hoboken, NJ 07030, USA
*
Author to whom correspondence should be addressed.
Submission received: 2 August 2018 / Revised: 24 August 2018 / Accepted: 27 August 2018 / Published: 27 August 2018
(This article belongs to the Special Issue Heavy-Tail Phenomena in Insurance, Finance, and Other Related Fields)

Abstract

:
Co-risk measures and risk contribution measures have been introduced to evaluate the degree of interaction between paired risks in actuarial risk management. This paper attempts to study the ordering behavior of measures on interaction between paired risks. For various co-risk measures and risk contribution measures, we investigate how the marginal distributions and the dependence structure impact on the level of interaction between paired risks. Also, several numerical examples based on Monte Carlo simulation are presented to illustrate the main findings.

1. Introduction

In risk analysis and actuarial science, measuring risks and their contribution is of both theoretical and practical interest. In the literature, there exists a large amount of research on risk measures for the univariate risk. For a random risk X, at level α ( 0 , 1 ) the popular value-at-risk (VaR)
VaR α ( X ) = inf { x : P ( X x ) α } ,
and expected shortfall (ES)
ES α ( X ) = 1 1 - α α 1 VaR t ( X ) d t .
See (Denuit et al. 2005; Jorion 2007) for comprehensive overviews on risk measures of univariate risks. However, in actuarial risk analysis, the interaction among risks within a portfolio has been found to play an important role in determining the measure of risk and thus should not be ignored. After the recent financial crisis, the regulatory and scientific discussion on interaction among risks have been significantly intensified. The interaction among risks of institutions in a market or components in a portfolio is closely related to the systemic risk, which usually occurs when a component in a portfolio or even the whole portfolio suffers from the loss of some other component in the same portfolio. Typically, there are two types of measures for the degree of interaction among risks–co-risk measures and risk contribution measures. For example, the conditional value-at-risk (CoVaR) due to (Adrian and Brunnermeier 2016) and the conditional expected shortfall (CoES) due to (Mainik and Schaanning 2014) are two typical co-risk measures. Specifically, both CoVaR and CoES of a random risk of a component employ its conditional distribution given that another component (or even the whole portfolio) is under stress. For more on different co-risk measures, one may refer to (Kylman 2011; Huang et al. 2012) and references therein. On the other hand, the risk contribution measure intends to quantify the effect of a stressful situation for one component on the risk of another component through comparing a co-risk measure and its univariate version. For example, the difference between CoVaR of a component and the VaR of the same component. For more examples of risk contribution measures, we refer readers to (Girardi and Ergün 2013; Mainik and Schaanning 2014; Adrian and Brunnermeier 2016) and reference therein. Also, see (Bisias et al. 2012) for a recent extensive overview on different systemic risk measures.
Paired risks are not uncommon in actuarial science and insurance practice. For example, in life insurance, when modelling probability behaviour of the insured’s mortality, it is important to divide the insureds according to their marriage status. It has been confirmed that there exists dependence between lifetimes of a couple. For more discussions on the dependence of couple’s mortalities, one may refer to (Denuit et al. 2001; Jagger and Sutton 1991; Jevtić and Hurd 2017). When measuring the contribution of a component to systemic risk, usually the risk of this institution and the risk of another institution or even the whole system will be jointly considered as paired risks. See for example, (Acharya et al. 2017; Karimalis and Nomikos 2018) and references therein.
The existing studies on co-risk measure and risk contribution measure mainly investigate how the change of marginal distributions and dependence structure affects the paired risks. For instance, (Mainik and Schaanning 2014) studied the dependence consistency of CoVaR and some other systemic risks measures, and they found that given the marginal distribution unchanged, a more concordant dependence structure results in a larger CoVaR. Recently, (Sordo et al. 2018) extended this finding to CoES, Δ CoVaR and Δ CoES under positive dependence structure through showing that a more concordant dependence structure tends to incur a larger CoES, Δ CoVaR and Δ CoES, respectively, whenever the marginal distribution possesses some specific stochastic orderings. Most research along this line investigates how the degree of dependence impacts on the risk level of a given component in a risk pair. However, risk interaction is two-way. The stressful event of one risk in a pair has its own effect on the other risk, and vice versa. To obtain a whole picture on the risk level of paired risks, it is reasonable to consider the co-risk measures or risk contribution measures of either components simultaneously. Furthermore, it has been shown that different components in the same financial system usually possess different degrees of systemic importance and hence different risk levels. In (Zhou 2010) several measures of risk interaction were estimated from market data to check whether a bank of larger size has a smaller risk, and it was found that a big bank did not necessarily take a smaller risk of bankruptcy and different banks did have different risk levels. This reveals that paired risks with different marginal may have different level of risk interaction, and thus several interesting questions concerned with risk interaction levels of paired risks arise naturally: for financial institutions in a portfolio or financial system, whether does the stressful situation of one component affect the risk of another component more heavily than the stressful situation of the latter affects the risk of the former one? Also, how do the risk distribution and dependence structure play a role in the interaction between two risks?
This paper aims to investigate how the marginal distributions and the symmetry of dependence structure of two paired risks impact on their co-risk measures and risk contribution measures. For example, the popular CoVaR, CoES, and Δ CoVaR, Δ CoES. In this study we employ the copula functional to formulate the dependence structure between the paired risks. For the symmetric copula, the stochastically higher marginal risk is found to have a larger CoVaR, and for asymmetric copula with more probability mass on the region above (beneath) the diagonal line, a stochastically larger marginal risk is proved to attain a higher CoVaR at the stress level above (below) a threshold. Concerning the CoES, Δ CoVaR and Δ CoES , the larger marginal risk in the sense of increasing convex order, dispersive order and excess wealth order is shown to have a larger CoES, Δ CoVaR and Δ CoES respectively in the context of a positively dependent and symmetric copula. As for the asymmetric copula, similar comparing results are developed at the stress level above a threshold, which is determined by the stress level of the conditioned risk and the copula.
The rest of the present paper is organised as follows: Section 2 recalls several related notions and technical results concerned with detailed discussions in the sequel. Section 3 presents the comparison results on CoVaR and CoES, and the comparison results on Δ CoVaR and Δ CoES are developed in Section 4. Finally, in Section 5 several numerical examples based on Monte Carlo simulation are employed to illustrate the main findings.

2. Some Preliminaries

For ease of reference, let us present several related concepts and technical results to be used during the discussion in the coming sections.
Stochastic orders
Denote F X , F Y and F X - 1 , F Y - 1 probability distribution functions and quantile functions of random variables X and Y, respectively.
Definition 1.
A random variable X is said to be smaller than another one Y in the
  • usual stochastic order, denoted as X st Y , if F X ( t ) F Y ( t ) for all t R ;
  • increasing convex order, denoted as X icx Y , if E [ φ ( X ) ] E [ φ ( Y ) ] for any increasing convex function φ, provided the expectations exit;
  • dispersive order, denoted as X disp Y , if F X - 1 ( p ) - F X - 1 ( q ) F Y - 1 ( p ) - F Y - 1 ( q ) for all p , q such that 0 < q < p < 1 ;
  • excess wealth order, denoted as X ew Y , if E [ max { X - F X - 1 ( p ) , 0 } ] E [ max { Y - F Y - 1 ( p ) , 0 } ] for all p ( 0 , 1 ) .
The excess wealth order due to (Shaked and Shanthikumar 1998), also called the right spread order in (Fernandez-Ponce et al. 1998), is found to be useful in economics and reliability. See for example (Kochar et al. 2002). Stochastic orders are found useful in financial and actuarial risk, reliability and operation research, etc. For comprehensive discussions on stochastic orders readers may refer to monographs (Belzunce et al. 2015; Müller and Stoyan 2002; Shaked and Shanthikumar 2006). As a useful characterization of the increasing convex order, the following lemma due to Theorem 2.1 of (Sordo and Ramos 2007) will be used in deriving some of the main results.
Lemma 1.
For two random variables X and Y, X icx Y if and only if
0 1 F X - 1 ( t ) d ϕ ( t ) 0 1 F Y - 1 ( t ) d ϕ ( t ) ,
for any increasing convex function ϕ : [ 0 , 1 ] R .
Co-risk measures
The systemic risk measure is usually utilized to quantify a financial institution’s contribution to the risk of other financial institutions or even the entire financial system. Typically, the co-risk measures can be viewed as the traditional risk measures conditioning on a certain risky scenario affecting another institution.
Definition 2.
For two random variables X and Y, and α , β ( 0 , 1 ) ,
  • the CoVaR of Y at stress level β given that X is under stress at level α is
    CoVaR α , β ( Y | X ) = VaR β ( Y | X > VaR α ( X ) ) ;
  • the CoES of Y at stress level β given that X is under stress at level α is
    CoES α , β ( Y | X ) = 1 1 - β β 1 CoVaR α , t ( Y | X ) d t .
One may refer to (Mainik and Schaanning 2014; Adrian and Brunnermeier 2016; Karimalis and Nomikos 2018) and references therein for more details on co-risk measures and their applications in quantitative risk management.
Risk contribution measures
Apart from the co-risk measures, risk contribution measures are also introduced to quantify the systemic risk in risk management. Unlike co-risk measures, this class of measures captures the marginal effect of an institution to the risk of another institution or the overall financial system.
Definition 3.
For two random variables X and Y, and α , β ( 0 , 1 ) ,
  • the ΔCoVaR of Y at stress level β given that X is under stress at level α is
    Δ CoVaR α , β ( Y | X ) = CoVaR α , β ( Y | X ) - VaR β ( Y ) ;
  • the ΔCoES of Y at stress level β given that X is under stress at level α is
    Δ CoES α , β ( Y | X ) = CoES α , β ( Y | X ) - ES β ( Y ) .
For more on risk contribution measures and their applications, one may refer to (Mainik and Schaanning 2014; Sordo et al. 2018) and references therein.
Statistical dependence
It is well-known that statistical dependence is crucial in studying joint probability behavior of multiple risks and hence of the portfolio aggregation. Ever since the novel work of (Sklar 1959), the copula function has become the most popular tool in modelling the statistical dependence among concerned random variables due to mathematical and statistical tractability.
Definition 4.
For a random vector ( X , Y ) with joint distribution function F, and univariate marginal distribution functions F X , F Y , if there exists some C : [ 0 , 1 ] 2 [ 0 , 1 ] such that
F ( x , y ) = C F X ( x ) , F Y ( y ) , f o r   a l l   x , y ( - , + ) ,
then C ( u , v ) : [ 0 , 1 ] 2 [ 0 , 1 ] is called the copula of ( X , Y ) .
Definition 5.
A random variable X is said to be stochastically increasing in Y, denoted as SI ( X | Y ) , if
[ X | Y = y 1 ] st [ X | Y = y 2 ] , f o r   a l l   y 1 y 2 ,
and a random vector ( X , Y ) is said to be positive dependent through the stochastic order (PDS) if both SI ( X | Y ) and SI ( Y | X ) .
According to (Cai and Wei 2012), a random vector ( X , Y ) is PDS if and only if its copula is PDS. This characterization will be utilized to derive some results in the sequel. Note that for any absolutely continuous copula C, the PDS property is equivalent to
2 C ( u , v ) u 2 0 and 2 C ( u , v ) v 2 0 , for   any   u , v ( 0 , 1 ) .
For more on copula theory one may refer to the monograph (Nelsen 1999). Also we recall one more technical lemma due to Theorem 3.1(b) in (Mainik and Schaanning 2014), serving as an implicit expression of the conditional value-at-risk in terms of copula.
Lemma 2.
For a vector ( U , V ) having distribution function C ( u , v ) , the copula of ( X , Y ) , if F X is continuous, then
CoVaR α , β ( Y | X ) = F Y - 1 F V | U α - 1 ( β ) a n d F V | U α ( v ) = v - C ( α , v ) 1 - α ,
for any α , β , v ( 0 , 1 ) .
In order to compare the degree of risk interaction of paired risks, we consider measures CoVaR, CoES, Δ CoVaR and Δ CoES. Given α ( 0 , 1 ) , for any β ( 0 , 1 ) , denote v α ( β ) and v ^ α ( β ) the solutions of equations
v - C ( α , v ) 1 - α = β and v - C ( v , α ) 1 - α = β , for   v   on   ( 0 , 1 ) ,
respectively. If ( U , V ) C , then U, V are both uniformly distributed, and then, v α ( β ) and v ^ α ( β ) are the β th quantiles of the conditional random variables V | U α and U | V α , respectively.
Arrangement monotonicity
A real function g ( x 1 , x 2 ) is said to be arrangement increasing (AI) if
g ( x 1 , x 2 ) g ( x 2 , x 1 ) , for   any x 1 < x 2 .
The function g is said to be arrangement decreasing (AD) when the inequality (1) is reversed.
Arrangement monotone functions are useful in risk management and operations research. For more on AI functions and various stochastic versions of AI, one may refer to (Boland and Proschan 1988; Li and Li 2018). In this study we will employ AI and AD copulas as two typical alternatives for the symmetric copula. Note that an AI (AD) copula is of more (less) dependence above the diagonal line on the plane.

3. Co-Risk Measures

Let us first investigate how the magnitude of marginal risks affects the corresponding CoVaR. In specific, we propose several sufficient conditions for
CoVaR α , β ( Y | X ) CoVaR α , β ( X | Y ) .
Theorem 1.
For a random vector ( X , Y ) with copula C ( u , v ) ,
  • if C ( u , v ) is symmetric, then X st Y is equivalent to (2) for α , β ( 0 , 1 ) ,
  • if C ( u , v ) is AI, then X st Y implies (2) for α , β ( 0 , 1 ) such that β [ α - C ( α , α ) ] / ( 1 - α ) , and
  • if C ( u , v ) is AD, then X st Y implies (2) for α , β ( 0 , 1 ) such that β [ α - C ( α , α ) ] / ( 1 - α ) .
Proof. 
Case 1: Let ( U , V ) = st ( F X ( X ) , F Y ( Y ) ) . Note that X st Y is equivalent to
F X - 1 ( p ) F Y - 1 ( p ) , for   any   p ( 0 , 1 ) .
Due to the symmetry of C, ( Y , X ) gets the copula
C ^ ( u , v ) = P ( F Y ( Y ) u , F X ( X ) v ) = P ( V u , U v ) = C ( v , u ) = C ( u , v ) ,
for u , v ( 0 , 1 ) . According to Lemma 2, we have, for α , β ( 0 , 1 ) ,
CoVaR α , β ( Y | X ) = F Y - 1 F V | U α - 1 ( β ) ,
and
CoVaR α , β ( X | Y ) = F X - 1 F U | V α - 1 ( β ) ,
where
F V | U α ( v ) = v - C ( α , v ) 1 - α = v - C ( v , α ) 1 - α = F U | V α ( v ) .
Owing to (3), X st Y is equivalent to
CoVaR α , β ( Y | X ) CoVaR α , β ( X | Y ) , for   α , β ( 0 , 1 ) .
Case 2: By Lemma 2 again we have
CoVaR α , β ( Y | X ) - CoVaR α , β ( X | Y ) = F Y - 1 ( v α ( β ) ) - F X - 1 ( v ^ α ( β ) ) , for   α , β ( 0 , 1 ) .
In light of (3), it suffices for us to verify that v α ( β ) v ^ α ( β ) for α , β ( 0 , 1 ) .
Note that, for any α , v ( 0 , 1 ) ,
F V | U α ( v ) = v - C ( α , v ) 1 - α and F U | V α ( v ) = v - C ( v , α ) 1 - α .
The AI property of C implies C ( u , v ) C ( v , u ) for 0 u v 1 , and hence it holds that
v ^ α ( β ) - C ( v ^ α ( β ) , α ) 1 - α = β = v α ( β ) - C ( α , v α ( β ) ) 1 - α v α ( β ) - C ( v α ( β ) , α ) 1 - α ,
for α , β ( 0 , 1 ) such that v α ( β ) α . This, in combination with the nondecreasing property of the distribution function, yields v α ( β ) v ^ α ( β ) . Similarly, for α , β ( 0 , 1 ) such that v ^ α ( β ) α , it holds that
v α ( β ) - C ( α , v α ( β ) ) 1 - α = β = v ^ α ( β ) - C ( v ^ α ( β ) , α ) 1 - α v ^ α ( β ) - C ( α , v ^ α ( β ) ) 1 - α ,
and hence we have v α ( β ) v ^ α ( β ) . Since v α ( β ) α and v ^ α ( β ) α are equivalent to
β = v α ( β ) - C ( α , v α ( β ) ) 1 - α α - C ( α , α ) 1 - α ,
and
β = v ^ α ( β ) - C ( v ^ α ( β ) , α ) 1 - α α - C ( α , α ) 1 - α ,
respectively, it holds that min { v α ( β ) , v ^ α ( β ) } α if and only if β [ α - C ( α , α ) ] / ( 1 - α ) . Therefore, we conclude that
CoVaR α , β ( Y | X ) CoVaR α , β ( X | Y ) ,
for α , β ( 0 , 1 ) such that β [ α - C ( α , α ) ] / ( 1 - α ) .
Case 3 can be verified in a similar manner to Case 2. □
According to Theorem 1, when paired risks have exchangeable copula, namely, the interaction between them is symmetric, the only factor determining the level of co-risk measure is the marginal distribution. When the interaction between risks is asymmetric, namely one risk may have more impact on the other one, both the marginal distribution and the dependence structure play a role in determining the level of co-risk measure. Also, under either one of the two scenarios a stochastically larger marginal risk intends to incur a larger co-risk.
As per Theorem 1, if one of the concerned paired risks is stochastically larger and the underlying copula of the two risks is AI, then its CoVaR conditioning on the other one may serve as a consistent upper bound on the systemic risk of this pair, especially in the situation that both risks are under high stress levels.
As a continued discussion on the effect of marginal distribution and dependence structure on the CoVaR, in what follows we have a parallel discussion on the CoES, which is actually the average of CoVaR on stress levels above a threshold. To be specific, we propose some conditions sufficient to
CoES α , β ( Y | X ) CoES α , β ( X | Y ) ,
for some typical α , β ( 0 , 1 ) .
Theorem 2.
For a random vector ( X , Y ) with a PDS copula C ( u , v ) ,
  • if C ( u , v ) is symmetric, then X icx Y implies (4) for all α , β ( 0 , 1 ) , and
  • if C ( u , v ) is AI, then X icx Y implies (4) for α , β ( 0 , 1 ) such that β [ α - C ( α , α ) ] / ( 1 - α ) .
Proof. 
According to the proof of Theorem 12 in (Sordo et al. 2018),
CoES α , β ( Y | X ) = 0 1 F Y - 1 ( t ) d ϕ α , β ( t ) ,
where
ϕ α , β ( t ) = 1 - 1 - α - t + C ( α , t ) ( 1 - α ) ( 1 - β ) I t > v α ( β ) .
is a convex function and the indicator function I ( A ) equals 1 if A is true and 0 if not. Similarly,
CoES α , β ( X | Y ) = 0 1 F X - 1 ( t ) d ϕ ^ α , β ( t ) ,
ϕ ^ α , β ( t ) = 1 - 1 - α - t + C ( t , α ) ( 1 - α ) ( 1 - β ) I t > v ^ α ( β ) .
Case 1: Due to the symmetric C ( u , v ) , it holds that ϕ α , β ( t ) = ϕ ^ α , β ( t ) for any α , β , t ( 0 , 1 ) . For a PDS C ( u , v ) , the ϕ α , β ( t ) is shown to be increasing and convex in the proof of Theorem 12 in (Sordo et al. 2018). Then, by Lemma 1, it follows from X icx Y that
CoES α , β ( Y | X ) - CoES α , β ( X | Y ) = 0 1 F Y - 1 ( t ) d ϕ α , β ( t ) - 0 1 F X - 1 ( t ) d ϕ ^ α , β ( t ) = 0 1 F Y - 1 ( t ) d ϕ α , β ( t ) - 0 1 F X - 1 ( t ) d ϕ α , β ( t ) 0 , for   all   α , β ( 0 , 1 ) .
Case 2: In a similar manner to the proof of Theorem 1, for α , β ( 0 , 1 ) , one can show that
β α - C ( α , α ) 1 - α v α ( β ) v ^ α ( β ) ,
and hence
I t > v α ( β ) I t > v ^ α ( β ) .
For any t > min { v α ( β ) , v ^ α ( β ) } , one has C ( α , t ) C ( t , α ) and hence ϕ α , β ( t ) ϕ ^ α , β ( t ) . Using that ϕ α , β ( 1 ) = ϕ ^ α , β ( 1 ) = 1 , we have, for any u ( 0 , 1 ) ,
u 1 d [ ϕ α , β ( t ) - ϕ ^ α , β ( t ) ] = ϕ ^ α , β ( u ) - ϕ α , β ( u ) 0 .
Noting that F X - 1 ( t ) is increasing in t, it follows from Lemma 7.1 of (Barlow and Proschan 1981) that
0 1 F X - 1 ( t ) d [ ϕ α , β ( t ) - ϕ ^ α , β ( t ) ] 0 .
On the other hand, similar to the proof of Case 1, we can verify that the PDS C ( u , v ) results in the increasing and convex ϕ α , β ( t ) . Thus, X icx Y implies
0 1 F Y - 1 ( t ) d ϕ α , β ( t ) - 0 1 F X - 1 ( t ) d ϕ α , β ( t ) 0 ,
and then it holds that,
0 1 F Y - 1 ( t ) d ϕ α , β ( t ) 0 1 F X - 1 ( t ) d ϕ α , β ( t ) 0 1 F X - 1 ( t ) d ϕ ^ α , β ( t ) .
That is,
CoES α , β ( Y | X ) CoES α , β ( X | Y ) ,
for α , β ( 0 , 1 ) such that β [ α - C ( α , α ) ] / ( 1 - α ) . □
Since CoES measures the averaged CoVaR at stress levels above some threshold, intuitively, a more divergent loss (in the sense of the increasing convex order) should result in a larger CoES when the concerned risk is under a sufficient high level of stress. This intuition is confirmed by the finding of Theorem 2.
It worths remarking that unlike the equivalence in the first assertion of Theorem 1, the sufficient condition of Theorem 2 under the symmetric copula case is not necessary. When the copula is symmetric, the inequality (4) is equivalent to that
0 1 F Y - 1 ( t ) d ϕ α , β ( t ) 0 1 F X - 1 ( t ) d ϕ α , β ( t ) ,
where ϕ α , β ( t ) defined in (5) is some specific increasing and convex function. But the equivalence of X icx Y provided in Lemma 1 requires that the above inequality holds for any increasing and convex function. Thus, the first assertion in Theorem 2 is not necessary.

4. Risk Contribution Measures

In this section, we switch our focus to the risk contribution measures Δ CoVaR and Δ CoES. In parallel to the previous study on co-risk measure, we bring forth several sufficient conditions for the following ordering result.
Δ CoVaR α , β ( Y | X ) Δ CoVaR α , β ( X | Y ) .
Theorem 3.
For a random vector ( X , Y ) with a PDS copula C ( u , v ) , if C ( u , v ) is
  • symmetric, then X disp Y implies (7) for α , β ( 0 , 1 ) ,
  • AI, then X disp Y implies (7) for α , β ( 0 , 1 ) such that β [ α - C ( α , α ) ] / ( 1 - α ) , and
  • AD, then X disp Y implies (7) for α , β ( 0 , 1 ) such that β [ α - C ( α , α ) ] / ( 1 - α ) .
Proof. 
Note that (7) is equivalent to
F Y h , α - 1 ( β ) - F Y - 1 ( β ) F X g , α - 1 ( β ) - F X - 1 ( β ) ,
where X g , α and Y h , α are distorted versions of X and Y induced by the concave distortion transforms
h α ( x ) = x - α + C ( α , 1 - x ) 1 - α and g α ( x ) = x - α + C ( 1 - x , α ) 1 - α ,
respectively. Moreover, it holds that, for α , β ( 0 , 1 ) ,
F Y h , α - 1 ( β ) = F Y - 1 ( 1 - h α - 1 ( 1 - β ) ) and F X g , α - 1 ( β ) = F X - 1 ( 1 - g α - 1 ( 1 - β ) ) .
Case 1: Note that the exchangeability of C implies that h α ( x ) = g α ( x ) . From Lemma 14 in (Sordo et al. 2018) the desired result follows immediately.
Case 2: Note that h α ( x ) = 1 - β and g α ( x ) = 1 - β can be rephrased as
1 - x - C ( α , 1 - x ) 1 - α = β and 1 - x - C ( 1 - x , α ) 1 - α = β ,
respectively. By the definitions of v α ( β ) and v ^ α ( β ) , we have
1 - h α - 1 ( 1 - β ) = v α ( β ) and 1 - g α - 1 ( 1 - β ) = v ^ α ( β ) .
When α , β are such that β α - C ( α , α ) 1 - α , it is shown in the proof of Theorem 1 that v ^ α ( β ) v α ( β ) .
Therefore, we can conclude that
F Y h , α - 1 ( β ) - F Y - 1 ( β ) F Y - 1 ( v ^ α ( β ) ) - F Y - 1 ( β ) F X - 1 ( v ^ α ( β ) ) - F X - 1 ( β ) = F X g , α - 1 ( β ) - F Y - 1 ( β ) .
Case 3: can be verified in a similar manner to Case 2. □
The risk contribution measure Δ CoVaR α , β ( Y | X ) quantifies the effect on Y brought forth by risk X under stress. As shown in (Sordo et al. 2018), under a stronger positive dependence structure, a more dispersive risk will result in a larger Δ CoVaR. Following the same logic, Theorem 3 reveals that for paired risks with symmetric or asymmetric dependence, a more dispersive marginal risk incurs a larger risk contribution measure at a higher or lower stress level, depending on the specific structure of the copula. It should be pointed out that for symmetric copula, (7) is equivalent to
F Y h , α - 1 ( β ) - F Y - 1 ( β ) F X h , α - 1 ( β ) - F X - 1 ( β ) ,
a special case of the necessary condition for X disp Y presented in Lemma 14 of (Sordo et al. 2018). Thus, the first assertion of Theorem 3 is sufficient but not necessary.
Next result presents sufficient conditions for the inequality
Δ CoES α , β ( Y | X ) Δ CoES α , β ( X | Y ) .
Theorem 4.
For a random vector ( X , Y ) with PDS copula C ( u , v ) ,
  • if C ( u , v ) is symmetric, then X ew Y implies (8) for α , β ( 0 , 1 ) , and
  • if C ( u , v ) is AI, then X ew Y implies (8) for α , β ( 0 , 1 ) such that β [ α - C ( α , α ) ] / ( 1 - α ) .
Proof. 
Similar to the proof of Theorem 23 in (Sordo et al. 2018), we have, for α , β ( 0 , 1 ) ,
Δ CoES α , β ( Y | X ) = 0 1 F Y - 1 ( t ) d ϕ α , β ( t ) - 0 1 F Y - 1 ( t ) d ψ β ( t ) ,
and
Δ CoES α , β ( X | Y ) = 0 1 F X - 1 ( t ) d ϕ ^ α , β ( t ) - 0 1 F X - 1 ( t ) d ψ β ( t ) ,
where ϕ α , β ( t ) is given in (5), ϕ ^ α , β ( t ) is given in (6), and ψ β ( t ) = max { 0 , ( t - β ) / ( 1 - β ) } . Note that copula C ( u , v ) is PDS and X ew Y implies
0 1 F Y - 1 ( t ) d ϕ α , β ( t ) - 0 1 F Y - 1 ( t ) d ψ β ( t ) 0 1 F X - 1 ( t ) d ϕ α , β ( t ) - 0 1 F X - 1 ( t ) d ψ β ( t ) .
Case 1: As is shown in the first assertion of Theorem 2, the exchangeability of C implies that ϕ α , β ( t ) = ϕ ^ α , β ( t ) , and hence the conclusion follows from (9).
Case 2: As is shown in the proof of the second assertion of Theorem 2, the AI copula C ( u , v ) along with β [ α - C ( α , α ) ] / ( 1 - α ) implies that ϕ α , β ( t ) ϕ ^ α , β ( t ) , and hence it holds that
0 1 F X - 1 ( t ) d ϕ α , β ( t ) - 0 1 F X - 1 ( t ) d ψ β ( t ) 0 1 F X - 1 ( t ) d ϕ ^ α , β ( t ) - 0 1 F X - 1 ( t ) d ψ β ( t ) .
Thus, we reach the conclusion by combining (9) and (10). □
Roughly speaking, the risk contribution measure based on CoES is an average of the risk contribution measure based on CoVaR. Naturally, similar property of Δ CoVaR may be inherited by the Δ CoES. Theorem 4 confirms this intuition and provides two scenarios, where a larger risk in the sense of excess wealth order leads to a larger Δ CoES. Note that for symmetric copula one can check that the inequality in (8) is actually equivalent to (9), is a special case of the general condition equivalent to X ew Y presented in Theorem 8 (ii) of (Sordo 2008). Therefore, the first assertion of Theorem 4 is not necessary.

5. Numerical Examples Based on Simulation

To close this study, we carry out Monte Carlo experiments to illustrate the main findings in the previous sections. In general, the joint distribution function of dependent paired risks is complex and the corresponding co-risk measures and risk contribution measures do not have explicit closed form. Lemma 2 facilitates us to obtain a sample of observations of the concerned risks. Moreover, it is well-known that for most regular parametric distributions, the error rate of using the empirical distribution and quantile to approximate the population distribution and quantile is bounded by the reciprocal squared root of the sample size. See for reference, (Reiss 1974; Van der Vaart 1998). Thus, we can illustrate the theoretical findings by using the sample version of co-risk and risk contribution measures based on large number of observations of concerned risks.
Throughout all experiments, we start with the Gumbel copula
C ˜ θ ( u , v ) = exp - ( ( - log u ) θ + ( - log v ) θ ) 1 / θ , θ 1 ,
which is clearly symmetric. As per (5.43) of (McNeil et al. 2005), the following modification
C θ ( u , v ) = u 1 - α v 1 - β C ˜ θ ( u α , v β )
is asymmetric for 0 α β 1 and θ 1 . For ease of reference, we call C θ ( u , v ) the asymmetric Gumbel copula with dependence parameters θ = ( θ , α , β ) .
In all simulation experiments, we set θ = ( 2 , 0 . 2 , 0 . 8 ) for the copula C θ ( u , v ) of the concerned paired risks. It is observed that C θ ( u , v ) C θ ( v , u ) for u v in Figure 1a and C θ ( u , v ) C θ ( v , u ) for v u in Figure 1b. Also, as is seen Figure 2 both second-order derivatives of C θ ( u , v ) with respect to u and v are non-positive, and this verifies that C θ ( u , v ) is PDS.
To produce a sample of the paired risks ( X , Y ) , we first generate a bivariate sample from the population having the copula of ( X , Y ) as the joint distribution, and then we calculate the marginal quantiles at this sample to obtain a desired sample of ( X , Y ) based on Definition 4.
Note that, if we directly simulate samples of paired risks following a given joint distribution and estimate the population CoVaR and other risk measures, there will be few suitable observations when the stress level is close to 1, the extreme level, which is often the concern of risk management. Without enough suitable observations, the estimated values could be very different from their true versions. To circumvent this rareness, we don’t directly draw observations from the joint population distribution of the paired risks. Instead, we first obtain a random sample of one margin U, and then generate a corresponding sample from the population of the conditional random variable V | U > α , where ( U , V ) has joint distribution C θ , the copula of ( X , Y ) , and α is a given stress level of U. The validness of this sample procedure is guaranteed by Lemma 2. Precisely, the simulation is carried out in accordance with the following steps.
  • For each stress level α = 1 101 , , 100 101 and i = 1 , , 10 5 , generate a sample of 10 5 observations v i , α of ( V | U > α ) and a sample of 10 5 observations u i , α of ( U | V > α ) , respectively.
  • For α = 1 101 , , 100 101 , based on v i , α ’s and u i , α ’s calculate respectively the adjusted empirical distribution functions
    F ^ ( t ) = 1 10 5 + 1 i = 1 10 5 I ( v i , α t ) and G ^ ( t ) = 1 10 5 + 1 i = 1 10 5 I ( u i , α t ) .
  • At each stress level α = 1 101 , , 100 101 , for each stress level β = 1 101 , , 100 101 , employ the sample β th-quantiles
    CoVaR ^ α , β ( Y | X ) = inf t : F ^ ( t ) α and CoVaR ^ α , β ( X | Y ) = inf t : G ^ ( t ) α
    to estimate the measures CoVaR α , β ( Y | X ) and CoVaR α , β ( X | Y ) , respectively. Also, calculate the sample versions
    CoES ^ α , β ( Y | X ) = 1 ( 1 + 10 5 ) ( 1 - β ) i = 1 10 5 I v i , α > CoVaR ^ α , β ( Y | X ) ,
    CoES ^ α , β ( X | Y ) = 1 ( 1 + 10 5 ) ( 1 - β ) i = 1 10 5 I u i , α > CoVaR ^ α , β ( X | Y ) ,
    for CoES α , β ( Y | X ) and CoES α , β ( X | Y ) , respectively.
  • As for the risk contribution measures Δ CoVaR and Δ CoES , the following empirical estimators are used.
    Δ CoVaR ^ α , β ( Y | X ) = CoVaR ^ α , β ( Y | X ) - VaR β ( Y ) ,
    Δ CoVaR ^ α , β ( X | Y ) = CoVaR ^ α , β ( X | Y ) - VaR β ( X ) ,
    Δ CoES ^ α , β ( Y | X ) = CoES ^ α , β ( Y | X ) - ES β ( Y ) ,
    Δ CoES ^ α , β ( X | Y ) = CoES ^ α , β ( X | Y ) - ES β ( X ) .
    It should be remarked here that to mitigate the approximation error we simply use the population version for the marginal VaR and ES when deriving these estimators.
In what follows, let us consider four scenarios of ( X , Y ) having different combinations of marginal distributions.
  • For X E ( 2 ) and Y E ( 1 ) , exponentially distributions with parameters 2 and 1, respectively, it is plain that X st Y . By the second assertion of Theorem 1, we have
    CoVaR α , β ( Y | X ) CoVaR α , β ( X | Y ) , for   all   β [ α - C θ ( α , α ) ] / ( 1 - α ) .
    Figure 3a plots CoVaR ^ α , β ( Y | X ) - CoVaR ^ α , β ( X | Y ) . The difference surface is seen to be always above the horizontal surface on the region { ( α , β ) [ 0 , 1 ] 2 : β [ α - C θ ( α , α ) ] / ( 1 - α ) } , and this confirms the finding of Theorem 1. On the other hand, although the copula C θ ( u , v ) fails to satisfy condition in the third assertion of Theorem 1, the surface above the horizontal surface on the region { ( α , β ) [ 0 , 1 ] 2 : β [ α - C θ ( α , α ) ] / ( 1 - α ) } in Figure 3a hints that the third assertion of Theorem 1 may still be true when the requirement on dependence structure is violated.
  • For X N ( 0 , 1 ) and Y N ( 0 , 2 ) , two normal distributions, according to Table 1.1 of (Müller and Stoyan 2002), X icx Y is valid. By the second assertion of Theorem 2, CoES α , β ( Y | X ) CoES α , β ( X | Y ) for β [ α - C θ ( α , α ) ] / ( 1 - α ) . This is illustrated by the surface of Figure 3b. Actually, the corresponding difference is still nonnegative when β [ α - C θ ( α , α ) ] / ( 1 - α ) .
  • For X W ( 1 , 1 ) and Y W ( 1 . 5 , 1 ) , two Weibull distributions, Example 16 of (Sordo et al. 2018) proves that X disp Y . Figure 3c plots the difference between Δ CoVaR ^ α , β ( Y | X ) and Δ CoVaR ^ α , β ( X | Y ) , and this difference surface confirms the finding of the second assertion of Theorem 3.
  • For X W ( 1 , 2 ) and Y W ( 1 , 1 ) , two Weibull distributions, as per Example 24 of (Sordo et al. 2018), we have X ew Y . Under this setting, the second assertion of Theorem 4 proves that Δ CoES α , β ( Y | X ) Δ CoES α , β ( X | Y ) for β [ α - C θ ( α , α ) ] / ( 1 - α ) . Also, this fact is justified by the difference surface of Figure 3d.
The findings in the previous sections assert that the effect of stressful situations of the larger marginal risk (in the sense of some stochastic order) on the other marginal risk is more significant than that of the smaller marginal risk on the larger one when the copula of the risk pair is AI, namely, the dependence structure possesses some asymmetry. However, it is not clear how large such a difference can be. As illustrated in the numerical examples, the difference varies quite differently and may be large at the extreme scenario, where both the stress levels α , β are close to 1. Note that larger value of the risk measures usually indicates higher risk, and a larger difference implies that one of the marginal risk may suffer more from the risky situation of the other one. Hence, in the practice of risk management concerning risk spread inside paired risks, one may treat the paired risks as a whole, and pay attention to one’s failure spreading to the other and vice versa, instead of just concerning one specific risk’s effect.

6. Concluding Remarks

Co-risk measures and risk contribution measures have been introduced to measure the degree of systemic risk ever since the financial crisis. Capturing the dependence of paired risks these measures reflect the interaction between the corresponding marginal risks. The present study contributes to the existing literature on risk interaction by showing that the underlying dependence structure plays a vital role in determining the level of risk interactions of the marginal risks between paired risks. Specifically, the larger marginal risk in the sense of some stochastic orders is found to have a larger CoVaR, CoES, Δ CoVaR and Δ CoES, respectively, at high stress levels (extreme scenarios usually concerned in risk management) in the context of AI copula. Also, the same phenomenon is observed at low stress levels for CoVaR and Δ CoVaR in the setting of AD copula. Therefore, the level of risk interaction of components in a risk pair is usually different as long as the marginal risks are of different distributions or the dependence structure is asymmetric.
This study reiterates the importance of accurately modelling the dependence structure of paired risks in practice. Under certain dependence structure, marginal risks in a pair could posses different levels of risk interaction even when they are homogenous. Furthermore, the effect of the stressful situation of one risk on the other one may be consistently larger than that vice versa at the extreme case. This indicates that in a portfolio the failure of one component may potentially spread faster or impact more heavily on another component than does the failure of other components. Our findings show the existence of a potential difference between interaction levels of marginal risks in a pair; However, several things are to be exploited: What is the exact difference? Are there any bounds on the difference under specific dependence structures or marginal distributions? Answers to these questions can provide more guidance to apply the co-risk measures and risk contribution measures into the practice of risk management. It is of both theoretical and practical interest to pursue these unsolve problems in the future study.

Author Contributions

Conceptualization, R.F. and X.L.; Methodology, R.F. and X.L.; Validation, R.F. and X.L.; Formal Analysis, R.F. and X.L.; Writing-Original Draft Preparation, R.F.; Writing-Review and Editing, R.F. and X.L.; Funding Acquisition, R.F.

Funding

This research was funded by National Natural Science Foundation of China grant number 11601303; Guangdong Natural Science Foundation grant number 2016A030310076.

Acknowledgments

The authors would like to thank the anonymous reviewers for their valuable comments, which helped improve the presentation of this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Acharya, Viral V., Lasse H. Pedersen, Thomas Philippon, and Matthew Richardson. 2017. Measuring systemic risk. The Review of Financial Studies 30: 2–47. [Google Scholar] [CrossRef]
  2. Adrian, Tobias, and Markus K. Brunnermeier. 2016. CoVaR. The American Economic Reviews 106: 1705–41. [Google Scholar] [CrossRef]
  3. Barlow, Richard E., and Frank Proschan. 1981. Statistical Theory of Reliability and Life Testing: Probability Models. New York: Holt, Rinehart and Winston. [Google Scholar]
  4. Belzunce, Felix, Carolina Martínez-Riquelme, and Julio Mulero. 2015. An Introduction to Stochastic Orders. London: Academic Press. [Google Scholar]
  5. Bisias, Dimitrios, Mark Flood, Andrew W. Lo, and Stavros Valavanis. 2012. A survey of systemic risk analytics. Annual Review of Financial Economics 4: 255–96. [Google Scholar] [CrossRef]
  6. Boland, Philip J., and Frank Proschan. 1988. Multivariate arrangement increasing functions with applications in probability and statistics. Journal of Multivariate Analysis 25: 286–98. [Google Scholar] [CrossRef]
  7. Cai, Jun, and Wei Wei. 2012. On the invariant properties of notions of positive dependence and copulas under increasing transformations. Insurance: Mathematics and Economics 50: 43–49. [Google Scholar] [CrossRef]
  8. Denuit, Michel, Jan Dhaene, Céline Le Bailly de Tilleghem, and Stéphanie Teghem. 2001. Measuring the impact of dependence among insured lifelengths. Belgian Actuarial Bulletin 1: 18–39. [Google Scholar]
  9. Denuit, Michel, Jan Dhaene, Mark Goovaerts, and Rob Kaas. 2005. Actuarial Theory for Dependent Risks—Measures, Orders and Models. Chichester: John Wiley & Sons. [Google Scholar]
  10. Fernandez-Ponce, Jose Maria, Subhash C. Kochar, and José Muñoz-Perez. 1998. Partial orderings of distributions based on right-spread function. Journal of Applied Probability 35: 221–28. [Google Scholar] [CrossRef]
  11. Girardi, Giulio, and A. Tolga Ergün. 2013. Systemic risk measurement: Multivariate GARCH estimation of CoVaR. Journal of Banking and Finance 37: 3169–80. [Google Scholar] [CrossRef]
  12. Huang, Xin, Hao Zhou, and Haibin Zhu. 2012. Systemic risk contributions. Journal of Financial Services Research 42: 55–83. [Google Scholar] [CrossRef]
  13. Jagger, Carol, and Christopher J. Sutton. 1991. Death after marital bereavement-is the risk increased? Statistics in Medicine 10: 395–404. [Google Scholar] [CrossRef] [PubMed]
  14. Jevtić, Petar, and Thomas R. Hurd. 2017. The joint mortality of couples in continuous time. Insurance: Mathematics and Economics 75: 90–97. [Google Scholar] [CrossRef]
  15. Jorion, Philippe. 2007. Value at Risk: The New Benchmark for Managing Financial Risk, 3rd ed.New York: McGraw-Hill Education. [Google Scholar]
  16. Karimalis, Emmanouil N., and Nikos Nomikos. 2018. Measuring systemic risk in the European banking sector: A copula CoVaR approach. The European Journal of Finance 24: 944–75. [Google Scholar] [CrossRef]
  17. Klyman, Jared. 2011. “Systemic Risk Measures: DistVaR and Other “Too Big To Fail” Risk Measures”. Ph.D. thesis, Princeton University, Princeton, NJ, USA. [Google Scholar]
  18. Kochar, Subhash C., Xiaohu Li, and Moshe Shak. 2002. The total time on test transform and the excess wealth stochastic orders of distributions. Advances in Applied Probability 34: 826–45. [Google Scholar] [CrossRef]
  19. Li, Chen, and Xiaohu Li. 2018. Stochastic arrangement increasing risks in financial engineering and actuarial science—A review. Quantitative Finance and Economics 2: 675–701. [Google Scholar] [CrossRef]
  20. Mainik, Georg, and Eric Schaanning. 2014. On dependence consistency of CoVaR and some other systemic risk measures. Statistics and Risk Modeling 31: 49–77. [Google Scholar] [CrossRef]
  21. McNeil, Alexander J., Rüdiger Frey, and Paul Embrechts. 2005. Quantitative Risk Management: Concepts, Techniques and Tools. Princeton: Princeton University Press. [Google Scholar]
  22. Müller, Alfred, and Dietrich Stoyan. 2002. Comparison Methods for Stochastic Models and Risks. New York: Springer. [Google Scholar]
  23. Nelsen, Roger B. 1999. An Introduction to Copulas. New York: Springer. [Google Scholar]
  24. Reiss, Rolf-Dieter. 1974. On the accuracy of the normal approximation for quantiles. Annals of Probability 2: 741–44. [Google Scholar] [CrossRef]
  25. Shaked, Moshe, and Jeyaveerasingam George Shanthikumar. 2006. Stochastic Orders. New York: Springer. [Google Scholar]
  26. Shaked, Moshe, and Jeyaveerasingam George Shanthikumar. 1998. Two variability orders. Probability in the Engineering and Informational Sciences 12: 1–23. [Google Scholar] [CrossRef]
  27. Sklar, Abe. 1959. Fonctions de répartition á n dimensions et leurs marges. Publications de l’Institut Statistique de l’Université de Paris 8: 229–31. [Google Scholar]
  28. Sordo, Miguel A. 2008. Characterizations of classes of risk measures by dispersive orders. Insurance: Mathematics and Economics 42: 1028–34. [Google Scholar] [CrossRef]
  29. Sordo, Miguel A., Alfonso J. Bello, and Alfonso Suárez-Llorens. 2018. Stochastic orders and co-risk measures under positive dependence. Insurance: Mathematics and Economics 78: 105–13. [Google Scholar] [CrossRef]
  30. Sordo, Miguel A., and Héctor A. Ramos. 2007. Characterizations of stochastic orders by L-functionals. Statistical Papers 48: 249–63. [Google Scholar] [CrossRef]
  31. Van der Vaart, Adrianus Willem. 1998. Asymptotic Statistics. Cambridge: Cambridge University Press. [Google Scholar]
  32. Zhou, Chen. 2010. Are banks too big to fail? Measuring systemic importance of financial institutions. International Journal of Central Banking 6: 205–50. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Surface of ratio C θ ( u , v ) / C θ ( v , u ) . (a) 0 u v 1 ; (b) 0 v u 1 .
Figure 1. Surface of ratio C θ ( u , v ) / C θ ( v , u ) . (a) 0 u v 1 ; (b) 0 v u 1 .
Risks 06 00088 g001
Figure 2. Surfaces of second-order derivatives. (a) 2 C θ ( u , v ) / u 2 ; (b) 2 C θ ( u , v ) / v 2 .
Figure 2. Surfaces of second-order derivatives. (a) 2 C θ ( u , v ) / u 2 ; (b) 2 C θ ( u , v ) / v 2 .
Risks 06 00088 g002
Figure 3. Difference surfaces (red) and β = [ α - C θ ( α , α ) ] / ( 1 - α ) (green). (a) CoVaR ^ α , β ( Y | X ) - CoVaR ^ α , β ( X | Y ) ; (b) CoES ^ α , β ( Y | X ) - CoES ^ α , β ( X | Y ) ; (c) Δ CoVaR ^ α , β ( Y | X ) - Δ CoVaR ^ α , β ( X | Y ) ; (d) Δ CoES ^ α , β ( Y | X ) - Δ CoES ^ α , β ( X | Y ) .
Figure 3. Difference surfaces (red) and β = [ α - C θ ( α , α ) ] / ( 1 - α ) (green). (a) CoVaR ^ α , β ( Y | X ) - CoVaR ^ α , β ( X | Y ) ; (b) CoES ^ α , β ( Y | X ) - CoES ^ α , β ( X | Y ) ; (c) Δ CoVaR ^ α , β ( Y | X ) - Δ CoVaR ^ α , β ( X | Y ) ; (d) Δ CoES ^ α , β ( Y | X ) - Δ CoES ^ α , β ( X | Y ) .
Risks 06 00088 g003

Share and Cite

MDPI and ACS Style

Fang, R.; Li, X. Some Results on Measures of Interaction between Paired Risks. Risks 2018, 6, 88. https://0-doi-org.brum.beds.ac.uk/10.3390/risks6030088

AMA Style

Fang R, Li X. Some Results on Measures of Interaction between Paired Risks. Risks. 2018; 6(3):88. https://0-doi-org.brum.beds.ac.uk/10.3390/risks6030088

Chicago/Turabian Style

Fang, Rui, and Xiaohu Li. 2018. "Some Results on Measures of Interaction between Paired Risks" Risks 6, no. 3: 88. https://0-doi-org.brum.beds.ac.uk/10.3390/risks6030088

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop