On Spillover Effect of Systemic Risk of Listed Securities Companies in China Based on Extended CoVaR Model

Based on the daily data from January 2, 2019, to September 30, 2020, this paper uses the extended CoVaR model to measure the spillover eﬀect of systemic risk among top 10 securities companies by market value in China, All Share Brokerage Index, All Share Financials Index, All Share Insurance Index, and CSI Banks Index. The conclusions are as follows: (1) there are risk spillover eﬀects among 10 securities companies, which are asymmetric and bidirectional and highly volatile in a short period of time; (2) the spillover eﬀect of systematic risk of securities companies is not necessarily related to the market value of securities companies but has a strong relationship with the stock market; (3) there are risk spillover eﬀects between the sample securities companies and the four major indexes, but there are signiﬁcant diﬀerences in the size of the spillover eﬀects; (4) the securities industry has a great risk spillover eﬀect on the ﬁnancial industry, but the risk spillover eﬀect of other ﬁnancial sectors on the securities industry is very small. Finally, we put forward countermeasures and suggestions.


Introduction
Once the systemic financial risk occurs, it will not only harm the financial field but also cause extremely serious losses to the macroeconomy and social wealth. Because of its complexity and strong infectivity, it is difficult to identify and measure systemic financial risk. In recent years, the contagion of financial risk in China has been increasing, and the possibility of systemic risk in financial field is increasing.
Since the 1990s, China has vigorously developed the direct financing market, and the securities industry has developed rapidly. As of June 30, 2020, China's 164 securities companies have total assets of 8.03 trillion yuan, net assets of 2.09 trillion yuan, net capital of 1.67 trillion yuan, customer transaction settlement fund balance of 1.64 trillion yuan, and total entrusted management capital of 1.183 billion yuan. However, the securities industry has the smallest size compared to the banking, insurance, and trust industries, and its financing channels are limited. At the same time, the market environment and policy constraints restrict the innovation of the securities industry. Traditional businesses such as securities brokerage and proprietary securities are still the main profit points of the securities industry, making the impact of the market easily spread to the securities industry. All these make the securities industry's antirisk ability relatively weak [1,2]. As time progresses, China's capital market has gradually become more open, which indicates that the development of securities industry will be faced with domestic and foreign risks. In addition, the deepening of the degree of mixed operation and promoting the registration system and other deepening reform measures in the securities industry have not only brought innovation and dividends but also increased the systematic risk of the securities industry.
Systemic risk can cause many market participants to suffer losses at the same time and spread rapidly in the market [3]. Compared with the risk of an individual risk, systemic risk has five basic characteristics: complexity, suddenness, rapid contagion, wide spread, and great harm [4]. erefore, many scholars are committed to studying how to scientifically monitor and measure systemic financial risks and prevent and resolve risk spillovers.
At present, the research methods of risk spillover effect can be summarized as follows: one is the financial network model proposed by Allen and Gale [5]. e core idea of the theory is to study the structural characteristics of financial network by collecting the data of mutual exposure or transaction between financial institutions, so as to study the contagion effect of financial risks in the network. e second is the systematic marginal expected loss (MES) method proposed by Acharya. is method calculates the risk contribution of a single participant to the market when the crisis breaks out. is method can measure the loss of a single financial institution and its proportion in the whole market [6]. e third is the CoVaR model proposed by Adriany and Brunnermeier. e model measures the loss of other institutions or markets when there is risk in one institution or market [7]. Since CoVaR method uses highfrequency data modeling, it can sensitively capture the contagion of financial risks and compute the risk spillover effect in the financial network more realistically, so it has gradually become an important tool for measuring financial risk contagion and spillover effect [8].
In recent years, some scholars have revised CoVaR model and applied the improved model to the estimation of systemic risk. Rungporn Roengpity and others used the CoVaR method to find that every large ai commercial bank has systemic risk spillovers to the banking system [9]. Li [13]. Girardi and Tolga Ergun used the CoVaR value to analyze the systemic risks of the US financial industry, and the results shows that there is no linear relationship between the size of financial enterprises and the size of systemic risk [14]. Reboredo and Ugolini used the CoVaR-Copula model to study the different characteristics of systemic risks in the European sovereign debt market before and after the Greek debt crisis broke out. ey found that the systemic risk of the crisis countries decreased after the crisis, while the systemic risk of Greece to the crisis countries increased, and the systemic risk to other countries decreased [15]. Härdle et al. applied generalized quantile regression to compute the CoVaR value of American financial companies, which organically combined financial tail events with network technology. According to the value of risk spillover, they ranked risk receivers and disseminators in the US financial market, and found that the risk of acceptance and transmission of depository institutions is greater, while insurance companies are less affected [16]. Xu et al. used the quantile regression method to estimate the CoVaR value of listed insurance companies in China. It shows that the intensity of risk spillover varies with different insurance companies; the higher the proportion of nontraditional insurance of insurance companies is, the greater the systemic risk is; solvency regulation has a significant effect on restraining the risk spillover level of insurance companies [17]. Yuan and Wang found that there were significant two-way risk spillovers among insurance companies and insurance industry, showing asymmetry and difference [18]. Long et al. calculated the CoVaR and ΔCoVaR of the CSI 300 secondary industry index by constructing a time-varying t-copula function and found that the systemic risk contribution of large industries such as banking and energy is not as high as that of small-scale industries like food and beverage [19]. Yang et al., respectively, measured the systemic risks of 56 Chinese listed enterprises, including financial institutions and real estate companies, using the methods of VaR, MES, CoVaR, and ΔCoVaR. e research shows that the level of systemic risk spillover in China is rising year by year, and the transmission center has changed correspondingly in the events of "bank money shortage" and "stock market circuit breaker mechanism" [20]. Li et al. used the ΔCoVaR to study and found that the risk spillover effect of China's banking and insurance sectors is higher than that of the securities industry [21]. Wang and Yang used EVT-Copula-CoVaR model to compute the risk spillover in Chinese stock market and found that it has a two-way and asymmetric risk spillover effect with other countries along "the Belt and Road" [22]. Xu et al. applied the systemic risk network based on LASSO-CoVaR to study the systemic risk of Chinese financial institutions from 2010 to 2017 and found that, for even the same institution, its risk has changed in different periods [23]. Zhang and Xu used ΔCoVaR to measure the systemic risk brought by China's listed banks participating in financial derivatives transactions. ey found that the risk has increased compared to before participating, and the impact of interest rate derivatives on systemic risk was greater than that of foreign exchange derivatives [24]. Tiwari et al. used ΔCoVaR and MES to capture system risk and found that the risk spillover of oil price changes in the securities markets of G7 countries increased with market turmoil [25].
To sum up, scholars' research on systemic risk measurement is more and more in-depth, and CoVaR model has been constantly expanded and improved. However, the existing research on risk spillover effect mainly focuses on commercial banks, insurance companies, or financial institutions as a whole, while there are few studies on securities companies and securities industry. In addition, most of the existing literatures using CoVaR model to measure systemic risk spillovers in the securities industry still remain on the basic CoVaR and ΔCoVaR methods, which will not capture the data characteristics of financial time series well.

