Next Article in Journal
Tax Shields, the Weighted Average Cost of Capital, and the Appropriate Discount Rate for a Project with a Finite Useful Life
Next Article in Special Issue
Does Fiscal Consolidation Affect Non-Performing Loans? Global Evidence from Heavily Indebted Countries (HICs)
Previous Article in Journal
Measuring the Performance of Private Pension Companies in Türkiye by Gray Relational Analysis Method
Previous Article in Special Issue
Properties of VaR and CVaR Risk Measures in High-Frequency Domain: Long–Short Asymmetry and Significance of the Power-Law Tail
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial Multivariate GARCH Models and Financial Spillovers

by
Rosella Giacometti
1,
Gabriele Torri
1,2,*,
Kamonchai Rujirarangsan
3 and
Michela Cameletti
3
1
Department of Management, University of Bergamo, Via dei Caniana 2, 24127 Bergamo, Italy
2
Department of Finance, VŠB-TU Ostrava, Sokolská třída 33, 701 21 Ostrava, Czech Republic
3
Department of Economics, University of Bergamo, Via dei Caniana 2, 24127 Bergamo, Italy
*
Author to whom correspondence should be addressed.
J. Risk Financial Manag. 2023, 16(9), 397; https://doi.org/10.3390/jrfm16090397
Submission received: 31 July 2023 / Revised: 31 August 2023 / Accepted: 4 September 2023 / Published: 6 September 2023
(This article belongs to the Special Issue Featured Papers in Mathematics and Finance)

Abstract

:
We estimate the risk spillover among European banks from equity log-return data via Conditional Value at Risk (CoVaR). The joint dynamic of returns is modeled with a spatial DCC-GARCH which allows the conditional variance of log-returns of each bank to depend on past volatility shocks to other banks and their past squared returns in a parsimonious way. The backtesting of the resulting risk measures provides evidence that (i) the multivariate GARCH model with Student’s t distribution is more accurate than both the standard multivariate Gaussian model and the Filtered Historical Simulation (FHS), and (ii) the introduction of a spatial component improves the assessment of risk profiles and the market risk spillovers.
JEL Classification:
C31; C32; G21; G32

1. Introduction

The interconnectedness of risk between banks is an increasingly hot topic. In the last decades, several countries have simultaneously faced severe economic conditions with spillover effects of risk across the EU. Due to the direct and indirect links among the banks, the stand-alone measurement of a Value-at-Risk (VaR) of each bank cannot provide a comprehensive representation of the risk (Adrian and Brunnermeier 2014; Billio et al. 2012; Rahman 2014).
Recently, multivariate GARCH models have been playing a crucial role in estimating risk interconnectedness. The constant conditional correlation GARCH model (CCC-GARCH) proposed by Bollerslev (1990) is computationally less complex than other multivariate models (see, among others, Bollerslev et al. 1988; Diebold and Nerlove 1989; Engle et al. 1990). However, it does not capture the dynamic interactions between the volatilities. The BEKK model proposed by Engle and Kroner (1995) (the acronym stands for Baba, Engle, Kraft, and Kroner) allows for the dependence of conditional variances and covariance of one variable on the lagged values of another variable, so that spillovers in variances can be modeled. However, it is highly computationally intensive due to the large number of parameters. A more parsimonious model is the Dynamic Conditional Correlation GARCH (DCC-GARCH), introduced by Engle (2002), that introduces an autoregressive process for the conditional correlation matrix, allowing us to model its dynamics with only two parameters in addition to the ones of the CCC-GARCH model.
An approach for modeling explicitly the volatility interactions is the introduction of a spatial component that accounts for the effect of direct bilateral exposures, closeness, or similarities between different financial institutions. Borovkova and Lopuhaa (2012) adopted a spatial GARCH approach to handle the spillover effects where the spatial weights are obtained from the GDP data and alternatively from the market capitalization of the US and European countries’ stock market and embedded in the Extended CCC-GARCH model (E-CCC), see Jeantheau (1998). As a result, they better capture the high kurtosis of squared returns. Keiler and Eder (2013) studied the systematic risk that integrates the interaction between the micro and macro stress situations as spatial econometrics parameters. Analogously, Chen (2017) showed that when the spatial weights are derived from credit rating downgrades, the multivariate spatial BEKK model can capture the spillover effects among the southern European stock index: Portugal, Italy, Ireland, Greece, and Spain (PIIGS). Zhang et al. (2018) applied the multivariate GARCH with a dynamic panel of spatial weight matrices based on the GDP. The work studies the countries’ interconnectedness of returns and uses the estimated parameters to forecast the portfolio risk of six stock indices.
Our contribution is to introduce a dynamic conditional correlation GARCH (DCC-GARCH) model with a spatial component based on the credit exposure similarity among banks derived from the EU-wide stress test data. We add the spatial components into a DCC-GARCH model to investigate whether we can better capture the spillover effects thanks to this additional information and alternative distributional assumptions of the DCC-GARCH model. We discuss and implement the estimation of the model using both Gaussian and Student’s t distributional assumptions. In particular, we estimate the individual risk via Value at Risk (VaR) and spillover risk via CoVaR. The results of the different models are compared, both amongst themselves and with the one obtained thanks to Filtered Historical Simulations (FHS). Finally, the results are backtested (Abad et al. 2014; Caporin 2008; Christoffersen and Pelletier 2004; Kupiec 1995) to evaluate the accuracy of the different VaR and CoVaR estimates, showing the superiority of the spatial DCC-GARCH model compared to the other models.
Finally, we point out that the applications of spatial DCC-GARCH models not only allow us to estimate CoVaR, but also permit investors to estimate more accurately the distribution of future returns of a set of assets to develop optimal portfolio strategies. The proposed framework has relevant applications for financial regulators interested in accurately measuring risk spillovers and systemic risk in financial systems, but it also caters to risk managers who want to measure the risk related to the interconnectedness among institutions.
The remainder of this paper is organized as follows. In Section 2, the spatial DCC-GARCH model and the estimation procedure are discussed, and the methodology for the financial application is presented. Section 3 presents the data and the empirical results, and in Section 4 we present our conclusions and discuss the results in relation to the literature.

2. Matherials and Methods

2.1. Modelling and Inference

A common feature of financial time series is the presence of volatility clustering (see, e.g., Cont 2001).1 Common tools used to address such features are Generalized Auto-Regressive Conditional Heteroscedasticity (GARCH) models (Bollerslev 1986), which generalize the ARCH models introduced by Engle (1982). Let r t be the return discrete-time process with zero mean. The standardized disturbances ε t are independent and identically distributed (i.i.d.) with zero mean, E ( ε t ε t 1 , ) = 0 , and unit variance, V a r ( ε t ε t 1 , ) = 1 . Then, the GARCH( p , q ) process for return r t is defined as
r t = h t ε t , t = 1 , , T
and
h t = ω + k = 1 q α k r t k 2 + k = 1 p β k h t k ,
where h t is the conditional variance, ω > 0 , α k 0 and β k 0 k .
When studying spillover risk, it is natural to look for multivariate extensions of the GARCH model to characterize the joint evolution of stock returns. Before presenting the multivariate model it is useful to define the following quantities of interest. Assuming a market with N assets, then at time t = 1 , , T we have:
  • r t is the vector of assets’ returns at time t,
  • H t is the conditional covariance matrix,
  • h t = d i a g ( H t ) is the vector of the univariate conditional variances,
  • D t is a squared matrix with the conditional standard deviations h t on the main diagonal and zero otherwise.
  • R t is the positive definite conditional correlation matrix,
  • Q t is the conditional covariance matrix of the standardized residuals,
  • Q ¯ is the unconditional covariance matrix of the standardized residuals,
