Multivariate Time-Varying GH Copula GARCH Model and Its Application in the Financial Market Risk Measurement

Taking full advantage of the strengths of G-H distribution, Copula function, and GARCH model in depicting the return distribution of financial asset, we construct the multivariate time-varyingG-HCopula GARCHmodel which can comprehensively describe “asymmetric, leptokurtic, and heavy-tail” characteristics, the time-varying volatility characteristics, and the extreme-tail dependence characteristics of financial asset return. Based on the conditional maximum likelihood estimator and IFM method, we propose the estimation algorithm of model parameters. Using the quantile function and simulation method, we propose the calculation algorithm of VaR on the basis of this model. To apply this model on studying a real financial market risk, we select the SSCI (China), HSI (Hong Kong, China), TAIEX (Taiwan, China), and SP500 (USA) from January 3, 2000, to June 18, 2010, as the samples to estimate the model parameters and to measure the VaRs of various index risk portfolios under different confidence levels empirically. The results of the application example are in line with the actual situation and the risk diversification theory of portfolio. To a certain extent, these results also justify the feasibility and effectiveness of the multivariate time-varyingG-H Copula GARCHmodel in depicting the return distribution of financial assets.


Introduction
Financial market risk has always been one of the hottest topics in the field of financial investment, and many financial researchers put forward many different financial market risk measurement methods.Among them, Value-at-Risk (VaR) management technology is an assessment and measurement method of financial risk that has risen in recent years, playing an increasingly important role in the risk management and investment decision.It has been widely adopted by the major banks, nonbank financial intermediaries, corporations, and financial regulators in the world and has become the standard of risk measurement and risk management in financial industry.Accurate calculation of VaR is one of the keys to estimate the probability distribution of future return on financial assets.Usually, it is assumed that financial asset returns are independent of each other and obey the normal distribution in the calculation of VaR, but the movement of financial asset return in the financial market is extremely complex.The return of all kinds of financial assets usually does not satisfy the normal distribution hypothesis.However, it often shows "asymmetric, leptokurtic, and heavy-tail" characteristics [1][2][3][4].At the same time, various financial asset returns do not satisfy the multivariate normal distribution hypothesis and present the extreme-tail dependence.On this occasion, a large error would be made by using normal distribution to fit financial asset return, and the estimation of the VaR may be overestimated or underestimated.To solve this problem, many scholars have proposed a lot of leptokurtic and heavy-tail distributions in recent years, such as the logistic distribution, Student's -distribution, and the - distribution.The logistic distribution and Student's distribution can comprehensively describe the leptokurtic characteristics of financial asset return series, but they could not make a good explanation for the heavy-tail characteristics of financial asset return series [5].The - distribution can comprehensively describe the asymmetric, leptokurtic, and heavy-tail characteristics of financial asset return series and it has a good fitting of the univariate unconditional return distribution of some financial assets; however, it could not 2 Mathematical Problems in Engineering reflect the time-varying volatility characteristics of financial asset return and the extreme-tail dependence characteristics of various financial assets return [6,7].Meanwhile, Copula function can connect the joint distribution and the marginal distribution of multiple random variables to construct flexible multivariate distribution functions, which can be used to measure the extreme-tail dependence of multiple financial assets return.GARCH model can comprehensively describe the time-varying volatility characteristics of financial asset return.Therefore, building the multivariate time-varying - Copula GARCH Model by combining the - distribution with Copula function and GARCH model can not only comprehensively describe the "asymmetric, leptokurtic, and heavy-tail" characteristics, the time-varying volatility, and extreme-tail dependence characteristics of the financial asset return and make measurement of VaR more accurately, but also enrich and expand the risk measurement theory and method of financial market theoretically and improve the risk control ability of investors, corporations, financial institutions, and policy authorities and reduce their unnecessary losses in practice.
The - distribution, Copula function, and GARCH model, as well as the financial risk measurement model based on them, have been researched in the existing relevant literatures.A lot of innovative research results with reference value have been brought out.Zhu and Pan [8] proposed three kinds of - VaR methods based on the portfolio gains, losses, and extreme losses according to the statistical characteristics of - distribution.Their empirical results showed that this method is superior to the commonly used delta-normal method.Kuester et al. [9] believed that the - distribution can describe skewness and kurtosis of the financial asset return simultaneously and it plays a very important role of VaR measurement of financial asset return.Degen et al. [10] discussed the application of - distribution in operational risk measurement.Jondeau and Rockinger [11], Rodriguez [12], Fischer et al. [13], and Sun et al. [14] combined time series model with various Copulas functions by using the Sklar theorem to build a lot of highly flexible multivariate time-varying models for risk measurement of portfolios.Liu et al. [15] proposed a GARCH- model with a time-varying coefficient of the risk premium.Their study indicated that the coefficient of the risk premium varies with the time, and even in a mature market the conditional skewness in the return distribution is negatively correlated with the time-varying coefficient of the risk premium.Wen et al. [16] built a -GARCH- model by separating investors' return into gains and losses on the basis of the characteristics of investors' risk preference.They found that investors become risk averse when they gain and risk-seeking when they lose, which effectively explains the inconsistent risk-return relationship.And the degrees of investors' risk aversion and risk-seeking are both in direct proportion to the value of gains and losses, respectively.Wen et al. [17] adopted aggregative indices of 14 representative stocks around the world as samples and established a TVRA-GARCH- model to investigate the influence of prior gains and losses on current risk attitude.The empirical results indicated that the prior gains increase people's current willingness to take risk asset at the whole market level.Huang et al. [18] combined Student's -marginal distribution with Archimedean Copula functions to build the conditional Copula GARCH model.They used this model to estimate the VaR of portfolios.Ghorbel and Trabelsi [19] built the conditional extremum Copula GARCH model by using extreme value theory (EVT) and measured the risk of financial asset according to this model.Chollete et al. [20] used multivariate regime-switching Copula function to build international financial asset return model and accordingly put forward the VaR calculation method.Huggenberger and Klett [21] proposed a measurement model of multivariate risk asset return VaR based on - distribution and Copula function.They used DAX30 (Germany), FTSE100 (UK), and CAC40 (French) from January 2000 to May 2010 as samples to test empirically.Wang et al. [22] applied the Gumbel Copula function in multivariate Archimedean Copula functions family to construct the joint distribution function which can describe the actual distribution and the correlation of various financial asset returns.They also used the Monte Carlo simulation technology to analyze the portfolios VaR and its composition under different confidence levels.The result showed that using the multidimensional Gumbel Copula function to construct the risk measurement model of financial asset can make the assets chosen by investors more robust, and it can also help investors to diversify and control the overall risk of the portfolios.Dai and Wen [23] proposed a computationally tractable robust optimization method for minimizing the CVaR of a portfolio under a general affine data perturbation uncertainty set.And they presented some numerical experiments with real market data to illustrate the behavior of robust optimization model.Liu et al. [24] proposed a pricing model for convertible bonds based on the utility-indifference method and got access to the empirical results by use of Information Technology.Furthermore, using the proposed theoretical model, they presented an empirical pricing study of China's market.They found that the theoretical prices are higher than the actual market prices 0.24-4.58%and the utilityindifference prices are better than the Black-Scholes (B-S) prices.
Based on the aforementioned analyses, the VaR is still the mainstream measurement method of financial market risk.In order to achieve the purpose of measuring VaR more precisely, it has been the hot issue of existing research literatures to construct the distribution functions as comprehensive as possible to describe the "asymmetric, leptokurtic, and heavy-tail" characteristics, the time-varying volatility characteristics, and the extreme-tail dependence characteristics of financial asset return through a variety of mathematical methods.However, in the process of constructing the return distribution model and measuring VaR of financial asset, existing results only grasp some characteristics of financial asset return distribution.They are not able to comprehensively describe the "asymmetric, leptokurtic, and heavy-tail" characteristics, the time-varying volatility characteristics, and the extreme-tail dependence characteristics of the financial asset return.The rationality and accuracy of VaR calculated based on the existing distribution models have a large space for further improvement.
In this paper, we would take full advantage of the strengths of - distribution, Copula function, and GARCH model in depicting the return distribution of financial asset to build multivariate time-varying - Copula GARCH model which can simultaneously describe "asymmetric, leptokurtic, and heavy-tail" characteristics, the time-varying volatility characteristics, and the extreme-tail dependence characteristics of financial asset return and propose the estimation method of model parameters and the calculation algorithm of VaR.Then, this paper selects the SSCI (China), HSI (Hong Kong, China), TAIEX (Taiwan, China), and SP500 (USA) from January 3, 2000, to June 18, 2010, as samples to estimate the parameters and calculate the VaRs of various index portfolios under different confidence levels.

