ORIGINAL RESEARCH article

Front. Phys., 18 November 2024

Sec. Social Physics

Volume 12 - 2024 | https://doi.org/10.3389/fphy.2024.1484589

Multi-quantile systemic financial risk based on a monotone composite quantile regression neural network

  • 1. School of Financial Technology, Anhui Business College, Wuhu, China

  • 2. School of Economics and Management, Southeast University, Nanjing, China

Abstract

This study proposes a novel perspective to calibrate the conditional value at risk (CoVaR) of countries based on the monotone composite quantile regression neural network (MCQRNN). MCQRNN can fix the “quantile crossing” problem, which is more robust in CoVaR estimating. In addition, we extend the MCQRNN method with quantile-on-quantile (QQ), which can avoid the bias in quantile regression. Building on the estimation results, we construct a systemic risk spillover network across countries in the Asia–Pacific region by considering the suffering and overflow effects. A comparison among MCQRNN, QRNN, and MCQRNN-QQ indicates the significance of monotone composite quantiles in modeling CoVaR. Additionally, the network analysis of composite risk spillovers illustrates the advantages of MCQRNN-QQ-CoVaR compared with QRNN-CoVaR. Moreover, the average composite systemic suffering index and the average composite systemic overflow index are introduced as country-specific measures that enable identifying systemically relevant countries during extreme events.

1 Introduction

The Sino–U.S. trade war and the COVID-19 epidemic have caused huge fluctuations in Asia–Pacific stock markets. Compared with other economic organizations, the Asia–Pacific Economic Cooperation (APEC) organization provides a diversified financial markets environment, including developed and developing countries. In addition, APEC’s organizational structure and cooperation mechanism are more flexible, which means member countries cooperate while maintaining autonomy. Moreover, the economic structure of APEC countries is highly complementary; for example, resource-rich countries tend to trade closely with countries with developed manufacturing industries. By establishing interconnectivity, APEC encourages deeper cooperation in infrastructure, trade, and investment among countries in the region. According to statistical data, APEC members account for more than 40% of global trade. Within the region, trade among members is higher than trade with non-members. Despite the large volume of intra-APEC trade, APEC trade relations may depend more on bilateral relationships of large countries such as China and the United States than other economic organizations such as the European Union. This means that the trade closeness of APEC is greatly affected by the policy changes of large countries. Consequently, in the stage of Sino–U.S. trade friction, the trade cooperation of member countries will undergo great changes.

In the COVID-19 phase, the economic conditions of China and the United States will directly affect the risk level of the organization’s members. Therefore, the research on systemic financial risks among APEC member countries in this paper is helpful for a deeper analysis of the risk contagion mechanisms between different economies and could provide a supplement to existing literature. Growing uncertainty results in countries facing cross-border risk shocks, making the issue of systemic risk a renewed focus of research by academics and regulators. Systemic risk caused by the bankruptcy of systemically important economies is primarily the failure of the financial system. From the aspect of international markets, when an important node is damaged by a shock, other markets may also be affected and could eventually be contagious to the entire financial system.

To measure the systemic risk, studies in recent years have begun to focus on the risk contagion or spillover effect []. The former is mainly from a theoretical modeling perspective []. Additionally, more focus on empirical measurement provides compositions for the “edge” of the financial network. For instance, Engle (2002) constructs the GARCH-DCC method to capture the risk spillover among market indexes []. Rodriguez (2007) measures the interdependence among East Asian stock via switching-parameter copulas []. Billio et al. (2012) prompt the empirical framework of a risk spillover network based on Granger causality []. Diebold and Yılmaz (2014) involve variance decomposition in risk spillover and analyze the vulnerability by financial networks []. Baruník and Křehlík (2018) study the risk spillover from aspects of heterogeneous frequency responses to shock [].

The risk modeling system of this paper is an addition to the conditional value at risk (), which is the systemic risk approach []. Adrian and Brunnermeier (2016) define as the change of system’s value at risk in the condition of one single institution’s loss, which has provided a new perspective for risk spillover effects []. Nevertheless, the original needs to assume a linear relationship between the return of market indexes and institutions’ stock prices. Thus, Hautsch et al. (2015) provide a paradigm of multivariate , which is based on the marginal effect of risk spillovers among financial institutions []. Furthermore, based on , Fan et al.(2018), Härdle et al. (2016) propose a tail event driven network technique (), where the is replaced as the partial differentiation of multivariate nonlinear [, ]. On this basis, Keilbar and Wang (2022) adapt the model approach based on a neural network method, which also uses partial differential to calculate the marginal effect among agents []. has already been employed for risk spillover among inter-industries and energy markets [, ]. Moreover, graph learning in attributed networks are used in risk spillover by different node-to-cluster distance functions [, ].

