Next Article in Journal
The Outperformance Probability of Mutual Funds
Next Article in Special Issue
Equity Option Pricing with Systematic and Idiosyncratic Volatility and Jump Risks
Previous Article in Journal
Bank Competition, Foreign Bank Entry, and Risk-Taking Behavior: Cross Country Evidence
Previous Article in Special Issue
Examination and Modification of Multi-Factor Model in Explaining Stock Excess Return with Hybrid Approach in Empirical Study of Chinese Stock Market
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

CVaR Regression Based on the Relation between CVaR and Mixed-Quantile Quadrangles

1
V.M. Glushkov Institute of Cybernetics, 40, pr. Akademika Glushkova, 03187 Kyiv, Ukraine
2
Applied Mathematics & Statistics, Stony Brook University, B-148 Math Tower, Stony Brook, NY 11794, USA
*
Author to whom correspondence should be addressed.
J. Risk Financial Manag. 2019, 12(3), 107; https://doi.org/10.3390/jrfm12030107
Submission received: 16 May 2019 / Revised: 19 June 2019 / Accepted: 20 June 2019 / Published: 26 June 2019
(This article belongs to the Special Issue Mathematical Finance with Applications)

Abstract

:
A popular risk measure, conditional value-at-risk (CVaR), is called expected shortfall (ES) in financial applications. The research presented involved developing algorithms for the implementation of linear regression for estimating CVaR as a function of some factors. Such regression is called CVaR (superquantile) regression. The main statement of this paper is: CVaR linear regression can be reduced to minimizing the Rockafellar error function with linear programming. The theoretical basis for the analysis is established with the quadrangle theory of risk functions. We derived relationships between elements of CVaR quadrangle and mixed-quantile quadrangle for discrete distributions with equally probable atoms. The deviation in the CVaR quadrangle is an integral. We present two equivalent variants of discretization of this integral, which resulted in two sets of parameters for the mixed-quantile quadrangle. For the first set of parameters, the minimization of error from the CVaR quadrangle is equivalent to the minimization of the Rockafellar error from the mixed-quantile quadrangle. Alternatively, a two-stage procedure based on the decomposition theorem can be used for CVaR linear regression with both sets of parameters. This procedure is valid because the deviation in the mixed-quantile quadrangle (called mixed CVaR deviation) coincides with the deviation in the CVaR quadrangle for both sets of parameters. We illustrated theoretical results with a case study demonstrating the numerical efficiency of the suggested approach. The case study codes, data, and results are posted on the website. The case study was done with the Portfolio Safeguard (PSG) optimization package, which has precoded risk, deviation, and error functions for the considered quadrangles.

1. Introduction

We start the introduction with a quick outline of the main result of this paper. The conditional value-at-risk (CVaR) is a popular risk measure. It is called expected shortfall (ES) in financial applications and it is included in financial regulations. This paper provides algorithms for the estimation of CVaR with linear regression as a function of factors. This task is of critical importance in practical applications involving low probability events.
By definition, CVaR is an integral of the value-at-risk (VaR) in the tail of a distribution. VaR can be estimated with the quantile regression by minimizing the Koenker–Bassett error function. This paper shows that CVaR can be estimated by minimizing a mixture of the Koenker–Bassett errors with an additional constraint. This mixture is called the Rockafellar error and it has been earlier used for CVaR estimation without a rigorous mathematical justification. One more equivalent variant of CVaR regression can be done by minimizing a mixture of CVaR deviations for finding all coefficients, except the intercept. In this case, the intercept is calculated using an analytical expression, which is the CVaR of the optimal residual without an intercept. The new mathematical result links quantile and CVaR regressions and shows that convex and linear programming methods can be straightforwardly used for CVaR estimation. Mathematical justification of the results involves a risk quadrangle concept combining regret, error, risk, deviation, and statistic notions.
Quantiles evaluating different parts of a distribution of a random value are quite popular in various applications. In particular, quantiles are used to estimate tail of a distribution (e.g., 90%, 95%, and 99% quantiles). This paper is motivated by finance applications, where a quantile is called VaR. Risk measure VaR is included in finance regulations for the estimation of market risk. VaR has several attractive properties, such as the simplicity of calculation, stability of estimation, and availability of quantile regression, for the estimation of VaR as a function of explanatory factors. The quantile regression (see Koenker and Bassett (1978), Koenker (2005)) is an important factor supporting the popularity of VaR. For instance, a quantile regression was used by Adrian and Brunnermeier (2016) to estimate institution’s contribution to systemic risk.
However, VaR also has some undesirable properties:
  • Lack of convexity: portfolio diversification may increase VaR.
  • VaR is not sensitive to outcomes exceeding VaR, which allows for stretching of the distribution without an increasing of the risk measured by VaR.
  • VaR has poor mathematical properties, such as discontinuity with respect to (w.r.t.) portfolio positions for discrete distributions based on historical data.
Shortcomings of VaR led financial regulators to use an alternative measure of risk, which is called conditional value-at-risk (CVaR) in this paper. This risk measure was introduced in Rockafellar and Uryasev (2000) and further studied in Rockafellar and Uryasev (2002) and many other papers. CVaR for continuous distributions equals the conditional expectation of losses exceeding VaR. An important mathematical fact is that CVaR is a coherent risk measure (see Acerbi and Tasche (2002), Rockafellar and Uryasev (2002)). Ziegel (2014) shows that CVaR is elicitable in a week sense. Fissler and Ziegel (2015) proved that (VaR, CVaR) is jointly elicitable, meaning elic(CVaR) ≤ 2, and more generally, that spectral risk measures have a low elicitation complexity. These results clarify the regression procedure of Rockafellar et al. (2014); their algorithm implicitly tracks the quantiles suggested by elicitation complexity.
Rockafellar and Uryasev (2000, 2002) have shown that CVaR of a convex function of variables is also a convex function. Due to this property, CVaR optimization problems can be reduced to convex and linear optimization problems.
This paper is based on risk quadrangle theory, which defines quadrangles (i.e., groups) of stochastic functionals Rockafellar and Uryasev (2013). Every quadrangle contains risk, deviation, error, and regret (negative utility). These elements of the quadrangle are linked by the statistic function.
The relation of quantile regression and CVaR optimization was explained using a quantile quadrangle (see Rockafellar and Uryasev (2013)). It was shown that the Koenker–Bassett error function and CVaR belong to the same quantile quadrangle. By minimizing the Koenker–Bassett error function with respect to one parameter, we obtain the CVaR deviation (which is the CVaR for the centered random value). The optimal value of the parameter, which is called statistic, equals VaR. Therefore, the linear regression with the Koenker–Bassett error estimates VaR as a function of factors. The fact that statistic equals VaR and is also used for building the optimization approach for CVaR (see Rockafellar and Uryasev (2000, 2002)).
Another important contribution that takes advantage of quadrangle theory is the regression decomposition theorem proved in Rockafellar et al. (2008). With this decomposition theorem, the regression problem is decomposed in two steps: (1) minimization of deviation from the corresponding quadrangle, and (2) calculation of the intercept by using statistic from this quadrangle. For instance, by applying the decomposition theorem to the quantile quadrangle, we can do quantile regression by minimizing CVaR deviation for finding all regression coefficients, except the intercept. Then, the intercept is calculated by using VaR statistic.
CVaR can be approximated using the weighted average of VaRs with different confidence levels, which is called the mixed VaR method. Rockafellar and Uryasev (2013) demonstrated that mixed VaR is a statistic in the mixed-quantile quadrangle. The error function, corresponding to this quadrangle (called the Rockafellar error) can be minimized for the estimation of the mixed VaR with linear regression. The Rockafellar error is a solution of a minimization problem with one linear constraint. Linear regression for estimating mixed VaR can be done by minimizing the Rockafellar error with convex and linear programming (Appendix A contains these formulations). Alternatively, this regression can be done in two steps with the decomposition theorem. The deviation in the mixed-quantile quadrangle is the mixed CVaR deviation, therefore all regression coefficients, except the intercept, can be found by minimizing this deviation. Further, the intercept can be found by using statistic, which is the mixed VaR.
Rockafellar et al. (2014) developed the CVaR quadrangle with the statistic equal to CVaR. Risk envelopes and identifiers for this quadrangle were calculated in Rockafellar and Royset (2018). This CVaR quadrangle is a theoretical basis for constructing the regression for estimating CVaR. Rockafellar et al. (2014) called the linear regression for estimation of CVaR using superquantile (CVaR) regression. The superquantile is an equivalent term for CVaR. Here we use the term “CVaR regression”. The CVaR regression plays a major role in various engineering areas, especially in financial applications. For instance, Huang and Uryasev (2018) used CVaR regression for the estimation of risk contributions of financial institutions and Beraldi et al. (2019) used CVaR for solving portfolio optimization problems with transaction costs.
This paper considers only discrete random values with a finite number of equally probable atoms. This special case is considered because it is needed for the implementation of the linear regression for the CVaR estimation. We have explained with an example how parameters of the optimization problems are calculated.
The equal probabilities property was used for calculating parameters of optimization problems. It is possible to calculate parameters with non-equal probabilities of atoms, but this is beyond the scope of the paper, which is focused on the linear regression.
We suggested two sets (Sets 1 and 2) of parameters for the mixed-quantile quadrangle. Set 1 corresponds to the two-step implementation of the CVaR regression in Rockafellar et al. (2014), and Set 2 is a new set of parameters. We proved that with Set 1, the statistic, risk, and deviation of the mixed-quantile and CVaR quadrangles coincide. Therefore, CVaR regression can be done by minimizing the Rockafellar error with convex and linear programming. For Set 2, the mixed-quantile and CVaR quadrangle share risk and deviation parameters. Also, the statistic of this mixed-quantile quadrangle (which may not be unique) includes statistic of the CVaR quadrangle. Therefore, minimizing the Rockafellar error correctly calculates all regression coefficients, but may provide an incorrect intercept. This is actually not a big concern because we know that the intercept is equal to the CVaR of an optimal residual without intercept.
Also, we demonstrated that the CVaR regression can be done in two steps with the decomposition theorem by using parameters from Sets 1 and 2 in the mixed-quantile deviation. A similar two-step procedure was used for CVaR regression in Rockafellar et al. (2014). Here we justify this two-step procedure through the equivalence of deviations in CVaR and mixed-quantile quadrangles with parameters from Sets 1 and 2.
This paper is organized as follows. Section 2 provides general results about quadrangles. In particular, we considered quantile, mixed-quantile, and CVaR quadrangles. Section 3 and Section 4 introduced and investigated the parameters from Sets 1 and 2, respectively. Section 5 provided optimization problem statements based on CVaR and mixed-quantile quadrangles and described the linear regression for CVaR estimation. Section 6 presented a case study and applied CVaR regression to the financial style classification problem. The case study is posted on the web with codes, data, and solutions. Appendix A provides convex and linear programming problems for minimization of the Rockafellar error; Appendix B provides Portfolio Safeguard (PSG) codes implementing regression optimization problems.

