Measuring systemic risk for financial time series: A dynamic bivariate Dvine model

Graphical abstract The framework of the bivariate time series model.


Introduction
Modeling multivariate time series is an important statistical application in different fields.There are two types of correlations embedded in multivariate time series: serial dependence, which refers to the dependence between different time points of a univariate time series, and cross-sectional dependence, which refers to the dependence between different dimensional time series at the same time point.Accurately measuring these two types of dependence is essential for evaluating the effectiveness of multivariate time series models.
Copula as a multivariate distribution with uniformly distributed marginals can flexibly model the dependence structure of multivariate variables, see Refs.[1-3] for details.The copula function can be used to model two types of dependence between multivariate variables.Currently, the application of the copula function to multivariate time series mainly focuses on cross-sectional dependence.When dealing with multivariate time series models, the popular method is to use the GARCH family models to analyze the serial dependence, and then employ the copula function to fit the correlation between the filtered data (including Refs.[4-7]).Moreover, Chen and Fan [8] summarized the class of multivariate time series models as SCOMDY models and studied the theoretical properties of such models.In the framework of SCOMDY, Patton [9] extended constant copula to time-varying copula by turning the copula parameters into functions of lagged variables and lagged exogenous variables.Giulio and Ergün [10] proposed a dynamic GAS copula model whose iteration over time-varying parameters relies on a scaled score function of predictive model density rather than simply moments of observations.Time-varying copulas have succeeded in describing the dynamic cross-sectional dependence between multidimensional distributions in the past few years.
Modeling serial dependence based on copulas has also attracted attention.Joe [11] proposed a class of stationary first-order Markov models based on bivariate copulas and marginal distributions.Chen and Fan [12] and Beare [13] generalized the class of models to semiparametric stationary Markov models and discussed asymptotic properties of the semiparametric Markov models.Zhao et al. [14] extended the univariate copulabased time series model to higher order by using the D-vine copula.(see Refs. [15-18] for more details of D-vine copula model).Bladt [19] constructed a new class of time series copula models using d-vines and v-transforms.Moreover, copula can also simultaneously capture the serial dependence and crosssectional dependence in a class of vine copula models.Brechmann and Czado [20] , Smith [21] , and Beare and Seo [22] introduced different vine copula models that can account for both types of dependence in multivariate time series.Nagler et al. [23] summarized a stationary vine copula model which can use a more general vine structure to fit multivariate time series including the above three models as special cases.
For financial time series, multivariate GARCH models focus on the linear dependence structure between variables; see, Refs.[24, 25].However, copula functions can capture the nonlinear and asymmetric dependence that exists in multivariate time series.Thus, we design a dynamic bivariate Dvine model in which serial dependence and cross-sectional dependence are modeled based on copula functions.Specific-ally, inspired by the work of Zhao et al [14] , we consider a class of univariate copula-based stationary Markov models (Dvine), in which D-vine copulas are used to capture the serial dependence in univariate time series and marginal distributions are suitable parametric distributions.Subsequently, we propose the dynamic bivariate Dvine model in which the timevarying copula is used to capture the dynamic cross-sectional dependence by connecting the two-dimensional Dvine model.Based on the complicated and dynamic nature of the financial market, the dynamic bivariate Dvine model can provide a flexible and nonlinear copula-based method compared to the existing time series models.

√ T
In the financial field, the ultimate goal of modeling time series is to measure the future financial situation through past observations, such as the conditional quantile of interested financial companies.Therefore, forecasting the future behavior of financial assets has become the main method of risk measurement (e.g., Refs.[10, 26-28]).Given past observations, the corresponding estimator of VaR of financial assets can be constructed via the univariate Dvine model.For our dynamic bivariate Dvine model, we can also provide forecasting methods for CoVaR of the financial system when the financial company is in distress.Moreover, given the asymptotic properties of the parameters estimated by the univariate Dvine model, we can obtain the -consistency and asymptotic normality of the estimators of the conditional quantiles by the Delta method.
Our proposed model enriches the existing literature in three ways.First and foremost, we generalize the existing method of modeling the dependence of time series by proposing the dynamic bivariate Dvine model.Second, we provide an alternative method to predict the conditional quantile based on our proposed model.In the empirical analysis, the backtesting results show that our method can rival and outperform the traditional GARCH method.Third, by taking the two-stage likelihood method, we can obtain the estimators of our model parameters and study the statistical properties of the estimators of the Dvine model.
The rest of our paper is organized as follows.In Section 2, we introduce the dynamic bivariate Dvine model.Section 3 introduces the parameter estimation for our dynamic bivariate Dvine model and estimators of the conditional quantiles.In Section 4, we establish the asymptotic properties of the estimators of the univariate Dvine model.In Section 5, real data applications on financial stock prices show good performance compared to the traditional GARCH method.Finally, we offer some concluding remarks in Section 6.The proofs of theorems are relegated to Supporting Information.

