Estimation of Conditional Weighted Expected Shortfall under Adjusted Extreme Quantile Autoregression ()
1. Introduction
Different regulators expect institutions under their jurisdiction to comply with the set guidelines on how to ensure capital adequacy. A major threat to capital adequacy is reserves for fluctuations in investment values due to market risk. The amount of this reserve is determined using risk measures that quantify the market risk from the downside distribution of returns. One of the basic risk measures is volatility commonly referred to as standard deviation. Popular measures of risk that have found immense use in financial risk management include Value at Risk (VaR) and Expected Shortfall (ES).
VaR is the maximum potential loss a portfolio can suffer at a certain confidence interval, say 99%, in a particular number of days called the holding period. In addition to comparing Value at Risk and the coherent Expected Shortfall proposed in [1], it was also noted in [2] that VaR improved the “what if?” sensitivity analysis in Greeks by expressing risk in terms of currency as one number with some probability attached. This eases risk reporting and aids in comparing different portfolios. However, VaR suffers from two major shortcomings that cripple its application in risk measurement. These are: VaR is not sub-additive with respect to subportfolios as well as risk variables and the fact that it underestimates true risk by failing to appreciate the severity of losses beyond the confidence threshold, [3]. Moreover, [4] showed that it is possible to construct two portfolios with different levels of tail risk but with the same VaR.
To remedy the above shortcomings the Basel Committee on Banking Supervision during the fundamental review of the trading book in the year 2013 adopted 97.5% Expected shortfall with a horizon of one day in place of 99% VaR in quantifying market risk for banks under the Basel III framework, [5]. This decision was followed by a lot of criticism based on earlier findings in [6] that ES is not elicitable and hence impossible to backtest. The fears were quelled by [7] who derived backtest procedures for ES that do not require elicitability. In fact elicitability makes backtesting easier but is not a necessary condition for backtesting procedures, [4]. However, following findings in [8] that ES is less robust compared to VaR, any misspecification of the loss/returns distribution greatly affects the risk estimates from ES. This can be attributed to the fact that ES averages quantiles in the extreme left tail of the returns distribution (right tail of the loss distribution) with equal weights meaning that the very extreme quantiles have the same effect on the final estimate as quantiles close to VaR.
The use of equal weights was remedied in [9] using a uniformly weighted sum of a systematic sample of quantiles below the value at risk as an estimator of the expected shortfall. The number of sample quantiles used in the estimation was chosen subjectively depending on the considered sample size of returns distribution. This approach was generalized in [10] through the use of nonuniform weights in the summation of the sample quantiles below the value at risk; where the left tail of the returns distribution was considered. The weights were defined using a continuously differentiable function that assigned more weight to quantiles near
thus improving the asymptotic efficiency of the ES estimator.
This paper, presents derivation of an estimator for Weighted Expected Shortfall (WES) based on extreme quantile regression in [11] to remedy the subjectivity in the WICQF1 estimator in [10]. To achieve this, dynamic weights are introduced in the structural form of expected shortfall. The weights depend on the distance between the extreme quantiles and the value at risk. This will allow us to use all the quantiles above VaR thus quashing the subjectivity in WICQF estimator.
2. Methodology
Consider a real valued financial time series
on a complete probability space
. Assume
is
-measurable where
is an increasing sequence of
-algebras representing information available up to time t. In particular, let
be the value of a portfolio at trading time t so that the return,
on the portfolio at time t is given by
(1)
For convenience, let the corresponding loss return be given by
(2)
The VaR at
of the return on this portfolio is given by the quantile at
of the loss distribution
(3)
where
. The corresponding ES of the return on this portfolio is the expected value of the losses beyond VaR which is given by
(4)
where
is the inverse of the loss distribution function
For a discrete ordered sample
the corresponding VaR and ES estimates are respectively given by
(5)
(6)
where appropriate approximations are made when
is not an integer.
Let us now define a conditional quantile autoregressive model on
of the form
(7)
where
is the central conditional
-quantile of
and
are heteroscedastic errors with zero
-quantile. Suppose we define
where
is the central conditional scale of
and
are assumed to be iid2 innovations with a common distribution F(.). To capture the ARCH effects in
we let
be a function of lagged values of
, that is
. This modification improves the QAR-QAR process in [11]
According to [11] the adjusted extreme conditional quantile of
is given by
(8)
defined such that if estimates of the quantiles of
at
and
are
and
respectively then we can estimate the function by
(9)
3. Estimation of Common Risk Measures
Using Equation (9) we obtain an estimator for the one step ahead VaR forecast as
(10)
where
and
are the corresponding one step
-quantile and scale estimators respectively from the linear conditional quantile process and;
and
are obtained by inverting the overall distribution of the iid errors. Note that the overall distribution of the iid errors is obtained by splicing GPD with empirical bulk distribution at the threshold as outlined in [11] to get
(11)
for any
.
and
are the estimated parameters of the Generalised Pareto Distribution (GPD) adopted in formulating the estimated overall distribution of the iid errors. m is the number of exceedances above a chosen threshold u from a sample of size N. Inverting
we get
(12)
Using Equation (4) the estimator of the expected shortfall from standardized errors is given by
(13)
where
is as defined in Equation (12). Similarly, using Equation (4) we have
(14)
where
is as defined in Equation (13).
Weighted Expected Shortfall
Suppose
, then we define the Weighted Expected Shortfall for a continuous random variable X as:
(15)
where
, are dynamic weights (risk spectrum) that vary from one quantile to another in the integrand and
is fixed.
is the quantile at level
.
Assumption 1 The dynamics of the weight function (risk spectrum) are such that for any risks X, Y;
and
(16)
where
is the quantile at level
for loss distribution of X,
is the quantile at level
for loss distribution of Y and
is the quantile at level
for loss distribution of
To reduce cumbersomeness in notations, we will drop the superscript on the quantile function,
and reintroduce it when referring to another loss distribution different from X.
Proposition 1 (Weights)
defined as
(17)
is an admissible risk spectrum.
Proof. We need to show that
satisfies the three conditions in [12] for an admissible risk spectrum in the definition of a spectral risk measure.
1) Expanding Equation (17) we obtain
(18)
as a result of non-negativity of the exponential function. Therefore
.
2) The first derivative of the weight function with respect to
is given by
(19)
for
which confirms that
is monotonically decreasing in the interval.
3) Observe that
(20)
Replacing
and
with their respective estimates in Equation (15) we deduce the corresponding WES estimate as:
(21)
where
is the estimated weighted expected shortfall of the standardized errors.
Theorem 1 (Consistency of WES estimator)
Proof. Observe that
is a continuous function of
. Since
by theorem 2 in [11] then by continuous mapping theorem in [13] we have the result.
For a discrete ordered sample
,
is given by
(22)
where
is as defined in Equation (6) and
(23)
Note that as indicated earlier if
is not an integer then appropriate approximations are made.
4. Properties of the Weighted Expected Shortfall
In this section we look at the fundamental properties, lemmas and theorems of the loss random variable as well as the weighted expected shortfall.
Lemma 1 For any random variable X, there exists a
random variable
such that
Proof. This has been proved through a distribution transform in [14] [Lemma 2.1].
Lemma 2 For
and X integrable we have
(24)
where
means
and
Proof. By definition
(25)
where
(26)
Lemma 3 For
and X integrable we have
(27)
Proof. From Lemma 2 we have
(28)
Definition 1 Let X be a bounded, integrable random variable. For
and
we define the generalized indicator function
(29)
Lemma 4 Let X be a bounded, integrable random variable. For
and
, the following holds
1)
2)
3)
4)
Proof. For proof of 1) and 2) see Lemma 3.5 in [14].
Prove of 3) follows trivially from admissibility of
in Proposition 1.
To prove 4) we note that if
, then
hence by lemma 3 we have
(30)
If
then
(31)
Theorem 2 (Subadditivity of WES) Given that
satisfies assumption 1 then
is subadditive. That is
(32)
Proof. By lemma 4 4) we have,
(33)
Let
We show that
. Note that
by lemma 4 1). When
then
since the highest value
can attain is 1 at which point
. When
, then
. Finally, if
then
since
. Therefore
By Lemma 4 3)
and so we obtain
Similarly, it can be shown that
Therefore
Remark Note that both assumption 1 and theorem 2 can be generalized for any number of losses.
Theorem 3 (Monotonicity of WES) For
is monotonic. That is if
always then
Proof. The result follows from Proposition 3.5(i) in [15] and monotonicity of both the quantile function and the Lebesgue integral.
Theorem 4 (Coherence of WES) Given that
satisfies assumption 1 then
is coherent.
Proof. We show that
satisfies the four axioms of definition 2.7 in [4]. We note that the proves of axiom 1) and 3) are trivial. Axiom 2) and 4) follows from theorems 2 and 3 respectively.
5. Application to Risk Measurement for NSE 20 Share Index
We compare risk estimates from VaR, ES and WES using NSE3 20 Share index data from January 2008 to March 2021. Returns and loss returns are calculated from the price data using Equations (1) and (2) respectively.
Table 1 reports summary statistics of the data which shows that the data is heavy tailed and skewed to the left. Moreso, the reported p-value from the ADF test implies that the data is stationary at 5% significance level.
From Figure 1 we observe some level of volatility clustering which is common in most financial data sets.
Q1, Lower quartile Q3, Upper quartile sd, standard deviation.
Figure 1. NSE 20 share, loss return time series.
The ACF and PACF plots of the data in Figure 2 suggests autocorrelation in the series of up to order two which informs the number of lags in our QAR process. Therefore using extreme quantile autoregression we obtain the following model for the loss return
(34)
where
follows the extreme value distribution given by Equation (11).
Based on p-values in Table 2 all the the estimated parameters of the quantile process are significant at 5% significance level except the constant of the central quantile process.
The ACF plot in Figure 3 confirms independence of the resulting standardized errors. Hence using the threshold
, to ensure that 10% of the errors are classified as extreme, we obtain the following estimates of the shape and scale parameters from the GPD fit.
As outlined in Table 3, the p-values of all the parameter estimates are less than 5% significance level implying that the distribution fits the data well. Figure 4 shows a plot of the corresponding risk estimates at 97.5% confidence level.
From Figure 4 we observe that risk estimates from ES are slightly higher than those from WES confirming that WES is less conservative compared to
Figure 2. Autocorrelation and Partial Autocorrelation Plot of the loss return series.
Figure 3. Autocorrelation Plot of the standardized errors.
Figure 4. Data plot superimposed with risk estimates from considered Risk measures.
Table 2. Parameter estimates of the central quantile and the scale.
Table 3. Parameter estimates of the Central quantile and the scale.
ES. However, risk estimates from both ES and WES are higher than those from VaR because VaR failed to appreciate severity of losses beyond the 97.5% confidence threshold.
6. Conclusion and Recommendations
We have improved the adjusted extreme conditional quantile estimator in [11] and used it to obtain the one-step-ahead risk estimators for Var, ES and WES. The three estimators were then used to quantify risk in the NSE 20 share index portfolio. Consistency and coherence of the proposed Weighted expected shortfall were also proved. It was observed that the WES is less conservative compared to ES.
Data Availability
The authors used NSE 20 share index data obtained from the Nairobi Stock Exchange (NSE).
Acknowledgements
The authors thank the Pan-African University of Basic Sciences, Technology and Innovation (PAUSTI) for funding this research.
NOTES
1WICQF, Weighted Conditional Integrated Quantile Function.
2Independent and identically distributed.
3Nairobi Stock Exchange.