Full generalizations of a univariate model, such as the VEC GARCH model (Bollerslev et al. 1988; Ling and McAleer 2003) or the BEKK model, (Engle and Kroner 1995) have been extensively discussed in the literature. Using matrix notation, it is possible to characterize a multivariate GARCH as follows:
r t = H t 1 / 2 ε t ,
and
H t = A 0 + k = 1 q A k r t k r t k + k = 1 p B k H t k ,
where H t , A k , B k are N × N matrices and ε t is an R N valued i.i.d. sequence of random variables with zero-mean and unit-variances (see Engle and Kroner (1995) for the restrictions required to ensure stationarity and positive semi-definiteness of the conditional covariance matrix).
Multivariate GARCH models have the drawback of having a large number of parameters, making the estimation complex and computationally challenging, hence these models are suitable only if the dimensionality N is very small. A solution to the dimensionality problem is to pose further restrictions on the multivariate process. A common restricted specification is the Constant Conditional Correlation model (CCC) proposed by Bollerslev (1990) that assumes that the conditional covariance matrix is constant over time, requiring focusing solely on the estimation of conditional variances. According to the CCC-GARCH model, Equation (1) is given by (3) and Equation (2) can be written as follows:
h t = ω + k = 1 q α k r t k 2 + k = 1 p β k h t k ,
where ω is the N × 1 dimensional vector of unconditional variances with ω R + , α k and β k are the N × 1 dimensional vector of ARCH and GARCH parameters of order q and p with α k , i R 0 + , β k , i R 0 + , and ⊙ is the Hadamard product.
The CCC-GARCH model assumes that the conditional covariance matrix, H t , can be factorized as
H t = D t R D t ,
where the correlation matrix is assumed to be constant throughout time R t = R , t and the conditional standard deviation matrix D t is a diagonal matrix given by
D t = d i a g h t .
The generic element of conditional covariance matrix H t is constructed as
H t i j = h i t ρ i j h j t , i j ; i , j = 1 , , N ,
where ρ i j = R i j is the constant conditional correlation coefficient between the ith and jth variables.
The multivariate GARCH model with a dynamic conditional correlation structure (DCC), introduced by Engle (2002), improves the dynamic relationship, assuming a time-varying correlation matrix as follows
H t = D t R t D t .
The dynamic correlation model allows R t to be time-varying, and its dynamics are modeled assuming a GARCH(1,1) process for the covariance of the standardized residuals. Hence R t is decomposed into
R t = d i a g ( Q t 1 ) Q t d i a g ( Q t 1 ) ,
where
Q t = Q ¯ ( 1 γ δ ) + γ ( ϵ t 1 ϵ t 1 ) + δ Q t 1 ,
where γ and δ are ARCH parameters and GARCH parameters of the DCC model, respectively. By following the GARCH model from Equation (2), the generic element of the time-varying conditional covariance matrix of the standardized residuals Q t i j = q i j , t can be expressed as
q i j , t = q ¯ i j 1 γ δ + γ ϵ i , t 1 ϵ j , t 1 + δ q i j , t 1 ,
where q ¯ i j = Q ¯ i j . The process is mean-reverting as long as 0 < δ < 1 and γ + δ < 1 . In the particular case of γ + δ = 1 , the process will follow the exponential smoother matrix of the standard residuals, as described in Engle (2002). Finally, the generic conditional correlation
ρ i j , t = q i j , t q i i , t q j j , t ,
can be written into matrix form as in Equation (10). Substituting the conditional correlation matrix into Equation (9), the DCC is given by
H t = D t R t D t = D t d i a g ( Q t 1 ) Q t d i a g ( Q t 1 ) D t .
Restricted GARCH models beyond CCC and DCC-GARCH have been discussed by Caporin (2008) and Billio et al. (2021), who introduce spatial matrices within BEKK models for measuring risk spillover. In these approaches, the interaction components of the model are based on spatial weight matrices provided exogenously (for instance on the basis of geographical distances among assets, or some similarity metrics). These models allow easier and more accurate estimation by effectively imposing restrictions on the parameter space. An alternative approach to improve the estimation of multivariate GARCH models is to introduce sparsity in the parameter estimates by using an L 1 penalization, as suggested by Dhaene et al. (2022).

2.1.1. Spatial DCC-GARCH

In this work, we introduce a spatial extension of the DCC-GARCH model. The model is based on the approach of Borovkova and Lopuhaa (2012). In particular, to enrich the DCC-GARCH model with a spatial component we introduce a spatial matrix W into the vector of the conditional variances h t . The resulting conditional variance is expressed as
h t = A 0 + k = 1 q ( A 1 , k + A 2 , k W ) r t k 2 + k = 1 p ( B 1 , k + B 2 , k W ) h t k ,
where A 0 = ( a 0 , 1 , , a 0 , N ) , A 1 , k , A 2 , k , B 1 , k , and B 2 , k are diagonal matrices. The term W = W i j is the weight matrix i , j = 1 , , N with j = 1 N w i j = 1 and w i i = 0 i , given by
W = 0 w 12 w 1 N w 21 0 w 2 N w N 1 w N 2 0 .
The i-th element of h t becomes
h t , i = a 0 , i + a 1 , i r t 1 , i 2 + a 2 , i X t 1 , i + b 1 , i h t 1 , i + b 2 , i Y t 1 , i ,
where X t 1 , i = j = 1 N w i j r t 1 , j 2 and Y t 1 , i = j = 1 N w i j h t 1 , j . The introduction of the spatial component results in two exogenous spatial variables in the conditional variance equation and two additional parameters a 2 , i and b 2 , i , which measure the influence of the aggregated lagged variances and squared returns of all the other assets. These two new variables measure the aggregated spillover effects. To complete the Spatial DCC-GARCH model, we then estimate the conditional correlation matrix following the two-step procedure described in Engle and Sheppard (2001), see Section 2.1.2.
The condition for the weak stationarity of the spatial GARCH model follows from the corresponding stationarity condition for E-CCC models, derived by Jeantheau (1998) and Conrad and Karanasos (2010) for E-CCC models. The positivity conditions on all GARCH coefficients are not necessary for the positivity of variance and in many empirical cases, these may be too restrictive, ruling out possible negative volatility feedback. One author Conrad and Karanasos (2010) studied the E-CCC models and stated necessary and sufficient conditions (in terms of the process parameters) for the positivity of variance; these conditions are summarized in Theorem 1 of their paper. It can be seen easily that our spatial DCC-GARCH(1,1) model is equivalent to the E-DCC model of order one, with the particular form of the parameter matrices A = A 1 + A 2 W and B = B 1 + B 2 W . So for the conditional variances to be positive, the conditions (C1)–(C3) of Theorem 1 of Conrad and Karanasos (2010) must apply. The proposed spatial DCC-GARCH(1,1) model is weakly stationary if the modulus of the largest eigenvalue of the matrix A 1 + B 1 + ( A 2 + B 2 ) W is less than 1. In that case, the unconditional variances are given by A 0 ( I ( A 1 + B 1 ) ( A 2 + B 2 ) W ) 1 . More specifically the unconditional variance of the ith bank is given by
σ i 2 = a 0 , i + ( a 2 , i + b 2 , i ) j = 1 n w i j σ j 2 1 ( a 1 , i + b 1 , i )
and it is positive for a 0 , i > 0 , a 2 , i + b 2 , i > 0 , ( a 1 , i + b 1 , i ) < 1 .

2.1.2. Estimation of the Multivariate Spatial GARCH(1,1) Model

We follow a two-step procedure for the DCC-GARCH estimation, as described in Engle and Sheppard (2001) and Engle (2002). The first step is devoted to the estimation of (16) where the exogenous variable Y t , i is not observable since it is a function of the conditional variance of the other assets. Hence, following Borovkova and Lopuhaa (2012), we start by estimating the standard univariate GARCH(1,1) models without the external regressors to obtain the initial parameters a 0 , i 0 , a 1 , i 0 , b 1 , i 0 and the estimated variances h 1 , i 0 , , h T , i 0 . Then we use an iterative procedure in which we alternate the following two steps:
  • Compute the exogenous variables Y t 1 , i given the weights w i j and the initially estimated variances h 1 , i 0 , , h T , i 0 ;
  • Estimate the complete set of parameters a 0 , i 1 , a 1 , i 1 , b 1 , i 1 , a 2 , i 1 , b 2 , i 1 and the new estimated variances h 1 , i 1 , , h T , i 1 according to Equation (16).
We iterate this procedure until the percentage variation of the estimate is less than a small threshold. For more details please refer to Borovkova and Lopuhaa (2012).
In the second step, as in Engle (2002), we maximize the quasi log-likelihood that, when the standardized error ϵ t follows a multivariate Gaussian distribution is
log L ( θ 2 | θ ^ 1 ; r 1 , , r T ) = 1 2 t = 1 T N log 2 π + 2 log | D t | + log ( | R t | ) + r t D t 1 R t 1 D t 1 r t = 1 2 t = 1 T N log 2 π + 2 log | D t | + log ( | R t | ) + ϵ t R t 1 ϵ t ,
where θ = ( θ 1 , θ 2 ) is the set of parameters of the multivariate distribution, with subsets θ 1 = ( A 0 , A 1 , A 2 , B 1 , B 2 ) being the spatial GARCH parameters estimated in the first step, and θ 2 = ( γ , δ ) the parameters of the time-varying conditional correlation that are estimated in the second step.
Excluding D t and other additive and multiplicative constants, we maximize the following function:
t = 1 T log ( | R t | ) + ϵ t R t 1 ϵ t .
The quasi-log-likelihood function under the Student’s t distribution is
log L ( θ 2 | θ ^ 1 ; r 1 , , r T ) = t = 1 T log Γ ν + N 2 log Γ ν 2 N 2 log π ν 2 + 1 2 log | D t R t D t | ν + N 2 log 1 + r t T D t 1 R t 1 D t 1 r t ν 2 ,
where ν is the degrees of freedom, θ 2 = ( γ , δ , ν ) is the set of parameters estimated in the second step, and Γ · is the Gamma function. The estimation of the model is implemented in R using the packages rugarch and rmgarch for the estimation of univariate GARCH models in the first step, and the DCC-GARCH model in the second step, respectively.
Concerning the complexity of estimation, we see that the spatial models add two parameters ( a i , 2 , b i , 2 ) for each asset. Hence the number of additional parameters scales linearly with the size of the dataset considered. We also see that Student’s t model has one extra parameter compared to the Gaussian model, and that in the limit for ν , the former converges to the latter. Moreover, the spatial model nests the non-spatial DCC-GARCH models, where the coefficients of the spatial components are restricted to zero. The Spatial DCC-GARCH model can therefore be considered parsimonious in terms of the number of parameters, especially compared to VEC GARCH or the BEKK model. One drawback of the proposed model is that it requires the exogenous identification of a spatial matrix.

2.1.3. Spatial Weight Matrix

To estimate the spatial DCC-GARCH described in Section 2.1.1, we need to specify the weight matrix W which incorporates the spatial structure defined a priori. The most intuitive way to compute the weights is to consider the geographical distance between the issuers’ market cities. However, according to Borovkova and Lopuhaa (2012), the obtained weights are not economically meaningful, and as an alternative, they consider a different set of information and compute distance in terms of GDP and market capitalization. In our work, we investigate whether the banks’ similarity of the structure of credit exposure provides some benefit in catching risk spillover effects. Hence we propose to consider the cosine similarity between exogenous information relative to the credit exposure of each bank derived from the EU-wide stress test under the European Banking Authority (EBA). The higher the cosine similarity the stronger the closeness of banks’s credit exposure. Suppose two attribute vectors of length L, U i , L = u i , 1 , u i , 2 , , u i , L and U j , L = u j , 1 , u j , 2 , , u j , L which describe the credit exposure information of bank i and j with i , j = 1 , , N . We define the cosine similarity as follows:
C i j = l = 1 L u i , l · u j , l l = 1 L u i , l 2 · l = 1 L u j , l 2 , i , j = 1 , , N , i j .
We set C i i = 0 i and we normalize the rows of C by dividing each element by the sum of the row. Doing so, we obtain the matrix W that is the spatial weight matrix used in Equation (15).