Methodology
In this section, we propose a copula-based time series model which can provide a new tool for capturing joint behavior in the financial market.As we know, the dependence structure that exists in financial time series is always time-varying and non-linear.We build the dynamic bivariate Dvine model to specify the dependence that exists in the bivariate time series .Let denote the joint distribution of conditional on the past information . The dynamic bivariate Dvine model for is where the marginal behavior of is captured by the univariate Dvine model with parameters , and the joint behavior of is characterized by the dynamic copula function with dependence parameter which means the cross-sectional dependence may change over time.We will introduce the specific forms in the following subsections.
The dynamic bivariate Dvine model is sensible and interpretable and links the univariate Dvine model by the dynamic copula function, which accurately reflects the nonlinear, timevarying dependence structure of financial time series.We emphasize the time-varying nature of the cross-sectional dependence in the dynamic bivariate Dvine model, which naturally extends the copula-based time series model in low-dimensional applications.

Univariate Dvine model
We develop a Dvine copula-based model for univariate time series which satisfies the requirement of strict stationarity and Markov properties, .Inspired by the work of Ref. [14], a class of univariate Dvine models is employed to construct the conditional marginal distribution at time conditional on the past observations , .
A time series is called a Dvine( ) copula process if, for all , the jont density of can be written as where is the parametric marginal density with parameter , is the bivariate copula in tree with parameter and are the set of all bivariate copulas parameters.The algorithm to obtain , the conditional distribution of , is provided in Supporting Information S.1.

p
The unconditional marginal distribution of our univariate time series model is estimated parametrically while Ref. [14]  use a non-parametric approach.Utilizing the Dvine( ) model to analyze the serial dependence present in univariate time series can provide a flexible way to analyze univariate time series in a non-linear form without relying on the first or second moment of observations as GARCH-type models do.

Modeling the parameter of cross-sectional dependence
The modeling process for the dependence parameter of the cross-sectional copula is pivotal, and we introduce the GAS method [29] to design time-varying mechanisms for the dependence parameter .Specifically, we assume to be the Student t copula.To capture the dynamic cross-sectional dependence, we update the correlation parameter based on past observations.We transform it as to ensure , where f t is the time-varying parameter which can be updated by the following autoregressive process, where denotes an appropriate function of the past information.Using the GAS framework, the driving mechanism is considered as the scaled score function of the observation density, and its specific form is given by where denotes the scaling matrix.In this paper, we define the scaling matrix as the square root matrix of the inverse of the Fisher information matrix with respect to the parameter and thus , where .The detailed formula of the driving mechanism is in Supporting Information S.2.

{X
t } Together, the dynamic bivariate Dvine model for the time series can be expressed as follows: In this model, the conditional marginal distribution of is established by the univariate Dvine model with the parameters , and the joint dynamic behavior of is characterized by the bivariate copula with the time-varying coefficient parameter .To our best knowledge, the dynamic bivariate Dvine model is the first copula-based time series model that can handle dynamic behavior.We emphasize the time-varying nature of the dynamic bivariate Dvine model, as it is built upon the assumptions of using a time-varying copula for the cross-sectional dependence.This makes it more suitable for financial time series modeling.

Estimation and prediction
In this section, we provide the parameter estimation for our proposed bivariate Dvine model and introduce the prediction method via the proposed dynamic bivariate Dvine model.