𝐺-𝐻 Distribution, Copula Function,
and Its Tail Dependence Index where  controls the skewness of  distribution.When  → 0,   () →  and  distribution tends to be symmetric.With the increase of the absolute value of , the degree of asymmetry increases.Changing the sign of  can change the asymmetric direction of  distribution, but it does not change its degree of asymmetry.

𝐻 Distribution.
Assuming that random variable  obeys the standard normal distribution and ℎ is a real number, then the random variable distribution stretches the tail of the standard normal distribution.ℎ controls the tail heaviness of  distribution.The larger the ℎ is, the heavier the tail is.Because  ℎ () is an even function,  distribution is symmetric.But the heaviness of its tail changes compared to the standard normal distribution.

𝐺-𝐻 Distribution.
The random variable  ,ℎ can be obtained by introducing both functions   () and  ℎ () to revise standard normal random variable .Consider Then,  ,ℎ can be obtained through linear transformation of  ,ℎ .Consider The distribution of the random variable  ,ℎ obeys the - distribution., , , and ℎ are real numbers. describes the asymmetry of - distribution, and ℎ describes the heavy-tail characteristics of - distribution.Obviously, (3) is a special form of (4).The random variable  ,ℎ in (3) is the random variable of - distribution after central standardization.

Copula Function and Its
Tail Dependence Index.Assuming that marginal distribution of random vector   =   (  ) ( = 1, 2, . . ., ) obeys uniform distribution (0, 1), according to the Sklar theorem, the joint distribution function of dimensional random vectors ( 1 , . . .,   ) can be represented as the following formula: where  is the Copula function of , which is a hypercube [0, 1]  multivariate density function defined on dimensional space R  .If the marginal distribution is continuous, there is a unique Copula function .Then On the contrary, given -dimensional Copula function ( 1 , . . .,   ) and its marginal distribution function  1 ( 1 ), . . .,   (  ), the density function of -dimensional joint distribution function is If   (  ) is the edge density, ( 1 , . . .,   ) denotes Copula density derived from (6).Thus, Since the joint distribution function of random variables defines the correlation among its components, Copula function determines the dependent structure among random variables uniquely.The upper tail index   and lower tail index   of tail dependence indicators can be defined as follows: According to Nelsen [25], Gauss Copula function generated by multivariate normal distribution function whose correlation matrix is R can be represented as follows: where u = ( 1 , . . .,   ) and Φ −1 is the inverse function of single normal distribution.Because the Gauss Copula function does not have the characteristics of tail dependence, we often use the -Copula function whose degree of freedom is  and correlation matrix is R to measure tail dependence structure of risk asset in empirical analysis; that is, where  −1  is the inverse function of simple standard Student's -distribution whose degree of freedom is .When  → ∞, -Copula function degenerates to Gauss Copula function.Its tail index   =   = 0; that is, the tail is independent.The tail index of -Copula function is where  +1 is simple standard Student's -distribution whose degree of freedom is  + 1. Considering that the innovation impacts on the price of risk asset in varying degrees at different times,  and  should have time-varying characteristics.For this reason, tail index also has the same characteristics.