2.2. Financial Application: CoVaR

Financial institutions use VaR to measure the standalone risk. However, the measurement of individual risk is not able to explain the linkages between other financial institutions and the financial system. Systemic risk is the possibility that an event at the institutional level could trigger severe instability or collapse of an entire industry or economy. The work Adrian and Brunnermeier (2014) introduces CoVaR to help regulators to measure risk spillovers.
The Value at Risk (VaR) at level q ( 0 , 1 ) of a random variable r with cumulative distribution function F r ( . ) is defined as
V a R q ( r ) = inf { x R : F r ( x ) q } ,
where 100 ( 1 q ) % denotes the confidence level of the VaR.2 Restricting our analysis to continuous probability distribution functions, VaR can be implicitly defined as the q-quantile of the probability distribution function
V a R q ( r ) = F r 1 ( q ) .
The Conditional Value-at-Risk (CoVaR) (see Adrian and Brunnermeier 2014), denoted by C o V a R q S | C ( r i ) , is implicitly defined by the q-quantile for a continuous probability distribution function of the financial system S conditional on some event related to C ( r i ) , where r i is the return of institution i such that
Pr r S CoVaR q S | C ( r i ) | C r i = q .
The CoVaR can capture the contribution of systemic risk by conditioning the VaR to a stressed situation. It captures the spillover of risk between a particular institution and the financial system, and it is commonly used to assess the systemic risk of a bank in a financial system. Inspired by this idea, we concentrate our attention on a CoVaR pairwise analysis between institutions in order to quantify the spillover between couples of banks.3
The conditioning event C r i in the original paper by Adrian and Brunnermeier (2014) is defined as the return of the conditioning asset i being equal to its negative VaR, that is C r i : = ( r i | r i = V a R ( r i ) ) . In this work, we follow the alternative approach of Girardi and Ergün (2013) that considers as a conditioning event the return r i being smaller or equal than the following quantity: C r i : = ( r i | r i V a R ( r i ) ) . This formulation allows us to consider more severe distress events and improves the consistency of the measure with respect to the dependence parameter, allowing for backtesting. Following Girardi and Ergün (2013), the redefined C o V a R q j | i is obtained solving
Pr r j C o V a R q j | i , r i V a R q i = q 2 .
Let f r j , r i be the bivariate probability distribution function of future returns, estimated using the DCC-GARCH model with either Gaussian or Student’s t innovations, C o V a R q j | i is implicitly defined as the quantity that solves
C o V a R q x | y V a R q y f ( x , y ) d y d x = q 2 .
We compute the integral (23) on a grid of 100 values for C o V a R q j | i to find the approximated solution under the different distributional assumptions.4

2.2.1. CoVaR Based on Filtered Historical Simulations (FHS)

In order to compare our result with a model-free approach, we consider the Filtered Historical Simulations (FHS).
FHS is a well-known tool for multivariate forecasting and simulation of time series that avoids the need for distributional assumptions for the returns’ joint dynamic, relying instead on past realizations. The main novelty of this approach compared to historical simulation is to rescale the innovation by the volatility that prevails on a specific day, allowing therefore to reflect the current market conditions (Barone-Adesi et al. 2002; Giannopoulos and Tunaru 2005; Gurrola-Perez and Murphy 2015). To provide a distribution-free benchmark model for the analysis, we compute the VaR and CoVaR via FHS. Consider a time window of length T and let r t be the series of historical returns with t [ 1 , T ] . The volatility weighted returns series can be computed as follows: z t = r t × σ ^ T + 1 / σ ^ t , where σ ^ t is the volatility estimated with an Exponentially Weighted Moving Average procedure (EWMA) with decay factor λ = 0.9 at time t and σ ^ T + 1 is the one-day-ahead estimate of volatility at the end of the estimation period. In practice, implementing FHS for the estimation of VaR and CoVaR requires the following steps:
  • compute the residual (or devol) time series, dividing the returns by EWMA estimated volatility σ ^ t . This allows us to sample from approximately serially independent and identically distributed data;
  • compute the estimated empirical distribution of r ^ T + 1 (revol), multiplying the devol time series z t by the latest estimate of volatility σ ^ T + 1 and assigning to each of the possible outcome a weight 1 / T ,
  • estimate V a R q , T + 1 i and C o V a R q , T + 1 j | i by computing the empirical quantile of r ^ i , T + 1 and r ^ j , T + 1 | r ^ i , T + 1 < V a R q , T + 1 i , respectively.
The FHS approach has the advantage of being non-parametric, although it has the drawback of requiring a large number of observations to accurately estimate risk, especially for the CoVaR. For this reason, it is not suitable for small values of q. For instance, with q = 0.01 the expected number of exceedances of the CoVaR for an estimation window of 10,000 daily observations (approx 40 years) is 1, while for q = 0.05 it is 25.

2.2.2. Backtesting VaR and CoVaR

In order to test the goodness of our VaR and CoVaR estimates we estimate a time series of length τ of one-day-ahead estimates, each computed on an estimation window of T = 1000 daily observations. We consider tests based on the number of violations and specifically unconditional and conditional coverage tests (Christoffersen and Pelletier 2004; Kupiec 1995), as well as tests based on asymmetric loss functions for the VaR and CoVaR (Caporin 2008). The model that provides estimates of VaR and CoVaR with the correct number and distribution of exceedances and/or lower loss function values will be considered the more accurate.

2.2.3. Tests Based on the Number of Violations

In order to determine the accuracy of the proposed model, we consider two tests based on the number of violations.
Denote by
  • r t i the ex-post realized returns of institution i with t = 1 , , τ ;
  • V a R q , t i the ex-ante Value-at-Risk forecasts at t 1 for time t, where q is the expected coverage;
  • I t i a sequence of violation for a given interval of the Value-at-Risk forecast:
    I t i = 1 , if r t i V a R q , t i 0 , if r t i > V a R q , t i .
The first test is the Kupiec test or unconditional coverage (UC) test (Kupiec 1995). The null hypothesis that the observed failure rate p is equal to the failure rate, suggested by the confidence level of VaR, q, is tested. Thus, the null hypothesis assumes that the observed violation rate is equal to the expected violation rate. If the null hypothesis is rejected, the model is considered inaccurate at the 95% confidence level.
The conditional coverage (CC) test proposed by Christoffersen and Pelletier (2004) indicates that the number of violations must be independently distributed along the testing period where the dependence can be described as a first-order Markov sequence with a transition probability matrix given by
Π = 1 π 01 π 01 1 π 11 π 11 ,
where π 01 is the probability that, conditional on today being a non-violation, the next period is a violation, and π 11 is the probability that, conditional on today being a violation, the next period is a violation. The hypothesis to test for the conditional coverage property is H 0 : π 01 = π 11 which assesses the independence of failures on consecutive time periods.
Girardi and Ergün (2013) proposed the backtesting of C o V a R q , t j | i via a straightforward application of the standard Kupiec and Christoffersen tests considering the violations (i.e., r t j C o V a R q , t j | i ) for those time periods in which r t i V a R q , t i . Having that in mind we compute a second hit sequence, I t j | i , on the sub-sample in which r t i V a R q t i as follows:
I t j | i = 1 , if r t j C o V a R q , t j | i 0 , if r t j > C o V a R q , t j | i ,
where the number of observations of the second hit sequence is equal to the number of violations of the first hit sequence. Hence for the tests on CoVaR, the sequence of violation I t j | i can be used instead of I t i .

2.2.4. Backtesting Based on Loss Functions

The backtesting based on the confidence level of VaR estimates shows the accuracy of an individual model; however, the comparison between the different models can be limited. To overcome the drawback, Lopez (1999) proposed backtesting based on a loss function. The method focuses on the magnitude of the failure when the violation occurs. Thus, the VaR estimates under the loss function can provide the model’s performance as a numerical score. The value of the loss function at time t can be given as
l t i = g r t , V a R q , t i if r t V a R q , t i h r t , V a R q , t i if r t > V a R q i , t ,
where g ( · ) and h ( · ) are the loss functions applied to exceedances and values within the VaR, respectively. Finally L i = t = 1 T l t i is defined as the total loss. The best model can be identified by the lowest total loss. Other works by Abad et al. (2014), Caporin (2008), and Cesarone and Colucci (2016) show several alternative specifications for the loss functions g ( · ) and h ( · ) , defined from the regulator and investor’s point of view. In the regulator’s view, we consider the size of the loss only if the violation occurs:
h r t , V a R q i = 0 if r t > V a R q i .
On the contrary, the investor is interested in both sides, as an overestimation of VaR may trigger limitations from the risk management, or lead to higher capital requirements imposed by the regulator. In particular we consider the functions
g ( r t , i , V a R q , t i ) = | r t + V a R q , t i |
and
h ( r t , i , V a R q , t i ) = q 1 q g ( r t , i , V a R q , t i ) .
We underline that the resulting loss function l t i is strictly related to the Koenker loss function used for the estimation of quantile regression, defined as
l ( X , ξ , q ) = ( 1 q ) ( X ξ ) + + q ( X ξ )
where ( · ) + = min ( X , 0 ) and ( · ) = min ( X , 0 ) . In case of independent and identically distributed returns the minimization arg min ξ l ( X , ξ , α ) is the value at risk. For further details we refer to Koenker and Bassett (1978), Rockafellar and Uryasev (2013), and Giacometti et al. (2021).
We extend the backtesting procedure to the case of the CoVaR as before, estimating the measure l t j | i on a sub-sample in which r t i V a R q i , t .

