Next Article in Journal
Can Recurrent Neural Networks Predict Inflation in Euro Zone as Good as Professional Forecasters?
Next Article in Special Issue
The Robust Efficiency Estimation in Lower Secondary Education: Cross-Country Evidence
Previous Article in Journal
Alzheimer Identification through DNA Methylation and Artificial Intelligence Techniques
Previous Article in Special Issue
Efficiency Evaluation of Regional Environmental Management Systems in Russia Using Data Envelopment Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Conditional Time-Varying General Dynamic Factor Models and Its Application to the Measurement of Volatility Spillovers across Russian Assets

by
Vladimir Balash
1,*,†,
Alexey Faizliev
1,†,
Sergei Sidorov
1,† and
Elena Chistopolskaya
2,†
1
Faculty of Mathematics and Mechanics, Saratov State University, 410012 Saratov, Russia
2
Faculty of Economics, Saratov State University, 410012 Saratov, Russia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2021, 9(19), 2484; https://doi.org/10.3390/math9192484
Submission received: 19 September 2021 / Revised: 29 September 2021 / Accepted: 30 September 2021 / Published: 4 October 2021
(This article belongs to the Special Issue Quantitative Analysis and DEA Modeling in Applied Economics)

Abstract

:
This study analyzes the spillover effects of volatility in the Russian stock market. The paper applies the Diebold–Yilmaz connectedness methodology to characterize volatility spillovers between Russian assets. The spectral representation of the forecast variance decomposition proposed by Baruník and Křehlik is used to describe the connectivity in short-term (up to 5 days), medium-term (6–20 days) and long-term (more than 20 days) time frequencies. Additionally, two new augmented models are developed and applied to evaluate conditional spillover effects in different sectors of the Russian economy for the period from January 2012 to June 2021. It is shown that spillover effects increase significantly during political and economic crises and decrease during periods of relative stability. The rising of the overall level of spillovers in the Russian stock market coincides in time with the political crisis of 2014, the intensification of anti-Russian sanctions in 2018 and the fall in oil prices and the start of the pandemic in 2020. With the consideration of the augmented models it can be argued that a significant part of the long-term spillover effects on the Russian stock market may be caused by the influence of external economic and political factors. However, volatility spillovers generated by internal Russian idiosyncratic shocks are short-term. Thus, the proposed approach provides new information on the impact of external factors on volatility spillovers in the Russian stock market.

1. Introduction

The purpose of this paper is to develop new models to analyze volatility spillover across the Russian financial market using an empirical technique for the time and frequency domains.
Over the past 10 years, many studies have investigated volatility spillover, including the interdependence relationships between global and regional markets [1,2], different segments of commodity markets [3,4,5,6,7,8], such as crude oil and the stock markets [9,10,11,12], bonds [13], banks and financial institutions [14], foreign exchange markets and cryptocurrency [15,16,17,18,19,20].
Such studies are aimed at assessing the relationship between the risk of investment in different asset classes, as well as the relationship between firms. The study of volatility is an important aspect of the estimation and modeling of risks, and the connectedness of volatilities can be used in the study of mechanisms and patterns of risk spreading, market, credit and systemic risks management, including fragility, contagion and system risk decomposition.
Recently, a variety of applicable approaches to defining and measuring interconnectedness and interdependence of financial time series have been proposed. For example, the measures based on pairwise correlations or linear Gaussian models have been widely used [21]. Many authors employ pairwise correlations [22,23], GARCH model [24], vector autoregression [25,26], CoVaR [27] and the marginal expected shortfall approach [28]. A large number of works have been devoted to the study of the cumulative effects of the interconnectedness of various market segments based on the analysis of the dynamics of volatility connectedness [2,25,29,30].
In this paper we use the approaches suggested by Diebold and Yilmaz [25] as well as by Baruník and Křehlik [30]. This methods are based on the use of a vector autoregression model for multivariate time series of volatility of a selected subset of assets. The use of the generalized forecast error variance decomposition makes it possible to construct a number of measures of both total and directional volatility spillovers, as well as to carry out decomposition for given time frequencies. Diebold and Yilmaz [29] took the nature of spillovers in a broader context, providing a rigorous framework for measuring spillovers of returns and volatility in financial markets and applying it to determine the effects of volatility propagation. The methodology was applied to measure the effects of connectivity in their follow-up work [31], in which they assessed systemic risk in the US financial sector.
As a characteristic of the connectivity of the system at various levels, from pairwise to system-wide, they proposed to use forecast error variance decomposition for the vector autoregression model. The variance decomposition shows how much of the future uncertainty of asset i is due to shocks in asset j. Variance decomposition is closely related to systemic risk measures such as the marginal expected shortfall [28,32] and CoVaR [27,33] and to various aspects of graph theory.
Representation of a complex system in the form of a graph is a convenient method for its analysis. The constructed graph can be used in various applications including portfolio selection [23], systemic risk analysis [34], stability analysis [35], detection of stock price manipulation [36]. Many studies apply and develop the market graph approach to analyze various structural properties and attributes of financial networks, such as maximum clique, maximum independent sets, degree distribution [37,38,39], clustering [40,41], graph dynamics [42]. The works [35,43] study the financial network’s structural properties and their topological stability, robustness against random vertex failures and fragilities to intentional attacks.
The analysis of network properties is also useful for portfolio investment and risk management, as well as for simplification of the portfolio selection process via targeting a group of stocks or cryptocurrencies belonging to certain regions of the stock market network, and the development of network-based algorithms for analyzing systemic risk [23,31,44,45,46].
Frequency-dependent connectivity measures are usually estimated based on the well-known Generalized Forecast Error Variance Decomposition (GFEVD) [25]. Identifying the interconnectedness in an economic system requires an understanding of its elements connectedness, that appears as responses to economics shocks, for different time-frequency (long, medium and short term). The Diebold and Yilmaz’s approach involves the decomposition of the forecast error over a sufficiently long forecasting period. Baruník and Křehlik [30] proposed a technique for decomposing the forecast variance on selected timeframe. This allows you to highlight the long-term, medium-term, and short-term impacts of a shock, which can conveniently be summed to a total aggregate effect, if needed. For this purpose, the authors defined the spectral representation of generalized forecast error variance decomposition [30].
Many authors touch upon the problem of the influence of macro variables on volatility and volatility jumps [47,48]. For example, the article [49] discusses the dependency of the density of graphs, built on a moving window, on the values of the macro variables.
Diebold’s approach and its extension [25,29,31] were applied to the US financial system, which is the most important part of the global market. In this case it can be assumed that shocks are mainly generated and propagated within the system itself. Presumably, in the case of the Russian stock market, a significant part of the shocks comes from the outside, and they affect companies in different sectors with various strength. If external shocks are not taken into account, then there is a possibility of exaggerating the volatility spillover values within the market. Therefore, we suggest considering augmented VAR models and compare the values of connectivity measures for standard V A R ( p ) and augmented models. The volatilities of commodity and stock markets were used as exogenous variables.
This approach also allows one to identify the starting moments of turbulence periods in the Russian financial market and to estimate their duration. According to the results of a retrospective analysis, two turbulence periods are well manifested, corresponding to the political crisis of 2014 and the combination of the oil shock and the beginning of the COVID epidemic in early 2020.
The number of studies examining the reaction of financial markets to news related to the development of the COVID-19 pandemic is rapidly increasing [17,50,51,52,53,54].
Our approach allows us to determine what part of the volatility spillovers was caused by the influence of external factors during the crisis. Moreover, we can decompose the total and net spillovers by the selected time frequencies using the proposed approach.
In addition, the cumulative spillover effects of volatility in the commodity and banking sectors on other industries are examined. As one might expect, volatility spillovers from the oil and metallurgical sectors outnumber backward impact from technology or service companies.
In this paper, the proposed tools are applied to describe the relationships between major Russian assets in dynamics before, after and during the financial crises of 2014 and 2020.
The work is organized as follows. Section 2 contains the methodology of spillover measurement. The methods for calculating spillover effects based on generalized forecast error variance decomposition are described. Two new ways to calculate conditional volatility spillovers are introduced. Section 3 describes the data used to measure the connectedness of companies in the RTS Index (RTS—“Russia Trading System”—is a free-float capitalization-weighted index of 50 Russian stocks traded on the Moscow Exchange, calculated in US dollars). In Section 4, a static analysis is carried out which describes the average or unconditional connectivity, while the dynamic conditional connectivity is analyzed using the moving window approach.

2. Methodology

2.1. Spillover Measurement

The methods proposed in this paper are based on the connectedness measurement framework developed by Diebold and Yilmaz [25,29,31]. This approach measures connectedness from three perspectives: total connectedness, total directional connectedness, and pairwise directional connectedness. To construct connectivity measures, the approach uses the decomposition matrices of the forecast variance for a vector autoregression (VAR) model of order p with N variables ( VAR ( p ) ):
x t = α + i = 1 p Φ i x t i + ϵ t , t = 1 , 2 , , T ,
α denotes a N × 1 intercept vector, x t = ( x 1 , t , , x N , t ) are dependent variables, Φ i is the matrix of coefficients, ϵ t ( 0 , ) is a vector of independently and identically distributed random variables (white noise with (possibly non-diagonal) covariance matrix).
This model assumes that the variables depend on their own lags of the order of p, as well as on the lags of the order p of each of the other variables in the system. Thus, the matrices of coefficients Φ i contain comprehensive information about the relationships between variables. It is convenient to work with a polynomial matrix of dimension N × N : Φ ( L ) = | I N Φ 1 L Φ p L P | . Then the model can be rewritten as Φ ( L ) x t = ϵ t . Assuming that the roots | Φ ( z ) | lie outside the unit circle, the VAR process has the following vector moving average representation M A ( ) :
x t = φ + i = 0 A i ϵ t i = A ( L ) ϵ t , t = 1 , 2 , , T ,
where the matrices A i N × N are calculated recursively from the relation Φ ( L ) = [ A ( L ) ] 1 :
A i = Φ 1 A i 1 + Φ 2 A i 2 + + Φ p A i p ,
where A 0 is the identity matrix of dimension N × N and A i = 0 for i < 0 , φ = ( I Φ 1 ) 1 α .
The approach of Koop et al. (1996) and Pesaran and Shin (1998) allows one to find the variance decomposition invariant to the order of variables in the model (GFEVD, generalized forecast error variance decomposition) [55,56]. Note that Pesaran and Shin’s approach for calculating GFEVD can be applyed to different linear and nonlinear multivariate models: VECM, GARCH, TVP-VAR, among many others [16,57,58,59,60,61].
It is possible to denote the expansion of the variance of the forecast error H-step-ahead by θ i j g ( H ) for H = 1 , 2 , , .
Then
θ i j g ( H ) = σ j j 1 h = 0 H 1 ( e i A h e j ) 2 h = 0 H 1 ( e i A h A h e i ) ,
where is the variance matrix for the error vector ϵ , σ j j is the standard deviation of the error for the j-th equation, and e i is the vector, the i-th of which is equal to one, the rest are zeros.
Following Diebold and Yilmaz, one can measure the connectedness of the system’s variables as follows. Since the sum of the elements in each row of the variance decomposition matrix is not equal to 1: j = 1 N θ i j g ( H ) 1 , then to calculate the pairwise spillover indices we should normalize the elements of the decomposition variance matrix:
θ ˜ i j g ( H ) = θ i j g ( H ) j = 1 N θ i j g ( H ) .
The overall spillover index is calculated as follows:
S i · g ( H ) = i , j = 1 ( i j ) N θ ˜ i j g ( H ) j = 1 N θ ˜ i j g ( H ) · 100 = j = 1 ( i j ) N θ ˜ i j g ( H ) N · 100 .
The generalized VAR approach allows one to estimate directional measures of volatility spillovers. The directional volatility spillover received by the i asset from all other j assets is calculated using the normalized elements of the generalized variance decomposition matrix:
S i · g ( H ) = i , j = 1 ( i j ) N θ ˜ i j g ( H ) j = 1 N θ ˜ i j g ( H ) · 100 = j = 1 ( i j ) N θ ˜ i j g ( H ) N · 100 .
The directional volatility spillover transmitted by the asset i to all other assets j is calculated in a similar way:
S · i g ( H ) = i , j = 1 ( i j ) N θ ˜ j i g ( H ) j = 1 N θ ˜ j i g ( H ) · 100 = j = 1 ( i j ) N θ ˜ j i g ( H ) N · 100 .
Finally, the net spillovers from the asset i to all other assets of j can be calculated as:
S i g ( H ) = S · i g ( H ) S i · g ( H ) .