This paper seeks to expand the research perspective on systemic financial risk by examining the composite risk spillover effects among financial markets to avoid possible errors in the setting of quantiles. We construct the composite risk spillover measure based on the multi-quantile by multi-quantile () and quantile-on-quantile regression (). has been widely involved in correlation and spillover effect [], although it has not been adopted by estimation. It is found that the concept of is suitable for systemic financial risk because there are two-sided quantile sets in the definition. However, in traditional methodology, both the stand-alone quantiles and exposed quantiles are set as a fixed small number (normally 5%). If is extended to multiple (stand-alone) quantiles, a tough problem, “quantile crossing,” will emerge and raise a paradox of “higher risk but less1 loss” [, ]. The problem of non-monotonicity of risk indicators arises when estimating CoVaR using single-quantile regression. However, financial risk indicators must ensure their monotonicity, so the issue of quantile crossover must be addressed [, ].

Some studies focus on this problem. Acharya et al. (2017) assessed the expected loss below one quantile in dealing with the problem of quantile crossing []. Catania and Luati (2023) used a semiparametric model to satisfy the condition of non-crossing quantiles []. By investigating the , this paper not only obtains a more robust risk spillover measure, but also facilitates the examination of the characteristics of nonlinear systemic risks in each market. In addition, although some studies have addressed the risk-related networks of financial markets in the Asia–Pacific region [], few visualize the systemic risk of quantile regression neural networks. The findings of this paper should make an important contribution to the field of capturing systemic risk spillovers among financial markets and recognition of risk sources.

This paper proposes a quantile-on-quantile regression to examine the two-sided quantile in CoVaR estimation. The effects of systemic risk are analyzed by three-dimensional surface plots in empirical research. We extend the quantile regression in the systemic risk approach with a monotone composite quantile regression neural network, which can not only be suitable for solving the nonlinear issue but also optimize the “quantile crossing” problem. Moreover, we introduce the composite systemic suffering indicator, the composite systemic overflow indicator, and the total composite overflow indicator as three country-specific measures to identify systemically relevant countries in the Asia–Pacific region during extreme events.

This article makes three main contributions. First, we estimate the systemic financial risk through a new perspective that adopts 3D surface plots. Second, multi-quantiles are adopted in the model to capture the multi-state characteristics of risk. Third, MCQRNN is used to relieve the quantile crossing problem.

The remainder of the paper will be organized as follows: Section 2 will (i) introduce the multi-quantile based on monotone composite neural network quantile regression () and (ii) describe the methodology of in details. After constituting the risk modeling system step by step, the empirical results based on Asia–Pacific stock markets and discussion will be presented in Section 3. Finally, a conclusion of this paper and suggestions for future study are drawn in Section 4.

2 Materials and methods

2.1 The monotone composite neural network of quantile regression () method

A quantile regression model based on a linear regression equation estimates parameters for quantiles under the variable by introducing an indicative function in the loss function [], defined as Equation 1.where is the loss function at the quantile level of (known as the pinball loss function). Where is the indicative function, the value is when the independent variable ; otherwise, the value is . However, this model only considers linear relationships between the variables, which cannot state the effect of non-linearity. For this reason, Taylor (2000) involved a neural network model called the quantile regression neural network approach () []. Cannon (2018) extended the to the monotone composite neural network of quantile regression (), which can mitigate the “quantile crossing” problem []. The comprehensive estimations of multi-quantile can be obtained by adjusted . Assuming the number of market indexes is and their price returns are , the conditional value at risk of market in level of quantile can be obtained by defined as2Equation 2.

The difference between and is introducing quantile as an input variable with a positive weight of . In addition, the nonlinear activation function is assumed to be invariant and known. The parameter of each node in the hidden layer consists of weights and the intercept , while the output layer parameters are and . These parameters are verified to be consistent and asymptotically normal under large sample and regularity conditions. Moreover, they converge to the true function at a certain rate [, ].

However, the loss function in the quantile regression is not differentiable everywhere. It limits the use of artificial neural networks’ regular algorithms (ANNs) in . Therefore, it is necessary to adjust the form of the loss function. One approach is to add a Huber norm [], which allows a smoothing approximation to the error term near the origin. Furthermore, is a hybrid , which makes it possible to use standard gradient-based optimization algorithms.

In addition, when the capacity of the neural network is large, it is prone to over-fitting problems. Choosing a modest neural network structure and hyperparameters is an effective approach often used in machine learning. In a single hidden layer network, the most important hyperparameter is the number of hidden nodes . Therefore, choosing the appropriate number of nodes can reduce the capacity of the neural network and avoid the over-fitting phenomenon. In addition, Bishop (1995) proposed alleviating this problem by introducing weight decay regularization []. Such regularization requires adding an additional penalty term to the weight parameter . Referring to model of Cannon [, ], the final estimator is set as Equation 3 included a quadratic penalty term.where is the parameter that regulates the weight of the quadratic penalty term in the loss function. When is , this regression is transformed to an ordinary . In this paper, simple sampling is used to train neural networks. According to the experience of [] and [], we selected 50% of the samples from the sample period as the training set to train the loss function, which is Equation 3, in the MCQRNN model. Except for and hyperparameters, the intercept and weights of the neural network are trained. Otherwise, too large will lead to the loss of non-linear characteristics of the model, when the transfer function is the sigmoidal hidden layer transfer function, such as hyperbolic tangent . To balance the degree of the over-fitting and the prediction accuracy, is set to equal to here. Furthermore, the optimization of the loss function as shown in Equation 3 can be achieved by using a quasi-Newton optimization algorithm, which is less complex and more appropriate for the computational complexity in this paper. The quantile regression neural network process is visualized in Figure 1.