Multivariate Time-Varying 𝐺-𝐻 Copula GARCH Model
Let r  = ( 1, , . . .,  , ) denote return time series of  risk assets.The prior information set before time  is where where the parameters satisfy the conditions   ,   ,   > 0 and   +   < 1.These parameters can ensure the stability of conditional volatility series.The innovation series {z  } obey - distribution whose parameter is (, ℎ) in ( 4).But in order to simplify the analysis, we only consider - distribution after central standardization given by (3) and its density function is written as    (  ).The Copula function (⋅ | I −1 ) is given by ( 11) and its density (⋅) can be represented as the following time-varying -Copula function whose degree of freedom is : where   ,  denotes the multivariate Student's -distribution whose degree of freedom is  and time-varying correlation matrix is   = ( ,, ) × , and The joint density function of  risk assets return is −1 where x  = ( 1, , . . .,  , ) and  , =  −1  ( V  ( , )).Then the likelihood function of overall samples is where  = {(  ,   ,   ,   , V  )  =1 , , ,   , } and x  = ( 1, , . . .,  , ).The value of correlation matrix   is similar to the timevarying correlation matrix of multivariate Copula GARCH model proposed by Jondeau and Rockinger [11].That is,   satisfies the following evolution equation where 0 ≤ ,  ≤