2.2. Time and Frequency Domains of Volatility Spillovers

Baruník and Křehlik used the Fourier transform to generalize the Diebold–Yilmaz approach. They proposed a methodology for decomposition of the spillover effect into different frequency bands [30].
Let us define the function of the frequency response function: A ( e i ω ) = h e i ω h A h , which is derived from the Fourier transform for the coefficients A h , for i = 1 .
Frequency connectivity measures based on spectral decomposition of variance are defined as follows:
θ i j g ( ) = 1 2 π π π Γ i ( ω ) ( f ( ω ) ) i , j d ω ,
where Γ i ( ω ) is the weight function with respect to the frequency fraction of the variance of the variable j:
Γ i ( ω ) = ( A ( e i ω ) A ( e + i ω ) ) i , i 1 2 π π π ( A ( e i λ ) A ( e + i λ ) ) j , j d λ ,
which represents the power of the i-th variable at a given frequency.
Decomposition of the variance into components corresponding to given time intervals is determined as follows:
θ i j g ( d ) = 1 2 π d Γ i ( ω ) ( f ( ω ) ) i , j d ω ,
where d = ( a , b ) : a , b ( π , π ) , a < b .
This allows one to find the values of the coefficients S · i g ( d ) , S i · g ( d ) , S i g ( d ) , S g ( d ) for the selected range frequencies d.

2.3. Conditional Spillover Effects

Previous studies have shown the existence of significant cross-country volatility spillovers [62,63,64,65]. In particular paper, ref. [2] has demonstrated that volatility spillovers between the USA and the BRICS countries increase during periods of crises. Moreover, the USA, Brazil, and China are net transmitters, while Russia, India, and South Africa are net receivers of volatility in their respective stock markets. For the net receiver countries, volatility spillovers within their domestic markets are largely due to external influences. Costola and Lorusso identify the booms and busts in the correspondence of political and war episodes that are related to spillover effects in the Russian economy. Their findings show that the Russian Oil and Gas and Metals and Mining sectors are net shock contributors of crude oil and have the highest spillovers to other Russian sectors [66]. Therefore, the Russian stock market should not be viewed as an isolated segment of the global stock market.
It is fair to assume that the movement of prices on world commodity markets and the dynamics of indices of the largest foreign stock markets contain sufficient information about shocks external to the Russian economy. In addition, the impact of foreign markets on the prices of Russian assets far outweighs the reverse impact due to differences in the size of the economy. If these assumptions are correct, then supplementing the econometric model with the appropriate exogenous variables should eliminate the impact of foreign shocks and estimate the conditional spillover effects.
We consider two ways of including foreign shocks for assessing spillover effects.
The first approach of assessing conditional spillover effects is based on a two-step procedure (TS). In the first step, the price volatility of commodities is used as control variables. At the second step, conditional directional spillover effects are calculated using the Equations (3)–(9).
1. First-stage regressions:
x k , t = z t β k 0 + η k , t 0 , x k , t 1 = z t β k 1 + η k , t 1 , · · · · · · · · x k , t p = z t β k p + η k , t p ,
where x k t and z t = ( z 1 , t , , z M , t ) are the values of the k th initial and exogenous variables for t = 1 , 2 , , T , β k is a vector of coefficients, η k t is an error term.
Let e k t i is residuals that we obtain of first stage:
e k , t 0 = x k , t z t b k 0 , t = 1 , 2 , , T e k , t 1 = x k , t 1 z t b k 1 , t = 2 , , T · · · · · · · · · · · · e k , t p = x k , t p z t b k p , t = p , , T
where b k i is OLS coefficient.
2. Second-stage regressions
e t 0 = α + i = 1 p Φ i e t + i i + ϵ t , t = 1 , , T p
where e t i = ( e 1 , t i , , e N , t i ) .
When using the rolling window approach, the coefficients b 0 , , b p are estimated once for the full sample. The matrices Φ are evaluated for each of the windows.
Conditional spillovers are found by Equations (2)–(9).
The second approach is based on the augmented VAR model (VARX). A vector auto-regressive model with exogenous variables VARX N , M ( p , s ) can be expressed as
x t = α + i = 1 p Φ i x t i + j = 1 s Ψ j z t j + ϵ t , t = 1 , 2 , , T ,
where x t = ( x 1 , t , , x N , t ) are dependent variables, z t = ( z 1 , t , , z M , t ) are exogenous variables, Φ i and Ψ j are coefficient matrices, ϵ t ( 0 , ) is a vector of independently and identically distributed random variables (white noise with (possibly non-diagonal) covariance matrix).
VARX model representation,
x t = φ + i = 0 A i ϵ t i + j = 1 s Ψ j i = 0 A i z t j i = φ + A ( L ) ϵ t + j = 1 s Ψ j A ( L ) z t j , t = 1 , 2 , , T ,
allows us to express the variance of the forecast through the lags of the vector of errors and exogenous variables.
Appendix A contains estimates of conditional volatility spillover for some data generation processes.
If the number of variables is large enough, then the number of estimated parameters grows rapidly depending on the number of lags, and the parameter estimates becomes unstable. In this case, regularization methods can be used for estimation. We use the Least Absolute Shrinkage and Selection Operator (LASSO) for the feature selection and regularization of data models [67].
To reduce the parameter space of the VARX, the VARX-L framework applies structured convex penalties to the least squares VARX problem, resulting in the objective
min α , Φ , Ψ t = 1 T x t α i = 1 p Φ i x t i j = 1 s Ψ j z t j 2 2 + i = 1 p Φ i 1 + j = 1 s Ψ j 1 ,
in which A 2 denotes L 2 -norm, A 1 denotes the L 1 norm.
Currently, there is already a sufficient number of works devoted to the problems of using sparse VAR models [68,69,70]. The penalty parameter estimated by sequential cross-validation using algorithms implemented in the BiGVAR package https://cran.r-project.org/web/packages/BigVAR, accessed on 1 July 2021.

3. Data and Preliminary Analysis

We investigate the problem of how volatility spillovers in the Russian stock market depend on different time frequencies. After a preliminary analysis of historical data, 27 large companies included in the RTS Index were selected.
Currently, the RTS index includes the most liquid stocks of the largest Russian issuers. However, some companies included in the RTS Index (for example, CBOM, DSKY, FIVE, LNTA, MOEX, POLY, RNFT, RUAL, SFIN) were not included in our dataset due to a short trading history (about 2–3 years), which is not enough to build an appropriate model. Price data for all companies for the period from 1 January 2012 to 1 June 2021 (2400 trading days) were downloaded from https://finam.ru, accessed on 1 September 2021. 5 min high frequency data were used. In total, the dataset included over 7,000,000 records. Information on the names, tickers of assets and their sector affiliations given in Table 1.
Daily volatilities were calculated using 5-min high frequency data.
Let P i t be price of asset i at time t, p t = l o g ( P t ) . Let r i t = p i , t 1 p i , t be its logprice.
Several methods have been proposed for calculating the realized volatility and the realized variance for the whole day for the optimal combination of price volatility in the trading period and overnight. It is well known that the information received during the non-trading period (e.g., overnight) may have a significant impact on volatility [71,72,73,74,75,76]. Hansen and Lunde [71] were the first to propose the method for calculating the whole day variance, as a solution to the problem of optimizing the weights of the square of the profitability of the trading time and night time. Further, Ahoniemi and Lann [74] confirmed the effectiveness of this method using empirical data.
The realized volatility can be calculated in different ways, depending on the data used and the assumptions about price movements during the trading session. Most often, for the computation of volatility, the authors restrict the analysis to daily logarithmic realized volatility, which is computed using 5-min high frequency returns during the trading hours. For the Moscow Exchange in different years, the durations of the trading period were different. Currently, the main trading session runs from 10:00 to 18:50, an additional session from 19:00 to 23:50 Moscow time.
Let [ a , b ] be the trading session interval. To calculate the realized volatility for a trading session, the following formula is usually used:
R V t = R V t [ a , b ] = j = 1 m ( p t j p t j 1 ) 2 ,
where m is the number of 5-min intervals.
In this case, price movements that occur between the periods of current closing and opening of the next trading session are not taken into account. Let r 1 t = p t p t 1 be the overnight return. The value of r 1 t 2 is usually used as a characteristic of volatility in the period between trading sessions.
Figure 1 shows an example of realized and overnight volatility charts for SBER. Some days overnight price changes significantly exceeded the fluctuations during the trading session. Considering only realized volatility in this case would lead to a significant underestimation of the overall price change per day.
For the selected companies, the overnight volatilities are about 1/3 of the intraday volatility, in average. These values roughly correspond to the results of the works [71,74]. Variance for the whole day estimate is based on a combination of realized volatility for the active part of the trades R V 2 , t and the square of the overnight yield. Consider a linear combination of r 1 , t and R V t :
R V t ( w ) = w 1 r 1 , t 2 + w 2 R V t [ a , b ] ,
where w = ( w 1 , w 2 ) .
The Appendix C shows the values of the weights estimates obtained for the daily volatility. For more than half of the assets, the weights w 1 , w 2 turned out to be approximately equal to each other. To avoid the overcomplicating the estimation process, we took w 1 = w 2 = 1 . As a result, we used an estimate of the realized variance for the whole day without adjustment.
Exogenous variables based on a time series of prices of world commodity markets and US stock indexes (Brent oil futures—ICE BRENT CRUDE OIL PRICES, Gold COMEX, US wheat futures—ZW, LME COPPER, LME NICKEL, SP5000) were constructed. Besides that, 5-min data from https://finam.ru (accessed on 1 July 2021) were used. For time series of commodity markets prices, information is available for the entire 24 h period without interruptions, which allows us to find realized volatility without additional adjustments.

4. Empirical Results

4.1. Spillover Results

Second-order vector autoregression is considered the main model. The method of structured penalized regression (LASSO) is used for its estimation. The forecasting horizon is chosen to be equal to 100 trading days ( H = 100 ), the frequency decomposition is carried out for three intervals of H: from 0 to 5 (short-term), from 6 to 20 (medium-term), and from 21 days to 100 trading days (long-term).
A full sample pairwise spillover matrix for VAR model is presented in Table 2. The i j -th entry corresponds to the i j -th directional connectivity, that is, the contribution to the variance of the 100-day volatility forecast error for i due to shocks from j.
The column (FROM) displays the average receiving directional connectivity, that is, the average of the off-diagonal elements of each row. The (TO) row provides to the average transferring directional connectivity, that is, the average value of off-diagonal elements of columns. The bottom line (NET = TO − FROM) shows the difference between the transferring and receiving spillover. The bottom-right element (in bold) is the total connectivity (sum over the FROM line or, equivalently, the sum over the TO column).
The lines TO (0; 5], TO (6; 20], TO (20; 100], NET (0; 5], NET (6; 20], NET (21; 100]) represent decomposition of the incoming volatility spillovers by the selected frequencies.
The largest elements of the table are located on the main diagonal (intrinsic connections). The measure of overall connectivity is 83%. Note that there are blocks of high pairwise directional connectivity, as in the case of Oil & Gas companies (ROSN, NVTK, LKOH, GAZP, TATN). In addition to this, the spread of the FROM degree distribution is noticeably smaller than that of the TO degree distribution. Companies of Energy and Financial Services sectors are characterized by positive NET ratios. This suggests that the uncertainty arising in these sectors of the economy is more of an influence on the fluctuations in asset prices of companies from other industries than the opposite effect of other sectors.
Using the TS and VARX approaches allows us to analyze the influence of exogenous variables on each of the original time series. In the case of the TS method, standard procedures based on the values of the adjusted R 2 or information criteria can be used to assess the statistical significance of the model and select an informative subset of predictors.
Figure 2 shows the values of the adjusted R 2 for the regression models for the full time range. The largest values of the R 2 coefficients were obtained for the volatility of shares of large exporters of oil, non-ferrous metals, as well as banks. This is consistent with the hypothesis that the volatility of the stock prices of exporters and financial companies depends on information from foreign markets. For companies focused on the domestic market (Consumer Cyclical, Real Estate, Communications Services), the share of risk that can be explained by external economic factors is noticeably less.
The volatility of oil prices was an important factor for the Energy companies, while for the Basic Materials Companies, the metals market volatility turned out to be significant. However, oil price volatility was not an important predictor for financial services companies. The volatility of the US stock markets and the gold market was more important.
The values of the conditional spillover measures, were obtained using the VARX and TS approaches, and us allowed to come to the following conclusions. The total value of NET-effects is not very sensitive to the model’s specification, but the frequency decomposition of volatility depends on whether exogenous variables are included in the model or not. A more detailed analysis of this phenomenon is presented in the next section.
Table 3 presents the decomposition of connectedness.