FIGURE 1

2.2 Calibrate systemic risk system

The calibration details of

are explained in this section. There are four steps involved in the systemic risk system calibration. The first step is the estimation of

based on

. Next, the results are used to estimate

with

and

for each country. In the third step, the composite risk spillover effects are calculated by resulting in an extreme risk spillover measure. Finally, the systemic risk measures are proposed based on the systemic risk network. The process of the risk modeling system is demonstrated in

Figure 2

as follows:

  • Step 1: Estimation of value at risk with

FIGURE 2

Because it is challenging to select common macro-state variables for all indexes, the linear quantile regression of value at risk is no longer suitable for measuring the tail risk of stock markets. Hence,

is adopted in the first step [

]. Given the asymmetric effects of the rise and fall of each index,

at the quantile level of

can be estimated in the following

Equation 4

[

]

3

:

where

is the absolute value of market index

’s return when the lagged log-return is larger than zero, and the rest is

.

is the absolute value when the lagged log-return is minus zero.

is the vector of the parameters. By utilizing the differential evolution algorithm [

] and the loss function similar to

Equation 1

, the

is an appropriate risk indicator for individual stock market

at the quantile level of

.

  • Step 2: Estimation of with

First, the method is adopted to estimate based on Equation 3. Additionally, given the quantiles of stand-alone and quantiles of exposed risk , the quantile-on-quantile conditional value at risk () can be obtained by the . is estimated similar to [] and []. To embed the dependency among financial markets, the estimation of by is introduced. Following [], the definition of is adjusted as Equation 5 to adapt to a multivariate model.

Assume the of market index is predictable via function , and other indexes’ current return , is the estimator of the function and can be trained by the algorithm []. Therefore, the in the condition of each risk state () can be obtained by Equation 6.

It should be noted that, in distinction to Keilbar and Wang, this paper estimates

for all quantiles of the whole range

4

. Thus, the

of each market index

can be expressed as a three-dimensional surface at each time point by setting the

and

axes to denote different quantile levels

and

, respectively. The

z

-axis values are

. There is an advantage of reflecting both the non-linear relationship between

varying with the risk condition

and the

at each of its own risk levels

.

  • Step 3: Calculation of composite risk spillover effects

To examine the margin impact of index on the , the partial derivative is taken, which is named . According to the format of the function estimated by , the risk spillover from to is expressed as Equation 7,

where the derivative of the transfer function is set as , and all parameters are trained by . This measure should be based on the accuracy of the estimation, but the , as a non-linear neural network model, is prone to over-fitting at a single quantile. Consequently, the adopted in this paper is able to reduce the impact of potential fitting error at a single quantile on the overall risk spillover. Considering the risk spillovers at different quantiles, the composite risk spillover can be defined as as Equation 8.

where the is the absolute value of risk spillover, and the risk-weighted composite risk spillover from index to is defined as . Because the absolute value of the risk indicators (, ) in the more extreme state is larger, the risk spillover in the extreme risk state couple accounts for more weight in the aggregate indicator.

To reflect the sensitivity of average-level composite risk spillover to the two-sided quantiles, we decompose the into two partial spillover indicators and , as Equation 9.

The partial spillover indicators are used to reflect the average impact of the two-sided risk state on the spillover from market to market in the sample period. The partial spillover indicator of side is the composite spillover effect at the exposed risk state of , which reflects the sensitivity to the risk state of the entire system. Correspondingly, the indicates the response of spillover effect to the change in stand-alone risk status . In addition, it is necessary to consider the spillover of market to under the most extreme conditions, that is, or . Hence, we consider the conditional partial spillover indicators under the fixed quantile or fixed as Equation 10.

The non-linear characteristics of inter-market risk spillover can be analyzed by comparing the two types of partial spillover indices under different quantiles, and the mutation quantile of spillover can be captured. In addition, to reflect the extreme risk spillover, the is defined as spillover at a single quantile. If this indicator is calculated by , it will be the same as the spillover defined by [].

Last but not least, similar to the partial spillover indicators, four partial terms can be defined as the average when the one-side quantile takes different values under the condition of the other fixed, that is, and . To reflect the tail risk in the extreme condition, the traditional estimated by and are also calculated as , . Taking the and as an example, Equation 11 is consistent with the partial spillover .