Estimation of the proposed model parameters
Based on the model we introduced in the previous section, the parameters that need to be determined are , where is the marginal distribution parameter of the Dvine( ) model, is the parameter set of the bivariate copulas in the Dvine( ) model, , and denotes the parameter for the timevarying copula of the bivariate time series.The estimation of is based on the likelihood approach.Denote as the observations of the bivariate time series, the complete log-likelihood function of is given by where the conditional cdf(cumulative distribution function) can be derived from the Dvine( ) model and .If we want to calculate the maximum likelihood estimates of all parameters of the model simultaneously, the calculation is complicated and inefficient.Joe [30] and Tsukahara [31] discussed the step-wise method in the copulabased models and proved its efficiency and theoretical properties.To improve our computational efficiency, we follow the step-wise maximum likelihood estimation method.The complete log-likelihood function in Eq. ( 6) is decomposed into several parts and optimized in successive steps.
In the first stage, for , the estimators of parameters contained in the Dvine( ) model are obtained using the univariate time series .In this paper, we mainly adopt parametric methods to fit the marginal distributions and transform the observed data into pseudo-observations.In general, we define an appropriate parametric marginal density for the -th univariate time series.The marginal parameters can be calculated by the following equation: Then the pseudo-observations can be calculated by the parametric distribution function , where .The vector can almost be seen as an i.i.d sample from the random vector and be used to estimate the parameters of bivariate copulas in Dvine models.

Given the pseudo-observations
, can be calculated by maximizing the joint log-likelihood of the Dvine( ) model.Let denote the bivariate copula in tree with parameter , the conditional distribution can be calculated by the bivariate copula function in tree .Then, the log-likelihood function of the Dvine( ) model for the pseudo-observations is given by where denotes the conditional distribution of given the intermediate variables and parameter set .Through the tree-by-tree sequential selection method proposed by Ref. [14], we determine the family and parameters of the copula function starting from the first tree, and obtain the estimated parameter of the copula in the second tree based on the result of the first tree.The process continues until an independent copula is selected in a certain tree.Consequently, for , we get where the vector denotes the estimated parameters from the tree.Then the conditional cdf are calculated by the partial derivative or -function of the copula .
θ c In the second stage, the parameter vector of the timevarying copula function can be obtained by the following equation: Compared to the full MLE, the two-stage MLE provides a more efficient way to obtain the parameter estimates of our bivariate time series model.

Prediction
In the financial time series field, tail risk measures have been widely employed to measure financial asset risk.Given the past observations , the dynamic bivariate Dvine model can make predictions for the future tail behavior of .At time , the conditional quantile or conditional VaR of at level given the past information can be expressed as, For exposition, we first give the explicit equation for theth conditional quantile of when using the Dvine(2) model, where denotes the conditional copula function in the tree , and denotes the inverse of the conditional copula function in tree .We provide explicit expressions for conditional quantiles, but this is based on the fact that the copula functions included in the Dvine model all have the inverse of the conditional distribution.Equations for theth conditional quantiles of when the Dvine model is of order 1 or higher are in Supporting Information S.3.
The CoVaR proposed by Adrian and Brunnermeier [27] is a conditional quantile that measures the contribution of financial institutions to systemic risk.To facilitate CoVaR backtesting, Giulio and Ergün [10] modified the A-B CoVaR to the G-E CoVaR, and measured the G-E CoVaR via the GARCH-copula method.Currently, the main methods used to estimate CoVaR include quantile regression and GARCHcopula models; see Refs.[27, 32].In this section, we can also use our model to give an explicit equation for estimating the CoVaR.We denote the -th quantile of conditional on proposed by Ref. [10] as , which can be seen in Eq. (11), As mentioned earlier, we use the time-varying GAS copula function to measure which can be obtained by solving Eq. ( 12), C t where denotes the conditional distribution of given the past information .Then, we can obtain the corresponding by calculating the inverse of the conditional distribution using Eq. ( 17).In Section 5, we present the performance of the prediction method based on our model in an empirical study.

The asymptotic properties of the proposed estimators
In this section, we discuss the consistency and asymptotic normality of our proposed estimators under sufficient conditions.
Assume is a strictly stationary Markov process.In the first stage, for , and are the estimators of the univariate time series model using the -dimensional univariate time series.In the second stage, is the estimator of the time-varying GAS copula given the conditional cumulative distribution function.In what follows, we provide the results on consistency and asymptotic normality of the proposed estimators and .In this paper, denotes the Euclidean norm.