Complexity
Compared with the simple GARCH model, on the one hand, the GJR-GARCH model can more deeply describe the asymmetry of volatility spillover between markets and is superior to the GARCH model and other related models in describing volatility and forecasting [26]; on the other hand, the DCC-GARCH model overcomes the shortcomings of the previous multivariate GARCH model. It not only has good computational advantages and can be used to estimate the correlation coefficient matrix that was originally difficult to calculate due to its scale but also can be used to study the dynamic correlation between multiple assets under the influence of market information and policy guidance in different periods [27].
e DCC-GJR-GARCH model combines the DCC-GARCH model and the GJR-GARCH model and combines the advantages of the two models, making the fitting of financial time series more realistic and convenient.
e DCC-GJR-GARCH-CoVaR model constructed on the basis of the DCC-GJR-GARCH model, uses the variance between the financial time series processed by the DCC-GARCH model, and calculates the CoVaR value representing the value of the risk spillover, which can well measure the size of the systemic risk spillover.
In view of this, we use the DCC-GJR-GARCH-CoVaR model to study the systemic risk spillover effect of listed securities companies and its changing rules, and the two-way risk spillover between them and the financial industry is comprehensively considered.

Construction of DCC-GJR-GARCH-CoVaR Model
is part first brings the basic definition of the CoVaR model into consideration, then explains how to use the DCC-GJR-GARCH model to estimate the conditional variance of yield series, and finally introduces the estimation steps of CoVaR value and ΔCoVaR value of yield series.