The quantile crossing problem in the estimation of can be described and counted as Equation 12,

For each market, the total number of “quantile crossing” problems can be accumulated by a monotonicity test. Because the

is more adaptable in financial markets, we employ

,

, and

to accumulate numbers of “quantile crossings” for each market.

  • Step 4: Network analysis for composite systemic risk

The average measures of composite systemic risk will be gained in the final step. First, because the risk spillover of market index on market index involves two quantiles and , the average indicators are composited and averaged. Referring to previous works [, ], the composite systemic suffering indicator should be aggregated as Equation 13. The composite index is weighted by the own risk of each spillover emitter. This is because if market has higher own risks, market will bear more systemic risk spillover.

Second, the composite systemic overflow indicator can be defined as Equation 14.

To reveal the trend of the spillover effect, the total composite overflow indicator can be calculated as Equation 15.

The total overflow indicator can be regarded as the weighted sum of , which reflects the changing trend of risk spillover time of each market in the sample.

Lastly, the adjusted adjacency matrix is defined as Equation 16.

The adjusted adjacency matrix accounts are the risk spillover indicators. Systemic spillover effects are thus determined by the marginal effects of the procedure, as well as by the and of the considered countries.

3 Results and discussion

3.1 Data and descriptions

We select the stock market indices of representative countries or regions from the Asia–Pacific Economic Cooperation (APEC) organization varying from January 2012 to December 2021 as sample data. Due to the weekday effect and excessive short-term volatility noise in daily frequency data, the daily returns of stock indexes taken from the WIND database are transformed into weekly data. The logarithmic return is calculated using the closing price on the last trading day of each week. Adopting weekly frequency data can effectively avoid the problem of time differences between the markets, and the data for each week can simply be treated as contemporaneous. The indexes and their abbreviations of regions for the 18 markets are shown in Table 1.

TABLE 1

RegionAbbrIndex
AustraliaAUS&P/ASX 200
CanadaCAS&P/TSX
ChileCLS&P CLX IPSA
ChinaCNShanghai&Shenzhen 300
Hong KongHKHengshen Index
IndonesiaIDJakarta Composite
JapanJPNikkei 225
MalaysiaMYKuala Lumpur KLCI
MexicoMXMMX
New ZealandNZNew Zealand NZ50
PhilippinesPHPhilippines Manila
Republic of KoreaKPKOSPI
Russian FederationRUMOEX
SingaporeSGFTSE Singapore STI
Tai WanTWTaiwan Weighted Index
ThailandTHSE THAI Index
United States of AmericaUSS&P500
Viet NamVNHo-Chi-Minh Index

Abbreviations of regions and indexes.

As shown in Table 2, the total number of observations is 9,396 because some markets have missing samples due to the holidays or other factors. The average weekly return of each market is all positive, which indicates that the indices prices of Asia–Pacific stock markets at the end of 2021 are higher than they were in 2011. It is worth noting that the minimum value of all samples is −20.13%, which is the weekly return of the Chile index in the fourth week of March 2020. The maximum value occurred in the next week, which is 15.81% in Japan. The World Health Organization recognized the COVID-19 outbreak as a global pandemic on 11 March 2020. After 2 weeks of declines, most global stock markets rebounded sharply in the last week of March. Except for the stock markets of Mexico, mainland China, and Hong Kong, all the financial markets have a kurtosis of 3 or more in their return distributions. The return curves show a sharp peak pattern, which indicates that the outliers are more dispersed. In addition, the skewness of all the markets in the sample is negative, indicating that the return series of each market is left skewness; that is, the probability of negative extreme value is higher than positive. The augmented Dickey–Fuller (ADF) value of each market return is negative and less than the test critical value at the 1% significant level, rejecting the null hypothesis of a unit root. All Jarque–Bera (JB) statistics are significant at the 1% level, which rejects the null hypothesis of Gaussian distribution for the market returns.

TABLE 2

abbrSampleMean (%)Std (%)Min (%)Median (%)Max (%)SkewnessKurtosisJBADF
AU5220.121.94−13.980.236.12−1.54088.9241900***−7.806***
CA5220.111.91−16.490.219.07−2.257718.9518,094***−8.064***
CL5220.012.54−20.140.0412.94−0.916010.6862,504***−7.68***
CN5220.153.00−14.020.2810.66−0.59342.303151***−7.024***
HK5220.052.40−9.970.267.60−0.34940.73922***−7.988***
ID5220.112.19−15.690.268.68−1.17498.0441,540***−7.154***
JP5220.242.87−17.430.2915.82−0.50164.980552***−7.958***
MY5220.001.42−9.790.005.49−0.32925.420633***−8.415***
MX5220.072.08−10.560.157.53−0.26212.324120***−8.051***
NZ5220.261.62−15.080.387.86−1.838416.6386,191***−8.393***
PH5220.092.53−19.260.1010.19−1.05999.9472,203***−7.523***
KP5220.092.19−14.130.259.26−0.87436.7791,049***−7.593***
RU5220.192.51−16.040.238.16−0.69884.424458***−8.843***
SG5220.031.88−11.700.149.16−0.25865.689693***−8.215***
TW5220.182.06−11.130.416.78−0.88003.443329***−8.126***
TH5220.092.12−18.960.257.54−1.699114.0184,435***−8.137***
US5220.262.13−16.230.3911.42−1.280411.5302,973***−8.483***
VN5220.282.56−15.720.548.25−1.00274.289489***−7.54***