2. Quantile, Mixed-Quantile, and CVaR Quadrangles

Rockafellar and Uryasev (2013) developed a new paradigm called the risk quadrangle, which linked risk management, reliability, statistics, and stochastic optimization theories. The risk quadrangle methodology united risk functions for a random value X in groups (quadrangles) consisting of five elements:
  • Risk ( X ) , which provides a numerical surrogate for the overall hazard in X .
  • Deviation D ( X ) , which measures the “nonconstancy” in X as its uncertainty.
  • Error ( X ) , which measures the “nonzeroness” in X .
  • Regret V ( X ) , which measures the “regret” in facing the mix of outcomes of X .
  • Statistic S ( X ) associated with X through and V .
These elements of a risk quadrangle are related as follows:
V ( X ) = ( X ) + E ( X )
( X ) = D ( X ) + E ( X )
( X ) = min C { C + V ( X C ) }
D ( X ) = min C { ( X C ) }
argmin C { C + V ( X C ) } = S ( X ) = arg min C { ( X C ) }
where E ( X ) denotes the mean of X and the statistic, S ( X ) , can be a set, if the minimum is achieved for multiple points.
Further, we use the following notations. The cumulative distribution function is denoted by F X ( x ) = p r o b { X x } . The positive and negative part of a number are denoted using:
[ t ] + = { t ,   for   t > 0 0 ,   for   t 0     and     [ t ] = { t ,   for   t < 0 0 ,   for   t 0
The lower and upper VaR (quantile) are defined as follows:
lower VaR:
V a R α ( X ) = { s u p { x ,   F X ( x ) < α } f o r   0 < α 1 i n f { x ,   F X ( x ) α }   f o r   α = 0
upper VaR:
V a R α + ( X ) = { i n f { x ,   F X ( x ) > α }   f o r   0 α < 1 s u p { x ,   F X ( x ) α }   f o r   α = 1
VaR (quantile) is a set if the lower and upper quantiles do not coincide:
V a R α ( X ) = [ V a R α ( X ) , V a R α + ( X ) ]
otherwise VaR is a singleton V a R α ( X ) = V a R α ( X ) = V a R α + ( X ) .
Conditional value-at-risk (CVaR) with the confidence level α ( 0 , 1 ) can be defined in many ways. We prefer the following constructive definition:
C V a R α ( X ) = min C { C + 1 1 α E [ X C ] + }
In financial applications, however, the most popular definition of CVaR is
C V a R α ( X ) = 1 1 α α 1 V a R β ( X ) d β .
For α = 0 , C V a R 0 ( X ) is defined as C V a R 0 ( X ) = lim ε 0   C V a R ε ( X ) = E ( X ) .
For α = 1 ,   C V a R 1 ( X ) is defined as C V a R 1 ( X ) = V a R 1 ( X ) if a finite value of V a R 1 ( X ) exists.
Quadrangles are named after statistic functions. The most famous quadrangle is the quantile quadrangle (see Rockafellar and Uryasev (2013)), named after the VaR (quantile) statistic. This quadrangle establishes relations between the CVaR optimization technique described in Rockafellar and Uryasev (2000, 2002) and quantile regression (see Koenker and Bassett (1978), Koenker (2005)). In particular, it was shown that CVaR minimization and the quantile regression are similar procedures based on the VaR statistic in the regret and error representation of risk and deviation.
Here is the definition of the quantile quadrangle for α ( 0 , 1 ) :
Quantile Quadrangle (Rockafellar and Uryasev (2013))
  • Statistic: S α ( X ) = V a R α ( X ) = VaR (quantile) statistic.
  • Risk: α ( X ) = C V a R α ( X ) = min C { C + V α ( X C ) } = CVaR risk.
  • Deviation: D α ( X ) = C V a R α ( X ) E [ X ] = min C { α ( X C ) } = CVaR deviation.
  • Regret: V α ( X ) = 1 1 α E [ X ] + = average absolute loss, scaled.
  • Error: α ( X ) = E [ α 1 α [ X ] + + [ X ] ] =   V α ( X ) E [ X ] = normalized Koenker–Bassett error.
The quantile quadrangle sets an example for development of more advances quadrangles. The following mixed-quantile quadrangle includes statistic, which is equal to the weighted average of VaRs (quantiles) with specified positive weights. Therefore, the error in this quadrangle can be used to build a regression for the weighted average of VaRs (quantiles). Since CVaR can be approximated by a weighted average of VaRs, the error function in this quadrangle can be used to build linear regression for the estimation of CVaR.
Mixed-Quantile Quadrangle (Rockafellar and Uryasev (2013)).
Confidence levels α k ( 0 , 1 ) ,   k = 1 , , r , and weights λ k > 0 , k = 1 r λ k = 1 . The error in this quadrangle is called the Rockafellar Error.
  • Statistic: S ( X ) = k = 1 r λ k V a R α k ( X ) =   mixed VaR (quantile).
  • Risk: ( X ) = k = 1 r λ k C V a R α k ( X ) = mixed CVaR.
  • Deviation: D ( X ) = k = 1 r λ k C V a R α k ( X E [ X ] ) = mixed CVaR deviation.
  • Regret: V ( X ) = min B 1 , , B r { k = 1 r λ k V α k ( X B k ) | k = 1 r λ k B k = 0 } = the minimal weighted average of regrets V α k ( X B k ) = 1 1 α k E [ X B k ] + satisfying the linear constraint on B 1 , , B r .
  • Error: ( X ) = min B 1 , , B r { k = 1 r λ k α k ( X B k ) | k = 1 r λ k B k = 0 } = Rockafellar error = the minimal weighted average of errors α k ( X B k ) = E [ α k 1 α k [ X B k ] + + [ X B k ] ] satisfying the linear constraint on B 1 , , B r .
The following CVaR quadrangle can be considered as the limiting case of the mixed-quantile quadrangle when the number of terms in this quadrangle tends to infinity. The statistic in this quadrangle is CVaR; therefore, the error in this quadrangle can be used for the estimation of CVaR with linear regression.
CVaR Quadrangle (Rockafellar et al. (2014)) for α ( 0 , 1 ) .
  • Statistic: S ¯ α ( X ) = C V a R α ( X ) = CVaR.
  • Risk: ¯ α ( X ) = 1 1 α α 1 C V a R β ( X ) d β =   CVaR2 risk.
  • Deviation: D ¯ α ( X ) = ¯ α ( X ) E [ X ] = 1 1 α α 1 C V a R β ( X ) d β E [ X ] = CVaR2 deviation.
  • Regret: V ¯ α ( X ) = 1 1 α 0 1 [ C V a R β ( X ) ] + d β = CVaR2 regret.
  • Error: ¯ α ( X ) = V ¯ α ( X ) E [ X ] = CVaR2 error.
The following section proves that for a discretely distributed random value with equally probable atoms, the CVaR quadrangle is “equivalent” to a mixed-quantile quadrangle with some parameters in the sense that statistic, risk, and deviation in these quadrangles coincide. This fact was proved for a set of random values with equal probabilities and variable locations of atoms.
The set of parameters considered in the following section is used in two-step CVaR regression in Rockafellar et al. (2014).

3. Set 1 of Parameters for Mixed-Quantile Quadrangle

Set 1 of parameters for the mixed-quantile quadrangle for a discrete uniformly distributed random value X consists of confidence levels α k ( 0 , 1 ) , k = 1 , , r and weights λ k > 0 such that k = 1 r λ k = 1 . Parameter r depends only on the number of atoms in X and the confidence level α of the CVaR quadrangle. We proved that statistic, risk, and deviation of the mixed-quantile quadrangle with the Set 1 of parameters coincide with the statistic, risk, and deviation of the CVaR quadrangle.
Let X be a discrete random value with support x i and Prob ( X = x i ) = 1 / ν for i = 1 , 2 , , ν , where ν is the number of atoms. Denote x m a x = max i = 1 , , ν x i . For this random value, C V a R 1 ( X ) = V a R 1 ( X ) = x m a x .
Set 1 of parameters:
  • partition of the interval [ α , 1 ] : β ν α 1 = α , and β i = i δ , for i = ν α ,   ν α + 1 , , ν , where δ = 1 / ν , ν α = ν α + 1 , with z being the largest integer less than or equal to z ; δ α = β ν α α .
  • weights: p ν α = δ α 1 α ,   p i = δ 1 α ,     i = ν α + 1 , , ν .
  • confidence levels: γ i = 1 β i β i 1 ln ( 1 β i 1 1 β i ) , i =   ν α , , ν 1 ; γ ν = 1 .
Lemma 1.
Let X be a discrete random value with ν equally probable atoms. Then, statistic, risk, and deviation of the CVaR quadrangle for X are given by the following expressions with parameters specified by Set 1:
1. 
CVaR statistic:
S ¯ α ( X ) = C V a R α ( X ) = i = ν α ν p i V a R γ i ( X )
2. 
CVaR2 risk:
¯ α ( X ) = 1 1 α α 1 C V a R β ( X ) d β = i = ν α ν p i C V a R γ i ( X )
3. 
CVaR2 deviation:
D ¯ α ( X ) = 1 1 α α 1 C V a R β ( X ) d β E [ X ] = i = ν α ν p i C V a R γ i ( X ) E [ X ]
Proof. 
Appendix C contains proof of the lemma. □
Note. 
Expression (1) is valid for arbitrary γ i ( β i 1 , β i ) , i = ν α , , ν . Equations (2) and (3) are valid for arbitrary γ ν [ β ν 1 , 1 ] .
We want to emphasize that the statement of Lemma 1 is valid for any discrete random value with equally probable atoms. The statement does not depend upon atom locations.
Corollary 1.
For the random value X defined in Lemma 1, statistic, risk, and deviation of the CVaR quadrangle coincide with statistic, risk, and deviation of the mixed-quantile quadrangle with r = ν ν α + 1 , λ k = p ν α 1 + k , α k = γ ν α 1 + k , k = 1 , , r .
Proof. 
Right hand sides in Equations (1)–(3) define statistic, risk, and deviation of the mixed-quantile quadrangle because p i > 0 ,   i = ν α , , ν , i = ν α ν p i = 1 and V a R γ i ( X ) ,   i = ν α , , ν , are singletons. □
Example 1.
Let X be a discrete random value with five atoms (−40; −10; 20; 60 100) and equal probabilities = 0.2.
Figure 1 explains how to calculate statistic in the CVaR quadrangle and mixed-quantile quadrangle with Set 1 of parameters for α = 0.5 . Bold lines show V a R α ( X ) as a function of α . C V a R 0.5 ( X ) equals the dark area under the V a R α ( X ) divided by 1 α . CVaR can be calculated as integral of V a R α ( X ) or as the sum of areas of rectangles. Figure 2 explains how to calculate risk in the CVaR quadrangle and mixed-quantile quadrangle with Set 1 of parameters for α = 0.5 . The bold continuous curve shows C V a R α ( X ) as a function of α . Risk ¯ 0.5 ( X ) is equal to the area under the CVaR curve divided by 1 α . This area can be calculated as the integral of CVaR or as the sum of areas of rectangles. The area of every rectangle is equal to the area under CVaR in the appropriate range of α . The equality of areas defines values of γ i . Parameters p i ,   γ i do not depend on the values of atoms.

4. Set 2 of Parameters for the Mixed-Quantile Quadrangle

This section gives an alternative expression for the risk ¯ α ( X ) and deviation D ¯ α ( X ) in the CVaR quadrangle for a discrete uniformly distributed random variable. This expression is based on the following Set 2 of parameters. This set of parameters has the same number of parameters as Set 1 but different values of weights and confidence levels. Similar to Section 3, let X be a discrete random value with support, x i , i = 1 , 2 , , ν , and Prob ( X = x i ) = 1 / ν for i = 1 , 2 , , ν . Denote x m a x =   max   i = 1 , , ν x i . For this random value C V a R 1 ( X ) = V a R 1 ( X ) = x m a x .
Set 2 of parameters:
  • partition of the interval [ β ν α 1 , 1 ] : δ = 1 / ν , β i = i δ , for i = ν α 1 ,   ν α , , ν , where ν α = ν α + 1 , with z being the largest integer less than or equal to z ; δ α = β ν α α .
  • confidence levels: β i , i = ν α 1 ,   ν α , , ν .
  • weights:
    q ν = 0 ;
    q ν 1 = δ 1 α × [ 2 ln ( 2 ) ] δ 1 α × 1.386294361 , ( if   ν 1 > ν α )
    q ν 2 = δ 1 α × 2 [ 3 ln ( 3 2 ) + ln ( 1 2 ) ] δ 1 α × 1.046496288 , ( if   ν 2 > ν α )
    q ν j = δ 1 α × j [ ( j + 1 ) ln ( j + 1 j ) + ( j 1 ) ln ( j 1 j ) ] , ( if   j > 2 , ν j > ν α )
    q ν α   = δ 1 α × j [ δ δ α + ( j + 1 ) ln ( 1 α δ j ) + ( j 1 ) ln ( j 1 j ) ] , ( if   ν α < ν 1 ,   j = ν ν α )
    q ν α 1 = δ 1 α × j [ δ α + ( j 1 ) ln ( δ ( j 1 ) 1 α ) ] , ( if   ν α 1 < ν 1 ,   j = ν ν α + 1 )
    if ν α = ν 1 , then q ν α = δ 1 α × 2 [ 1 + ln ( 1 α δ ) ] 1 , q ν α 1 = δ 1 α × 2 [ ln ( δ 1 α ) 1 ] + 2
    if ν α = ν , then q ν α 1 = 1 .
Lemma 2.
Let X be a discrete random value with ν equally probable atoms. Then, risk and deviation of the CVaR quadrangle for X are given by the following expressions with parameters from Set 2.
1. 
CVaR2 Risk:
¯ α ( X ) = 1 1 α α 1 C V a R β ( X ) d β = i = ν α 1 ν 1 q i C V a R β i ( X )
2. 
CVaR2 Deviation:
D ¯ α ( X ) = 1 1 α α 1 C V a R β ( X ) d β E [ X ] = i = ν α 1 ν 1 q i C V a R β i ( X ) E [ X ]
Proof. 
Appendix D contains proof of the lemma. □
Note. 
Equations (4) and (5) are valid if C V a R β ν 1 is replaced by C V a R γ with an arbitrary γ [ β ν 1 , 1 ] .
Corollary 2.
For the random value X defined in Lemma 2, risk and deviation of the CVaR quadrangle coincide with risk and deviation of the mixed-quantile quadrangle with r = ν ν α + 1 , λ k = q ν α 2 + k , α k = β ν α 2 + k , k = 1 , , r .
Proof. 
Right hand sides in Equations (4) and (5) define risk and deviation of the mixed-quantile quadrangle because q i > 0 ,   i = ν α 1 , , ν 1 , and i = ν α 1 ν 1 q i = 1 . □
Lemma 3.
Let X be a discrete random value with equally probable atoms x i , i = 1 , 2 , , ν , P r o b ( X = x i ) = 1 / ν . Then, statistic of the mixed-quantile quadrangle defined by the Set 2 of parameters is a range containing statistic of the CVaR quadrangle.
Proof. 
Appendix E contains proof of the lemma. □

5. On the Estimation of CVaR with Mixed-Quantile Linear Regression

This section formulates regression problems using the CVaR quadrangle and mixed-quantile quadrangle. For discrete final distributions with equally probable atoms, we prove some equivalence statements for the CVaR and mixed-quantile quadrangles. Further, we demonstrate how to estimate CVaR by using the linear regression with error and deviation from the mixed-quantile quadrangle.
We want to estimate variable V using a linear function f ( Y ) = C 0 + C T Y of the explanatory factors Y = ( Y 1 , , Y n ) . Let ˜ be an error from some quadrangle (further we consider, mixed-quantile and CVaR quadrangles), and D ˜ and S ˜ be a deviation and a statistic, respectively, corresponding to this quadrangle. Below we consider optimization statements for solving regression problems.
General Optimization Problem 1
Minimize error ˜ and find optimal C 0 * ,   C *
m i n C 0 ,   C m ˜   ( Z ( C 0 , C ) )
where Z ( C 0 , C ) = V C 0 C T Y .
General Optimization Problem 2
  • Step 1. Find an optimal vector C * by minimizing deviation:
    m i n   C m D ˜ ( Z 0 ( C ) )
    where Z 0 ( C ) = V C T Y .
  • Step 2. Assign C 0 * :
    C 0 * S ˜ ( Z 0 ( C * ) )
Error ˜   ( X ) is called nondegenerate if:
  i n f X : E X = D ˜   ( X ) > 0   for   constants   D 0 .
Rockafellar et al. (2008), p. 722, proved the following decomposition theorem.
Theorem 1.
(Error-Shaping Decomposition of Regression). Let ˜ be a nondegenerate error and D ˜ = min C { ˜ ( X C ) } be the corresponding deviation, and let S ˜ be the associated statistic. Point ( C 0 * , C * ) is a solution of the General Optimization Problem 1 if and only if C * is a solution of the General Optimization Problem 2, Step 1 and C 0 * S ˜ ( Z 0 ( C * ) ) with Step 2.
According to the decomposition theorem, when ˜ , D ˜ , and S ˜ are elements of the CVaR quadrangle, the following Optimization Problems 1 and 2 are equivalent.
Optimization Problem 1
Minimize error from the CVaR quadrangle:
  m i n C 0 ,   C m ¯ α ( Z ( C 0 , C ) )
where Z ( C 0 , C ) = V C 0 C T Y .
Optimization Problem 2
  • Step 1. Find an optimal vector C * by minimizing deviation from the CVaR quadrangle:
    m i n C m   D ¯ α ( Z 0 ( C ) )
    where Z 0 ( C ) = V C T Y .
  • Step 2. Calculate:
    C 0 * = C V a R α ( Z 0 ( C * ) )
In Optimization Problem 2 in Step 2 the statistic equals S ¯ ( Z 0 ( C * ) ) = C V a R α ( Z 0 ( C * ) ) , which is the specification of the inclusion operation in the Optimization Problem General 2 in Step 2.
Optimization Problem 2 is used in Rockafellar et al. (2014) for the construction of the linear regression algorithms for estimating CVaR.
According to the decomposition theorem, when ˜ , D ˜ , and S ˜ are elements of a mixed-quantile quadrangle, the following Optimization Problems 3 and 4 are equivalent.
Optimization Problem 3
Minimize error from the mixed-quantile quadrangle:
m i n C 0 ,   C m ( Z ( C 0 , C ) )
Optimization Problem 4
  • Step 1. Find an optimal vector C * by minimizing deviation from the mixed-quantile quadrangle:
    m i n   C m   D ( Z 0 ( C ) )
  • Step 2. Assign:
    C 0 * S ( Z 0 ( C * ) )
Corollaries 1 and 2 can be used for constructing the linear regression for estimating CVaR. Let Y be a random vector of factors for estimating the random value V . We consider that the linear regression function f ( Y ) = C 0 + C T Y approximates CVaR of V , where C 0 ,   C m are variables in the linear regression. The residual is denoted by Z ( C 0 , C ) = V ( C 0 + C T Y ) and Z 0 ( C ) = V C T Y .
Further we provide a lemma about linear regression problems based on Corollary 1 with the Set 1 of parameters. The main statement here is that the Optimization Problems 3 and 4 for the mixed-quantile quadrangle can be used to solve linear regression problems for estimating CVaR. This is the case because CVaR and mixed-quantile quadrangles have the same Statistic and Deviation.
Lemma 4.
Let the residual random value Z ( C 0 , C ) = V ( C 0 + C T Y ) be discretely distributed with ν equally probable atoms. Let us consider the CVaR quadrangle with error ¯ α ( Z ( C 0 , C ) ) , deviation D ¯ α ( Z 0 ( C ) ) , and statistic S ¯ α ( Z 0 ( C * ) ) . Let us also consider the mixed-quantile quadrangle with the error ( Z ( C 0 , C ) ) , deviation D ( Z 0 ( C ) ) , and statistic S ( Z 0 ( C * ) ) with parameters defined by Set 1 and r = ν ν α + 1 , λ k = p ν α 1 + k , α k = γ ν α 1 + k , k = 1 , , r . Then, the Optimization Problems 1–4 are equivalent, i.e., the sets of optimal vectors of these optimization problems coincide. Moreover, let ( C 0 * , C * ) be a solution vector of the equivalent Optimization Problems 1–4. Then:
C 0 * = C V a R α ( Z 0 ( C * ) )
¯ α ( C 0 * , C * ) =   ( C 0 * , C * ) =   D ( Z 0 ( C * ) ) = D ¯ α ( Z 0 ( C * ) )
Proof. 
This lemma is a direct corollary of the decomposition Theorem 1 and Corollary 1 of Lemma 1. Indeed, Corollary 1 implies that the Optimization Problems 2 and 4 are equivalent. Further, the decomposition theorem implies that Optimization Problems 1 and 2 and the Optimization Problems 3 and 4 are equivalent. □
Further we provide a lemma about linear regression problems based on Corollary 2 with the Set 2 of parameters. The main statement is that the Optimization Problem 4, Step 1 for the Mixed-Quantile Quadrangle can be used to solve linear regression problem for estimating CVaR. This is the case because CVaR and Mixed-Quantile Quadrangles have the same Deviation. After obtaining vector of coefficients C * , intercept is calculated, C 0 * = C V a R α ( Z 0 ( C * ) ) .
Lemma 5.
Let the residual random value Z ( C 0 , C ) = V ( C 0 + C T Y ) be discretely distributed with ν equally probable atoms. Let the mixed-quantile quadrangle with deviation D ( Z 0 ( C ) ) be defined by parameters of Set 2 and r = ν ν α + 1 , λ k = q ν α 2 + k , α k = β ν α 2 + k , k = 1 , , r . Then, ( C 0 * , C * ) is a solution of Optimization Problem 1, if and only if, ( C 0 * , C * ) is a solution of the following two-step procedure:
  • Step 1. Find an optimal vector C * by minimizing deviation from the mixed-quantile quadrangle:
    m i n   C m   D ( Z 0 ( C ) )
  • Step 2. Calculate C 0 * = C V a R α ( Z 0 ( C * ) ) .
Proof. 
This lemma is a direct corollary of the decomposition Theorem 1 and Corollary 2 or Lemma 2. Indeed, the Corollary 2 implies that the Optimization Problems 2 and 4 are equivalent. Further, since deviations of the CVaR and mixed-quantile quadrangles coincide, we can use Step 1 to calculate optimal coefficients C * . Further, the intercept is calculated with C 0 * = C V a R α ( Z 0 ( C * ) ) because CVaR is the statistic in the CVaR quadrangle. □
For the Set 2 of parameters, the deviations in the CVaR and mixed-quintile quadrangles coincide. The two-step procedure in Optimization Problem 4 can be used to solve linear regression problems with Set 2 parameters for the mixed-quantile deviation. Also, the minimization of the Rockafellar error with the Set 2 of parameters may result in a correct C * . The statistic of CVaR belongs to the statistic of the mixed-quantile quadrangle. Therefore, the optimization of the Rockafellar error with the Set 2 of parameters may lead to a wrong value of intercept C 0 * . This potential incorrectness can be fixed by assigning C 0 * = C V a R α ( Z 0 ( C * ) ) .

6. Case Study: Estimation of CVaR with Linear Regression and Style Classification of Funds

The case study described in this section is posted online (see Case Study (2016)). The codes and data are available for downloading and verification. Every optimization problem is presented in three formats: Text, MATLAB, and R. Calculations were done with a PC with a 3.14 GHz processor.
We have applied CVaR regression to the return-based style classification of a mutual fund. We regress a fund return by several indices as explanatory factors. The estimated coefficients represent the fund’s style with respect to each of the indices.
A similar problem with a standard regression based on the mean squared error was considered by Carhart (1997) and Sharpe (1992). They estimated the conditional expectation of a fund return distribution (under the condition that a realization of explanatory factors is observed). Basset and Chen (2001) extended this approach and conducted the style analyses of quantiles of the return distribution. This extension is based on the quantile regression suggested by Koenker and Bassett (1978). The Case Study (2014), “Style Classification with Quantile Regression” implemented this approach and applied quantile regression to the return-based style classification of a mutual fund.
For the numerical implementation of CVaR linear regression, we used the Portfolio Safeguard (2018) package. Portfolio Safeguard (PSG) can solve nonlinear and mixed-integer nonlinear optimization problems. A special feature of PSG is that it includes precoded nonlinear functions: CVaR2 error (cvar2_err) and CVaR2 deviation (cvar2_dev) from the CVaR quadrangle, Rockafellar error (ro_err) from the mixed-quantile quadrangle, and CVaR deviation (cvar_dev) from the quantile quadrangle.
We implemented the following equivalent variants of CVaR regression:
  • Minimization of the CVaR2 error (PSG function cvar2_err).
  • Two-step procedure with CVaR2 deviation (PSG function cvar2_dev).
  • Minimization of the Rockafellar error (PSG function ro_err) with the Set 1 of parameters.
  • The two-step procedure using mixed CVaR deviation with the Set 1 and Set 2 of parameters. This is calculated as a weighted sum of CVaR deviations (PSG function cvar_dev) from the quantile quadrangle.
PSG automatically converts the analytic problem formulations to the mathematical programming codes and solves them. We included in Appendix A convex and linear programming problems for the minimization of the Rockafellar error with the Set 1 of parameters. These formulations are provided for verification purposes. They can be implemented with standard commercial software. For instance, the linear programming formulation can be implemented with the Gurobi optimization package. If Gurobi is installed in the computer, PSG can use Gurobi code as a subsolver. With the CARGRB solver in PSG, by setting the linearize option to 1, it is possible to solve the linear programming problem with Gurobi. However, this conversion will deteriorate the performance, compared to the default PSG solver VAN. For small problems it will not be noticeable. However, for problems with a large number of scenarios (e.g., with 108 observations), the standard PSG solver VAN will dramatically outperform the Gurobi linear programming implementation. In this case, Gurobi may not even start on a small PC because of a shortage of memory. Nevertheless, if the number of observations is small (e.g., 103) and the number of factors is very large (e.g., 107), it is recommended that the linear programming formulation is used.
We regressed the CVaR of the return distribution of the Fidelity Magellan Fund on the explanatory variables: Russell 1000 Growth Index (RLG), Russell 1000 Value Index (RLV), Russell Value Index (RUJ), and Russell 2000 Growth Index (RUO). The dataset includes 1264 historical daily returns of the Magellan Fund and the indices, which were downloaded from the Yahoo Finance website. The data (design matrix for the regression) is posted on the Case Study (2016) website.
The CVaR regression was done with the confidence levels α = 0.75 and α = 0.9. Calculation results are in Table 1 and Table 2, respectively. Here is the description of the columns of the tables:
  • Optimization Problem #: Optimization Problem number, as denoted in Section 5; also it is the problem number in the case study posted online, see, Case Study (2016).
  • Set #: Set of parameters for the mixed-quantile quadrangle.
  • Objective: Optimal value of the objective function.
  • RLG: coefficient for the Russell Value Index.
  • RLV: coefficient for the Russell 1000 Value Index.
  • RUJ: coefficient for the Russell Value Index.
  • RUO: coefficient for the Russell 2000 Growth Index.
  • Intercept: regression intercept.
  • Solving Time: solver optimization time.
Table 1 and Table 2 show calculation results for the considered equivalent problems. We observe that regression coefficients coincide for all problems in Table 1 and Table 2. This confirms the correctness of theoretical results and the numerical implementation. Also, we want to point out that the regression coefficients are quite similar for α = 0.75 (Table 1) and α = 0.9 (Table 2).
The calculation time in majority of cases was around 0.02–0.04 s, except for the case with the mixed CVaR deviation for Set 1, which took 0.11 s. The PSG calculation times were quite low because the solver “knows” analytical expressions for the functions and can take advantage of this knowledge.

7. Conclusions

The quadrangle risk theory Rockafellar and Uryasev (2013) and the decomposition theorem Rockafellar et al. (2008) provided a framework for building a regression with relevant deviations. Solution of a regression problem is split in two steps: (1) minimization of deviation from the corresponding quadrangle, and (2) determining of intercept by using statistic from this quadrangle. For CVaR regression, Rockafellar et al. (2014) reduced the optimization problem at Step 1 to a high-dimension linear programming problem. We suggested two sets of parameters for the mixed-quantile quadrangle and investigated its relationship with the CVaR quadrangle. The Set 1 of parameters corresponds to CVaR regression in Rockafellar et al. (2014), where the Set 2 is a new set of parameters.
For the Set 1 of parameters, the minimization of error from CVaR Quadrangle was reduced to the minimization of the Rockafellar error from the mixed-quantile quadrangle. For both sets of parameters, the minimization of deviation in CVaR quadrangle is equivalent to the minimization of deviation in mixed-quantile quadrangle.
We presented optimization problem statements for CVaR regression problems using CVaR and mixed-quantile quadrangles. Linear regression problem for estimating CVaR were efficiently implemented in Portfolio Safeguard (2018) with convex and linear programming. We have done a case study for the return-based style classification of a mutual fund with CVaR regression. We regressed the fund return by several indices as explanatory factors. Numerical results validating the theoretical statements are placed to the web (see Case Study (2016)).

Supplementary Materials

Data and codes used in the case study can be downloaded: 1. Case Study (2016): Estimation of CVaR through Explanatory Factors with CVaR (Superquantile) Regression. http://www.ise.ufl.edu/uryasev/research/testproblems/financial_engineering/on-implementation-of-cvar-regression/. 2. Case Study (2014): Style Classification with Quantile Regression. http://www.ise.ufl.edu/uryasev/research/testproblems/financial_engineering/style-classification-with-quantile-regression/.

Author Contributions

Conceptualization, S.U.; Formal analysis, V.K.; Investigation, A.G.; Methodology, S.U.; Software, V.K.; Supervision, S.U.; Writing—original draft, A.G.

Funding

Research of Stan Uryasev was partially funded by the AFOSR grant FA9550-18-1-0391 on Massively Parallel Approaches for Buffered Probability Optimization and Applications.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. CVaR Regression with Rockafellar Error: Convex and Linear Programming

The value of the Rockafellar error with given set of parameters λ k , α k ( α k ( 0 , 1 ) , k = 1 , r ,   k = 1 r λ k = 1 ) for a random value X is a minimum w.r.t. a set of variables B 1 , , B r of a mixture of Koenker–Bassett Error functions with one linear constraint on these variables:
R o c k a f e l l a r _ E r r o r   ( X ) λ 1 , α 1 , , λ r , α r = min B 1 , , B r { k = 1 r λ k α k ( X B k ) | k = 1 r λ k B k = 0 }
where α k ( X B k ) = E [ α k 1 α k [ X B k ] + + [ X B k ] ] is the normalized Koenker–Bassett error.
By using regret from the mixed-quantile quadrangle, we express the Rockafellar error as follows:
R o c k a f e l l a r _ E r r o r   ( X ) λ 1 , α 1 , , λ r , α r = min B 1 , , B r { k = 1 r λ k 1 α k E [ X B k ] + | k = 1 r λ k B k = 0 } E [ X ] .
For the linear regression problem, the random variable X is defined by a set of differences between observed values V i and linear functions C 0 + C T Y i , where Y i is a vector of explanatory factors, i = 1 ,   2 , , ν . Vectors C and Y have m components, C = ( C 1 , ,   C m ) , Y = ( Y 1 , , Y m ) , and C 0 is a scalar. Residuals X i = V i C 0 C T Y i are values (scenarios) of atoms of the random value X . We consider that all atoms have equal probabilities. The estimation of V with factors Y is done by minimizing the error w.r.t. variables C , C 0 . Further we use the Set 1 of parameters. Let us denote:
E [ V ] = 1 v i = 1 v V i ,   E [ Y ] =   1 v i = 1 v Y i .

Appendix A.1. Convex Programming Formulation for CVaR Regression

Minimize the Rockafellar error:
min B 1 , , B r , C 0 , C 1 , ,   C m , { k = 1 r λ k ( 1 α k ) v i = 1 v [ V i C 0 C T Y i B k ] + E [ V ] + C 0 + C T E [ Y ] }
subject to the constraint:
k = 1 r λ k B k = 0 .
This optimization problem has a convex objective and one linear constraint.

Appendix A.2. Linear Programming Formulation for CVaR Regression

Equations (A1) and (A2) are reduced to linear programming with additional variables and constraints:
min A 11 , , A r v B 1 , , B r , C 0 , C 1 , ,   C m , { k = 1 r λ k ( 1 α k ) v i = 1 v   A k i E [ V ] + C 0 + C T E [ Y ] }
subject to constraints:
k = 1 r λ k B k = 0
A k i V i C 0 C T Y i B k ,   k = 1 , , r ,   i = 1 , , v
A k i 0 ,   k = 1 , , r ,   i = 1 , , v
The linear function C 0 * + C * T Y estimates C V a R α ( V ) as a function of explanatory factors Y , where C 0 * and C * are optimal values of variables for Equations (A1) and (A2) or (A3)–(A6).

Appendix B. Codes Implementing Regression Optimization Problems

This appendix contains codes implementing Optimization Problems 1–4 described in Section 5. Codes and solution results are posted at internet link: Case Study (2016). Codes are written in Portfolio Safeguard (PSG) Text, MATLAB, and R environments. Here are codes in the Text environment.
Optimization Problem 1
Code in PSG Text format:
  • minimize
  • cvar2_err(0.75,matrix_s)
The keyword “minimize” indicates that the objective function is minimized. The objective function cvar2_err(0.75,matrix_s) calculates error ¯ α ( Z ( C 0 , C ) ) in CVaR quadrangle with confidence level α = 0.75. The matrix_s contains scenarios of residual of the regression Z ( C 0 , C ) = V C 0 C T Y .
Optimization Problem 2
Code in PSG Text format:
  • minimize
  • cvar2_dev(0.75,matrix_s)
  • value:
  • cvar_risk(0.75,matrix_s)
The code includes two parts. The first part begins with the keyword “minimize” indicating that the objective function is minimized. It implements Step 1 of the Optimization Problem 2, which minimizes deviation from the CVaR quadrangle (for determining the optimal vector C * of regression coefficients without an intercept). The PSG function cvar2_dev(0.75,matrix_s) calculates deviation D ¯ α ( Z 0 ( C ) ) with α = 0.75. The matrix_s contains scenarios of the residual of the regression without an intercept Z 0 ( C ) = V C T Y .
The second part of the code begins with the keyword “value.” This part implements Step 2 of the Optimization Problem 2 for calculating the optimal value of intercept C 0 * . The PSG function cvar_risk(0.75,matrix_s) calculates C V a R α ( Z 0 ( C * ) ) with α = 0.75 at the optimal point C * .
Optimization Problem 3
Code in PSG Text format with parameters from Set 1.
  • minimize
  • ro_err(matrix_s, matrix_coeff)
The keyword “minimize” indicates that the objective function is minimized. The objective function ro_err(matrix_s, matrix_coeff) calculates the Rockafellar error ( Z ( C 0 , C ) ) in the mixed-quantile quadrangle. The matrix_s contains scenarios of the residual of the regression Z ( C 0 , C ) = V C 0 C T Y . The matrix_coeff includes vectors of weights and confidence levels for Set 1 with α = 0.75.
Optimization Problem 4
Code in PSG Text format with parameters from Set 1 for α = 0.75.
  • minimize
  • vector_c*cvar_dev(vector_a, matrix_s)
  • value:
  • cvar_risk(0.75,matrix_s)
The code includes two parts. The first part begins with the keyword “minimize” indicating that the objective function is minimized. It implements Step 1 of the Optimization Problem 4, which minimizes deviation from the mixed-quantile quadrangle (for determining optimal vector C * of regression coefficients without intercept). The inner product vector_c*cvar_dev(vector_a, matrix_s) calculates the mixed CVaR deviation D ( Z 0 ( C ) ) . The function cvar_dev corresponds to the CVaR deviation from the quantile quadrangle. Vector vector_c contains weights for CVaR Deviation mix corresponding to Set 1. Vector vector_a contains confidence levels defined by Set 1. The matrix_s contains scenarios of the residual of the regression Z ( C 0 , C ) = V C 0 C T Y .
The second part of the code begins with the keyword “value.” This part implements Step 2 of the Optimization Problem 4, calculating the optimal value of intercept C 0 * . The PSG function cvar_risk(0.75,matrix_s) calculates C V a R α ( Z 0 ( C * ) ) with α = 0.75 at the optimal point C * .

Appendix C. Proof of the Lemma 1

Proof. 
According to the definition, β v = 1 . However, while proving this lemma, we consider β v a bit smaller than 1, i.e., β v = 1 ε > β v 1 ,   ε > 0 , to avoid division by 0. Then, we will consider limit ε 0 to finish the proof. Note that for α < 1 and for the considered partition, β i 1 < β i for all i = ν α ,   , ν .
First, let us prove that γ i ( β i 1 , β i ) for β i < 1 . Consider three functions of σ ,   δ for σ > 0 ,   δ 0 :
f 1 ( δ ) = δ σ + δ ,   f 2 ( δ ) = ln ( 1 + δ σ ) ,   f 3 ( δ ) = δ σ
When δ = 0 , all functions equal to 0 and have equal derivatives. When δ > 0 , it is true for derivatives that:
f 1 ( δ ) < f 2 ( δ ) < f 3 ( δ )
Hence, when > 0 , similar inequalities are valid for functions:
δ σ + δ < ln ( 1 + δ σ ) < δ σ .
There exist some δ γ such that 0 < δ γ < δ and δ σ + δ γ = ln ( 1 + δ σ ) . If δ = β i β i 1 , σ = 1 β i , γ = 1 σ δ γ , then γ = γ i and β i 1 < γ < β i . Therefore, γ i ( β i 1 , β i ) .
Further, we show how to calculate the integral α 1 ε C V a R β ( X ) d β as a sum of integrals over the partition:
¯ α ( X ) = lim ε 0 1 1 α α 1 ε C V a R β ( X ) d β = lim β v 1 1 1 α i = ν α ν β i 1 β i C V a R β ( X ) d β .
Let us denote C β i = C V a R β i ( X ) and V i = V a R γ i ( X ) . Note that V a R γ ( X ) is a singleton for every γ ( β i 1 , β i ) and it equals V i , because γ i ( β i 1 , β i ) .
Below, we use value V i for the closed interval [ β i 1 , β i ] while calculating the integral over this interval because the value of the integral does not depend on the finite values of V a R γ ( X ) at the boundary points β i 1 , β i .
Using the definition of CVaR (Rockafellar and Uryasev (2002), Proposition 8 CVaR for scenario models) we write:
β i 1 β i C V a R β ( X ) d β = β i 1 β i 1 1 β ( C β i ( 1 β i ) + V i ( β i β ) ) d β = C β i ( 1 β i ) β i 1 β i 1 1 β d β + V i β i 1 β i β i β 1 β d β = C β i ( 1 β i ) ln ( 1 β i 1 1 β i ) + V i [ ( β i β i 1 ) ( 1 β i ) ln ( 1 β i 1 1 β i ) ]
Let us make the transformation of Equation (A7) using the expression for γ i in the Set 1 definition:
β i 1 β i C V a R β ( X ) d β = ( β i β i 1 ) [ V i + 1 β i β i 1 ln ( 1 β i 1 1 β i ) ( C β i V i ) ( 1 β i ) ] = ( β i β i 1 ) [ V i + 1 1 γ i ( C β i V i ) ( 1 β i ) ] = ( β i β i 1 ) 1 1 γ i [ C β i ( 1 β i ) + V i ( β i γ i ) ] = ( β i β i 1 ) C V a R γ i ( X )
The last equality is valid because γ i ( β i 1 , β i ) .
Taking into account that lim ε 0   γ ν = 1 , C V a R 1 ε ( X ) = C V a R 1 ( X ) , and V a R 1 ε ( X ) = V a R 1 ( X ) , we obtain:
¯ α ( X ) = lim ε 0 i = ν α ν β i β i 1 1 α C V a R γ i ( X ) = i = ν α ν p i C V a R γ i ( X )
Deviation is calculated as follows:
D ¯ α ( X ) = i = ν α ν p i C V a R γ i ( X ) E [ X ]
By the definition of CVaR (Rockafellar and Uryasev (2002)):
C V a R α ( X ) = i = ν α ν p i V a R γ i ( X ) = S ¯ α ( X )
Lemma 1 is proved. □

Appendix D. Proof of Lemma 2

Proof. 
Similar to proof of Lemma 1, we consider β v as a bit smaller than 1, i.e., β v = 1 ε > β v 1 ,   ε > 0 , to avoid division by 0. Then, we consider limit ε 0 to finish the proof. We denote C β i = C V a R β i ( X ) and V i = V a R γ ( X ) for any γ ( β i 1 , β i ) because V a R γ ( X ) does not change on this interval.
Additionally we denote δ i = β i β i 1 , i = ν α ,   , ν , and σ i = 1 β i , i = ν α 1 ,   , ν . Note that all δ i > 0 , all σ i > 0 , and δ i = δ j = δ for all i , j < ν .
¯ α ( X ) = lim ε 0 1 1 α α 1 ε C V a R β ( X ) d β = lim β v 1 1 1 α i = ν α ν β i 1 β i C V a R β ( X ) d β .
Equation (A7) in Lemma 1 is valid for any interval [ θ , β i ] such that β i 1   θ β i , therefore:
θ β i C V a R β ( X ) d β = C β i ( 1 β i ) ln ( 1 θ 1 β i ) + V i [ ( β i θ ) ( 1 β i ) ln ( 1 θ 1 β i ) ]
Let us express V i from C V a R β i 1 ( X ) using the definition of CVaR from Rockafellar and Uryasev (2002) (Equation (25)), then insert this expression into Equation (A9):
C β i 1 = 1 1 β i 1 ( C β i ( 1 β i ) + V i ( β i β i 1 ) ) = 1 σ i 1 ( C β i σ i + V i δ i ) , V i = 1 δ i ( C β i 1 σ i 1 C β i σ i )
By substituting V i into Equation (A9), we obtain:
θ β i C V a R β ( X ) d β = C β i σ i ln ( 1 θ σ i ) + 1 δ i ( C β i 1 σ i 1 C β i σ i ) [ ( β i θ ) σ i ln ( 1 θ σ i ) ]
The last equation contains two CVaRs, C β i and C β i 1 . Let us express the coefficients of these CVaRs.
C β i in Equation (A10) has the following coefficient:
q 1 i = σ i [ ln ( 1 θ σ i ) ( 1 + σ i δ i ) β i θ δ i ]
and C β i 1 has the coefficient:
q 2 i 1 = σ i 1 δ i [ ( β i θ ) σ i ln ( 1 θ σ i ) ]
Coefficient q 2 i 1 0 because the value in square brackets is the same as in Equations (A7) and (A9), and is a result of the integration of a non-negative function over the interval [ θ , β i ] . The coefficient q 1 i 0 because it is an integral of the non-negative function σ i δ i [ δ i β i + β 1 β ] over the same interval.
When summing up integrals to obtain 1 1 α α 1 ε C V a R β ( X ) d β , every C β i , aside from C β ν α 1 and C β ν , enters the sum two times with coefficients depending on ν , α , β i 1 ,   β i , and β i + 1 . Once in Equation (A11) for i , and the second time in Equation (A12) for i + 1 . All coefficients are non-negative.
Let us explain this in more detail.
If i is such that ν α < i < ν , then β ν α < β i < β ν and θ = β i 1 for i in Equation (A11) and θ = β i for i + 1 in Equation (A12). Then, the coefficient for C β i in Equation (A8) equals:
q i = 1 1 α ( q 1 i + q 2 i ) = 1 1 α ( σ i [ ln ( σ i 1 σ i ) ( σ i 1 δ i ) 1 ] + σ i δ i + 1 [ δ i + 1 σ i + 1 ln ( σ i σ i + 1 ) ] ) = σ i 1 α [ σ i 1 δ i ln ( σ i 1 σ i ) σ i + 1 δ i + 1 ln ( σ i σ i + 1 ) ]
If i = ν α < ν , then θ = α in Equation (A11) and θ = β i in Equation (A12). Then, the coefficient for C β i in Equation (A8) equals:
q i = 1 1 α ( q 1 i + q 2 i ) = 1 1 α ( σ i [ ln ( 1 α σ i ) ( σ i 1 δ i ) δ α δ i ] + σ i δ i + 1 [ δ i + 1 σ i + 1 ln ( σ i σ i + 1 ) ] ) = σ i 1 α [ 1 δ α δ i + σ i 1 δ i ln ( 1 α σ i ) σ i + 1 δ i + 1 ln ( σ i σ i + 1 ) ]
If i = ν α 1 , then C β i enters in the sum only in Equation (A12) with θ = α . Then:
q i = 1 1 α q 2 i = σ i 1 α [ δ α δ i + 1 σ i + 1 δ i + 1 ln ( 1 α σ i + 1 ) ]
If i = ν > ν α , then C β i enters in the sum only in Equation (A11) with θ = β i 1 . Then:
q i = 1 1 α q 1 i = σ i 1 α [ ln ( σ i 1 σ i ) ( σ i 1 δ i ) 1 ]
Also, in the case when i = ν = ν α , C β i enters in the sum only in Equation (A11) with θ = α . Then:
q i = 1 1 α q 1 i = σ i 1 α [ ln ( 1 α σ i ) ( 1 + σ i δ i ) δ α δ i ]
Then, the risk in the CVaR quadrangle equals:
¯ α ( X ) = lim ε 0 i = ν α 1 ν q i C V a R β i ( X ) .
Deviation is calculated as follows: D ¯ α ( X ) =   ¯ α ( X ) E [ X ] .
Taking into account that all β i are fixed for i < ν , the limit operation affects only coefficients for i = ν . Namely, lim ε 0   β ν = 1 , lim ε 0   σ ν = 0 , lim ε 0   σ ν ln ( σ ν ) = 0 , lim ε 0   δ ν = δ . Then, the limit values of coefficients q i are equal to:
  • q ν = 0 in both cases, when ν > ν α and ν = ν α .
  • q i = 1 1 α × σ i δ [ σ i 1 ln ( σ i 1 σ i ) ] for i such that ν α < i = ν 1 .
  • q i = 1 1 α × σ i δ [ σ i 1 ln ( σ i 1 σ i ) + σ i + 1 ln ( σ i + 1 σ i ) ] for i such that ν α < i < ν 1 .
  • q i = 1 1 α × σ i δ [ δ δ α + σ i 1 ln ( 1 α σ i ) + σ i + 1 ln ( σ i + 1 σ i ) ] for i = ν α < ν 1 .
  • q i = 1 1 α × σ i δ [ δ α + σ i + 1 ln ( σ i + 1 1 α ) ] for i = ν α 1 < ν 1 .
If the number of atoms used in calculating ¯ α ( X ) is 3 or 2, then the coefficients have values:
  • q i = 1 1 α × σ i δ [ δ δ α + σ i 1 ln ( 1 α σ i ) ] for i = ν α = ν 1 .
  • q i = 1 1 α × σ i δ [ δ α ] = 1 for i = ν α 1 = ν 1 .
It can be shown that i = ν α 1 ν q i = 1 by sequentially summing up coefficients.
By recalling that σ i = δ ( v i ) , we can rewrite the equations for coefficients q i using δ and j = v i .
Lemma 2 is proved. □

Appendix E. Proof of the Lemma 3

Proof. 
Distribution of random value X defines the partition of the interval [ 0 , 1 ] : δ = 1 / ν , β i = i δ , for i = 0 ,   1 , , ν .
CDF F X ( x ) = p r o b { X x } is a non-decreasing and right-continuous function and it is constant for every right-open interval [ β i 1 , β i ) , i = 1 , , ν . Therefore, from the definitions of V a R γ ( X )   and   V a R γ + ( X ) , we have:
V a R γ ( X ) = V a R γ + ( X ) = V a R γ ( X ) = c o n s t   for   γ ( β i 1 , β i ) .
V a R β i + ( X ) = V a R γ ( X )   for   γ ( β i , β i + 1 )   and   i = 0 ,   , ν 1 .
V a R β i ( X ) = V a R γ ( X )   for   γ ( β i 1 , β i )   and   i = 1 ,   , ν .
For i = 0 we have β i = 0 and V a R 0 ( X ) = .
Thus, according to the definition of statistic for the mixed-quantile quadrangle, for Set 2:
S α I I ( X ) = i = ν α 1 ν 1 q i V a R β i ( X )
where V a R β i ( X ) are intervals (that may have zero lengths for some i ), and q i and ν α are defined above for the Set 2. Therefore, S α I I ( X ) is also an interval and it has a non-zero length if not all V a R β i ( X ) ,   i = ν α 1 , , ν 1 are equal.
If ν α 1 = 0 , then S α I I ( X ) is a left-open interval with the lower bound .
Let ν α 1 > 0 . For the Set 1, we defined confidence levels as internal points in intervals γ i ( β i 1 , β i ) . Using these definitions of γ i , we can express S α I I ( X ) as the interval:
S α I I ( X ) = [ i = ν α 1 ν 1 q i V a R γ i ( X ) ,   i = ν α 1 ν 1 q i V a R γ i + 1 ( X )   ]
To simplify notations, let us denote: V i = V a R γ i ( X ) ; C γ = C V a R γ ( X ) ; and L , U are bounds for S α I I ( X ) such that S α I I ( X ) = [ L , U ] .
Statistic S ¯ α ( X ) of the CVaR quadrangle is defined in Lemma 1 (see Equation (1)) as S ¯ α ( X ) = i = ν α ν p i V i .
Therefore, we wish to prove that:
L = i = ν α 1 ν 1 q i V i   S ¯ α ( X ) i = ν α 1 ν 1 q i V i + 1 = U
Let us prove the right inequality for the upper bound.
According to Lemmas 1 and 2, it is valid that:
¯ α ( X ) = i = ν α ν p i C γ i = i = ν α 1 ν 1 q i C β i
Let d i = q i 1 p i . Then, from the last equality:
i = ν α ν p i C γ i = i = ν α ν q i 1 C β i 1 = i = ν α ν p i C β i 1 + i = ν α ν d i C β i 1
i = ν α ν d i C β i 1 = i = ν α ν p i ( C γ i C β i 1 ) 0
The last inequality is valid because p i > 0 , γ i > β i 1 , and C γ i C β i 1 .
Because C β i 1 may have arbitrary but ordered values ( C β i 1 C β i   due   to   ordered   β i ), i = ν α , ,   ν , the last inequality is valid for any ordered sequence of values z i 1 z i . Therefore, it is valid for VaRs V i that:
i = ν α ν d i V i = i = ν α ν ( q i 1 p i ) V i 0
This leads to the inequality for the upper bound:
U = i = ν α ν q i 1 V i i = ν α ν p i V i = S ¯ α ( X ) .
Let us prove the inequality for the lower bound L :
L = i = ν α ν q i 1 V i 1   = i = ν α ν ( p i + d i ) ( V i Δ V i )   = S ¯ α ( X ) i = ν α ν p i Δ V i   + i = ν α ν d i V i 1
where Δ V i = V i V i 1 .
Let us calculate the upper estimate of L . Let us set V ν α 1 = V ν α , therefore Δ V ν α = 0 . This increases the right hand side of Equation (A13). By recalling that p i = δ / ( 1 α ) = 1 ν ( 1 α ) ,   i = ν α + 1 , ,   ν , for the Set 1, we have:
i = ν α ν p i Δ V i   = δ ( V ν V ν α ) / ( 1 α )
Because i = ν α ν p i = i = ν α ν q i 1 , then i = ν α ν d i = 0 and i = ν α ν d i V i 1 = i = ν α ν d i ( V i 1 D ) for any D .
Because i = ν α ν d i z i 0 for any increasing ordered z i , and we can set z ν α = 0 ,   z i = 1 ,   i = ν α + 1 , ,   ν , then:
i = ν α ν d i z i = i = ν α + 1 ν d i = d ν α = p ν α q ν α 1 0
Taking into account that p ν α , q ν α 1 0 and p ν α δ 1 α , we have i = ν α + 1 ν d i δ 1 α . Then:
i = ν α ν d i V i 1 = i = ν α ν d i ( V i 1 V ν α 1 ) i = ν α + 1 ν d i ( V ν 1 V ν α 1 ) ( V ν 1 V ν α 1 ) δ 1 α
Let us return to estimation L by taking into account that V ν α 1 = V ν α :
L S ¯ α ( X ) i = ν α ν p i Δ V i   + i = ν α ν d i V i 1 S ¯ α ( X ) δ ( V ν V ν α ) 1 α + δ ( V ν 1 V ν α 1 ) 1 α = S ¯ α ( X ) δ ( V ν V ν 1 ) 1 α S ¯ α ( X )
Therefore, we have proved that L S ¯ α ( X ) U .
Lemma 3 is proved. □

References

  1. Acerbi, Carlo, and Dirk Tasche. 2002. On the Coherence of Expected Shortfall. Journal of Banking and Finance 26: 1487–503. [Google Scholar] [CrossRef]
  2. Adrian, Tobias, and Markus K. Brunnermeier. 2016. CoVaR. American Economic Review 106: 1705–41. [Google Scholar] [CrossRef]
  3. Basset, Gilbert W., and Hsiu-Lang Chen. 2001. Portfolio Style: Return-based Attribution Using Quantile Regression. Empirical Economics 26: 293–305. [Google Scholar] [CrossRef]
  4. Beraldi, Patrizia, Antonio Violi, Massimiliano Ferrara, Claudio Ciancio, and Bruno Antonio Pansera. 2019. Dealing with complex transaction costs in portfolio management. Annals of Operations Research, 1–16. [Google Scholar] [CrossRef]
  5. Carhart, Mark M. 1997. On Persistence in Mutual Fund Performance. Journal of Finance 52: 57–82. [Google Scholar] [CrossRef]
  6. Case Study. 2014. Style Classification with Quantile Regression. Available online: http://www.ise.ufl.edu/uryasev/research/testproblems/financial_engineering/style-classification-with-quantile-regression/ (accessed on 24 June 2019).
  7. Case Study. 2016. Estimation of CVaR through Explanatory Factors with CVaR (Superquantile) Regression. Available online: http://www.ise.ufl.edu/uryasev/research/testproblems/financial_engineering/on-implementation-of-cvar-regression/ (accessed on 24 June 2019).
  8. Fissler, Tobias, and Johanna F. Ziegel. 2015. Higher order elicitability and Osband’s principle. arXiv. [Google Scholar] [CrossRef]
  9. Huang, Wei-Qiang, and Stan Uryasev. 2018. The CoCVaR Approach: Systemic Risk Contribution Measurement. Journal of Risk 20: 75–93. [Google Scholar] [CrossRef]
  10. Koenker, R. 2005. Quantile Regression. Cambridge: Cambridge University Press. [Google Scholar]
  11. Koenker, Roger, and Gilbert Bassett. 1978. Regression Quantiles. Econometrica 46: 33–50. [Google Scholar] [CrossRef]
  12. Portfolio Safeguard. 2018. American Optimal Decisions, USA. Available online: http://www.aorda.com (accessed on 24 June 2019).
  13. Rockafellar, R. Tyrrell, and Johannes O. Royset. 2018. Superquantile/CVaR Risk Measures: Second-order Theory. Annals of Operations Research 262: 3–29. [Google Scholar] [CrossRef]
  14. Rockafellar, R. Tyrrell, and Stan Uryasev. 2000. Optimization of Conditional Value-At-Risk. Journal of Risk 2: 21–41. [Google Scholar] [CrossRef]
  15. Rockafellar, R. Tyrrell, and Stan Uryasev. 2002. Conditional Value-at-Risk for General Loss Distributions. Journal of Banking and Finance 26: 1443–71. [Google Scholar] [CrossRef]
  16. Rockafellar, R. Tyrrell, and Stan Uryasev. 2013. The Fundamental Risk Quadrangle in Risk Management, Optimization and Statistical Estimation. Surveys in Operations Research and Management Science 18: 33–53. [Google Scholar] [CrossRef]
  17. Rockafellar, R. Tyrrell, Stan Uryasev, and Michael Zabarankin. 2008. Risk Tuning with Generalized Linear Regression. Mathematics of Operations Research 33: 712–29. [Google Scholar] [CrossRef]
  18. Rockafellar, R. Terry, Johannes O. Royset, and Sofia I. Miranda. 2014. Superquantile Regression with Applications to Buffered Reliability, Uncertainty Quantification and Conditional Value-at-Risk. European Journal Operations Research 234: 140–54. [Google Scholar] [CrossRef]
  19. Sharpe, William F. 1992. Asset Allocation: Management Style and Performance Measurement. Journal of Portfolio Management (Winter) 18: 7–19. [Google Scholar] [CrossRef]
  20. Ziegel, Johanna F. 2014. Coherence and elicitability. arXiv. [Google Scholar] [CrossRef]
Figure 1. Five equally probable atoms. VaR and CVaR with Set 1 of parameters for α = 0.5 .
Figure 1. Five equally probable atoms. VaR and CVaR with Set 1 of parameters for α = 0.5 .
Jrfm 12 00107 g001
Figure 2. Five equally probable atoms. Risk in the CVaR quadrangle and mixed-quantile quadrangle with Set 1 of parameters for α = 0.5 .
Figure 2. Five equally probable atoms. Risk in the CVaR quadrangle and mixed-quantile quadrangle with Set 1 of parameters for α = 0.5 .
Jrfm 12 00107 g002
Table 1. Optimization outputs: estimating CVaR with the linear regression, α = 0.75.
Table 1. Optimization outputs: estimating CVaR with the linear regression, α = 0.75.
Optimization Problem #Set #ObjectiveRLGRLVRUJRUOInterceptSolving Time (s)
1N/A0.012480.4860.581−0.0753−6.22 × 10−36.98 × 10−30.02
2N/A0.012480.4860.582−0.0753−6.22 × 10−36.98 × 10−30.02
3Set 10.012470.4860.582−0.0753−6.22 × 10−36.96 × 10−30.03
4Set 10.012510.4860.582−0.0752−6.22 × 10−36.98 × 10−30.11
4Set 20.012480.4860.582−0.0753−6.23 × 10−36.98 × 10−30.03
Table 2. Optimization outputs: estimating CVaR with the linear regression, α = 0.9.
Table 2. Optimization outputs: estimating CVaR with the linear regression, α = 0.9.
Optimization Problem #Set #ObjectiveRLGRLVRUJRUOInterceptSolving Time (s)
1N/A0.0166560.4720.606−0.078−7.052 × 10−31.05 × 10−20.02
2N/A0.0166560.4720.606−0.078−7.052 × 10−31.05 × 10−20.02
3Set 10.0166560.4720.606−0.078−7.052 × 10−31.05 × 10−20.02
4Set 10.0166560.4720.606−0.078−7.052 × 10−31.05 × 10−20.11
4Set 20.0166560.4720.606−0.078−7.052 × 10−31.05 × 10−20.04

Share and Cite

MDPI and ACS Style

Golodnikov, A.; Kuzmenko, V.; Uryasev, S. CVaR Regression Based on the Relation between CVaR and Mixed-Quantile Quadrangles. J. Risk Financial Manag. 2019, 12, 107. https://doi.org/10.3390/jrfm12030107

AMA Style

Golodnikov A, Kuzmenko V, Uryasev S. CVaR Regression Based on the Relation between CVaR and Mixed-Quantile Quadrangles. Journal of Risk and Financial Management. 2019; 12(3):107. https://doi.org/10.3390/jrfm12030107

Chicago/Turabian Style

Golodnikov, Alex, Viktor Kuzmenko, and Stan Uryasev. 2019. "CVaR Regression Based on the Relation between CVaR and Mixed-Quantile Quadrangles" Journal of Risk and Financial Management 12, no. 3: 107. https://doi.org/10.3390/jrfm12030107

APA Style

Golodnikov, A., Kuzmenko, V., & Uryasev, S. (2019). CVaR Regression Based on the Relation between CVaR and Mixed-Quantile Quadrangles. Journal of Risk and Financial Management, 12(3), 107. https://doi.org/10.3390/jrfm12030107

Article Metrics

Back to TopTop