Introduction to CoVaR Model.
Adrian and Brunnermeier [28] proposed CoVaR based on VaR. CoVaR is defined as when the loss X i of financial institutions is VaR i q , under the confidence level of q, the VaR is faced by financial institution j, that is, us, the systemic risk spillover value (ΔCoVaR) of financial institution i to financial institution j or the whole system can be expressed as follows: Among them, VaR j q is the value at risk faced by financial institution i when it is not under extreme risk conditions.

Fitting Conditional Variance of DCC-GJR-GARCH
Model. GARCH model can explain heteroscedasticity effect and volatility aggregation of financial time series. In addition, Engel's research found that the DCC-GARCH model can effectively resolve the difficulty in too many parameters in the multivariate GARCH model by allowing variable conditional correlation coefficients and meet the characteristics that the influence degree between financial time series is not fixed but changes with time [28].
e GJR-GARCH model first introduced by Glosten et al., through bringing in the dummy variables, makes the parameters in the basic GARCH model can be negative [29]. erefore, this article uses the DCC-GJR-GARCH model introduced by Glosten and other scholars to estimate the conditional variance.
We assume that the residual series e t of financial institutions' return series is white noise series, the covariance matrix is H t , and ] is the degree of freedom of multivariate joint t distribution: where Ω t−1 represents the previous information set and H t is the covariance matrix, which can be expressed as in which D t of equation (4) is the variance matrix of the residual error e t , and it can be shown as follows: Among them, h ii,t is the conditional variance of the residual error e t , and it is characterized by GARCH (1, 1), that is, e expression of schematic function I t−1 of equation (6) is as follows: When there is good news in the market (e i,t−1 ≥ 0), the variance h ii,t is affected by α i . When the market appears in bad news (e i,t−1 < 0), the variance h ii,t is affected by (α i + g).R t of equation (4) is the dynamic correlation coefficient matrix. e expression is as follows: where q nn,t is the conditional variance of standardized residuals. erefore, the dynamic correlation coefficient (ρ ij,t ) in R t is as follows:ρ ij,t � q ij,t / ������ q ii,t q jj,t Q t of equation (8) can be expressed as where Q is the unconditional variance matrix of standardized residuals ε t and ε t � D −1 t e t is the vector normalized residuals.

Estimation Steps of CoVaR. Based on the basic definition of CoVaR and the conditional variance fitted by the DCC-GJR-GARCH model, we can estimate the CoVaR value and
ΔCoVaR value of the research object. e specific estimation steps are as follows: the first step is to collect the time series data of financial institutions and establish GARCH (p, q) model to measure the volatility of their returns: Among them, r i,t is the yield of financial institution i in period t; h 2 i,t is the variance of the residual error e i,t in the t period; e 2 i,t−1 is the variance of the residual error e i,t in t−1 period, namely, ARCH term; h 2 t−1 denotes the GARCH term; and dummy variable I t−1 represents the asymmetric effect of new information.β i and c i are constant terms, and A p (L) and B p (L) are lagging operators.
In the second step, the GARCH model is used to estimate the VaR i value of financial institutions. e formula is as follows: Among them, r i,t is the one-step forward estimate of the return rate of financial institutions i obtained by the GARCH model, and h i,t is the one-step forward estimation value of h ii,t . F − 1 is the inverse function of the distribution function of t distribution, and q is the confidence level [31].
In the third step, using the DCC-GJR-GARCH model and the results of the previous step, the CoVaR value of the mechanism j under the influence of the mechanism i and macrostate is constructed. e formula is as follows: where r j,t is the yield of financial institution j and r i,t and r m,t−1 represent the impact of financial institution i and macrosituation on financial institution j, respectively. A p (L) and B p (L) are lag operators, and θ j|i is the regression coefficient of financial institution i in the formula. Introducing r i,t � VaR i,t into the mean value equation, the following formula is obtained: According to the mean value r j,t and the conditional variance h j,t , the CoVaR and ΔCoVaR of financial institution j are obtained. e formula is as follows: where VaR i,q,t is the same as the VaR i value of financial institution i calculated by formula (12).