Asymptotic properties of the estimators in univariate time series model
Based on the parametric sequential maximum likelihood estimation method, we can estimate all parameters more simply and efficiently.Yet the theoretical results on the fully parametric MLE seem to be few, and we extend our results in this area.To establish the asymptotic properties of the estimators of the univariate time series model proposed by our paper, we need to introduce some notations.
) p where denotes the unconditional marginal distribution with parameter , and denotes the parameter set of bivariate copulas in the Dvine( ) model.
Under regularity conditions, the parametric sequential MLE is defined as the solution of the following equation: where can be seen as the true parameters when the model is correctly specified.Then we need to impose the following conditions: is the parametric space of the Dvine( ) model, is the compact set of ; if and only if .
is well defined for ; For all , the probability that is Lipschitz continuous at is 1.
is well defined and continuous in .
The function is continuously differentiable with respect to and satisfies H × Θ for any compact set .
is well defined and continuous in .
(C7) β The -mixing coefficients of satisfy , where is the inverse survival function of .
Theorem 1.If the above conditions (C1)-(C3) hold for the univariate time series model, then we have: .
) Theorem 2. If the above conditions (C1)-(C7) hold for the univariate time series model, then we have: and

Asymptotic properties of the conditional quantile estimators
Similarly, one can show the asymptotic results of the estimator of the conditional quantile of under some appropriate conditions.The conditional quantile estimator of given can be evaluated explicitly by the univariate Dvine model.More specifically, we let be the conditional distribution function parametrized by , and we denote the value of our interest for some functional .Given the estimator of , we define .In this paper, is defined as the conditional quantile.
Assume the maps are Frechet differentiable at and with gradient .
(ii) T → ∞ If and the above conditions (C1)-(C7) hold, then The estimator of the conditional quantile converges to the pseudo-true at a rate .Theorem 3 exhibits the consistency and asymptotic normality of such predictions.

&
In this section, we illustrate the methods introduced in the previous sections with real financial datasets.The financial datasets are retrieved from Yahoo Finance and cover the period from January 1, 2015, to May 17, 2022, containing the turbulent period in financial markets caused by the outbreak of COVID-19.The first five datasets are mainly the daily closing prices of Amazon (AMAZ), Boeing (BA), Coca-Cola (COCA), General Motors (GM), and IBM companies stocks which belong to different industry branches.The last dataset is the daily closing price of the S P 500 index (SP500) which can be seen as a proxy for the financial system.There are 1856 trading days in each dataset and we take the log difference of the daily closing price as the empirical data in our analysis.p p p In the process of modeling the whole dataset, we first choose the Dvine( ) model to analyze the serial dependence in the univariate time series.Table 1 shows the descriptive statistics for daily returns.The results indicate that all corresponding empirical distributions have asymmetric and heavytailed characteristics.By comparing different families of distributions, we employ the skewed-t distribution to fit the marginal distribution.After obtaining the pseudo-observations of individual time series, we select the order and the specific copula function for the parametric Dvine( ) model according to the tree-by-tree selection method proposed in Ref. [14].Table 2 exhibits the estimates for the Dvine( ) models fitted to the financial datasets.
Certainly, the conditional cdf can be evaluated by the Dv- Jarque-Bera < 0.001 < 0.001 < 0.001 < 0.001 < 0.001 < 0.001 Mean refers to the average value; Skew refers to the skewness of the empirical distribution; Kurt refers to the kurtosis of the empirical distribution; Std refers to the standard deviation; Jarque-Bera refers to the p-value of the Jarque-Bera test.
should be iid uniform if the fitted models are reasonable.Traditional modeling methods for univariate time series include GARCH-type models, we also fit the standard GARCH model proposed by Bollerslev [4] , the exponential GARCH model proposed by Nelson [33] and the GJR-GARCH model proposed by Glosten et al [34] .Hence, we can compare the performance of the parametric Dvine(p) model with GARCH-type models by testing the uniformity and independence of the conditional cdf.Table 3 shows the p-values of the Kolmogorov-Smirnov test and the Ljung-Box test using 10 lags.From the results of Table 3, we can see that the uniformity and independence of the conditional cdf are not rejected for those financial datasets, and the p-values of our model perform slightly better than the p-values of the GARCH-type models on both types of tests.