Descriptive statistics.

3.2 Estimation of multi-quantile and

3.2.1 Estimation of multi-quantile

The series for each market return can be calculated based on the methodology described in Section 2.2. Different from the conventional calculation, a list of quantiles [5%, 10%, , 95%] is selected to calculate the multi-quantile () at each quantile.

The line diagrams of the United States and China are shown in Figure 35. Scatter points in gray represent the weekly returns, while the line represents the at each quantile level. As can be seen from the diagrams, the lines approximately envelop the points of return. From the perspective of the whole-time sequence, the Chinese stock market fluctuates more. In contrast, the US stock market returns are more concentrated overall, except for two extreme values: late 2011 and early 2020. Since the second half of 2011, the US stock market has suffered severe shocks, especially due to the downgrade of the US credit rating and the continued deterioration of the European debt crisis, which triggered investor panic and increased US stock volatility. In addition, 2020 was the time of the worldwide outbreak of the COVID-19 pandemic. These two unexpected events hugely impacted the US stock market, while during the remaining time, the US stock market was relatively stable compared to the Chinese stock market. The patterns of the charts reflect the differences in the maturity of the two markets. In terms of general trends, the Chinese stock market is more volatile, with larger absolute market returns under extreme conditions because of the frequent policy intervention and excessive proportion of individual investors. Meanwhile, the of the US market is smaller in absolute terms, indicating that its market is able to return to a steady state relatively quickly after a short-term shock. Furthermore, the lines of China at different quantiles are asymmetrical in the vertical dimension. Its envelope area is larger below the zero value; in other words, its value stays more in the negative zone. In contrast, the values of the United States are more concentrated in the positive zone. According to the distribution of returns, the Chinese stock market tends to suffer losses, while the US stock market tends to make profits.

FIGURE 3

3.2.2 Comparison of estimated by and

According to the count method described in Equation 11, the number of cross-quantile occurrences in each market is shown in Table 3. Because 18 pairs of adjacent quantiles and 522 periods in each market are compared, the first column of the table shows that the quantile crossing is a common problem in the estimation of . Even when considering the average level, the quantile crossing is still obvious when comparing 18 pairs of quantiles. However, completely eliminates this problem. To further analyze the level of risk in the market under different conditions, as described in Section 2.2, the three-dimensional mesh-surface graphs were plotted to present the at each quantile of a stand-alone state and quantile of an exposed state. Figure 4 presents the MCQRNN-QQ-CoVaR plots for other markets indexes6.

TABLE 3

IDQRNNMCQRNN
AU1,39621000
CA1,52920000
CL2,20230000
CN1,36130000
HK88010000
ID1,14920000
JP86010000
MY1,54130000
MX1,39920000
NZ1,47710000
PH1,75320000
KP84110000
RU1,63640000
SG1,15330000
TW82510000
TH38100000
US1,77010000
VN78700000
Total22,940331000

Cross-quantile count by the monotonicity test.

FIGURE 4

We select China and the United States as outstanding examples. Figure 5 represents the level of and non-linear characteristics in the United States and mainland China7.

FIGURE 5

By comparing the surfaces in Figures 5A, B and in Figures 5E, F, it is found that the surface calculated by is smoother than that calculated by . That means less bias in calculations for an extreme quantile, such as by . On the same axis scale, the surface of the United States is flatter, which means fewer sensitivities to the quantile of stand-alone risk.

The reasons for the differences between China and the United States are market structure, investor base, and risk distribution. The American market is the largest and most diversified in the world, with companies representing a wide array of industries and sectors. This diversification helps to mitigate stand-alone risks associated with individual companies, sectors, or events. Additionally, American institutional investors like mutual funds, pension funds, and hedge funds play a significant role. These institutions usually employ sophisticated risk management strategies, including diversification and hedging, which further diminish sensitivity to stand-alone risks. Furthermore, the US market offers a wide array of financial instruments, such as options, futures, and swaps, that allow for the hedging of specific risks. This availability of hedging tools enables market participants to isolate and manage stand-alone risks effectively. Because the estimated is based on each quantile separately, there is a non-monotonic trend with the change of quantile in Figure 5E. In order to compare the differences between and in different conditions, we analyze the of single quantiles and multi-quantiles separately.