Data Selection and Processing.
is paper uses the stock price index of the top 10 securities companies with total market value which listed in China's A-share market before January 1, 2019. ey are GF Securities Co., Ltd. representing the securities industry, All Share Financials Index (000992. SH) representing the financial industry, All Share Insurance Index (h30186. CSI) representing the insurance industry, and CSI Banks Index (399986. SZ) representing the banking industry as a research sample, use the CSI 500 (000905. SH) index as a macrovariable, and choose the closing price on each trading day as the original data. e data selection time limit is from January 2, 2019, to September 30, 2020. All data are from Wind database. e formula for calculating the daily return of sample companies stocks and the indexes is as follows: Among them, P t is the closing price of the stocks or indexes on day t.

Descriptive Statistics of Data.
We conducted descriptive statistics on the daily yield series of samples (see Table 1 for details).
According to Table 1, the average returns of 10 securities companies stocks are all greater than 0. Among them, the largest average yield is the stock of China Securities Co., Ltd., reaching 0.4123; the smallest is the stock of Guotai Junan Securities Co., Ltd., which is 0.0497. e standard deviation of the stock return sequences is generally more than 2 units larger than the mean and much larger for the mean between 0.05 and 0.4, and the difference in the yield of the same stock can reach 20%, which all shows that the returns of stocks have changed drastically. We can observe that the kurtosis of the stock return sequences of all securities companies are greater than 3, showing a certain degree of right deviation.
e JB statistics are all very large, and the P values are all less than 0.01, showing the characteristic of the financial time series with peaks and thick tails. erefore, this paper chooses t distribution to describe the stock return sequences, which is more consistent with the characteristics of the data itself.
Next, we conduct descriptive statistics on the daily yield series of the CSI 500 Index, All Share Financials Index, All Share Insurance Index, CSI Banks Index, and CSI 500 Index in the inspection period (see Table 2 for details).
According to Table 2, the average returns of 5 indexes are all greater than 0, which is much smaller than the mean of the stocks of securities companies. Among them, the largest average yield is the All Share Brokerage Index, reaching 0.1196; the smallest is CSI Banks Index, which is 0.0158. Compared with the mean, the standard deviation of the yield series of indexes is still very large and is generally more than 1 unit higher than the mean between 0.01 and 0.1, indicating that the yield of indexes has also changed drastically. e kurtosis of all yield series of indexes is all greater than 3, and there is a certain degree of left or right deviation. e JB statistics are all large, and the P values do not exceed 0.01. e original hypothesis that the yield series is a Gaussian distribution is rejected at the 1% significance level, indicating that the t distribution is better for describing the yield series of indexes.

Stability Test.
When the variable is time series data, if it is not tested for stationarity, there may be pseudoregression in regression analysis. Here, we conduct the ADF test on the daily yield series of ten listed securities company stocks, CSI 500 Index, All Share Financials Index, All Share Insurance Index, CSI Bank Index, and All Share Brokerage Index during the inspection period (see Tables 3 and 4 for details).
According to Tables 3 and 4, whether it is the stock yield series of securities companies or the index yield series, their data during the inspection period have passed the ADF test, and the original hypothesis that the yield series have unit root is rejected by all results at the 1% significance level.

Autocorrelation Test.
When the autocorrelation of random error term appears in time series data, the real variance of parameter estimation value will be underestimated, the value of t statistic will be overestimated, and the original unimportant explanatory variables may be mistakenly considered important and retained, which will eventually lead to the reduction in the accuracy of the model prediction.
We use Eviews software to calculate the ACF value and PACF value of ten securities company stocks and five index yield series, and the original hypothesis that the series having autocorrelation is rejected by all results at the 5% significance level.

ARCH Effect Test.
Engle [32] pointed out that heteroscedasticity may occur in time series data and thus proposed a method for observing variance changes in time series, namely, ARCH test method [33]. We assume the number of lag periods in the ARCH test is one and conduct the ARCH effect test on the daily yield series of ten securities company stocks and five indexes during the inspection period (see Tables 5 and 6 for details).
According to Table 5, the P value of the n * R-squared statistic of the stock yield series of ten securities companies is all below 0.1; that is, the original hypothesis that the random error term of the yield series has no heteroscedasticity is rejected at the significance level of 10%. e series have the ARCH effect under the assumption that the random error term lags for one period.
According to Table 6, in addition to the CSI 500 index, which is a macrovariable, in the test results, the P values of the Obs. n * R-squared of the All Share Financials Index, All Share Insurance Index, CSI Banks Index, and All Share Brokerage Index yield series are all below 0.1, and the original hypothesis that the random error term of the yield series has no heteroscedasticity is rejected at the significance level of 10%. e series have the ARCH effect.