4.2. Moving-Window Analysis

Analyzing connectivity across the entire sample gives a good estimate of the averaged aspects of each connectivity measure, but does not track changes over time. Therefore, it is also of interest to study interconnectedness in dynamics. One of the most common solutions to this problem is based on the Moving-Window approach. The advantages of this method are the simplicity and consistency with the mechanisms underlying the time-varying parameters. Overlapping windows of 252 calendar days (1 year) width and 1 day increments were taken to track time-varying connectivity in real time.
Diebold–Yilmaz connectivity measures provide aggregated information on how the systemic risk has changed over the period under study. However, it does not show whether shocks affect the system in the short or long term. Therefore, it is of interest to study the frequency sources of connectivity, since it is obvious that the response to shocks from different transmitters is heterogeneous and has different effects on short-term, medium-term and long-term systemic risks.
Figure 3 shows the dynamics of the total volatility connectivity index S g of the system for the VAR model (purple line), and also its decomposition S g = d s S g ( d s ) into short-term ( d 1 [ 1 ; 5 ] ) , medium-term ( d 2 ( 6 ; 20 ] ) and long-term ( d 3 ( 21 ; 100 ] ) components. The time on the scale corresponds to the end of the window.
The total connectivity in the considered period ranged from 35% to 90%. Large fluctuations in index values can be explained by the different strength and speed of the spread of shocks through the system, for example, weakening of spillover effects during safe periods, or strengthening during periods of political, financial crises, pandemic. During periods of low volatility, the overall connectedness decreased and in crisis periods, on the contrary, it increased.
For example, during the political crisis in Ukraine and the decision to join Crimea to the Russian Federation in 2014, there was an increase in volatility spillovers between Russian assets. The jump in connectivity in the second quarter of 2018 also roughly corresponds in time to the tightening of the sanctions regime against some Russian companies. Peaks in connectivity in 2020 coincide with the beginning of several adverse factors, such as the start of the COVID epidemic and a sharp drop in oil prices. High connectivity persisted for some time after the acute phases of crises. This is typical for both the 2014 political crisis and the 2020 pandemic.
It is important to note that periods of high full connectivity are mainly explained by high-frequency components, with the exception of the period of a pandemic (lockdown, economic closure, a strong drop in oil prices), which is due to the low-frequency component. After periods of crisis, markets begin to stabilize, which reduces uncertainty in low volatility periods and growing markets. This leads to the fact that shocks that create uncertainty in the system in the future will be transmitted much faster, and their impact on the system will be of a short-term nature.
The volatility of volatility has been documented by many authors. To explain this phenomenon, it is suggested that periods, when connectedness is created at high frequencies, are characteristic of financial markets in states of their growth or stability [30]. In this case, markets are able to process information quickly and a shock to one asset affects short-term cyclical behavior. If connectivity appears at lower frequencies, then shocks are transmitted over longer periods. This can be explained by the expectations of investors influencing systemic risk in the long term, which are transferred to other assets in the portfolio. This explanation does not take into account external influences and is correct if the analyzed system is large enough or isolated.
For the past two decades, raw material exporting companies have played a leading role in the Russian economy. The volatility of their shares is significantly influenced by the conjuncture of the world markets for oil, metals, etc. That is, the Russian economy is not closed or isolated, since it is significantly involved in the global exchange of goods through the export of energy resources, raw materials and metals; therefore, volatility shocks in the prices of assets traded in the Russian financial market can be generated by external factors that simultaneously affect more than one asset.
Thus, some volatility spillovers between the Russian assets may be caused by foreign commodity and stock markets fluctuation. The following mechanism of volatility shocks propagation can be assumed. Price shocks in the oil, metals and raw materials markets directly affect the asset prices of Russian exporting companies and indirectly influence the asset prices of all other Russian companies. This leads to a significant overestimation of volatility spillovers for the Energy and Basic Materials sectors. To identify the spillover effects, cleared of the influence of external markets, it is possible either to include in the VAR model exogenous variables reflecting the price volatility of commodity markets, or to try to use indirect methods of estimation. We compared two methods of estimation: VARX and 2-step procedure. Both volatilities and commodity price levels were considered exogenous variables. It is possible that in the short term, price volatility in commodity markets directly affects the price volatility of stocks in exporting companies. Higher or lower prices should lead to a revaluation of the future earnings of exporters.
The following most important explanatory variables were selected: volatility in the prices of oil, gold, aluminum, wheat, volatility of the SP500 index. Price levels turned out to be less important predictors compared to volatility.
The estimation results for TS and VARX methods are presented below. As the final specification of the VARX model, a vector autoregressive model with two lags is chosen.
We use the ratio of total mean squared errors to compare the accuracy of long-term forecasts of different models ( h = 20 or h = 100 ):
M S E ( x j , t ( H ) ) = h = 0 H 1 ( e i A h A h e i ) j j .
We denote
M S E r a t i o = M S E V A R X / M S E V A R
We also compare the mean square forecast errors for the volatility of each of the assets for a given horizon. The Table A2 shows the MSE ratio for the DY and VARX models.
The MSE ratios for the DY and VARX models for different time horizons, as well as their rolling estimates for h = 20 are given in Table A2 and Figure A1 (see Appendix B).
Exogenous factors can explain about 40% of the mean squared errors of long-term forecasts of volatility of oil companies’ stocks, while for companies focused on the domestic market the impact of commodity prices is much weaker. The forecast error using the VARX model is reduced by only 10–20% compared to the VAR model.
For the full sample, we obtain the following estimates: j M S E j V A R X ( H = 20 ) = 111.39 , M S E j V A R ( H = 20 ) = 140.95 , M S E r a t i o = 0.79 , that is, exogenous variables explained about 21% of the long-term forecast error.
Figure A1 shows a rolling estimate of the MSE-ratio for the 20-day horizon for the DY and VARX models. The impact of exogenous factors on volatility spillovers in the Russian financial market in 2016–2019 was relatively small. But during the crises of 2014 and 2020, more than half of the fluctuations could be attributed to external factors. Thus, it can be argued that, during stable periods, volatility flows are mainly generated within Russia. During crises, a significant part of volatility flows are due to external factors. As can be noted from the Figure 3, Figure 4 and Figure 5, the values of the total spillover coefficient obtained by the VAR, TS and VARX methods are quite close. However, adding exogenous variables to the model significantly reduces the estimates of conditional medium-term and long-term spillovers during crisis periods. The addition of exogenous variables is most noticeable in the case of VARX.
Estimates of conditional medium and long-term spillovers declined both during the economic and political crisis associated with the imposition of sanctions and the fall in oil prices in 2014, and the beginning of the spread of the COVID pandemic in 2020. This testifies in favor of the assumption that the spillover effects of volatility on the Russian stock market can be divided according to the types of shocks that generate them. With the exclusion of external sources of uncertainty, the systemic risk generated directly within the Russian economy spreads rather quickly. More than half of the adaptation takes place within a week, and within a month the influence of almost 85% of overflows is exhausted. In the Figure 4 and Figure 5, the estimate of the long-term component of the total overflow index changes relatively weakly. Generally, the long-term component accounted for less than 15% of the total spillover effect, while in the low volatility periods period of 2016–2018 for less than 5%. Possibly higher estimates of the long-term component of 2014–2015 associated with the fact that the exogenous variables we used not take into account all the features of a given period.
Table 4 shows the average values of the NET coefficients in total and for the frequency (21–100) for all moving windows. The overall NET ratios obtained by all three methods are in good agreement with each other. The positive spillover values for companies in the Energy, Utilities and Financial Services sectors suggest that the risks arising in these sectors spread further to the economy as a whole and outweigh the return flows in terms of their impact. Companies in the Telecommunications and Basic Materials sectors have a negative net ratio. It means that, for these companies, the influence of external sources on risk exceeds in magnitude the outflow of their own risks to the outside.
For instance, the coefficient for AFKS is equal to −1.68 and corresponds to the average value of the difference between transmitted and received flows from all other companies (TO = 33.3%, FROM = 78.9%, (33.3 − 78.9)/27 = −1.68). Note that a FROM ratio of around 80% is typical for most companies. But the transmitted spillover is the lowest in our sample. For LKOH, the outgoing effects are much higher than the incoming ones, which are reflected in the positive value of the coefficient (TO = 119.6%, FROM = 85.2%, (119.6 − 85.2)/27 = 1.27). For PIKK, the incoming effects are much higher than the outgoing ones, which is reflected in the negative value of the coefficient (TO = 43.6%, FROM = 72.1%, (43.6 − 72.1)/27 = −1.06). For Sberbank, we get (TO = 118.1%, FROM = 89.2%, (118.1 − 89.2)/27 = 1.07). For all evaluation periods, the net effects for SBER and LKOH remain positive, while AFKS and PIKK are negative. Financial services and Energy sectors have been a source of instability, while Real Estate and Communications companies have had a negative spillover effect.
NET coefficients obtained by the moving window method for AFKS and LKOH, as well as PIKK and SBER are presented in Figure 6 and Figure 7.
Figure 8 and Figure 9 present the results of frequency decomposition for the periods of crises of 2014–2015 and 2020, respectively. In February and March of 2020, the most important external shocks affecting economic processes in the Russian Federation were a sharp drop in oil prices in March 2020, the spread of the COVID pandemic and the imposition of restrictions on population mobility and economic activities. In March 2020, several simultaneous surges in volatility were observed, attributable to the breakdowns of OPEC agreements and a sharp drop in oil prices at the beginning of the month, as well as the official introduction of quarantine measures at the end of the month. The VAR model estimates lead to the conclusion of a sharp and long-term increase in long-term connectivity after the beginning of March 2020. The shocks at the start of the pandemic, as well as the oil crisis, created long-term ties above all. The diminishing of the volatility took several months.
In contrast, the use of the TS and VARX approaches with the inclusion of exogenous variables (oil, gold, etc.) indicate a fairly quick reaction of the markets. The oil and pandemic shocks were reflected in an increase in the overall interconnectedness by more than 30%, primarily due to the short-term component, and were significantly less pronounced in the medium and long-term relationships. Indeed, for most Russian companies, there was a jump in stock price volatility in March 2020 that was largely abated in a few weeks. The graphs of the realized volatility values for several companies at the beginning of 2020 are given in Figure 8.
Next, we will consider the evolution of the frequency decomposition for the period of the political and economic crisis of 2014 in Russia (Figure 9). The addition of oil and gold price level to exogenous variables did not practically change the interconnection in this time period. One may note that Black Monday (3 March 2014, marked with a dotted line), associated with the political situation in the country, led to an increase in connectivity. At the same time, a more significant jump in medium and long-term connectivity was observed during the period of falling oil prices (from 20 June 2014 to 1 April 2015, marked with dashed lines).

4.3. Directional Volatility Spillovers