3. Results

3.1. Data

We consider ten years of weekly data from seven representative banks in Italy, France, Germany, the United Kingdom, the Netherlands, Spain, and Belgium. The data span from 20 September 2010 to 18 September 2020, including 2566 daily equity log-returns. The data are downloaded from Refinitiv Eikon. We perform a rolling analysis with an estimation window of T = 1000 daily observations (approximately 4 years), forecasting one-day-ahead VaR and CoVaR, for a total of τ = 1566 out-of-sample daily observations. We use the same windows for both the DCC-GARCH models and the FHS estimation.
Table 1 reports the descriptive statistics and tests the output of the log returns for the out-of-sample period. We see that all the banks in the sample with the exception of KBC Group had negative average returns. The series have typically negative skewness and excess kurtosis, as expected from equity time series. The results of the Engle ARCH test (Engle 1982) indicate that the null hypothesis of homoscedastic returns is rejected, suggesting the need for GARCH models. The autocorrelograms and partial autocorrelograms of the returns, not reported for brevity, do not highlight relevant serial correlation structure, while the autocorrelogram of squared residuals (also omitted for brevity) show significant and persistent autocorrelations, confirming the heteroscedasticity of the data. Next, we consider the correlation between the banks. Figure 1 shows that correlations are positive and high. Figure 2 studies the evolution of correlations, computed using 6-month rolling windows. The right panel represents the dynamics of the 21 pair-wise correlations over time, while the left panel shows the Frobenius norm of the correlation matrices to provide a synthetic representation. We see that the Frobenius norm changes over time, suggesting that a CCC-GARCH model is not appropriate for the dataset. On the contrary, the time-varying correlation matrix is consistent with the assumptions of a DCC-GARCH, and the high variability in the individual correlations leaves space for spatial models, which could better characterize the multivariate stochastic process.

Spatial Weight Data

For the construction of the spatial matrix, we analyze the data from the EU-wide stress test under the European Banking Authority (EBA). The EBA stress test aims to evaluate financial institutions’ resilience to adverse market conditions. It also provides the overall assessment of systematic risk in the European banking system. In the EU-wide stress test analysis report, we consider the base scenarios for each bank and the relative credit exposure information: exposure values, risk exposure amounts, stock of provision, and leverage ratio under the internal ratings-based (IRB) approach or Standardized approach (STA) referred to credit exposure in specific asset classes,5 as presented on the EBA’s website (EBA 2021). We compute for each bank the vector of percentage exposure in each class with respect to the total exposure and the similarity between couples of vectors for the different banks, as in Equation (21). This indicator provides a broad view of the similarity between the banks’ credit structures and exposures.
We then rescale the values such that each row sums to 1, as shown in Table 2. The spatial matrix is based on the EU-wide stress test of 2018 and is kept fixed for the entire analysis. To ascertain the matrix weight’s consistency over time, we perform the inequality test for couples of matrices by Jennrich (1970) on matrices computed in different years. In particular, we compare the normalized cosine similarity matrix weight from the EU-wide stress test of 2014 vs. 2016, 2016 vs. 2018, and 2018 vs. 2014. We do not reject the null hypotheses at a 1 % significance level, suggesting that the spatial components of the EU-wide stress test do not change significantly over time.

3.2. Empirical Results

We estimate the DCC-GARCH(1,1) models using the equity log returns, considering four specifications that differ in terms of distribution and inclusion of the spatial component: Gaussian DCC (GaussDCC), spatial Gaussian DCC (GaussSpDCC), Student’s t DCC (tDCC), and spatial Student’s t (tSpDCC).6 The procedure is numerically stable, and only for a small percentage of the estimation window, do the DCC-GARCH models fail to converge. In such cases, we carry over the result from the previous estimation window.
Table 3 reports the Akaike, Bayesian, Shibata, and Hannan-Quinn information criteria, averaged across the 1566 estimation windows. Information criteria allow us to assess the quality of the model in relation to the data, controlling for both the quality of the fit and the number of parameters. Therefore, we use them to assess whether the inclusion of a distribution with more parameters (Student’s t, compared to the Gaussian) and the use of the spatial component actually improve the quality of the fit or, on the contrary, the added complexity of the model affects negatively the estimation. We see that the models based on the Student’s t distribution have lower information measures according to all the measures considered (denoting a better model for the data), and that the introduction of the spatial component improves the performance of the model.
To further illustrate the characteristics of the model, Table 4 shows the average value of the coefficients computed on the 1566 rolling windows (for brevity we report the averages of the parameters related to each bank). We notice that a relatively large part of the spatial parameters ( a 2 , i , b 2 , i ) are statistically significant at the 95% significance level, suggesting that the spatial components have some explanatory power, motivating their introduction in the model. We also see that the DCC parameters γ and δ are almost always statistically significant (confirming the presence of time-varying correlations). The ν parameter is also significant and has values close to 6, suggesting that the innovations have fat-tails and that the Student’s t model is more suitable than the Gaussian one (that is the limit of the Student’s t model with ν ) Table A1 in the Appendix A reports the coefficients of the tSpDCC model for the first estimation window as an example.

3.2.1. VaR and CoVaR

In this section, we use the spatial DCC-GARCH model to compute pair-wise CoVaRs and study risk spillover in the European banking system. We report in Table 5 the estimates of the average V a R 5 % i (diagonal elements) and C o V a R 5 % j | i (off-diagonal elements) for the four DCC-GARCH models and for the Filtered Historical Simulations (FHS) estimates. The one-day-ahead forecast of V a R 5 % i is computed using the conditional variance estimate and the parametric distribution of the model. The corresponding C o V a R 5 % is computed numerically, according to (23) using the time-varying covariance matrices.
We observe that, as expected, the CoVaR is always greater in absolute values than the VaR figures on the diagonal. Furthermore, we see that the estimates of the Value at Risk are similar across the five considered models. On the contrary, the estimates of CoVaR are significantly larger for the DCC-GARCH models based on the Student’s t distribution, suggesting that the Gaussian model may potentially underestimate the risk of joint distress and risk spillover, as expected by the stylized fact of fat tails and high tail correlations in financial time series (see, e.g., Cont 2001). Finally, we see that the estimates of CoVaR of the tSpDCC model are slightly smaller than the tDCC model. The FHS estimates yield similar results to the other models in terms of VaR, and lie in the middle between the Gaussian and Student’s t models in terms of CoVaR. Figure 3 shows the out-of-sample equity log returns and the estimate of the Value at Risk. We see that the dynamics are similar for all the models, and that the two models with a spatial component share some similar dynamics in specific time periods.
Figure 4 shows the V a R 5 % j of each bank j together with the average C o V a R 5 % j | i for i j (for brevity we report only the results estimated using the tSpDCC model). We see that for all the banks, the CoVaR is always higher than the VaR. Focusing on the dynamics, we see two main spikes in the series that affected all banks: one in mid-2016 (corresponding to the Brexit Referendum), and one when the COVID-19 crisis started in March 2020. The former shock was short-lived, and risk measures returned to normal levels quickly, while the Covid crisis had more long-lasting effects, with CoVaR decreasing slowly in the following months (although with differences across banks, for instance, Intesa San Paolo risk levels returned to a normal level quicker than ING Group). In other periods, the dynamics of CoVaR are diversified across banks, with some institutions (in particular Intesa San Paolo) characterized by several spikes (likely related to idiosyncratic or regional shocks) while other institutions such as Credit Agricole or Barclays characterized by a more stable risk profile.
The last panel of Figure 4 shows the average ratio between C o V a R 5 % j | i and V a R 5 % j for each couple of institutions, and the rescaled average log-returns of the seven banks included as reference. We see that the shocks associated with Brexit and Covid had the effect of increasing the ratio, thus increasing the risk spillover in the system. The effect is persistent, as the ratio remained higher for the periods after the shocks, despite the levels of risks reduced more quickly after the event, suggesting the presence of long-lasting spillover effects.

3.2.2. Backtesting Results