Analysis of Systemic Risk Spillover Effect among Securities
Companies.
e analysis steps of systemic risk spillover effect among securities companies are as follows: first, use   (17) to calculate the stock returns of 10 sample companies. Second, use formula (10) and equation (11) to fit the mean and variance of the 10 stocks' return series, so as to obtain the degree of freedom of t distribution. irdly, predict the mean forecast value r i,t and variance forecast value h i,t by the estimation formula one-step forward and substitute the degree of freedom of t distribution and the 95% confidence level into equation (12) to calculate the value of 10 stock return series at 95% confidence level. Fourth, combine the return series of 10 stocks in pairs to get 45 combinations and use the DCC-GJR-GARCH (1, 1) model of Student t distribution to fit the mean and variance of the return series of these 45 portfolios. At the same time, the return series of the CSI 500 Index and the stock return series as the influencing items are added into the mean estimation formula to obtain the following two formulas: Fifthly, predict the mean r j,t and the conditional variance h j,t by the estimation formula one-step forward, substitute the degree of freedom of t distribution and 95% confidence level into equation (15), and calculate the CoVaR j|i value of 182 combinations at 95% confidence level.
Sixth, substitute CoVaR n and VaR n into equation (16) to get 182 ΔCoVaR j|i values. e results are shown in Table 7.
It can be seen from Table 7 that there is a certain risk spillover effect among securities companies, and the risk spillover effects among different companies are quite different.
In terms of the spread of systemic risks, we can observe that the two largest companies in the industry: Citic Securities Co., Ltd. and China Securities Co., Ltd., have played an important role in spreading risks, and their spreading ΔCoVaR values are generally 2-3 units higher than normal level; on the contrary, GF Securities Co., Ltd. and Guotai Junan Securities Co., Ltd. are less likely to spread the risk, and the ΔCoVaR value of their spread is about 2-3 units lower than the normal level. e ΔCoVaR value spread by other securities companies is at a normal level, and the possibility of spreading systemic risks is moderate. is shows that securities companies with a larger volume have more business volume than other securities companies and have closer business ties with other securities companies.    (1) represents the degrees of freedom of the chi-square distribution and also represents the number of lag periods assumed in the ARCH test.
6 Complexity ey are very important to the securities system and spread systemic risks. e possibility is greater.
In addition, in terms of systemic risk acceptance, we can observe that Shen Wan Hong Yuan Securities Co., Ltd. is more likely to accept the risk and its accepted ΔCoVaR value is about 1 unit higher than the normal level; on the contrary, the probability of China Securities Co., Ltd. accepting the risk is very small, much lower than the normal level, and the accepted ΔCoVaR value is about 3 units lower than the normal level. e ΔCoVaR value accepted by the remaining securities companies is at a normal level, and the possibility of accepting systemic risks is moderate. is shows that the level of risk management of China Securities Co., Ltd. is relatively high, while the level of risk management of Shen Wan Hong Yuan Securities Co., Ltd. is relatively low.
By comparing the differences in risk spread and risk acceptance of different securities companies, it is not difficult to find that there is an asymmetric two-way spillover effect among securities companies; for example, the ΔCoVaR value of China Galaxy Securities Co., Ltd. for other securities companies is significantly higher than its accepted ΔCoVaR value, generally 3-5 units higher. On the contrary, the risk spillover value accepted by Guotai Junan Securities Co., Ltd. is higher, which is about 2 units higher than the risk spillover value spread by Guotai Junan Securities Co., Ltd. We also find that there is no direct relationship between the size of the company's market value and the risk spillover effect between securities companies; for example, as the two securities companies with the largest market value in the industry, China Securities Co., Ltd. and Citic Securities Co., Ltd. receive the ΔCoVaR value from other securities companies at a normal level, or even smaller.
In order to further study the change trend of risk spillover value among securities companies, we draw the time series diagram of systemic risk spillover effect of securities companies (see Figures 1 to 10).
It can be seen from Figures 1 to 10 that the systemic risk spillover effects of different securities companies are quite different, which verifies the analysis results in Table 3. e change trend of risk spillover among securities companies is basically the same, and the abnormal risk spillover mainly occurs in three periods, namely, March 2019 to May 2019, December 2019 to January 2020, and end June to early July 2020.
We compare the trend of the SSE Index (000001. SH), which represents the stock market in China. We can discover that the SSE Index has experienced significant fluctuations in these three periods. From March to May 2019, the market trading sentiment was hot, and the SSE Index rose from about 2600 points to about 3000 points. From December 2019 to January 2020, affected by the epidemic situation, investors did not have confidence in the next market, and the SSE Index plummeted on the first day after the Spring Festival holiday. During this period, the systemic risk spillover phenomenon among securities companies has risen dramatically but is decreased rapidly and maintained at a normal level. From the end of June to the beginning of July, 2020, the stock market entered a hot market again. e SSE index rose from about 3000 points to about 3300 points. e systemic risk spillover effect between securities companies increased rapidly but then quickly dropped to a normal level.