Parameter Estimation Algorithm of the Multivariate Time-Varying 𝐺-𝐻 Copula GARCH Model
On the basis of Huggenberger and Klett [21], this section will use dynamic correlation matrix   instead of static correlation matrix in the multidimensional discrete-time stochastic process to estimate the parameters of multivariate timevarying - Copula GARCH model established in Section 3.
Assuming that Θ denotes the parameter space defined by the model and (r 1 , r 2 , . . ., r  ) denotes the log return samples of -dimensional risk asset which is generated by multivariate conditional density function    |I −1 (r  | I −1 ,  0 ), where  0 ∈ Θ, I −1 is  algebra of time  − 1 and before, the maximum likelihood estimation of parameter vector  can be calculated by the following equation: where The probability density function and distribution function can be obtained in the process of model built in Section 3. Using the IFM method proposed by Joe [26], we can convert (22) into an optimization problem.Therefore, we need to divide the parameter vector  into two subparameter vectors   and   : that is  = (  ,   ), where   = (  1 , . . .,    ),    is the parameter vector of th marginal distribution, and   is the parameter vector of Copula function.Because IFM method is a two-step likelihood estimation method, the model parameters should be estimated through the following two steps.
Step 1. Solving the maximum likelihood estimator of the parameter vector of each risk asset return, This means that we need to estimate parameters vector θ  of  distributions continuously.
Step 2. Taking each θ  into the likelihood equation ( 22 In the maximum likelihood estimation, we need to use the derivative function of the density function of - marginal distribution with respect to the component of parameter vector.Because the density function of - marginal distribution is very complex, this paper uses the implicit function differentiation rule to take its derivative.The estimator θ2 of parameter vector  obtained by the above-mentioned two-step method obeys normal distribution consistently and asymptotically under the standard regularity conditions proposed in Huggenberger and Klett [21], Joe [26], and Patton [27]; that is, where Because the matrixes Σ and Ω can be estimated by the estimated parameter vector consistently, Thus ( 26) can be used to calculate the standard deviation of the estimator θ2, .

VaR Algorithm Based on the Multivariate Time-Varying 𝐺-𝐻 Copula GARCH Model
After estimating the parameters of multivariate time-varying - Copula GARCH model, VaR of the risk portfolio can be measured.VaR of risk portfolio indicates the expected maximum losses of risk portfolio held by investors within a given confidence level and a certain period of time.Assuming that r  = ( 1, , . . .,  , ) ( = 1, 2, . . ., ) are the return samples of  risk assets which satisfy the multivariate time-varying - Copula GARCH (1, 1) model in Section 3 and ∑  =1    , is the portfolio of  risk assets in which the weight of the risk asset  is   ( = 1, 2, . . ., ) that can be less than 0 because of permitting short-purchasing and short-selling the risk assets and meet ∑  =1   = 1, the VaR of risk portfolio under confidence level  at time  should satisfy Pr(∑  =1    , ≤ VaR  ) = .The confidence level  can reflect the different risk preferences of investors or financial institutions to a certain extent.Choosing a larger confidence level means that investors or financial institutions have greater risk aversion, and they hope to get a forecast result with larger probability.
Although the conditional distributions of  1,+1 ,  2,+1 , . . .,  ,+1 can be calculated through the known marginal distributions, it is very difficult to calculate quantile from time-varying Copula density function, and it is adverse to measure and calculate the VaR of risk portfolio.Therefore, this paper measures the dynamic risk of portfolio and its estimation value approximately through simulating - Copula GARCH model.Based on the parameters  () of the sample, the return series of risk assets {[ (,) 1,1+ , . . .,  (,) ,1+ ],  = 1, . . ., } and the one-step measurement and estimation values of VaR of their portfolios can be obtained through estimating (ℎ ()  1,+1 , . . ., ℎ (+1) ,+1 ) according to the volatility equation of - Copula GARCH model and calculating  () +1 by using Copula dynamic evolution equation and then repeating the following algorithm for  times ( ≥ 3).

Application of the Multivariate Time-
Varying - Copula GARCH Model The moment estimation results of the daily log return of SSCI, HSI, TAIEX, and SP500 are shown in Table 1.
Table 1 shows that the skewness of daily log returns of SSCI, HSI, TAIEX, and SP500 is less than 0 and their kurtosis is much larger than that of standard normal distribution which is equal to 3.These results demonstrate that the daily log returns of these indices have the right skew and leptokurtic characteristics.Therefore, it is appropriate to fit the daily log return of SSCI, HSI, TAIEX, and SP500 by applying - distribution which has leptokurtic, heavy-tail characteristics, and it is reasonable to apply the multivariate     7. From Table 7, the following results can be obtained: (1) The inequalities VaR (SSCI) < VaR (HSI) < VaR (SP500) < VaR (TAIEX) can be satisfied for any confidence level.It shows that the risk of extreme losses in Shanghai stock market is higher than that in Hong Kong stock market, Taiwan stock market, and USA stock market.This measurement result is in line with the actual situation that the maturity of Shanghai stock market is far lower than that of Hong Kong stock market, Taiwan stock market, and USA stock market.
(2) For any confidence level, the extreme losses risk of the investors who equally allocate their total assets among SSCI, HSI, TAIEX, and SP500 is lower than that of the investors who put their total assets into one index asset.The extreme losses risk of the investors increases with the concentration of risk asset in the index portfolios.This measurement result is consistent with the risk diversification theory of portfolio.

Conclusion
Considering the "asymmetric, leptokurtic, and heavy-tail" characteristics, the time-varying volatility characteristics, and extreme-tail dependence characteristics of financial asset return, this paper combined the - distribution, Copula function, and GARCH model to construct a multivariate time-varying - Copula GARCH model which can comprehensively describe the "asymmetric, leptokurtic,   Note: the ratio of index portfolio is ranked by the sequence of SSCI : HSI : TAIEX : SP500 in Table 7.
and heavy-tail" characteristics, the time-varying volatility characteristics, and extreme-tail dependence characteristics of financial asset return.It proposed the parameter estimation algorithm of the multivariate time-varying - Copula GARCH model by using condition maximum likelihood method and IFM two-step method.An algorithm was constructed to calculate VaR by using the quantile function and the simulation method based on - Copula GARCH model.In addition, this paper selected the daily log return of SSCI (China), HSI (Hong Kong, China), TAIEX (Taiwan, China), and SP500 (USA) from January 3, 2000, to June 18, 2010, as samples to estimate the parameters of the multivariate time-varying - Copula GARCH model, and it also estimated the VaR for various index risk asset portfolios under different confidence levels.The research results showed that the multivariate time-varying - Copula GARCH model constructed in this paper could reasonably estimate and measure the extreme losses of risk portfolios in financial market, and the measurement results were in line with the actual situation of stock market and the risk diversification theory of portfolio.The achievement of this paper provided a practical and effective method for measuring the extreme losses of financial market.

Table 1 :
Moment estimation results of the daily log return of SSCI, HSI, TAIEX, and SP500.
USA and China, as the most developed capitalism country and the largest developing country in the world, respectively, rank top two of the world economy.Their stock markets should have strong representation in the world.At the same time, due to the historical reasons, there exist several regions with different political systems such as Mainland China, Hong Kong, Taiwan, and Macau in Greater China.Macau is similar to Hong Kong on the whole.For the above-mentioned reasons, this paper selects the SSCI (China), HSI (Hong Kong, China), TAIEX (Taiwan, China), and SP500 (USA) from January 3, 2000, to June 18, 2010, as data samples to estimate the VaR of various index portfolios under different confidence levels by using the multivariate time-varying - Copula GARCH model.The data comes from Yahoo Finance website: http://finance.yahoo.com/.

Table 2 :
Parameter estimates of the four-variate time-varying G-H Copula GARCH (1, 1) model based on SSCI index risk asset.: * * * in the table denotes that the parameter is significant at 1% level. Note

Table 3 :
Parameter estimates of the four-variate time-varying G-H Copula GARCH (1, 1) model based on HSI index risk asset.: * * * in the table denotes that the parameter is significant at 1% level. Note

Table 4 :
Parameter estimates of the four-variate time-varying G-H Copula GARCH (1, 1) model based on TAIEX index risk asset.* * in the table denotes that the parameter is significant at 1% level.

Table 5 :
Parameter estimates of the four-variate time-varying G-H Copula GARCH (1, 1) model based on SP500 index risk asset.Consider  1 +  1 = 0.92,  2 +  2 = 0.9573,  3 +  3 = 0.9571, and  4 +  4 = 0.9599.This shows that the volatility persistence of Shanghai stock market is the strongest, Taiwan and Hong Kong stock market rank second and third, and the volatility persistence of USA stock market is minimum.It indicates that the investors' expectation of risk compensation in the emerging markets represented by China's stock market is stronger than that in the mature markets represented by the USA's stock market and the price discovery efficiency of innovation in the emerging markets represented by China's stock market is lower than that in the mature markets represented by USA's stock market.Inaddition, the sum of the coefficients and  is very close to 1, which indicates that the impact and shock of innovation on the index volatility of each stock market has a long memory.(2)Thedegree of freedom  = 14.57and the correlation coefficients   of -Copula in Table6show that there exists the strongest correlation between Hong Kong stock market and Taiwan stock market, and the correlation between Shanghai stock market and Hong Kong stock market is also relatively large.
6.3.VaR Measurement Based on the Multivariate Time-Varying - Copula GARCH Model.Based on the multivariate time-varying - Copula GARCH model with SSCI, HSI, TAIEX, and SP500 whose parameters have been estimated, the VaRs of various index portfolios under different confidence levels can be measured.The measurement results are shown in Table

Table 6 :
Estimates of correlation coefficients, time-varying parameters, and degree of freedom of four-variate time-varying G-H Copula GARCH (1, 1) model based on SSCI, HSI, TAIEX, and SP500.: * * * and * * in the table denote that the parameter is significant at 1% and 5% level, respectively. Note

Table 7 :
VaR estimation results based on the four-variate time-varying G-H Copula GARCH model with SSCI, HSI, TAIEX, and SP500.