The dynamic analysis of the overall connectivity provided an opportunity to elicit the main factors affecting the correlations of the volatility of the pivotal Russian shares in 2012–2021. Thus, it is possible to examine the dynamics of the directional connectivity over time. Figure 10 shows the full directional connectivity (TO and FROM) plots for each company.
The fact that FROM connectivity plots are smoother compared to TO connectivity plots is consistent with previous studies for other [25,77] markets. This is due to how volatility shocks propagate through the system. Idiosyncratic shocks are a significant component of volatility, but some of them have little effect on the volatility of other assets, while others significantly affect the volatility of many stocks in the short or long term. For example, for large companies with strong ties to other companies, the shock can be expected to be passed on to a large number of assets. As a result, individual shocks are more pronounced in the TO coefficient. The FROM indicator reflects the aggregated impact of shocks of all assets, and due to averaging, its fluctuations are weaker.
For most companies, FROM connectivity reached its highest level during the crises of 2014 and 2020. At the same time, a similar level shift in the TO connectivity indicators is not observed. That is, shocks caused by crises first affect individual companies, and then these shocks are transmitted further. During crises, these shocks became more frequent and affected more companies each time than during low volatility periods. None of the companies in our sample is immune to risk spillovers from other organizations and is not exposed to systemic risk.

4.4. Pairwise Directional Spillovers

This section deals with the differences between directed connectivity graphs constructed from the results of evaluating the VAR, TS, VARX model for some periods. The vertices of the graph correspond to the companies, and the edges correspond to the values of the indicators θ i j g ( d ) from Equation (9). We simultaneously analyze the indicators of the pairwise-spillover effect as a whole and with a breakdown into short-term, medium-term and long-term components. The first column shows graphs based on short-term relationships; the second column is based on long-term and the third column is based on overall relationships. The graphs corresponding to the frequency from 5 to 20 days are not presented, since they are quite similar to the short-term ones. The lines correspond to different models: VAR, TS, VARX, respectively. Green, orange and gray directional links correspond to the first, fifth and tenth percentile of all net directional links for the periods under consideration. The blue node color corresponds to companies that have a negative NET effect (companies that take shocks), and the red color of the node corresponds to companies with a positive NET effect (companies that distribute shocks). The node size indicates the magnitude of the NET effect in absolute terms.
The structures of graphs built during stable and crisis periods are also compared. Figure 11 shows an example of graphs built from data for the period from 1 March 2015 to1 March 2017. This period of time can be characterized as a period of stable external conditions, a stable level of prices in the markets for raw materials and oil, relatively weak volatility of most of the analyzed time series.
Note that during the low volatility periods, the graphs constructed according to various models are very similar, both for short-term, long-term frequencies, and overall. At the same time, the bulk of communications falls on short-term and medium-term frequencies. The positive NET effect is typical for Energy and Financial Services companies; for companies from other industries, incoming volatility spillovers exceed outgoing ones. Overflows with a frequency of 20 days or more account for less than 10 percent of the total flow. This time period can be considered an example of the behavior of the Russian market in the absence of significant external shocks, when most of the volatility spillovers arise under the influence of events taking place within the Russian economy. Shocks are exhausted in a relatively short period of time (up to 5 or 20 days).
During crises, the frequency of volatility shocks increases. If shocks are caused by external factors, then volatility jumps appear in several or most of the assets at the same time. Figure 12 and Figure 13 show the graphs of directional connections of crisis periods (political crisis of 2014 and oil crisis and COVID in the spring of 2020, respectively).
Note that during crisis periods, the graphs based on overall effects are very close for different model specifications.
The lists of companies that were transmitters and receivers of volatility during two crises of 2014 and 2020 are quite different. Unlike the 2014 crisis, in 2020, in addition to Real Estate and Communications Services, the receivers of volatility were added Consumer Defensive, Consumer Cyclical, Industrials and Basic Materials. In the case of using the VAR model, most of the flows occur on a frequency of a month or more. Those shocks are long term. The addition of exogenous variables (TS and VARX models) to the model leads to a redistribution of a part of the long-term spillover effects in favor of the medium-term and short-term ones.
Volatility spillovers for oil exporting countries have been the subject of analysis by many researchers. For example, Mikhailov [78] analyze the features of the bidirectional volatility spillover effect between the stock and foreign exchange markets for four emerging markets. The results indicate that bidirectional spillover effect turned out to be the most visible in the Russian stock market, because investors in Russian stocks are the most careful and react sharply to bad news and changes in the exchange rate. Pavlova and et.al. [79] study the dynamic spillover of crude oil prices and volatilities on sovereign risk premia of ten oil-exporting countries. The results indicate that Venezuela, Colombia, Russia, and Mexico are the top recipients of crude oil shocks. The effect of political variables, and aggregate demand and supply shocks are relatively less than the oil-specific shocks. The study [80] examines the volatility spillovers between gold, oil futures, and stock markets. Authors find evidence of time-varying volatility spillovers, which are intensified under major events. A portfolio management analysis reveals that a mixed portfolio (commodity and stock markets) provides a higher level of hedging effectiveness for both emerging and developed markets. Kai Shi [81] attempts to comprehensively decode the connectedness among the stock markets of BRICS countries between 1 August 2002 and 31 December 2019 not only in the time domain but also in the frequency domain. The paper shows that China’s and Russia’s stock markets play the influential role for return spillover and volatility spillover across BRICS markets, respectively. Paper [82] examines the direction and extent of the asymmetric volatility connectedness among international equity markets. The authors analyze asymmetric volatility connectedness using realized volatility and identify the magnitude of the volatility spillover and the connectedness in networks. It is shown that macroeconomic shocks increase volatility asymmetry.
As far as we know, no one has previously used the method we have developed in this paper for evaluating conditional spillover effects. In most of the previous studies, time varying spillovers are firstly estimated by the moving window method. Then, the dependence of the spillovers on selected factors [66] is analyzed. This is easy to do for a general spillover. However, if one would like to estimate directional and pairwise effects, then the number of equations becomes too large, which complicates the analysis and the formulation of conclusions. Our approach suggests a uniform way to include in the model the variables that reflect the impact of shocks external to the system under study.
The estimates of conditional spillovers significantly depend on the choice of exogenous variables. The huge number of factors can potentially influence volatility spillovers between Russian assets. Perhaps the use of the FAVAR approach could solve the problem of factors selection. In our study, we choose the exogenous variables based on generally accepted opinions about the most important foreign factors affecting the prices of Russian export goods and cross-border cash flows. However, we did not consider the impact of the government and the Central Bank monetary and fiscal policy. The comparison of the impact of external and internal shocks can be the subject of additional research. To do this, the conditional spillover coefficients should be multiplied by the ratio of the forecast variance of the conditional and unconditional models. This could reveal how the absolute values of spillovers change when exogenous variables are added to the model.

5. Conclusions

This article examines the volatility spillovers between the stocks of the 27 largest Russian assets. The approach proposed in this paper extends:
  • the estimation method developed by Francis X. Diebold and Kamil Yilmaz [29] for estimation both total and directional spillovers;
  • the Frequency Domains Decomposition technique proposed by Jozef Baruník and Tomas Křehlik [30].
The results of estimating volatility spillovers for the model without exogenous and with exogenous variables were compared. We analyzed 5-minute data on the stock prices of the 27 largest Russian companies from 1 Jaunary 2012 to 1 June 2021. The volatilities of commodity prices and foreign stock markets were used as exogenous variables.
We estimated the parameters of the models, both for the sample as a whole, and using the sliding window method, lasting 1 year with a step of 1 day. Frequency decomposition was performed for the intervals of 0–5 days, 6–20 days, 21 days or more.
In 2012–2021, the transmitters of spillover effects were the companies of Oil and Gas and Financial Services sectors, which, on average, were characterized by positive values of NET-coefficients. Companies from Real Estate, Consumer Cyclical and Communications Services are characterized by negative values of NET-effects. This is consistent with the leading role of the Energy sector as the main provider of export revenues. During periods of stable external economic and political conditions, a large share of the spillovers accounted for the short-term component.
The Russian stock market is closely integrated into the global financial system and, as shown in previous studies, is the net receivers of volatility. Moreover, such spillovers are especially intensified during times of crises in the world commodity and stock markets, as well as periods of increasing political risks. To determine the net effects generated by the domestic risk factors impacts of Russian companies the conditional spillovers were calculated. For this, the basic VAR model was extended by exogenous variables.
In times of crisis, methods that take into account and do not take into account external factors led to significantly different results of frequency decomposition. We present the results of the analysis of the directed relationship graphs for the periods of crises of 2014 and 2020. The crisis of 2014 coincides in time with a period of external political uncertainty, and 2020 with a sharp drop in oil prices and the onset of the COVID pandemic.
Application of the Baruník and Krehlík method without exogenous variables led to the conclusion that the long-term component predominated. However, in the case of adding exogenous variables to the model, conditional overflows were redistributed in favor of short-term ones. This suggests that volatility shocks generated within the Russian economy are absorbed rather quickly (more than 50% per week, 90% per month). However, it takes much longer to overcome external shocks associated with changes in the conjuncture of commodity markets or political crises.
The clusters of companies that are transmitters and receivers of volatility during the crises of 2014 and 2020 do not coincide. In 2020, unlike in 2014, Consumer Defensive, Consumer Cyclical, Industrials and Basic Materials were added to the volatility receivers.
The results of the study provide new information on the impact of external factors on volatility flows in the Russian stock market, the significance of shifts in external and internal transmitters of volatility, switching between low and high volatility modes.
The research results can also be used to assess the potential selection of risk-adjusted portfolios and thus provide new information for investor portfolio decisions.
We believe that comparisons of unconditional and conditional estimates of volatility spillovers can be useful in the assessment and management of systemic risk. Unconditional measures depict the connectivity of asset price volatility without considering the impact of external political or economic factors. Conditional measures allow us, due to the inclusion of exogenous variables in the model, to obtain measures of the connectivity of the system under study, cleared from the influence of the selected factors. In addition, comparing unconditional and conditional networks of financial connectivity allows one to trace the most important links in the sense of systemic risk propagation.
Of course, in future studies, it is advisable to expand the composition of exogenous variables to include a larger number of both external and internal factors. Spillover shocks can result from decisions of the central bank or government, for example, changes in key rate, taxes, tariffs, or measures aimed at supporting specific industries. The inclusion of such variables would allow one to examine the response of the spillovers around financial markets to monetary or fiscal policy shocks.

Author Contributions

Conceptualization, V.B. and A.F.; methodology, V.B.; software, V.B. and A.F.; validation, S.S., V.B. and E.C.; formal analysis, V.B., A.F. and S.S.; writing—original draft preparation, V.B. and A.F.; writing—review and editing, S.S. and E.C.; visualization, A.F. and V.B.; supervision, V.B. All authors have read and agreed to the published version of the manuscript.

Funding

The work was supported by the Russian Science Foundation, Project 19-18-00199.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data are freely available and was taken from www.finam.ru.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Simulation Study

To show the importance of conditional measures, we examine the processes that generate volatility spillovers by simulation. Suppose the data is generated from the following equations:
y 1 , t = β 1 y 1 , t 1 + γ 1 y 2 , t 1 + s 1 z t + ϵ 1 , t y 2 , t = γ 2 y 1 , t 1 + β 2 y 2 , t 1 + s 2 z t + ϵ 2 , t z t = ρ z t 1 + η t
where ( ϵ 1 , t , ϵ 2 , t , η 1 , t ) N ( 0 , I ) .
We are considering the following cases:
  • Bivariate VAR(1) model ( s 1 = s 2 = 0 , ρ = 0 );
  • Augmented bivariate VAR(1) model ( s 1 = s 2 > 0 , ρ = 0 );
  • Augmented bivariate VAR(1) model with autocorrelated exogenous variable ( s 1 = s 2 > 0 , 0 < ρ < 1 );
  • Three-variable vector autoregressive model:
    y 1 , t = β 1 y 1 , t 1 + γ 1 y 2 , t 1 + s 1 z t 1 + ϵ 1 , t y 2 , t = γ 2 y 1 , t 1 + β 2 y 2 , t 1 + s 2 z t 1 + ϵ 2 , t z t = α 1 y 1 , t 1 + α 2 y 2 , t 1 + ρ z t 1 + η t ,
    where ( ρ > 0 , s 1 = s 2 > 0 , α 1 0 , α 2 0 ).