Considering that in this paper takes into account two quantiles, partial in two directions are shown in Figures 5C, D according to Equation 11. On the one hand, as shown in Figure 5C, the red line and the blue line respectively represent the calculated by and , when the quantile is 0.01. Meanwhile, the green and gray line respectively represent the average level of calculated by and at each quantile. In extreme conditions, the blue line tends to be a straight line, while the red line fluctuates around it. The phenomenon of “quantile crossing” occurs when the of is higher than that of , which indicates that is more robust than . The “quantile crossing” problem is weakened at the average level as shown in green line. On the other hand, as shown in Figure 5D, four lines respectively represent partial exposed to the entire risk, which states nonlinear characteristics. The results of and are similar at the average level, as shown by the green and gray lines in Figure 5. However, under extreme conditions, the red lines are always below the blue lines, and the phenomenon of quantile crossing still appears. Therefore, the improved calculated by is more robust when considering multiple and two-sided quantiles to avoid the “quantile crossing” problem.

3.3 Analysis of composite risk spillovers

3.3.1 Network analysis of composite risk spillovers

According to Equation 8, after the estimation of , the composite risk spillover can be obtained as the adjacency matrix shown in Figure 6 and the network graph shown in Figure 7. In order to draw the risk spillover network maps across markets, the adjacency matrix can be obtained by the average level of sample period . A weighted directed network can be plotted on the basis of this adjacency matrix.

FIGURE 6

FIGURE 7

As shown in Figure 6, each cell of the matrix represents a risk spillover correlation between the two markets. Where the color is darker, the level of risk spillover represented by the cell is higher. The cells with relatively dark colors in the graph are, respectively, the risk spillovers of HK →CN, CA →US, and US →CA. The mutual risk spillover between the United States and Canada can be explained by their geographical location, economic connections, and political policies. Canada is adjacent to the United States, and the two countries have a very close political relationship and an active trade association. Similar to the US and CA, the geographical location and economic connections between CN and HK are quite tight. However, the risk spillover level of HK →CN is the highest, while that of CN →HK is much lower. It is obvious that the high level of risk spillover effect from Hong Kong to the mainland China is due to the economic linkage between the two countries and the effect of Shanghai-Hong Kong Stock Connect program. In contrast, the economic policy of the Chinese system is different from that of western systems. The financial institutions in mainland China are not aggressive in investing. Moreover, the trade from mainland China to Hong Kong concentrates on domestic goods, which are at low prices. Those goods are why the risk spillover from CN to HK is relatively low. Therefore, the Hong Kong stock market is more mature and less susceptible to shocks.

In the following, we present a two-way weighted network to analyze systemic risks with a clearer visual structure. First, the adjacency matrix needs to be read via the NetworkX package in Python. Figure 7 shows a network map using the mean value of the samples. The arrow indicates the direction of the risk spillover. Both the size of the arrow and the width of the line segment indicate the intensity of risk spillover. Note that the width of the line segment states the level of the spillover of the larger one in the two-way relationship, in which the thinner one is covered. Therefore, the level of the risk spillover can only be judged based on the arrows in the comparison of two-way relationships. As can be seen, the most prominent line segment in this map is from HK to CN because of the Shanghai-Hong Kong Stock Connect and Shenzhen-Hong Kong Stock Connect programs. The level of risk spillover between the United States and Canada is also high, but the two-way relationship is symmetrical with almost equal size arrows. Similar two-way relationships also exist between TW and KP, JP and KP, ID and PH, and AU and CA. Such two-way relationships can be explained by the frequent trade interactions. It implies that stock markets are not only barometers of the economy but also effective reflections of the economic trades and global value chains through risk spillovers among financial markets.

Figure 8 is a frequency histogram of the risk spillover relationships. Four colors represent four algorithms adopted in calculating the risk spillover. The intensity of risk spillovers can be seen to exhibit a right-skewed spike with a thick tail. This indicates that while most risk spillovers are at low levels, the risk in extreme conditions is substantially outside the average range. In addition, the distribution of the red line is relatively flat, which means the risk spillover may be overestimated by at a single quantile. The results of the remaining three algorithms are similar, though the peak of the green line is slightly skewed to the right.

FIGURE 8

3.3.2 Comparison of composite risk spillover calculated by QRNN and MCQRNN

A 3D-mesh surface is also employed to illustrate the one-way spillover from Hong Kong to mainland China, which is the most significant correlation in the Asia–Pacific region. As can be seen from the result from Figure 9A, the from Hong Kong to the mainland is fluctuating and outstanding at the extreme quantile . However, as shown in Figure 9B, can show relatively gentle overflow changes, especially showing no mutation characteristics at the extreme level. Similar to the algorithm comparison diagram of in Figure 5, it can be seen from Figure 9C that spillover levels represented by the blue and gray lines are more stable than those of the red and green lines. Meanwhile, the extreme situation of does not appear. This suggests that the estimators at the extreme quantile may be very sensitive to outliers under the partial differential spillover method. In contrast, can both make the estimation results of more robust and obtain a more accurate assessment of the overflow level. In addition, from Figures 9C, D, the partial spillover increases with the decline of quantile and ; that is, during the challenging period, the risk spillover from Hong Kong to Chinese mainland is higher.