We now present two backtesting analyses in order to assess the quality of the estimation of VaR and CoVaR based on different models. First, we report the results for the backtesting of V a R 0.05 i with the Kupiec (unconditional) and Christoffersen (conditional) coverage tests. Table 6 provides the p-value of the unconditional coverage (UC) and conditional coverage (CC) tests. We see that, regarding the UC test, in most cases, the null hypothesis of correct exceedances is not rejected at the 95% confidence level, meaning that the models identify correctly the expected number of exceedances. The spatial component, according to this statistic does not provide benefits, leading typically to lower p-values and higher rejection rates. Concerning the CC test, the null hypothesis of correct and independent exceedances is rejected in a higher number of cases, indicating the presence of some residual clustering of the exceedances. Summing up, according to the CC tests, the introduction of the spatial component improves the estimation of Value at Risk, leading to higher p-values and lower rejection rates. The FHS estimation has good performance in terms of the UC test, while it is the worst one considered in terms of the CC test (the null hypothesis is rejected for all banks except for Banco de Sabadell).
Table 7 reports a summary of the UC and CC tests for the CoVaR. Panels A and B show the average number of exceedances, the average p-value across all the combinations of assets, and the percentage of p-values that are smaller than 5% (i.e., the null hypothesis is rejected at the 95% confidence level) for the UC test and CC test, respectively.7 We see that for the GaussDCC and GaussSpDCC based on the Gaussian distribution, the null hypotheses for both the UC and CC tests are rejected for all the bilateral CoVaRs. On the contrary, for the Student’s t models the null hypotheses are never rejected for both the UC and CC tests. Finally, we observe that the spatial model tSpDCC performances are aligned with those of tDCC. The results of the FHS are slightly better than the Gaussian DCC-GARCH models but are worse than the Student’s t DCC-GARCH models. The results, although they do not highlight relevant differences between the spatial and non-spatial models, do confirm that the models based on the Gaussian distribution are not suitable to properly measure spillover risk, while Student’s t models have much better results when compared to FHS estimation of CoVaR.
Finally, we assess the quality of the estimation of VaR and CoVaR using the loss function methodology described in Section 2.2.4. Table 8 reports the value of the loss function for the investor and the regulator (lower values are better). We see that, consistently with the results for CC and UC tests, the introduction of the spatial component has a positive effect for both the models based on Gaussian and Student’s t distribution. Contrary to expectations, the Student’s t models do not outperform the Gaussian models according to this metric, showing similar or slightly worse performance. The similarity between Gaussian and Student’s t distribution may be due to the fact that we are considering a low confidence level ( 95 % for q = 0.05 ): for a higher confidence level (e.g., 99 % , q = 0.01 ) the shape of the tails may matter more, and the Student’s t model may perform better. In this analysis, we do not consider higher confidence levels as the backtesting of the CoVaR would become impossible.
Looking at Table 9, we see the average value of the loss functions for the CoVaR estimations. The results highlight once again a beneficial effect of the spatial component, that leads to improvements regardless of the distribution and the loss function used. We also see that the average loss is clearly smaller for the Student’s t models, highlighting their ability to better estimate the spillover risk compared to the Gaussian models. The results for the FHS show that such a model is not as accurate as tSpDCC and tDCC, although it performs better than Gaussian DCC-GARCH models.
Overall, the results suggest that the introduction of a spatial component improves the performance of the DCC-GARCH model, both in terms of the information criteria, and in terms of the out-of-sample estimation of CoVaR, as confirmed by the backtesting. The results also confirm the better fit of the Student’s t models compared to the Gaussian models.

4. Discussion and Conclusions

As the spillover effects of risk become a problem in interconnected banking systems, this study introduces a spatial DCC-GARCH(1,1) to provide a more accurate measurement of joint tail risk in a parsimonious way. The model aims to improve the standard DCC-GARCH model (Engle 2002), without introducing the estimation complexity typical of VEC GARCH and BEKK model (Bollerslev et al. 1988; Engle and Kroner 1995; Ling and McAleer 2003). After discussing the estimation of the model, we perform an empirical analysis on the equity log returns of seven large European banks, using a matrix that reflects the similarity in credit structure and exposures. We compare four multivariate GARCH specifications (with and without spatial components, and with two alternative distributions for the innovations). The comparison of information criteria suggests that the proposed spatial model with Student’s t innovation provides a better fit compared to the alternatives.
A common approach to measuring the spillover risk is the usage of pairwise CoVaR Adrian and Brunnermeier (2014), which measures the tail risk of an institution conditional to the distress of another institution. We estimate pairwise CoVaRs using our spatial DCC-GARCH(1,1) model. Compared to other GARCH-based estimation procedures for CoVaR (see, e.g., Girardi and Ergün 2013), our framework allows us to consider a time-varying correlation matrix thanks to the DCC component, and a network dimension thanks to the spatial component. Indeed, as highlighted in the literature, the network component of risk is more and more relevant (see, e.g., Billio et al. 2012; Diebold and Yılmaz 2014), and spatial GARCH models allow us to include it while maintaining the estimation feasible.
To test the reliability of the VaR and pairwise CoVaR, we first examine their accuracy via the UC and CC tests considering different GARCH specifications and a non-parametric FHS approach. The results show that the Student’s t spatial DCC-GARCH(1,1) model (tSpDCC) provides the lowest rejection rate for CoVaR 5 % compared to other models. Second, we investigate the models via backtesting based on loss functions. The analysis confirms that the Student’s t spatial DCC GARCH(1,1) model outperforms the other DCC-GARCH specifications, as well as the filtered historical simulations model in terms of the estimation of CoVaR 5 % . Overall, from a methodological perspective, we conclude that the multivariate GARCH model with the Student’s t and the spatial component obtained thanks to the proposed similarity matrix can improve the assessment of credit risk profiles and the credit risk market’s spillover.
Concerning the economic analysis and interpretation of the empirical results, the spillover analysis shows that the dynamics of CoVaR were diversified across European banks, and that in the out-of-sample period (2014–2020) there were two main shocks common to all the institutions: the Brexit Referendum (mid-2016) and the COVID-19 Crisis (started in the first half of 2020). Both periods were associated with spikes in VaR and CoVaR, and persistent increases in the ratio of CoVaR over VaR (denoting therefore a long-lasting increase in the interconnectedness and risk spillovers). The persistent increase in CoVaRs after Brexit is consistent with Li (2020), which studies the behavior of European stock markets in a multivariate time-varying setting, finding that market co-volatility continues to be substantial and persists after Brexit despite the fact that the market adjusted quickly to the shock. Similarly, the increase of CoVaR during the Covid Crisis confirms results from the literature that show how spillover and interconnectedness increased in the first part of the Covid period. In particular, Aslam et al. (2021) studied twelve European markets using the methodology from Diebold and Yılmaz (2014) on high-frequency data and found more stable spillovers in the Covid period compared to the previous, and Foglia et al. (2022) show an increase of volatility connectedness during the Covid period across 30 major Eurozone banks.
Finally, we point out that the proposed framework not only has relevant applications for financial regulators, but it is also relevant for the asset management industry. Indeed, the proposed model allows us to estimate the joint distribution of future returns, providing asset managers with reliable inputs for optimal portfolio strategies (Meucci 2005), and risk managers with data useful to measure the risk of investment funds and to conduct stress-testing analyses based on hypothetical scenarios (see, e.g., Alexander and Sheedy 2008; Koliai 2016).
Future works may further extend the current work by testing alternative spatial matrices and applying the model to other markets and larger datasets.

Author Contributions

Conceptualization, R.G. and M.C.; methodology, R.G., G.T. and M.C.; software, G.T., R.G. and K.R.; validation, M.C.; data curation, R.G., G.T. and K.R.; writing—original draft preparation, R.G., G.T. and K.R.; writing—review and editing, M.C.; supervision, R.G. All authors have read and agreed to the published version of the manuscript.

Funding

Gabriele Torri acknowledges the financial support from the Czech Scientific Foundation under project 23-07128S, and an SGS research project of VSB–TUO (SP2023/19).

Data Availability Statement

Restrictions apply to the availability of these data. Data were obtained from the financial data provider Refinitiv Eikon and are available with the permission of Refinitiv.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Estimated parameters of Student’s t spatial DCC-GARCH(1,1) model (September 2010–September 2014).
Table A1. Estimated parameters of Student’s t spatial DCC-GARCH(1,1) model (September 2010–September 2014).
h 1 , t h 2 , t h 3 , t h 4 , t h 5 , t h 6 , t h 7 , t = 1 . 39 × 10 5 ( 0.000 ) 8 . 31 × 10 6 ( 0.000 ) 2.22 × 10 16 ( 1.000 ) 4.37 × 10 6 ( 0.460 ) 4.92 × 10 6 ( 0.847 ) 3.49 × 10 5 ( 0.523 ) 6.07 × 10 6 ( 0.341 ) + 0 . 033 0 0 0 0 0 0 ( 0.008 ) 0 0 0 0 0 0 0 0.019 0 0 0 0 0 0 ( 0.110 ) 0 0 0 0 0 0 0 0 . 032 0 0 0 0 0 0 ( 0.045 ) 0 0 0 0 0 0 0 0.015 0 0 0 0 0 0 ( 0.545 ) 0 0 0 0 0 0 0 0.039 0 0 0 0 0 0 ( 0.829 ) 0 0 0 0 0 0 0 0.078 0 0 0 0 0 0 ( 0.250 ) 0 0 0 0 0 0 0 0.056 0 0 0 0 0 0 ( 0.099 ) + 0 . 067 0 0 0 0 0 0 ( 0.026 ) 0 0 0 0 0 0 0 0 . 087 0 0 0 0 0 0 ( 0.000 ) 0 0 0 0 0 0 0 0.009 0 0 0 0 0 0 ( 0.500 ) 0 0 0 0 0 0 0 0 . 061 0 0 0 0 0 0 ( 0.006 ) 0 0 0 0 0 0 0 0.041 0 0 0 0 0 0 ( 0.565 ) 0 0 0 0 0 0 0 0.013 0 0 0 0 0 0 ( 0.534 ) 0 0 0 0 0 0 0 0.035 0 0 0 0 0 0 ( 0.241 ) W r 1 , t 1 2 r 2 , t 1 2 r 3 , t 1 2 r 4 , t 1 2 r 5 , t 1 2 r 6 , t 1 2 r 7 , t 1 2 +
0 . 936 0 0 0 0 0 0 ( 0.000 ) 0 0 0 0 0 0 0 0 . 951 0 0 0 0 0 0 ( 0.000 ) 0 0 0 0 0 0 0 0 . 927 0 0 0 0 0 0 ( 0.000 ) 0 0 0 0 0 0 0 0 . 945 0 0 0 0 0 0 ( 0.000 ) 0 0 0 0 0 0 0 0 . 927 0 0 0 0 0 0 ( 0.033 ) 0 0 0 0 0 0 0 0 . 871 0 0 0 0 0 0 ( 0.000 ) 0 0 0 0 0 0 0 0 . 918 0 0 0 0 0 0 ( 0.000 ) + 0 . 047 0 0 0 0 0 0 ( 0.015 ) 0 0 0 0 0 0 0 0 . 061 0 0 0 0 0 0 ( 0.000 ) 0 0 0 0 0 0 0 0.016 0 0 0 0 0 0 ( 0.258 ) 0 0 0 0 0 0 0 0.031 0 0 0 0 0 0 ( 0.326 ) 0 0 0 0 0 0 0 0.016 0 0 0 0 0 0 ( 0.951 ) 0 0 0 0 0 0 0 0.017 0 0 0 0 0 0 ( 0.494 ) 0 0 0 0 0 0 0 0.006 0 0 0 0 0 0 ( 0.924 ) W h 1 , t 1 h 2 , t 1 h 3 , t 1 h 4 , t 1 h 5 , t 1 h 6 , t 1 h 7 , t 1
Q 1 , 1 , t Q 1 , 7 , t Q 2 , 1 , t Q 2 , 7 , t Q 7 , 1 , t Q 7 , 7 , t = 1 0 . 008 0.000 0 . 976 0.000 Q ¯ + 0 . 008 0.000 ϵ 1 , t 1 ϵ 2 , t 1 ϵ 3 , t 1 ϵ 4 , t 1 ϵ 5 , t 1 ϵ 6 , t 1 ϵ 7 , t 1 ϵ 1 , t 1 ϵ 2 , t 1 ϵ 3 , t 1 ϵ 4 , t 1 ϵ 5 , t 1 ϵ 6 , t 1 ϵ 7 , t 1 + 0 . 976 0.000 Q 1 , 1 , t 1 Q 1 , 7 , t 1 Q 2 , 1 , t 1 Q 2 , 7 , t 1 Q 7 , 1 , t 1 Q 7 , 7 , t 1
ν = 7 . 181 0.000
Parameters of the Student-t Spatial DCC-GARCH(1,1) model estimated on the first rolling window (September 2010–September 2014). p-value in brackets, values significant at 95% confidence level in bold.