The length of the simulated time series is T = 2500 . The estimates are computed averaging over the 100 simulations, and the standard error is the sample standard deviation. The simulations are carried out with the given values of the following parameters: β 1 = β 2 = 0.6 , γ 1 = 0.4 , γ 2 = 0.1 .
We start with the bivariate autoregressive case with uncorrelated errors. Further, we assume that the exogenous variable is z i i d ( 0 , 1 ) . Then the data generation process can be represented as V A R ( 1 ) with error covariance matrix:
= 1 + s 1 s 2 s 1 s 2 s 1 s 2 1 + s 1 s 2
For this model, all assumptions of the Diebold–Yilmaz approach are fulfilled, and the Baruník-Křehlik method can be used for spectral representation for variance decompositions.
We then consider the case where the exogenous variable is generated by AR(1) processe with ρ = 0.5 . Finally, we assume that the variable z is endogenous ( s 1 = s 2 = 0.4 , ρ = 0.5 , α 1 = α 2 = 0.1 ).
The purpose of the simulation example is to demonstrate the consequences of restrictions when choosing the composition of time series for assessing volatility spillovers. Variables y 1 , y 2 correspond to the subset of indicators selected for analysis. The variable z reflects the influence of factors external to the system under consideration. We analyze 4 cases: there are no external factors (model (i)); exogenous non-autocorrelated external variables (model (ii)); exogenous variables with time dependency (model (iii)); endogenous external factors (model (iv)).
For cases (iii) and (iv), assumption DY does not hold, so we use the methods described in Section 2.3 for estimating conditional spillovers. As exogenous variables, we used z t , z t 1 . Thus, the regression equations used in the first step of the TS procedure were as follows:
y 1 , t = β 1 y 1 , t 1 + γ 1 y 2 , t 1 + s 11 z t + s 12 z t 1 + ϵ 1 , t y 1 , t = β 1 y 1 , t 1 + γ 1 y 2 , t 1 + s 21 z t + s 22 z t 1 + ϵ 1 , t
Also, in the augmented vector autoregression model, the variable z and its first lag were used. The results obtained by simulations are presented in Table A1.
Table A1. Simulation results.
Table A1. Simulation results.
ModelTotal(0, 5)(6, 20)(21, 100)
(i) s 1 = s 2 = 0 , ρ = 0
DY-BK20.3222.4358.4759.413
(1.239)(0.142)(0.38)(0.896)
(ii) s 1 = s 2 = 0.4 , ρ = 0
DY-BK20.2952.4138.459.431
(1.353)(0.173)(0.475)(0.917)
TS20.3192.4158.4629.441
(1.139)(0.161)(0.401)(0.776)
VARX20.3192.4158.4629.441
(1.14)(0.161)(0.402)(0.776)
(iii) s 1 = s 2 = 0.4 , ρ = 0.5
DY-BK27.1442.149.32415.68
(1.318)(0.152)(0.406)(1.164)
TS20.4682.4158.4969.557
(1.255)(0.171)(0.423)(0.816)
VARX20.4682.4158.4969.557
(1.256)(0.171)(0.424)(0.816)
(iv) s 1 = s 2 = 0.4 , ρ = 0.5 , α 1 = α 2 = 0.1
DY-BK45.0490.4792.07142.5
(0.758)(0.062)(0.481)(1.277)
TS20.4192.4188.4959.506
(0.955)(0.147)(0.352)(0.604)
VARX20.4182.4188.4949.505
(0.954)(0.147)(0.352)(0.604)
Under the assumption of the Diebold and Yilmaz model unconditional and conditional volatility spillovers are close to theoretical values.
However, if external factors are generated by the autoregression process ( ρ > 0 ), then the spillovers found by the DY-BK method change significantly. The value of the long-term component changes especially noticeably. For the four considered models of data generation processes, the average values of the long-term component were respectively 9.413; 9.431; 15.68; 42.50.
The methods we have proposed in Section 2.3 for assessing conditional spillover effects (TS, VARX) make it possible to obtain estimates that are close to theoretical ones.
Note that the estimates obtained by the TS and VARX methods practically coincide.
However, when analyzing real data, we performed the first step of the TS procedure for the entire sample. The second step is either for the sample as a whole, or for sliding windows with a step-by-step procedure for selecting predictors. This made it possible to guarantee the participation of exogenous variables in the estimation process. In the case of the rolling window and VARX procedure, the same weights were used for both endogenous and exogenous variables of the LASSO method. Therefore, for some windows, not all exogenous variables in each of the equations of the VARX model could be used.

Appendix B. The Mean Squared Error Ratio of the h-Step Forecast

Table A2. Mean square error ratio of the h-step forecast for full sample.
Table A2. Mean square error ratio of the h-step forecast for full sample.
Ticker(0;5](0;20](0;100]
AFKS0.970.920.91
MTSS0.950.90.89
MVID0.920.880.88
PHOR0.930.880.87
PIKK0.920.880.87
FEES0.910.840.83
MGNT0.880.830.82
LSRG0.890.820.82
PLZL0.890.810.81
RTKM0.880.810.81
ALRS0.880.80.79
AFLT0.850.780.78
HYDR0.870.780.77
CHMF0.850.780.77
MAGN0.860.770.77
GAZP0.810.770.76
VTBR0.840.760.76
NLMK0.830.760.75
IRAO0.870.730.71
SNGS0.80.710.71
SBER0.780.710.7
TRNFP0.80.70.69
GMKN0.780.680.67
NVTK0.750.670.66
ROSN0.70.610.61
TATN0.680.590.59
LKOH0.670.580.58
Figure A1. Dynamic M S E V A R / M S E D Y .
Figure A1. Dynamic M S E V A R / M S E D Y .
Mathematics 09 02484 g0a1

Appendix C. Optimal Measure of Daily Volatility

We shall employ the estimator that Hansen & Lunde (2005) [71] refer to as the Newey-West modified realized variance:
R V N W [ a , b ] = j = 1 m r t j 2 + 2 h = 1 q ( 1 h q + 1 ) j = 1 m h r t j r t j + m
The optimal weights were calculated as follows
  • w 1 = ( 1 ϕ ) μ 0 μ 1 ,
  • w 2 = ϕ μ 0 μ 2 ,
  • where
  • μ 0 = E ( I V ) , μ 1 = E ( r 1 ) , μ 2 = E ( R V N W [ a , b ] ) , η 1 2 = v a r ( r 1 ) , η 2 2 = v a r ( R V N W [ a , b ] ) , η 12 = c o v ( r 1 , t , R V N W [ a , b ] ) ,
  • ϕ = μ 2 2 η 1 2 μ 1 μ 2 η 12 μ 2 2 η 1 2 + μ 1 2 η 2 2 2 μ 1 μ 2 η 12 .
To find the parameter η 12 , robust estimates were used.
Table A3. Optimal measure of daily volatility: Empirical estimates of components.
Table A3. Optimal measure of daily volatility: Empirical estimates of components.
TICKER μ 0 μ 1 μ 2 η 1 η 2 η 12 ϕ w 1 w 2
1AFKS10.4422.14210.37116.58839.8980.2140.8020.9650.808
2AFLT5.5160.9895.4584.58111.6010.1520.8270.9660.836
3ALRS5.3930.7165.2782.7099.0060.1290.8321.2640.850
4CHMF3.8600.7163.8234.0325.7730.0750.9340.3580.943
5FEES6.4411.0356.3594.34213.5050.1580.7971.2650.807
6GAZP2.7450.6132.7254.2326.3180.0380.8990.4520.906
7GMKN3.4470.6773.4032.5336.7700.0630.7801.1180.791
8HYDR4.4640.6774.4423.2228.8440.0780.8520.9770.856
9IRAO6.3530.9206.2344.92214.9550.1430.8331.1530.849
10LKOH3.0670.6373.0256.3358.5370.0380.9260.3580.939
11LSRG8.2891.0838.1285.39113.5900.2240.8990.7700.917
12MAGN5.1020.8115.0453.4377.7110.1190.8860.7170.896
13MGNT4.5990.7564.5504.7908.2850.0760.9240.4620.934
14MTSS3.7670.7473.7253.8408.3780.0520.8400.8080.849
15MVID8.1451.4028.09810.79519.5950.2570.9100.5210.916
16NLMK4.4230.7324.3632.6436.3760.0960.8610.8420.873
17NVTK4.6260.7644.5845.41911.1140.0670.8960.6320.904
18PHOR8.3410.8248.1243.05520.7930.1590.6783.2620.696
19PIKK6.4080.6536.3552.85114.4620.0950.7872.0920.793
20PLZL8.2681.0548.1744.52319.3260.1240.7681.8240.776
21ROSN3.2380.7423.2126.2028.4340.0610.9110.3910.918
22RTKM3.5990.6113.5703.7448.9500.0520.8570.8430.864
23SBER3.5010.8383.4764.2868.0220.0860.8320.7040.838
24SNGS4.1240.6804.0796.15110.4140.0530.9260.4470.937
25TATN4.6540.7884.5864.9158.3140.0890.9230.4570.936
26TRNFP4.7780.5644.7203.2779.7570.0840.8880.9460.899
27VTBR3.6010.6393.5623.4716.2630.0630.9060.5300.916