&
As for cross-sectional dependence, dynamic and nonlinear correlations between financial companies and the S P 500 index in financial markets can be inferred.Hence, we focus on the time-varying GAS copula functions including the Student's t GAS copula and normal GAS copula.To show the strengths of the time-varying GAS copula functions over other types of copula functions, the time-varying DCC copula and the constant copula are fitted for the cross-sectional dependence.In Table 4, the corresponding statistics of the copula functions are listed.The minimum Akaike information criterion (AIC) values can be obtained by fitting the Student's t GAS copula function.The results are reasonable because there exist certain upper and lower tail correlations between financial companies such as IBM and the financial system SP500.Moreover, the results show that the time-varying GAS functions are better than the constant copula and time-varying DCC copula functions.The Student's t GAS copula function can effectively fit the dependence structure between our financial dataset.For demonstration, the time-varying avercorrelation among the five financial companies and the S&P500 index estimated by the Student's t GAS copula function is plotted in Fig. 1, which reflects a high positive correlation in financial markets.

VaR and CoVaR
VaR and CoVaR as the mainstream approaches for measuring risks in financial markets are evaluated by our proposed models.In our paper, we aim to evaluate the one-day-ahead VaR and CoVaR at the confidence level.We first calculate the for each financial dataset at time .For comparison, we also consider GARCH-type models to calculate VaR for our univariate time series, including the standard GARCH model, the exponential GARCH model, and the p  The bold denotes the optimal AIC, BIC and log-likelihood.
Measuring systemic for financial time series: A dynamic bivariate Dvine model Chen et al.
GJR-GARCH model.Fig. 2 shows the estimation results of VaR for financial time series where a black line represents the original dataset, the VaR calculated by our model is marked with a red line, and lines in other colors represent the estimated VaR based on GARCH-type models.It shows that the estimated VaR calculated by our method can have excellent fitting results when the financial markets are in turbulence.Almost all financial datasets fluctuate sharply during the period from 2020 to 2022 due to the outbreak of COVID-19, yet our method accurately measures the risk value of financial data during periods of high volatility.As a result, the VaR building on our Dvine model is more sensitive to the occurrence of financial risk events.Thus, our method can be attractive to financial risk managers.
n e P B P uc P cc α = 0.05 p Moreover, we evaluate our estimated VaR using the twosided binomial score test, Kupiec [35] test and Christoffersen [36] test.Table 5 shows the results of the backtesting of VaR and compares it with other models.There are 1853 samples for backtesting and the actual number of violations where the true loss exceeds the estimated VaR are listed in column 3 (labeled ).The p-values of the two-sided binomial score test, Kupiec [35] test and Christoffersen [36] test are given in columns 4, 6 and 8 (labeled , and ) and the largest ones are listed in bold.From the results of Table 5, the unconditional coverage property and conditional coverage property are not rejected at the level of significance for all datasets.The number of violations is in line with the expected number of violations under the null hypothesis that the Dvine( ) model is the correct model.In addition, but the pvalues obtained by our method are also slightly better than the results calculated by the GARCH method.
To demonstrate the performance of our model in predicting VaR, the period of the training dataset is set from January 1, 2015 to February 1, 2022, while the period of the testing dataset is set from February 2, 2022 to May 17, 2022.Table 6 shows the results of the predicted VaR obtained by our model.The p-values of the two-sided binomial score test, Kupiec [35] test and Christoffersen [36] test demonstrate the accuracy of the predicted VaR.When we calculate the mean square error (MSE) and correlation between the predicted VaR and the calculated VaR, the results show a small difference and high correlation between the predicted VaR and calculated VaR.The backtesting results prove that our predicted VaR is reliable.
CoVaR s| j q,t t j 5% Subsequently, we estimate the of the financial system at time conditional on the financial company being in financial distress (when the company's return falls below its VaR).To compare the performance of our model, we also fit the eGARCH-t GAS copula model to evaluate CoVaR.The results of the estimated CoVaR of SP500 conditional on the financial company being in financial distress are shown in Fig. 3.
By analyzing the CoVaR in Fig. 3, we find that the CoVaR sharply plunge during the COVID-19 pandemic period.The values of VaR and CoVaR obtained by our model are sensitive to the high volatility period.Its risk contagion is more significant in crisis periods using our method compared to the GARCH method.The finding reflects that companies in different industries are closely linked in the international market and that different companies will have a risk contagion effect on the financial system.The magnitude is different but the   trend is generally similar.To better interpret the risk spillover effect on the financial system conditional on financial company being financial distress, we define as Eq.(15), ∆CoVaR s| j q,t = CoVaR s| j q|0.05,t − CoVaR s| j q|0.5,t .
CoVaR s| j q|0.5,t q where denotes the -th quantile of the financial j system conditional on financial company being its median state.
By evaluating the measure, we can more intuitively determine the risk spillover effect of the financial company on the financial system and its magnitude.The estimation results of the value of are shown in Fig. 4. The estimated can reflect the impact of the financial company on the S P500 index in a similar way to the previ- ous volatility.Meanwhile, the results suggest that IBM and GM companies can have a greater influence on the financial system while other companies have a relatively small influence on the US financial market.The estimation results demonstrate that some industrial companies may be more prone to contribute to financial risk contagion among financial markets than Internet companies in crisis periods.
To ensure the prediction accuracy under our model, we also backtest our CoVaR according to the method proposed in Ref. [10].For each financial company , we obtain its observations when the actual loss at time is greater than its estimated .Then, in those observations, we can construct a hit sequence again which compares the estimated at time with the actual loss of the financial system , CoVaR s| j q,t t j

5%
The hit sequence returns a 1 if the actual loss of the financial system is greater than its estimated on day in which financial company is in financial distress.Then we can backtest our CoVaR measure by using the Kupiec [35] and Christoffersen [36] tests.Table 7 reports the test results for Co-VaR.As show in Table 7, our CoVaR measure can satisfy the unconditional and conditional coverage properties at the confidence level.The backtest results of our method are also quite good compared to traditional eGARCH methods.This shows that the calculation of our model is accurate and efficient.Therefore, using the model constructed in this paper enables accurate measurement of risk spillovers in extreme risk situations.

5%
Table 8 shows the results of the predicted CoVaR obtained by our model.As we can see from the table, our predicted CoVaR can satisfy the unconditional coverage property and conditional coverage property at the confidence level.The results of MSE and correlation also show the similarity between the predicted CoVaR and calculated CoVaR.Therefore, we can also obtain reliable results by using our model to estimate the conditional quantile.

Conclusions
Our paper proposes a dynamic bivariate Dvine model to deal p p with the dependence that exists in bivariate time series.We apply a class of stationary Dvine( ) models to fit serial dependence in univariate time series and then model cross-sectional dependence via a time-varying copula model.To improve the computational efficiency, we use a step-wise estimation method to obtain the parameter estimates we need and provide the theoretical proof of the asymptotic properties of the parameter estimates of the Dvine( ) model under sufficient conditions.Our method is a new extension to the modeling of time series, and we consider the copula-based models to forecast the VaR and CoVaR of financial time series.
Using applications to the financial stock prices, our VaR and CoVaR can provide predictions of tail risk for financial time series, especially during periods of high volatility caused by the outbreak of the novel coronavirus.Backtesting the conditional quantiles calculated by our method proves the reliability of the results and the improvement over the previous methods.Further research is needed to produce better algorithms to improve the efficiency of models on higher dimensional data.

Fig. 1 .
Fig. 1.The time-varying average correlation among the financial companies and system.

Fig. 3 .
Fig. 3.The CoVaR of SP500 conditional on the financial company AMAZ, BA, COCA, GM or IBM being in financial distress.
Measuring systemic risk for financial time series: A dynamic bivariate Dvine model Chen et al.

Table 1 .
Descriptive statistics for daily returns.

Table 2 .
Summary of parameter estimates for the Dvine( ) models. p

Table 3 .
-values for tests for uniformity and independence of u-PIT val- LBrefers to the p-value of the Kolmogorov-Smirnov test and refer to the p-value of the Ljung-Box test using 10 lags.

Table 4 .
Estimation results of time-varying copula function.

Table 5 .
Backtesting of the calculated VaR at level .

Table 6 .
The performance of the predicted VaR at level .test of conditional coverage and independence (cc), MSE denotes the mean square error between predicted VaR and calculated VaR, and Corr denotes the correlation between the predicted VaR and the calculated VaR.

Table 8 .
The performance of the predicted CoVaR at level .values and statistics of the LR test of unconditional coverage (uc), the subscript denotes the LR test of conditional coverage and independence (cc), MSE denotes the mean square error between predicted CoVaR and calculated CoVaR, and Corr denotes the correlation between the predicted CoVaR and the calculated CoVaR.Measuring systemic risk for financial time series: A dynamic bivariate Dvine model Chen et al.