Notes

1
On the contrary the return are not typically serially correlated, consistently with the efficient market hypothesis Fama (1970).
2
We use the convention that a higher value of r is preferable to a lower value, as in the cases when r represents returns or wealth. Other works consider a random variable x such that lower values of x are preferable, as in the case of losses (e.g., Artzner et al. 1999). The signs in the definitions then need to be changed accordingly.
3
The construction of a network of bilateral CoVaR is inspired by Adrian and Brunnermeier (2014). An alternative approach to extend CoVaR to a network dimension is proposed by Torri et al. (2021) and uses penalized multivariate quantile regression.
4
Under both Student t and Gaussian distribution it can be shown that for positive correlation C o V a R q x | y is always greater than V a R q x . Hence we consider a grid of values ranging from V a R q x to infinity, such that the corresponding quantiles are equally spaced. In case of null or negative correlations, we use 0 as the lower value of the grid.
5
Central governments Institutions, Corporates, Retail, Equity, Securitization, and Other non-credit obligation assets.
6
In preliminary analyses we tested different orders of GARCH and we also considered an ARMA component, finding no relevant differences. We omit them for brevity.
7
We do not report the complete results for all the bivariate CoVaR estimations for brevity. The results are available upon request.

References

  1. Abad, Pilar, Sonia Benito, and Carmen López. 2014. A comprehensive review of value at risk methodologies. The Spanish Review of Financial Economics 12: 15–32. [Google Scholar] [CrossRef]
  2. Adrian, Tobias, and Markus Brunnermeier. 2014. CoVaR. Staff Reports 348. New York: Federal Reserve Bank of New York. [Google Scholar] [CrossRef]
  3. Alexander, Carol, and Elizabeth Sheedy. 2008. Developing a stress testing framework based on market risk models. Journal of Banking & Finance 32: 2220–36. [Google Scholar]
  4. Artzner, Philippe, Freddy Delbaen, Jean-Marc Eber, and David Heath. 1999. Coherent measures of risk. Mathematical Finance 9: 203–28. [Google Scholar] [CrossRef]
  5. Aslam, Faheem, Paulo Ferreira, Khurrum Shahzad Mughal, and Beenish Bashir. 2021. Intraday volatility spillovers among european financial markets during COVID-19. International Journal of Financial Studies 9: 5. [Google Scholar] [CrossRef]
  6. Barone-Adesi, Giovanni, Kostas Giannopoulos, and Les Vosper. 2002. Backtesting derivative portfolios with filtered historical simulation (fhs). European Financial Management 8: 31–58. [Google Scholar] [CrossRef]
  7. Billio, Monica, Massimiliano Caporin, Lorenzo Frattarolo, and Loriana Pelizzon. 2021. Networks in risk spillovers: A multivariate garch perspective. Econometrics and Statistics, in press. [Google Scholar] [CrossRef]
  8. Billio, Monica, Mila Getmansky, Andrew W. Lo, and Loriana Pelizzon. 2012. Econometric measures of connectedness and systemic risk in the finance and insurance sectors. Journal of Financial Economics 104: 535–59. [Google Scholar] [CrossRef]
  9. Bollerslev, Tim. 1986. Generalized autoregressive conditional heteroskedasticity. Journal of Econometrics 31: 307–27. [Google Scholar] [CrossRef]
  10. Bollerslev, Tim. 1990. Modelling the coherence in short-run nominal exchange rates: A multivariate generalized arch model. The Review of Economics and Statistics 72: 498–505. [Google Scholar] [CrossRef]
  11. Bollerslev, Tim, Robert F. Engle, and Jeffrey M. Wooldridge. 1988. A capital asset pricing model with time-varying covariances. Journal of Political Economy 96: 116–31. [Google Scholar] [CrossRef]
  12. Borovkova, Svetlana, and Rik Lopuhaa. 2012. Spatial GARCH: A spatial approach to multivariate volatility modeling. SSRN Electronic Journal. [Google Scholar] [CrossRef]
  13. Caporin, Massimiliano. 2008. Evaluating value-at-risk measures in the presence of long memory conditional volatility. The Journal of Risk 10: 79–110. [Google Scholar] [CrossRef]
  14. Cesarone, Francesco, and Stefano Colucci. 2016. A quick tool to forecast value-at-risk using implied and realized volatilities. Journal of Risk Model Validation 10: 71–101. [Google Scholar] [CrossRef]
  15. Chen, Xiurong. 2017. Impact effects and spatial volatility spillover effects of sovereign credit rating downgrades–empirical analysis of multivariate spatial-bekk-garch model based on symbolic transfer entropy. Boletín Técnico 55: 614–23. [Google Scholar]
  16. Christoffersen, Peter, and Denis Pelletier. 2004. Backtesting Value-at-Risk: A Duration-Based Approach. Journal of Financial Econometrics 2: 84–108. [Google Scholar] [CrossRef]
  17. Conrad, Christian, and Menelaos Karanasos. 2010. Negative volatility spillovers in the unrestricted eccc-garch model. Econometric Theory 26: 838–62. [Google Scholar] [CrossRef]
  18. Cont, Rama. 2001. Empirical properties of asset returns: Stylized facts and statistical issues. Quantitative Finance 1: 223. [Google Scholar] [CrossRef]
  19. Dhaene, Geert, Piet Sercu, and Jianbin Wu. 2022. Volatility spillovers: A sparse multivariate garch approach with an application to commodity markets. Journal of Futures Markets 42: 868–87. [Google Scholar] [CrossRef]
  20. Diebold, Francis X., and Marc Nerlove. 1989. The dynamics of exchange rate volatility: A multivariate latent factor arch model. Journal of Applied Econometrics 4: 1–21. [Google Scholar] [CrossRef]
  21. Diebold, Francis X., and Kamil Yılmaz. 2014. On the network topology of variance decompositions: Measuring the connectedness of financial firms. Journal of Econometrics 182: 119–34. [Google Scholar] [CrossRef]
  22. EBA. 2021. 2021 EU-Wide Stress Test—Methodological Note. Technical Report. La Défense: European Banking Authority. Available online: https://www.eba.europa.eu/sites/default/documents/files/document_library/RiskAnalysisandData/EU-wideStressTesting/2021/LaunchoftheST/962559/2021EU-widestresstest-MethodologicalNote.pdf (accessed on 3 September 2023).
  23. Engle, Robert. 2002. Dynamic conditional correlation: A simple class of multivariate generalized autoregressive conditional heteroskedasticity models. Journal of Business & Economic Statistics 20: 339–50. [Google Scholar]
  24. Engle, Robert F. 1982. Autoregressive conditional heteroscedasticity with estimates of the variance of United Kingdom inflation. Econometrica 50: 987–1007. [Google Scholar] [CrossRef]
  25. Engle, Robert F., and Kenneth F. Kroner. 1995. Multivariate simultaneous generalized arch. Econometric Theory 11: 122–50. [Google Scholar] [CrossRef]
  26. Engle, Robert F., and Kevin Sheppard. 2001. Theoretical and Empirical Properties of Dynamic Conditional Correlation Multivariate Garch. Working Paper 8554. Cambridge: National Bureau of Economic Research. [Google Scholar] [CrossRef]
  27. Engle, Robert F., Victor K. Ng, and Michael Rothschild. 1990. Asset pricing with a factor-arch covariance structure: Empirical estimates for treasury bills. Journal of Econometrics 45: 213–37. [Google Scholar] [CrossRef]
  28. Fama, Eugene F. 1970. Efficient capital markets: A review of theory and empirical work. The Journal of Finance 25: 383–417. [Google Scholar] [CrossRef]
  29. Foglia, Matteo, Abdelhamid Addi, and Eliana Angelini. 2022. The eurozone banking sector in the time of COVID-19: Measuring volatility connectedness. Global Finance Journal 51: 100677. [Google Scholar] [CrossRef]
  30. Giacometti, Rosella, Gabriele Torri, and Sandra Paterlini. 2021. Tail risks in large portfolio selection: Penalized quantile and expectile minimum deviation models. Quantitative Finance 21: 243–61. [Google Scholar] [CrossRef]
  31. Giannopoulos, Kostas, and Radu Tunaru. 2005. Coherent risk measures under filtered historical simulation. Journal of Banking & Finance 29: 979–96. [Google Scholar]
  32. Girardi, Giulio, and A. Tolga Ergün. 2013. Systemic risk measurement: Multivariate garch estimation of CoVaR. Journal of Banking & Finance 37: 3169–80. [Google Scholar] [CrossRef]
  33. Gurrola-Perez, Pedro, and David Murphy. 2015. Filtered Historical Simulation Value-at-Risk Models and Their Competitors. London: Bank of England. [Google Scholar]
  34. Jeantheau, Thierry. 1998. Strong consistency of estimators for multivariate arch models. Econometric Theory 14: 70–86. [Google Scholar] [CrossRef]
  35. Jennrich, Robert I. 1970. An asymptotic χ2 test for the equality of two correlation matrices. Journal of the American Statistical Association 65: 904–12. [Google Scholar]
  36. Keiler, Sebastian, and Armin Eder. 2013. CDS Spreads and Systemic Risk: A Spatial Econometric Approach. Discussion Papers 01/2013. Frankfurt: Deutsche Bundesbank. [Google Scholar] [CrossRef]
  37. Koenker, Roger, and Gilbert Bassett Jr. 1978. Regression quantiles. Econometrica: Journal of the Econometric Society 46: 33–50. [Google Scholar] [CrossRef]
  38. Koliai, Lyes. 2016. Extreme risk modeling: An evt–pair-copulas approach for financial stress tests. Journal of Banking & Finance 70: 1–22. [Google Scholar]
  39. Kupiec, Paul H. 1995. Techniques for verifying the accuracy of risk measurement models. The Journal of Derivatives 3: 73–84. [Google Scholar] [CrossRef]
  40. Li, Hong. 2020. Volatility spillovers across european stock markets under the uncertainty of brexit. Economic Modelling 84: 1–12. [Google Scholar] [CrossRef]
  41. Ling, Shiqing, and Michael McAleer. 2003. Asymptotic theory for a vector arma-garch model. Econometric Theory 19: 280–310. [Google Scholar] [CrossRef]
  42. Lopez, Jose Antonio. 1999. Methods for Evaluating Value-at-Risk Estimates. Economic Policy Review. San Francisco: Federal Reserve Bank of San Francisco, pp. 3–17. [Google Scholar] [CrossRef]
  43. Meucci, Attilio. 2005. Risk and Asset Allocation. Berlin: Springer, vol. 1. [Google Scholar]
  44. Rahman, Dima. 2014. Are banking systems increasingly fragile? Investigating financial institutions? CDS returns extreme co-movements. Quantitative Finance 14: 805–30. [Google Scholar] [CrossRef]
  45. 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]
  46. Torri, Gabriele, Rosella Giacometti, and Tomáš Tich ỳ. 2021. Network tail risk estimation in the european banking system. Journal of Economic Dynamics and Control 127: 104125. [Google Scholar] [CrossRef]
  47. Zhang, Wei-Guo, Guo-Li Mo, Fang Liu, and Yong-Jun Liu. 2018. Value-at-risk forecasts by dynamic spatial panel GJR-GARCH model for international stock indices portfolio. Soft Computing 22: 5279–97. [Google Scholar] [CrossRef]