FIGURE 9

3.3.3 Trend of total composite overflow indicator

To compare the overflow dynamic throughout the sample period, the time series diagram is drawn in Figure 10. Four lines in various colors represent the overall overflow levels of the two algorithms at the extreme level or the average level, based on the computational method of in Section 2.4.

FIGURE 10

Compared with the composite overflow indicators calculated by the multi-quantile algorithm (in green and gray), the peak of overflow levels in the extreme condition represented by the red and blue lines are relatively higher because the peaks of overflow levels at the multi-quantile are flattened by averaging. Moreover, although spillover instability under extreme quantile conditions is reduced, the fluctuations of overflow calculated by the method reveal more significance during periods of higher systemic risk. In other words, also presents more robust and significant results even when calculating overflow under extreme conditions. The most prominent period was during the COVID-19 pandemic, which showed a higher level of spillover than other periods. Under the influence of this extreme event, the global real economy has stagnated, and production has been interrupted, leading to investor panic and insufficient investment confidence. Therefore, risk accumulates, and global asset prices fall. In contrast, the peak value calculated by at a single quantile in this period is not different from that in other periods. Therefore, the single-quantile method is worse than others.

Although the fluctuations of each line are different, there are four significant periods with high overall composite overflow levels. The first period is from May 2012 to September 2013, which corresponds to the EU debt crisis and the US stock market crash. The second period is in the second half of 2015 before the Chinese stock market crash occurred. The third period begins in 2018, which corresponds to the Sino-American trade war. The last fluctuant period is from March 2020 to February 2021, which is caused by the outbreak of COVID-19.

3.4 Comparison of systemic risk models

In this part, we analyze the average overflow of each market to the systemic , and their suffering , relying on the methodology in8 Section 2.4.

The suffering indicators and the overflow indicators , which are respectively calculated by single-, multi-, single-, and multi-, are drawn as bar charts in Figure 11. As the legend shows, the suffering indicators are in light color and on the left side of each market, while the overflow indicators are darker and on the right side. In addition, the calculated by and are illustrated as light and dark orange in the secondary axis. It is obvious that the light red bars stand out, indicating that China suffers the highest risk overflow when the calculated by at a single quantile. However, when the algorithm is substituted by with multiple quantiles, China is no longer the highest suffering market. In addition, the obtained by and are relatively close, and the of each market states no significant correlation with both and indicators.

FIGURE 11

Although Figure 11 shows differences in index calculations under the four algorithms, the comprehensive index is smaller than that calculated by a single sub-site. Whether under single or multi-component sites, the exponents obtained by the algorithm are all smaller than those obtained by the algorithm. In order to investigate whether the results of algorithms affect the ranking of systemic risk in each market, we also list the top ten markets of systemic risk index under various algorithms. As shown in Table 4, markets with higher risk overflow are HK, CA, US, and SG. This result is consistent with the global financial markets’ practical experiences. In contrast, JP, CN, HK, and PH suffer more systemic risk because of the higher . On the one hand, this may be related to the fact that these countries are more dependent on trade exports and have poor economic resilience, resulting in suffering more risk overflow. On the other hand, higher may also indicate that investors in these stock markets react strongly to the shocks. Compared with , the of China is driven higher than it calculated by , which means that the method weakens the impact of the extreme condition.

TABLE 4

IDSystemic suffering indicatorSystemic overflow indicator
Single quantileMultiple quantileSingle quantileMultiple quantile
QRNNMCQRNNQRNNMCQRNNQRNNMCQRNNQRNNMCQRNN
AU0.07420.06580.05680.05070.06810.06100.06530.0577
CA0.06560.04260.05260.04700.11870.07100.09920.0830
CL0.11320.03920.07390.05920.09810.05640.05590.0458
CN0.17130.10190.09690.07630.06980.03480.04520.0343
HK0.08740.05390.07800.06720.13740.10040.09890.0882
ID0.11760.05850.06550.05200.07350.05750.06330.0545
JP0.13920.09620.09240.08090.10220.05430.06200.0574
KP0.07100.06790.05980.05300.08500.07100.08440.0717
MX0.07470.04440.06180.05380.07170.04170.05890.0483
MY0.06330.01890.04330.03390.09190.05310.05550.0464
NZ0.06050.02820.04560.03310.07850.04640.04800.0357
PH0.10440.07920.07680.06340.09680.05770.06150.0513
RU0.08070.04070.07200.05830.07160.03560.04530.0357
SG0.05800.04210.05290.05030.13850.07870.08870.0769
TH0.08090.03460.05490.04520.09910.04950.05940.0468
TW0.06250.05040.05730.04820.09850.04980.05960.0473
US0.08650.07030.06730.06160.11410.08040.09140.0776
VN0.10280.07590.06630.04840.03750.02850.03430.0259