References

  1. Ng, A. Volatility spillover effects from Japan and the US to the Pacific–Basin. J. Int. Money Financ. 2000, 19, 207–233. [Google Scholar] [CrossRef]
  2. McIver, R.P.; Kang, S.H. Financial crises and the dynamics of the spillovers between the U.S. and BRICS stock markets. Res. Int. Bus. Financ. 2020, 54, 101276. [Google Scholar] [CrossRef]
  3. Wu, F.; Guan, Z.; Myers, R.J. Volatility spillover effects and cross hedging in corn and crude oil futures. J. Futur. Mark. 2011, 31, 1052–1075. [Google Scholar] [CrossRef]
  4. Roy, R.P.; Sinha Roy, S. Financial contagion and volatility spillover: An exploration into Indian commodity derivative market. Econ. Model. 2017, 67, 368–380. [Google Scholar] [CrossRef]
  5. Yip, P.S.; Brooks, R.; Do, H.X.; Nguyen, D.K. Dynamic volatility spillover effects between oil and agricultural products. Int. Rev. Financ. Anal. 2020, 69, 101465. [Google Scholar] [CrossRef]
  6. Yoon, S.M.; Al Mamun, M.; Uddin, G.S.; Kang, S.H. Network connectedness and net spillover between financial and commodity markets. N. Am. J. Econ. Financ. 2019, 48, 801–818. [Google Scholar] [CrossRef]
  7. Al-Yahyaee, K.H.; Mensi, W.; Sensoy, A.; Kang, S.H. Energy, precious metals, and GCC stock markets: Is there any risk spillover? Pac.-Basin Financ. J. 2019, 56, 45–70. [Google Scholar] [CrossRef]
  8. Jiang, Y.; Fu, Y.; Ruan, W. Risk spillovers and portfolio management between precious metal and BRICS stock markets. Phys. A Stat. Mech. Its Appl. 2019, 534, 120993. [Google Scholar] [CrossRef]
  9. Geng, J.B.; Du, Y.J.; Ji, Q.; Zhang, D. Modeling return and volatility spillover networks of global new energy companies. Renew. Sustain. Energy Rev. 2021, 135, 110214. [Google Scholar] [CrossRef]
  10. An, S.; Gao, X.; An, H.; An, F.; Sun, Q.; Liu, S. Windowed volatility spillover effects among crude oil prices. Energy 2020, 200, 117521. [Google Scholar] [CrossRef]
  11. Karali, B.; Ramirez, O.A. Macro determinants of volatility and volatility spillover in energy markets. Energy Econ. 2014, 46, 413–421. [Google Scholar] [CrossRef]
  12. Chen, Y.; Li, W.; Qu, F. Dynamic asymmetric spillovers and volatility interdependence on China’s stock market. Phys. A Stat. Mech. Its Appl. 2019, 523, 825–838. [Google Scholar] [CrossRef]
  13. Fernandez-Rodriguez, F.; Gomez-Puig, M.; Sosvilla-Rivero, S. Volatility spillovers in EMU sovereign bond markets. Int. Rev. Econ. Financ. 2015, 39, 337–352. [Google Scholar] [CrossRef]
  14. Alemany, A.; Ballester, L.; González-Urteaga, A. Volatility spillovers in the European bank CDS market. Financ. Res. Lett. 2015, 13, 137–147. [Google Scholar] [CrossRef]
  15. Liow, K.H. Volatility spillover dynamics and relationship across G7 financial markets. N. Am. J. Econ. Financ. 2015, 33, 328–365. [Google Scholar] [CrossRef]
  16. Giudici, P.; Leach, T.; Pagnottoni, P. Libra or Librae? Basket based stablecoins to mitigate foreign exchange volatility spillovers. Financ. Res. Lett. 2021, 102054. [Google Scholar] [CrossRef]
  17. Pagnottoni, P.; Spelta, A.; Pecora, N.; Flori, A.; Pammolli, F. Financial earthquakes: SARS-CoV-2 news shock propagation in stock and sovereign bond markets. Phys. A Stat. Mech. Its Appl. 2021, 582, 126240. [Google Scholar] [CrossRef]
  18. Giudici, P.; Pagnottoni, P. High Frequency Price Change Spillovers in Bitcoin Markets. Risks 2019, 7, 111. [Google Scholar] [CrossRef] [Green Version]
  19. Yi, S.; Xu, Z.; Wang, G.J. Volatility connectedness in the cryptocurrency market: Is Bitcoin a dominant cryptocurrency? Int. Rev. Financ. Anal. 2018, 60, 98–114. [Google Scholar] [CrossRef]
  20. Zieba, D.; Kokoszczynski, R.; Sledziewska, K. Shock transmission in the cryptocurrency market. Is Bitcoin the most influential? Int. Rev. Financ. Anal. 2019, 64, 102–125. [Google Scholar] [CrossRef]
  21. Engle, R.; Kelly, B. Dynamic Equicorrelation. J. Bus. Econ. Stat. 2011, 30, 212–228. [Google Scholar] [CrossRef]
  22. Millington, T.; Niranjan, M. Partial correlation financial networks. Appl. Netw. Sci. 2020, 5. [Google Scholar] [CrossRef] [Green Version]
  23. Peralta, G.; Zareei, A. A network approach to portfolio selection. J. Empir. Financ. 2016, 38, 157–180. [Google Scholar] [CrossRef]
  24. Kumar, A.S.; Anandarao, S. Volatility spillover in crypto-currency markets: Some evidences from GARCH and wavelet analysis. Phys. A Stat. Mech. Its Appl. 2019, 524, 448–458. [Google Scholar] [CrossRef]
  25. Diebold, F.X.; Yilmaz, K. Better to give than to receive: Predictive directional measurement of volatility spillovers. Int. J. Forecast. 2012, 28, 57–66. [Google Scholar] [CrossRef] [Green Version]
  26. BenSaïda, A.; Litimi, H.; Abdallah, O. Volatility spillover shifts in global financial markets. Econ. Model. 2018, 73, 343–353. [Google Scholar] [CrossRef]
  27. Adrian, T.; Brunnermeier, M.K. CoVaR. Am. Econ. Rev. 2016, 106, 1705–1741. [Google Scholar] [CrossRef]
  28. Acharya, V.; Pedersen, L.; Philippon, T.; Richardson, M. Measuring Systemic Risk. Rev. Financ. Stud. 2017, 30, 2–47. [Google Scholar] [CrossRef]
  29. Diebold, F.X.; Yilmaz, K. Measuring Financial Asset Return and Volatility Spillovers, with Application to Global Equity Markets*. Econ. J. 2009, 119, 158–171. [Google Scholar] [CrossRef] [Green Version]
  30. Barunik, J.; Krehlík, T. Measuring the Frequency Dynamics of Financial Connectedness and Systemic Risk. J. Financ. Econom. 2018, 16, 271–296. [Google Scholar] [CrossRef]
  31. Diebold, F.X.; Yilmaz, K. On the network topology of variance decompositions: Measuring the connectedness of financial firms. J. Econom. 2014, 182, 119–134. [Google Scholar] [CrossRef] [Green Version]
  32. Acharya, V.V.; Pedersen, L.H.; Philippon, T.; Richardson, M. Measuring Systemic Risk; Working Papers (Old Series) 1002; Federal Reserve Bank of Cleveland: Cleveland, OH, USA, 2010.
  33. Adrian, T.; Brunnermeier, M.K. CoVaR; Staff Reports 348; Federal Reserve Bank of New York: New York, NY, USA, 2008.
  34. Hu, D.; Zhao, J.L.; Hua, Z.; Wong, M.C.S. Network-based Modeling and Analysis of Systemic Risk in Banking Systems. MIS Q. 2012, 36, 1269–1291. [Google Scholar] [CrossRef] [Green Version]
  35. Huang, W.Q.; Zhuang, X.T.; Yao, S. A network analysis of the Chinese stock market. Phys. A Stat. Mech. Its Appl. 2009, 388, 2956–2964. [Google Scholar] [CrossRef]
  36. Shi, F.B.; Sun, X.Q.; Shen, H.W.; Cheng, X.Q. Detect colluded stock manipulation via clique in trading network. Phys. A Stat. Mech. Its Appl. 2019, 513, 565–571. [Google Scholar] [CrossRef]
  37. Boginski, V.; Butenko, S.; Pardalos, P.M. On Structural Properties of the Market Graph. Innov. Financ. Econ. Netw. 2003, 48, 29–35. [Google Scholar]
  38. Boginski, V.; Butenko, S.; Pardalos, P. Network models of massive datasets. Comput. Sci. Inf. Syst. 2004, 1, 75–89. [Google Scholar] [CrossRef]
  39. Xiao, F.; Liu, X.F.; Tse, C. Dynamics of Network of Global Stock Markets. Account. Financ. Res. 2012, 1. [Google Scholar] [CrossRef]
  40. Onnela, J.P.; Saramäki, J.; Kertész, J.; Kaski, K. Intensity and coherence of motifs in weighted complex networks. Phys. Rev. E 2005, 71, 065103. [Google Scholar] [CrossRef] [Green Version]
  41. Emmert-Streib, F.; Dehmer, M. Identifying critical financial networks of the DJIA: Toward a network-based index. Complexity 2010, 16, 24–33. [Google Scholar] [CrossRef]
  42. Boginski, V.; Butenko, S.; Pardalos, P. Mining market data: A network approach. Comput. Oper. Res. 2006, 33, 3171–3184. [Google Scholar] [CrossRef]
  43. Namaki, A.; Shirazi, A.; Raei, R.; Jafari, G. Network analysis of a financial market based on genuine correlation and threshold method. Phys. A Stat. Mech. Its Appl. 2011, 390, 3835–3841. [Google Scholar] [CrossRef]
  44. Lee, T.K.; Cho, J.H.; Kwon, D.S.; Sohn, S.Y. Global stock market investment strategies based on financial network indicators using machine learning techniques. Expert Syst. Appl. 2019, 117, 228–242. [Google Scholar] [CrossRef]
  45. Pagnottoni, P. Neural Network Models for Bitcoin Option Pricing. Front. Artif. Intell. 2019, 2, 5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Balash, V.; Chekmareva, A.; Faizliev, A.; Grigoriev, A.; Sidorov, S. Analysis of Financial Network Topological Dynamics of the Russian Stock Market from 2012 to 2019. J. Phys. Conf. Ser. 2020, 1564, 012030. [Google Scholar] [CrossRef]
  47. Lahaye, J.; Laurent, S.; Neely, C.J. Jumps, cojumps and macro announcements. J. Appl. Econom. 2011, 26, 893–921. [Google Scholar] [CrossRef] [Green Version]
  48. Sidorov, S.; Revutskiy, A.; Faizliev, A.; Korobov, E.; Balash, V. Stock Volatility Modelling with Augmented GARCH Model with Jumps. IAENG Int. J. Appl. Math. 2014, 44, 212–220. [Google Scholar]
  49. Hu, M.; Zhang, D.; Ji, Q.; Wei, L. Macro factors and the realized volatility of commodities: A dynamic network analysis. Resour. Policy 2020, 68, 101813. [Google Scholar] [CrossRef]
  50. Spelta, A.; Pecora, N.; Flori, A.; Giudici, P. The impact of the SARS-CoV-2 pandemic on financial markets: A seismologic approach. Ann. Oper. Res. 2021. [Google Scholar] [CrossRef]
  51. Zhang, D.; Hu, M.; Ji, Q. Financial markets under the global pandemic of COVID-19. Financ. Res. Lett. 2020, 36, 101528. [Google Scholar] [CrossRef]
  52. Corbet, S.; Hou, Y.G.; Hu, Y.; Oxley, L.; Xu, D. Pandemic-related financial market volatility spillovers: Evidence from the Chinese COVID-19 epicentre. Int. Rev. Econ. Financ. 2021, 71, 55–81. [Google Scholar] [CrossRef]
  53. Corbet, S.; Larkin, C.; Lucey, B. The contagion effects of the COVID-19 pandemic: Evidence from gold and cryptocurrencies. Financ. Res. Lett. 2020, 35, 101554. [Google Scholar] [CrossRef]
  54. Albulescu, C.T. COVID-19 and the United States financial markets’ volatility. Financ. Res. Lett. 2021, 38, 101699. [Google Scholar] [CrossRef]
  55. Koop, G.; Pesaran, M.; Potter, S.M. Impulse response analysis in nonlinear multivariate models. J. Econom. 1996, 74, 119–147. [Google Scholar] [CrossRef]
  56. Pesaran, H.; Shin, Y. Generalized impulse response analysis in linear multivariate models. Econ. Lett. 1998, 58, 17–29. [Google Scholar] [CrossRef]
  57. Giudici, P.; Pagnottoni, P. Vector error correction models to measure connectedness of Bitcoin exchange markets. Appl. Stoch. Model. Bus. Ind. 2020, 36, 95–109. [Google Scholar] [CrossRef] [Green Version]
  58. Antonakakis, N.; Gabauer, D.; Gupta, R.; Plakandaras, V. Dynamic connectedness of uncertainty across developed economies: A time-varying approach. Econ. Lett. 2018, 166, 63–75. [Google Scholar] [CrossRef] [Green Version]
  59. Gabauer, D.; Gupta, R. On the transmission mechanism of country-specific and international economic uncertainty spillovers: Evidence from a TVP-VAR connectedness decomposition approach. Econ. Lett. 2018, 171, 63–71. [Google Scholar] [CrossRef] [Green Version]
  60. Gamba-Santamaria, S.; Gomez-Gonzalez, J.E.; Hurtado-Guarin, J.L.; Melo-Velandia, L.F. Stock market volatility spillovers: Evidence for Latin America. Financ. Res. Lett. 2017, 20, 207–216. [Google Scholar] [CrossRef] [Green Version]
  61. He, F.; Liu, Z.; Chen, S. Industries Return and Volatility Spillover in Chinese Stock Market: An Early Warning Signal of Systemic Risk. IEEE Access 2019, 7, 9046–9056. [Google Scholar] [CrossRef]
  62. Yilmaz, K. Return and volatility spillovers among the East Asian equity markets. J. Asian Econ. 2010, 21, 304–313. [Google Scholar] [CrossRef] [Green Version]
  63. Alter, A.; Beyer, A. The dynamics of spillover effects during the European sovereign debt turmoil. J. Bank. Financ. 2014, 42, 134–153. [Google Scholar] [CrossRef] [Green Version]
  64. Akhtaruzzaman, M.; Boubaker, S.; Sensoy, A. Financial contagion during COVID–19 crisis. Financ. Res. Lett. 2021, 38, 101604. [Google Scholar] [CrossRef] [PubMed]
  65. Zhou, X.; Zhang, W.; Zhang, J. Volatility spillovers between the Chinese and world equity markets. Pac.-Basin Financ. J. 2012, 20, 247–270. [Google Scholar] [CrossRef]
  66. Costola, M.; Lorusso, M. Spillovers among Energy Commodities and the Russian Stock Market; MPRA Paper 108990; University Library of Munich: Munich, Germany, 2021. [Google Scholar]
  67. Nicholson, W.B.; Matteson, D.; Bien, J. BigVAR: Tools for Modeling Sparse High-Dimensional Multivariate Time Series. arXiv 2017, arXiv:1702.07094. [Google Scholar]
  68. Davis, R.A.; Zang, P.; Zheng, T. Sparse Vector Autoregressive Modeling. J. Comput. Graph. Stat. 2016, 25, 1077–1096. [Google Scholar] [CrossRef] [Green Version]
  69. Basu, S.; Michailidis, G. Regularized estimation in sparse high-dimensional time series models. Ann. Stat. 2015, 43, 1535–1567. [Google Scholar] [CrossRef] [Green Version]
  70. Souto, M.; Moreira, A.; Veiga, A.; Street, A.; Garcia, J.D.; Epprecht, C. A high-dimensional VARX model to simulate monthly renewable energy supply. In Proceedings of the 2014 Power Systems Computation Conference, Wrocław, Poland, 18–22 August 2014; pp. 1–7. [Google Scholar] [CrossRef]
  71. Hansen, P.R.; Lunde, A. A Realized Variance for the Whole Day Based on Intermittent High-Frequency Data. J. Financ. Econom. 2005, 3, 525–554. [Google Scholar] [CrossRef]
  72. Hansen, P.; Lunde, A. An Optimal and Unbiased Measure of Realized Variance Based on Intermittent High-Frequency Data; Department of Economics, Stanford University: Stanford, CA, USA, 2003. [Google Scholar]
  73. Taylor, N. A note on the importance of overnight information in risk management models. J. Bank. Financ. 2007, 31, 161–180. [Google Scholar] [CrossRef]
  74. Ahoniemi, K.; Lanne, M. Overnight stock returns and realized volatility. Int. J. Forecast. 2013, 29, 592–604. [Google Scholar] [CrossRef]
  75. Oldfield, G., Jr.; Rogalski, R.J. A Theory of Common Stock Returns Over Trading and Non-Trading Periods. J. Financ. 1980, 35, 729–751. [Google Scholar] [CrossRef]
  76. Tsiakas, I. Overnight information and stochastic volatility: A study of European and US stock exchanges. J. Bank. Financ. 2008, 32, 251–268. [Google Scholar] [CrossRef]
  77. Aslam, F.; Ferreira, P.; Mughal, K.S.; Bashir, B. Intraday Volatility Spillovers among European Financial Markets during COVID-19. Int. J. Financ. Stud. 2021, 9, 5. [Google Scholar] [CrossRef]
  78. Mikhaylov, A. Volatility Spillover Effect between Stock and Exchange Rate in Oil Exporting Countries. Int. J. Energy Econ. Policy 2018, 8, 321–326. [Google Scholar]
  79. Pavlova, I.; de Boyrie, M.E.; Parhizgari, A.M. A dynamic spillover analysis of crude oil effects on the sovereign credit risk of exporting countries. Q. Rev. Econ. Financ. 2018, 68, 10–22. [Google Scholar] [CrossRef]
  80. Mensi, W.; Shafiullah, M.; Vo, X.V.; Kang, S.H. Volatility spillovers between strategic commodity futures and stock markets and portfolio implications: Evidence from developed and emerging economies. Resour. Policy 2021, 71, 102002. [Google Scholar] [CrossRef]
  81. Shi, K. Spillovers of Stock Markets among the BRICS: New Evidence in Time and Frequency Domains before the Outbreak of COVID-19 Pandemic. J. Risk Financ. Manag. 2021, 14, 112. [Google Scholar] [CrossRef]
  82. Mensi, W.; Maitra, D.; Vo, X.V.; Kang, S.H. Asymmetric volatility connectedness among main international stock markets: A high frequency analysis. Borsa Istanb. Rev. 2021, 21, 291–306. [Google Scholar] [CrossRef]