Figure 1. Correlations, marginal distributions, and bivariate scatterplots with LOESS local regression lines of equity daily log-returns.
Figure 1. Correlations, marginal distributions, and bivariate scatterplots with LOESS local regression lines of equity daily log-returns.
Jrfm 16 00397 g001
Figure 2. Frobenius norm of rolling correlation matrix (left panel), and 21 pair-wise rolling correlations (right panel). Rolling correlations are computed using a 6-month window of daily observations.
Figure 2. Frobenius norm of rolling correlation matrix (left panel), and 21 pair-wise rolling correlations (right panel). Rolling correlations are computed using a 6-month window of daily observations.
Jrfm 16 00397 g002
Figure 3. Out-of-sample equity log-returns of the seven banks considered in the study and VaR estimates.
Figure 3. Out-of-sample equity log-returns of the seven banks considered in the study and VaR estimates.
Jrfm 16 00397 g003
Figure 4. Out-of-sample equity log-returns, VaR estimated using the tSpDCC model (blue lines), and the average of the CoVaRs for each bank (red lines). The last panel reports the average ratio between the CoVaR and the VaR (moving average over 10 days), and the (rescaled) average log-return of the seven banks for reference.
Figure 4. Out-of-sample equity log-returns, VaR estimated using the tSpDCC model (blue lines), and the average of the CoVaRs for each bank (red lines). The last panel reports the average ratio between the CoVaR and the VaR (moving average over 10 days), and the (rescaled) average log-return of the seven banks for reference.
Jrfm 16 00397 g004
Table 1. Descriptive statistics of daily equity log-returns. The table reports statistics on the univariate distributions (mean, standard deviation, skewness, and kurtosis), and the p-values of the Engle ARCH test (null hypothesis: the process is homoscedastic). We use the following abbreviations for each bank: Intesa Sanpaolo S.p.A.–Turin, Italy (ISP), Crédit Agricole Group–Montrouge, France (ACA), Deutsche Bank AG–Frankfurt am Main, Germany (DB), Barclays Plc–London, United Kingdom (BCS), ING Groep NV–Amsterdam, Netherlands (ING), Banco de Sabadell S.A.–Alicante, Spain (SAB), KBC Group NV–Bruxelles, Belgium (KBC).
Table 1. Descriptive statistics of daily equity log-returns. The table reports statistics on the univariate distributions (mean, standard deviation, skewness, and kurtosis), and the p-values of the Engle ARCH test (null hypothesis: the process is homoscedastic). We use the following abbreviations for each bank: Intesa Sanpaolo S.p.A.–Turin, Italy (ISP), Crédit Agricole Group–Montrouge, France (ACA), Deutsche Bank AG–Frankfurt am Main, Germany (DB), Barclays Plc–London, United Kingdom (BCS), ING Groep NV–Amsterdam, Netherlands (ING), Banco de Sabadell S.A.–Alicante, Spain (SAB), KBC Group NV–Bruxelles, Belgium (KBC).
BankMeanStDevSkewnessKurtosisARCH Test (p-Value)
ISP 1 × 10 4 0.0256−0.810111.6198< 2.2 × 10 16
ACA 1 × 10 4 0.0253−0.381011.2140< 2.2 × 10 16
DB 6 × 10 4 0.02440.12038.1958< 2.2 × 10 16
BCS 4 × 10 4 0.0235−0.656812.9321< 2.2 × 10 16
ING 1 × 10 4 0.0236−0.470511.4810< 2.2 × 10 16
SAB 8 × 10 4 0.0248−0.343510.9513< 2.2 × 10 16
KBC 1 × 10 4 0.0254−0.300210.0297< 2.2 × 10 16
Table 2. Normalized cosine similarity matrix.
Table 2. Normalized cosine similarity matrix.
BankISPACADBBCSINGSABKBC
ISP00.14820.17630.12210.18750.18730.1786
ACA0.145300.18070.17450.14800.15620.1953
DB0.16890.176600.13470.13060.21400.1753
BCS0.13880.20250.160000.16680.15200.1799
ING0.21300.17150.15490.166700.15960.1343
SAB0.18340.15600.21880.13080.137500.1734
KBC0.17610.19640.18040.15590.11650.17460
Table 3. Akaike, Bayesian, Shibata, and Hannan-Quinn information criteria. The table reports the average across estimation windows. The smallest value for each measure is highlighted in bold.
Table 3. Akaike, Bayesian, Shibata, and Hannan-Quinn information criteria. The table reports the average across estimation windows. The smallest value for each measure is highlighted in bold.
AkaikeBayesShibataHannan-Quinn
GaussDCC−39.18−38.96−39.18−39.09
GaussSpDCC−39.34−39.06−39.35−39.23
tDCC−39.96−39.71−39.98−39.87
tSpDCC−40.04−39.71−40.04−39.91
Table 4. The table reports the average values of the coefficients of the different specifications across the estimation windows of the multivariate GARCH. The percentages of estimation windows in which the coefficients are significantly different from 0 at the 95% confidence level are reported in italics. For brevity we report the values of a 0 , i , a 1 , i , a 2 , i , b 1 , i , b 2 , i averages across the seven banks considered in the analysis.
Table 4. The table reports the average values of the coefficients of the different specifications across the estimation windows of the multivariate GARCH. The percentages of estimation windows in which the coefficients are significantly different from 0 at the 95% confidence level are reported in italics. For brevity we report the values of a 0 , i , a 1 , i , a 2 , i , b 1 , i , b 2 , i averages across the seven banks considered in the analysis.
a 0 a 1 b 1 a 2 b 2 γ δ ν
GaussDCC00.0650.907--0.0160.885-
19%57%97%--74%97%-
GaussSpDCC00.0590.7220.0620.1090.0240.858-
37%29%81%31%12%93%99%-
tDCC00.0660.893--0.0140.9155.90
27%69%96%--95%99%100%
tSpDCC00.0480.7200.0510.1150.0150.9246.15
27%32%74%25%11%92%100%100%
Table 5. Estimated V a R 5 % i (diagonal elements) and C o V a R 5 % j | i (off-diagonal elements) for the spatial and non-spatial DCC GARCH models with Student’s t and Gaussian innovations, and for FHS. The reported values are the average of the out-of-sample estimates computed across 1566 daily rolling windows.
Table 5. Estimated V a R 5 % i (diagonal elements) and C o V a R 5 % j | i (off-diagonal elements) for the spatial and non-spatial DCC GARCH models with Student’s t and Gaussian innovations, and for FHS. The reported values are the average of the out-of-sample estimates computed across 1566 daily rolling windows.
VaR 5 % i and CoVaR 5 % j | i —FHS
BankISPACADBBCSINGSABKBC
ISP0.0340.0650.0760.0600.0640.0830.068
ACA0.0640.0330.0810.0570.0640.0780.064
DB0.0660.0650.0410.0720.0620.0770.061
BCS0.0650.0690.0900.0340.0590.0720.063
ING0.0650.0710.0760.0610.0320.0780.066
SAB0.0670.0660.0770.0580.0650.0420.063
KBC0.0650.0670.0760.0610.0660.0790.028
VaR 5 % i and CoVaR 5 % j | i —GaussDCC
BankISPACADBBCSINGSABKBC
ISP0.0350.0550.0630.0500.0520.0630.048
ACA0.0560.0340.0620.0500.0510.0620.048
DB0.0560.0550.0390.0510.0510.0620.047
BCS0.0550.0540.0630.0320.0510.0610.047
ING0.0570.0560.0640.0510.0310.0630.050
SAB0.0550.0540.0610.0490.0500.0400.047
KBC0.0560.0550.0620.0500.0520.0620.030
VaR 5 % i and CoVaR 5 % j | i —GaussSpDCC
BankISPACADBBCSINGSABKBC
ISP0.0350.0570.0640.0500.0520.0630.048
ACA0.0560.0350.0640.0500.0520.0620.048
DB0.0560.0560.0400.0510.0520.0610.047
BCS0.0550.0560.0640.0320.0510.0600.047
ING0.0580.0580.0650.0510.0320.0630.050
SAB0.0560.0560.0620.0490.0510.0390.047
KBC0.0570.0570.0630.0500.0530.0610.030
VaR 5 % i and CoVaR 5 % j | i —tDCC
BankISPACADBBCSINGSABKBC
ISP0.0330.0880.1000.0780.0830.1040.078
ACA0.0870.0330.0990.0790.0830.1020.077
DB0.0870.0870.0380.0810.0830.1010.075
BCS0.0850.0850.1010.0310.0820.1000.075
ING0.0890.0890.1020.0820.0310.1050.078
SAB0.0870.0850.0970.0770.0820.0400.075
KBC0.0880.0880.0990.0790.0830.1020.029
VaR 5 % i and CoVaR 5 % j | i —tSpDCC
BankISPACADBBCSINGSABKBC
ISP0.0330.0870.0980.0750.0810.1000.075
ACA0.0850.0340.0970.0760.0810.0980.075
DB0.0850.0860.0380.0770.0810.0970.073
BCS0.0830.0850.0980.0300.0810.0960.073
ING0.0870.0880.1000.0780.0310.1000.076
SAB0.0850.0840.0950.0740.0800.0390.073
KBC0.0860.0870.0970.0760.0810.0980.029
Table 6. p-Values of the unconditional (Panel A) and conditional (Panel B) for V a R 5 % i . The null hypotheses are “Correct Exceedances” (unconditional test) and “correct & independent exceedances” (conditional test). p-values lower than 5% are highlighted in bold.
Table 6. p-Values of the unconditional (Panel A) and conditional (Panel B) for V a R 5 % i . The null hypotheses are “Correct Exceedances” (unconditional test) and “correct & independent exceedances” (conditional test). p-values lower than 5% are highlighted in bold.
Panel A: Unconditional Coverage Test for VaR 5 % i , p-Value.
BankFHSGaussDCCGaussSpDCCtDCCtSpDCC
ISP0.7000.1130.0670.4590.113
ACA0.7560.1430.0040.5340.037
DB0.5340.6700.5890.2700.122
BCS0.6700.7890.3270.7560.589
ING0.5130.5340.3900.8800.789
SAB0.3790.9350.5890.8440.756
KBC0.2240.1800.2720.6150.700
Panel B: Conditional Coverage Test for VaR 5 % i , p -Value.
BankFHSGaussDCCGaussSpDCCtDCCtSpDCC
ISP0.0010.0020.1250.0170.119
ACA0.0330.0020.0170.0070.069
DB0.0220.0130.0420.0080.009
BCS0.0380.2550.3660.2000.413
ING0.0010.0020.0120.0010.013
SAB0.4550.9970.6550.8820.774
KBC0.0230.0950.1560.0790.035
Table 7. Average exceedances, p-values, and percentage of p-values smaller than 0.05 of the unconditional (Panel A) and conditional (Panel B) for C o V a R 5 % j | i . The null hypotheses are “Correct Exceedances” (unconditional test) and “correct & independent exceedances” (conditional test). Averages and percentages are computed across all the couples of assets. For q = 5 % the expected number of exceedances is 3.9.
Table 7. Average exceedances, p-values, and percentage of p-values smaller than 0.05 of the unconditional (Panel A) and conditional (Panel B) for C o V a R 5 % j | i . The null hypotheses are “Correct Exceedances” (unconditional test) and “correct & independent exceedances” (conditional test). Averages and percentages are computed across all the couples of assets. For q = 5 % the expected number of exceedances is 3.9.
Panel A: Unconditional Coverage Test for CoVaR 5 % j | i .
FHSGaussDCCGaussSpDCCtDCCtSpDCC
Average exceedances12.018.917.14.45.0
Average p-value0.0470.0000.0000.6390.472
(% p-values < 0.05)(88.1%)(100%)(100%)(0%)(0%)
Panel B: Conditional Coverage Test for CoVaR 5 % j | i .
FHSGaussDCCGaussSpDCCtDCCtSpDCC
Average p-value0.0610.0000.0000.4460.508
(% p-values < 0.05)(81.0%)(100%)(100%)(0%)(0%)
Table 8. Backtesting based on loss function for V a R 5 % i . Panel A reports the investor’s point of view and Panel B reports the regulator’s point of view (see Section 2.2.4). For each line the best value is highlighted in bold.
Table 8. Backtesting based on loss function for V a R 5 % i . Panel A reports the investor’s point of view and Panel B reports the regulator’s point of view (see Section 2.2.4). For each line the best value is highlighted in bold.
Panel A: Backtesting Based on Loss Function. VaR 5 % i (Investor’s Point of View).
BankFHSGaussDCCGaussSpDCCtDCCtSpDCC
ISP4.114.063.974.083.96
ACA4.214.244.144.244.18
DB4.704.744.764.714.77
BCS3.953.873.753.913.80
ING3.743.743.703.773.72
SAB4.924.854.834.894.87
KBC3.543.613.483.553.47
mean4.174.164.094.164.11
Panel B: Backtesting Based on Loss Function. VaR 5 % i (Regulator’s Point of View).
BankFHSGaussDCCGaussSpDCCtDCCtSpDCC
ISP1.241.151.051.291.16
ACA1.441.351.171.441.33
DB1.341.471.441.561.59
BCS1.151.211.111.351.28
ING1.061.111.051.191.12
SAB1.451.581.601.631.64
KBC1.141.090.971.101.02
mean1.261.281.201.371.31
Table 9. Backtesting based on loss function for C o V a R 5 % j | i . The null hypotheses are “Correct Exceedances” (unconditional test) and “correct & independent exceedances” (conditional test). Panel A reports the investor’s point of view and Panel B reports the regulator’s point of view (see Section 2.2.4). For each line, the best value is highlighted in bold.
Table 9. Backtesting based on loss function for C o V a R 5 % j | i . The null hypotheses are “Correct Exceedances” (unconditional test) and “correct & independent exceedances” (conditional test). Panel A reports the investor’s point of view and Panel B reports the regulator’s point of view (see Section 2.2.4). For each line, the best value is highlighted in bold.
Panel A: Backtesting Based on Loss Function. CoVaR 5 % j | i (Investor’s Point of View).
FHSGaussDCCGaussSpDCCtDCCtSpDCC
average loss0.3920.4650.4140.3930.371
Panel B: Backtesting Based on Loss Function. CoVaR 5 % j | i (Regulator’s Point of View).
FHSGaussDCCGaussSpDCCtDCCtSpDCC
average loss0.2330.3590.3070.1530.142
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Giacometti, R.; Torri, G.; Rujirarangsan, K.; Cameletti, M. Spatial Multivariate GARCH Models and Financial Spillovers. J. Risk Financial Manag. 2023, 16, 397. https://doi.org/10.3390/jrfm16090397

AMA Style

Giacometti R, Torri G, Rujirarangsan K, Cameletti M. Spatial Multivariate GARCH Models and Financial Spillovers. Journal of Risk and Financial Management. 2023; 16(9):397. https://doi.org/10.3390/jrfm16090397

Chicago/Turabian Style

Giacometti, Rosella, Gabriele Torri, Kamonchai Rujirarangsan, and Michela Cameletti. 2023. "Spatial Multivariate GARCH Models and Financial Spillovers" Journal of Risk and Financial Management 16, no. 9: 397. https://doi.org/10.3390/jrfm16090397

APA Style

Giacometti, R., Torri, G., Rujirarangsan, K., & Cameletti, M. (2023). Spatial Multivariate GARCH Models and Financial Spillovers. Journal of Risk and Financial Management, 16(9), 397. https://doi.org/10.3390/jrfm16090397

Article Metrics

Back to TopTop