Analysis of Systemic Risk Spillover Effect between Securities
Companies and Financial Industry. Systemic risk spillover effect exists not only between securities companies but also between securities companies and the whole financial industry [32]. Using the calculation method of ΔCoVaR between securities companies, we calculate the risk spillover effect between 10 sample securities companies and securities industry (Tables 8 and 9).
It can be seen from Tables 8 and 9 that there is a great asymmetry in the risk spillover between 10 securities

Analysis of Systemic Risk Spillover Effect between Securities Industry and Financial Industry.
e lack of liquidity in the financial sector will enhance the contagion of the crisis in financial institutions so that the systemic financial risk may have a significant cross sectoral effect [32,34].
Using the calculation method of systemic risk spillover effect among securities companies, we calculate the ΔCoVaR between All Share Brokerage Index, All Share Financials Index, All Share Insurance Index, and CSI Banks Index. e results are shown in Table 10.
According to Table 10, we can see that the value of risk they absorb and release is very different.      Complexity e securities industry receives relatively small risk spillovers from the other three indexes, generally around −0.7. Among them, the banking industry has minimal risk spillovers to the securities industry, which is −0.1638, followed by the entire financial industry, and the insurance industry has the largest risk spillover, which is −0.8112. However, the risk spillover value released by the securities industry is very high, generally around −4. Among them, the securities industry releases the largest risk spillover value to the entire financial industry, which is −4.3827, the risk spillover value to the banking industry is the second, and the risk spillover value to the insurance industry is the smallest [35].
In order to directly depict the risk spillover effects between industries, this paper draws the time series diagrams of the two-way risk spillover effects between the securities industry and the financial subsectors (Figures 11 and 12).
It can be seen from Figures 11 and 12 that the change trend of risk spillover effect between indexes is similar to that between securities companies, and there are also large fluctuations in the following three periods, namely, March 2019 to May 2019, December 2019 to January 2020, and end June to early July 2020.
In addition, we found that the securities department received less risk but spread greater risk to the system.

Conclusions and Suggestions
is paper selects the top 10 securities companies by market value in China, All Share Brokerage Index, All Share Financials Index, All Share Insurance Index, and CSI Banks Index as the research objects, uses the DCC-GJR-GARCH model to calculate the risk spillover values among the research objects, and draws the sequence chart of systemic risk spillover according to the results, so as to explore the change rules and reasons of systemic risk spillover effect of Chinese securities companies. e conclusions are as follows: (1) there are risk spillover effects among 10 securities companies which are asymmetric and bidirectional and highly volatile in a short period of time; (2) the spillover effect of systematic risk of securities companies has no directly relationship with the market value of securities companies but has a strong relationship with the stock market; (3) there are risk spillover effects between the sample securities companies and the four major indexes, and the difference of these effect values is obvious; (4) the risk spillover effect of securities industry on other financial sectors is great, but on the contrary, it is very small.
Based on this, we put forward the following suggestions: (1) implement differentiated risk supervision measures for different types of securities companies to improve the flexibility of risk supervision; (2) steadily advance the reform of the securities industry and stabilize the fluctuation of market trading sentiment; (3) control the scale and quantity of securities companies and reduce the possibility of systemic risk infection; (4) promote the accuracy of systematic risk measurement and enhance the ability of systematic risk identification.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.