Figure 1. Realized and overnight volatility for SBER.
Figure 1. Realized and overnight volatility for SBER.
Mathematics 09 02484 g001
Figure 2. The value of the R 2 for regression models.
Figure 2. The value of the R 2 for regression models.
Mathematics 09 02484 g002
Figure 3. Dynamic frequency connectedness. DY-BK approach.
Figure 3. Dynamic frequency connectedness. DY-BK approach.
Mathematics 09 02484 g003
Figure 4. Dynamic frequency conditional connectedness. Two step approach.
Figure 4. Dynamic frequency conditional connectedness. Two step approach.
Mathematics 09 02484 g004
Figure 5. Dynamic frequency conditional connectedness. VARX approach.
Figure 5. Dynamic frequency conditional connectedness. VARX approach.
Mathematics 09 02484 g005
Figure 6. Net volatility spillover for LKOH and AFKS.
Figure 6. Net volatility spillover for LKOH and AFKS.
Mathematics 09 02484 g006
Figure 7. Net volatility spillover for SBER and PIKK.
Figure 7. Net volatility spillover for SBER and PIKK.
Mathematics 09 02484 g007
Figure 8. Dynamic frequency connectedness during period of crisis 2020.
Figure 8. Dynamic frequency connectedness during period of crisis 2020.
Mathematics 09 02484 g008
Figure 9. Dynamic frequency connectedness during 2014–2016.
Figure 9. Dynamic frequency connectedness during 2014–2016.
Mathematics 09 02484 g009
Figure 10. Overall Directional Connectedness.
Figure 10. Overall Directional Connectedness.
Mathematics 09 02484 g010
Figure 11. Frequency decomposition of pairwise volatility spillovers networks (window: from 1 March 2015 to 1 March 2017).
Figure 11. Frequency decomposition of pairwise volatility spillovers networks (window: from 1 March 2015 to 1 March 2017).
Mathematics 09 02484 g011
Figure 12. Net Pairwise Directional Connectedness by pandemic crisis period (from 6 January 2020 to 6 January 2021).
Figure 12. Net Pairwise Directional Connectedness by pandemic crisis period (from 6 January 2020 to 6 January 2021).
Mathematics 09 02484 g012
Figure 13. Net Pairwise Directional Connectedness by politic crisis period (from 1 January 2014 to 1 March 2015).
Figure 13. Net Pairwise Directional Connectedness by politic crisis period (from 1 January 2014 to 1 March 2015).
Mathematics 09 02484 g013
Table 1. Companies with tickers classified by industry.
Table 1. Companies with tickers classified by industry.
SectorTickerName
ALRSAlrosa
CHMFSeverstal
GMKNNornickel
Basic MaterialsMAGNMagnitogorsk Iron and Steel Works
NLMKNLMK Group
PHORCompany PhosAgro
PLZLPolyus
AFKSAFK Sistema
Communications ServicesMTSSMTS
RTKMRostelecom
Consumer DefensiveMGNTMagnit
and Consumer CyclicalMVIDM. Video
GAZPGazprom
LKOHLukoil
NVTKNovatek
Energy (Oil & Gas )ROSNRosneft
SNGSSurgutneftegas
TATNTatneft
TRNFPTransneft
Financial ServicesSBERSberbank of Russia
VTBRVTB Bank
Industrials (Airlines)AFLTAeroflot
Real EstateLSRGLSR Group
PIKKPIK Group
FEESFGC UES
UtilitiesHYDRRusHydro
IRAOInter RAO
Table 2. Full-sample connectedness table.
Table 2. Full-sample connectedness table.
AFKSAFLTALRSCHMFFEESGAZPGMKNHYDRIRAOLKOHLSRGMAGNMGNTMTSSMVIDNLMKNVTKPHORPIKKPLZLROSNRTKMSBERSNGSTATNTRNFPVTBRFROM
AFKS21.012.631.931.963.033.292.574.002.604.002.062.793.467.551.973.144.062.511.871.363.542.853.963.543.031.863.442.93
AFLT1.1815.022.502.964.233.503.304.182.854.802.163.643.602.221.914.245.231.901.521.114.473.935.093.694.441.944.403.15
ALRS1.093.3115.933.073.883.253.343.703.964.382.793.643.372.201.964.004.462.172.091.703.943.624.093.584.092.553.853.11
CHMF0.863.412.6013.95.333.463.404.113.293.842.895.253.342.281.735.774.182.081.741.583.523.544.543.323.602.014.433.19
FEES0.943.672.343.5014.983.332.985.735.263.762.674.603.562.251.854.534.132.411.810.903.403.934.123.663.332.334.043.15
GAZP1.152.892.072.763.4513.353.393.592.716.341.592.953.772.401.164.046.091.471.110.865.863.666.025.035.212.284.793.21
GMKN1.183.232.813.103.394.1211.953.902.935.811.513.253.782.401.714.305.631.751.351.315.073.365.324.805.362.414.273.26
HYDR1.363.452.643.115.543.543.2813.083.664.762.153.633.433.021.614.424.931.801.651.074.183.644.664.074.272.624.443.22
IRAO1.652.782.572.655.803.153.054.4817.054.211.814.073.412.491.953.814.232.182.231.233.163.654.203.303.483.404.003.07
LKOH1.163.162.362.483.005.013.753.842.9910.881.493.003.802.321.013.806.861.301.221.186.623.355.565.686.543.084.523.30
LSRG1.242.842.473.073.893.112.213.512.883.6719.933.263.642.062.974.104.313.262.801.723.533.583.723.593.631.493.552.97
MAGN1.193.312.584.275.203.303.104.033.764.112.3614.223.642.361.635.743.882.001.572.083.503.594.803.253.632.464.453.18
MGNT1.383.022.312.733.664.003.393.583.164.982.353.6414.942.392.214.475.302.071.520.974.243.274.654.424.402.464.493.15
MTSS4.792.612.032.663.163.602.944.442.564.271.513.053.3819.211.903.584.692.301.971.123.533.074.603.573.312.223.912.99
MVID1.142.421.872.403.252.441.892.612.612.764.393.253.722.0830.613.543.743.091.631.302.713.073.082.823.071.253.252.57
NLMK1.003.472.564.524.843.753.504.263.164.352.505.293.882.301.6512.094.782.131.431.483.953.515.053.734.022.154.633.26
NVTK1.173.392.342.713.334.553.533.943.016.362.013.073.962.521.774.1710.751.691.421.175.823.595.275.145.952.864.533.31
PHOR1.092.732.262.614.012.652.443.003.262.974.593.183.402.624.084.033.6523.742.601.632.513.723.582.792.431.442.972.82
PIKK0.842.612.412.683.522.402.253.313.633.123.832.853.032.282.053.023.603.0827.892.022.903.412.872.933.091.692.692.67
PLZL0.872.302.642.982.512.302.872.172.463.012.173.982.401.921.633.392.982.402.2136.282.402.373.092.192.271.093.122.36
ROSN1.133.242.442.522.975.103.783.722.517.401.602.833.622.131.113.816.821.191.220.9911.443.215.605.426.772.784.653.28
RTKM0.983.602.562.954.133.963.104.003.244.602.643.603.442.311.994.094.702.361.871.073.9815.644.703.844.282.154.203.12
SBER1.233.632.263.043.375.004.073.933.245.831.583.683.662.661.374.555.751.681.171.265.293.5311.494.014.542.555.613.28
SNGS1.152.662.032.453.334.733.383.842.516.621.942.734.082.291.553.836.401.501.310.935.703.274.5714.285.622.984.323.17
TATN1.053.112.672.673.104.563.563.762.796.721.803.123.792.101.433.946.501.281.371.096.263.484.915.1712.253.104.433.25
TRNFP1.142.432.662.573.863.342.924.083.785.361.483.783.552.451.403.685.242.051.591.444.303.164.534.404.8316.083.893.11
VTBR1.323.442.363.233.734.333.284.123.705.122.093.773.882.531.694.525.211.661.321.374.733.556.024.304.512.3911.833.27
TO1.232.942.312.803.683.553.013.703.064.562.223.403.432.451.753.954.721.971.611.264.043.294.393.794.062.203.9683.34
NET−1.69−0.21−0.81−0.390.540.34−0.250.48−0.021.26−0.740.230.28−0.54−0.820.691.41−0.85−1.06−1.100.760.171.110.610.81−0.900.69
TO (0;5]0.300.800.620.780.951.050.881.040.781.370.380.910.930.670.241.131.350.460.350.221.220.901.311.081.170.641.13
TO (6;20]0.360.890.690.801.101.070.881.100.911.380.671.001.020.720.541.181.440.570.450.351.220.981.291.131.250.661.16
TO (20;100]0.581.241.011.121.651.441.181.561.401.821.191.441.451.080.941.651.930.960.780.701.601.411.751.571.650.921.58
NET (0;5]−0.14−0.13−0.12−0.170.15−0.250.060.210.380.29−0.13−0.01−0.21−0.02−0.260.000.34−0.15−0.45−0.200.130.030.200.140.270.010.03
NET (6;20]−0.52−0.06−0.22−0.130.170.12−0.190.110.220.35−0.240.060.12−0.16−0.250.250.42−0.29−0.32−0.310.210.020.310.100.25−0.270.26
NET (20;100]−1.02−0.01−0.45−0.170.240.49−0.160.16−0.580.63−0.350.140.35−0.35−0.320.460.66−0.39−0.29−0.580.430.140.570.370.31−0.620.34
Table 3. The decomposition of connectedness.
Table 3. The decomposition of connectedness.
Spillover “From”
DY-BKTwo StepVARX
TickerAll123All123All123
AFKS2.9215.0%30.4%54.6%2.8821.3%35.6%43.1%2.8325.5%40.2%34.3%
AFLT3.1429.6%30.5%39.9%3.0039.4%32.2%28.4%3.1039.5%34.9%25.6%
ALRS3.1023.8%29.4%46.8%2.9835.9%33.0%31.1%3.0631.2%37.2%31.6%
CHMF3.1629.9%29.3%40.8%3.0741.4%31.7%27.0%3.1536.8%35.8%27.4%
FEES3.1425.4%29.6%45.0%3.1033.1%33.8%33.1%3.1432.5%37.0%30.6%
GAZP3.2040.8%29.7%29.5%3.1151.7%28.7%19.6%3.1749.0%31.0%20.1%
GMKN3.2425.3%33.2%41.5%3.1040.5%32.9%26.6%3.1736.8%36.8%26.4%
HYDR3.2125.8%30.7%43.4%3.1536.2%33.9%29.9%3.2033.1%37.4%29.6%
IRAO3.0713.0%22.4%64.6%2.8827.8%29.6%42.6%2.9722.9%36.1%40.9%
LKOH3.2932.7%31.1%36.2%3.2351.4%28.4%20.2%3.2846.0%32.4%21.7%
LSRG2.9617.3%30.8%51.9%2.8628.5%35.4%36.1%2.9725.0%40.1%34.9%
MAGN3.1629.0%29.7%41.3%3.0840.8%31.7%27.5%3.1436.9%36.0%27.1%
MGNT3.1436.3%28.6%35.1%3.0844.8%29.9%25.3%3.1343.0%32.6%24.4%
MTSS2.9922.8%29.2%48.0%2.9227.5%33.1%39.4%2.9831.1%37.9%30.9%
MVID2.5519.6%30.8%49.6%2.4628.9%35.4%35.8%2.5425.4%39.3%35.3%
NLMK3.2534.6%28.7%36.7%3.1845.2%30.3%24.6%3.2440.8%34.1%25.1%
NVTK3.3030.7%30.9%38.5%3.2543.8%31.2%25.0%3.2940.6%34.5%24.9%
PHOR2.8221.6%30.3%48.0%2.8127.2%35.3%37.5%2.8325.4%38.9%35.7%
PIKK2.6430.2%29.3%40.5%2.5646.2%29.9%23.9%2.6138.9%34.5%26.6%
PLZL2.3618.1%27.7%54.2%2.0440.8%32.5%26.7%2.2730.2%38.0%31.9%
ROSN3.2733.1%31.0%35.9%3.1951.0%28.5%20.6%3.2546.5%32.2%21.3%
RTKM3.1228.0%31.0%41.0%3.0336.5%34.0%29.5%3.1136.1%36.4%27.5%
SBER3.2733.9%30.1%36.0%3.2045.2%30.7%24.0%3.2542.7%33.6%23.7%
SNGS3.1729.7%32.5%37.7%3.1141.4%34.0%24.6%3.1639.1%35.8%25.1%
TATN3.2427.7%30.7%41.6%3.1845.0%30.3%24.7%3.2441.2%34.3%24.5%
TRNFP3.1120.4%29.8%49.8%2.9635.3%33.4%31.3%3.0332.8%37.1%30.1%
VTBR3.2533.9%27.7%38.4%3.1744.4%29.8%25.8%3.2242.3%32.8%24.9%
Spillover “To”
DY-BKTwo stepVARX
Tickerall123all123all123
AFKS1.2424.2%29.2%46.6%1.4934.2%32.9%32.9%1.3533.3%36.6%30.1%
AFLT2.9327.1%30.4%42.4%2.6937.8%33.1%29.1%2.9335.8%36.1%28.1%
ALRS2.3226.6%30.0%43.4%2.3538.8%32.2%29.0%2.3934.5%36.5%28.9%
CHMF2.7028.8%29.7%41.6%2.8440.0%31.9%28.2%2.9037.0%35.5%27.5%
FEES3.7025.6%29.7%44.7%3.7636.9%32.7%30.4%3.7033.8%36.4%29.8%
GAZP3.5629.7%30.0%40.3%3.2142.3%31.2%26.4%3.3839.5%34.7%25.8%
GMKN2.9430.0%29.9%40.2%2.7642.5%31.1%26.4%3.0138.4%35.1%26.5%
HYDR3.6928.1%29.7%42.2%3.7739.4%31.9%28.8%3.7236.9%35.4%27.6%
IRAO3.0925.2%29.4%45.5%3.0237.3%32.4%30.3%3.1032.6%36.1%31.3%
LKOH4.5629.9%30.1%39.9%4.2642.8%31.2%26.1%4.4839.7%34.6%25.6%
LSRG2.2417.0%30.0%53.0%2.2627.9%34.9%37.3%2.1127.2%39.0%33.8%
MAGN3.3627.2%29.8%43.0%3.3038.7%32.3%29.0%3.3735.6%36.0%28.4%
MGNT3.4127.4%30.0%42.6%3.3139.3%32.1%28.6%3.3236.4%35.7%27.9%
MTSS2.4727.0%29.0%44.0%2.5836.6%31.6%31.8%2.5335.8%35.7%28.6%
MVID1.7213.8%31.2%55.0%1.3827.6%35.8%36.7%1.6421.2%41.2%37.6%
NLMK3.9628.5%29.8%41.7%3.8639.9%31.9%28.2%3.9737.1%35.4%27.4%
NVTK4.7128.6%30.5%41.0%4.4841.2%31.7%27.1%4.5938.5%35.1%26.4%
PHOR1.9923.1%28.5%48.4%2.1133.7%32.8%33.5%1.9030.5%36.9%32.6%
PIKK1.5822.1%28.5%49.4%1.6734.8%32.2%33.0%1.4732.1%36.7%31.2%
Spillover “From”
DY-BKTwo StepVARX
TickerAll123All123All123
PLZL1.2717.6%27.3%55.1%1.1232.8%31.6%35.7%1.2027.3%38.3%34.4%
ROSN4.0430.1%30.2%39.7%3.8042.8%31.2%26.0%3.9639.9%34.7%25.4%
RTKM3.3027.3%29.8%42.9%3.2338.4%32.1%29.5%3.3835.8%35.8%28.4%
SBER4.3430.0%29.8%40.2%4.0542.1%31.3%26.6%4.2139.2%34.7%26.1%
SNGS3.7928.6%29.9%41.5%3.5341.5%31.3%27.1%3.6738.6%35.0%26.4%
TATN4.0728.7%30.7%40.7%3.8441.8%31.5%26.8%4.0038.6%35.2%26.2%
TRNFP2.2229.0%29.5%41.5%2.1741.2%31.3%27.5%2.1838.8%34.8%26.3%
VTBR3.8829.2%29.9%40.9%3.7441.1%31.6%27.3%3.9038.0%35.2%26.8%
Table 4. Overall and long term NET spillovers.
Table 4. Overall and long term NET spillovers.
OverallLong term
TickerVARTwo StepVARXTickerVARTwo StepVARX
AFKS−1.68−1.38−1.48AFKS−1.02−0.75−0.56
PIKK−1.06−0.89−1.14PHOR−0.39−0.35−0.39
PLZL−1.09−0.93−1.07TRNFP−0.62−0.33−0.34
PHOR−0.83−0.69−0.92LSRG−0.35−0.19−0.32
MVID−0.83−1.08−0.9PLZL−0.58−0.15−0.31
LSRG−0.72−0.6−0.86MVID−0.32−0.37−0.28
TRNFP−0.88−0.79−0.85ALRS−0.45−0.25−0.27
ALRS−0.79−0.63−0.67IRAO−0.58−0.31−0.24
MTSS−0.52−0.34−0.45PIKK−0.29−0.06−0.23
CHMF−0.46−0.23−0.25MTSS−0.35−0.33−0.2
AFLT−0.2−0.31−0.17CHMF−0.17−0.03−0.07
GMKN−0.29−0.34−0.17GMKN−0.16−0.09−0.04
IRAO0.020.140.13AFLT−0.01−0.070.03
MGNT0.260.230.19HYDR0.160.140.08
GAZP0.350.10.21MAGN0.140.110.1
MAGN0.20.230.23RTKM0.140.060.11
RTKM0.180.20.27FEES0.240.110.14
SNGS0.620.420.51MGNT0.350.170.16
HYDR0.480.620.52SNGS0.370.190.17
FEES0.560.660.56GAZP0.490.240.24
VTBR0.630.570.68VTBR0.340.20.25
ROSN0.770.610.71TATN0.310.240.25
NLMK0.710.680.73NLMK0.460.310.27
TATN0.820.650.76ROSN0.430.330.31
SBER1.070.840.95SBER0.570.30.33
LKOH1.271.031.2NVTK0.660.40.39
NVTK1.421.231.3LKOH0.630.460.44
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Balash, V.; Faizliev, A.; Sidorov, S.; Chistopolskaya, E. Conditional Time-Varying General Dynamic Factor Models and Its Application to the Measurement of Volatility Spillovers across Russian Assets. Mathematics 2021, 9, 2484. https://doi.org/10.3390/math9192484

AMA Style

Balash V, Faizliev A, Sidorov S, Chistopolskaya E. Conditional Time-Varying General Dynamic Factor Models and Its Application to the Measurement of Volatility Spillovers across Russian Assets. Mathematics. 2021; 9(19):2484. https://doi.org/10.3390/math9192484

Chicago/Turabian Style

Balash, Vladimir, Alexey Faizliev, Sergei Sidorov, and Elena Chistopolskaya. 2021. "Conditional Time-Varying General Dynamic Factor Models and Its Application to the Measurement of Volatility Spillovers across Russian Assets" Mathematics 9, no. 19: 2484. https://doi.org/10.3390/math9192484

APA Style

Balash, V., Faizliev, A., Sidorov, S., & Chistopolskaya, E. (2021). Conditional Time-Varying General Dynamic Factor Models and Its Application to the Measurement of Volatility Spillovers across Russian Assets. Mathematics, 9(19), 2484. https://doi.org/10.3390/math9192484

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

Article Metrics

Back to TopTop