Composite systemic risk indicators of each market.

4 Discussion and conclusion

To improve the traditional paradigm of risk spillovers among financial markets, this paper has calculated a multi-quantile based on . This study broadens the perspective of risk spillover research to financial market indexes and takes into account both the analysis of tail risk spillovers and risk spillovers under normality. The following conclusions were drawn based on empirical analyses of the Asia–Pacific region:

First, by visualizing the partial , the “quantile crossing” problem is found on the estimation of but can be relieved by . This issue can be described as “a worse condition may cause less risk loss.” It is not only a logical problem but also reveals that may be sensitive to quantile selection. Fortunately, this problem rarely occurs on at the quantile of exposure risk state and can be relieved by substituting for . On the other hand, through the estimation of partial at the other side quantile, the non-linear characteristic of each market is visualized. Different from the at the stand-alone risk state, the value at risk declines rapidly when exposed risk rises to an extreme state. This concludes that the non-linear algorithm is suitable for estimating precisely without over-fitting.

Second, the overestimation of spillover may occur when calculated by at a single quantile, compared with or at multiple quantiles. Based on the comparison of two 3D-mesh surfaces and the line charts of partial of four algorithms, the over-fitting in extreme conditions may contribute to the overestimation of composite spillover. The estimation of at a single quantile shows less robustness than other three algorithms in the trend chart of total composite overflow.

Third, the stock market in mainland China is highly exposed to the risk spillovers from the Hong Kong stock market. In addition to the short geographical distance between them, another reason for this relatively asymmetric risk spillover may be that investors in mainland China are more concerned about the opposite, but investors in Hong Kong are more independent and have more complete information. In addition, the mutual spillovers between the United States and Canada are also significant, which may be due to the special geographic relationship between the United States and Canada as well as tight trade cooperation and economic dependency between two markets. Cross-market comparisons show that the model supports the traditional view that Hong Kong, Canada, United States, and Singapore are more important markets in the Asia–Pacific region. In contrast, the Chinese mainland and Japanese markets received the most spillovers during the sample period.

This paper studies the systemic risk and risk spillover under multiple quantiles, providing a reference for stock investment and risk regulation in the Asia–Pacific market. This method can not only be applied to the study of inter-institutional risk spillover but can also be helpful in capturing the nonlinear characteristics of individuals’ systemic risk. However, this paper still has some shortcomings. Limited by the time and space complexity of the algorithms, it is impossible to use the rolling window to estimate and calculate the daily samples. Therefore, the out-of-sample prediction effect of the model cannot be investigated. Further study is needed to improve the efficiency of the model and expand the sample.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.

Author contributions

CR: conceptualization, investigation, methodology, software, supervision, writing–original draft, and writing–review and editing. ZZ: conceptualization, formal analysis, investigation, resources, and writing–original draft. DZ: formal analysis, investigation, resources, and writing–original draft.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This research is funded by the National Nature Science Foundation of China (NSFC), grant number 72173018; the Humanity and Social Science Research Project of Anhui Educational Committee, grant number 2024AH052470; and the High-level Talent Research Initiation Project at Anhui Business College, grant number 2024KYQD05.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphy.2024.1484589/full#supplementary-material

Footnotes

1.^“Higher risk but less loss” means the greater the risk, the smaller the loss because of the quantile crossing.

2.^ is the stand-alone level; of is a vector.

3.^ represents the conditional exposed risk.

4.^in is a vector, which is the same as in Equation 2.

5.^MQ-CAViaR charts for other stock markets are shown in Supplementary Appendix Figure SA1.

6.^Supplementary Appendix Figure SA2 presents the QRNN-QQ-CoVaR 3D-surface plots for other market indexes.

7.^Supplementary Appendix Figure SA3 represents the stand-alone plots for other market indexes. Supplementary Appendix Figure SA4 represents the exposed partial plots for other market indexes.

8.^In Section 3.4, we calculate the average values for and , respectively. The equations are and .

References

Summary

Keywords

multiple quantile risk spillover, MCQRNN, QQ-CoVaR, systemic financial risk, quantile crossing

Citation

Ren C, Zhu Z and Zhou D (2024) Multi-quantile systemic financial risk based on a monotone composite quantile regression neural network. Front. Phys. 12:1484589. doi: 10.3389/fphy.2024.1484589

Received

24 August 2024

Accepted

09 October 2024

Published

18 November 2024

Volume

12 - 2024

Edited by

Hui-Jia Li, Nankai University, China

Reviewed by

Peican Zhu, Northwestern Polytechnical University, China

Ge Gao, Beijing Sport University, China

Updates

Copyright

*Correspondence: Chao Ren,

ORCID ID: Ziyan Zhu, orcid.org/ 0000-0002-2658-9867; Donghai Zhou, orcid.org/ 0000-0001-9160-8713

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics