logo资料库

论文研究 - 基于GARCH模型和极值理论的中国股票市场风险度量。.pdf

第1页 / 共13页
第2页 / 共13页
第3页 / 共13页
第4页 / 共13页
第5页 / 共13页
第6页 / 共13页
第7页 / 共13页
第8页 / 共13页
资料共13页,剩余部分请下载后查看
Risk Measurement of Chinese Stock Market Based on GARCH Model and Extreme Value Theory
Abstract
Keywords
1. Introduction
2. Empirical Methodology
2.1. VaR
2.2. GARCH Model
2.3. POT Model
3. Empirical Study
3.1. Data Description
3.2. GARCH Model Establishment
3.3. POT Model Establishment
3.4. VaR Calculation and Failure Rate Test
4. Conclusions
Fund
Conflicts of Interest
References
Open Journal of Business and Management, 2019, 7, 963-975 http://www.scirp.org/journal/ojbm ISSN Online: 2329-3292 ISSN Print: 2329-3284 Risk Measurement of Chinese Stock Market Based on GARCH Model and Extreme Value Theory Shixue Du*, Guoqiang Tang, Shijun Li School of Science, Guilin University of Technology, Guilin, China How to cite this paper: Du, S.X., Tang, G.Q. and Li, S.J. (2019) Risk Measurement of Chinese Stock Market Based on GARCH Model and Extreme Value Theory. Open Journal of Business and Management, 7, 963-975. https://doi.org/10.4236/ojbm.2019.72065 Received: March 26, 2019 Accepted: April 22, 2019 Published: April 25, 2019 Copyright © 2019 by author(s) and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/ Open Access Abstract The risk of the stock market has always been a hot issue for investors. This paper selects the daily closing price data of the Shanghai and Shenzhen 300 Index of the Shanghai stock exchange, the Shenzhen 300 Index of Shenzhen stock exchange, the Hang Seng Index of the Hong Kong stock exchange mar- ket and Taiwan Weighted Index of the Taiwan stock market and calculates logarithm and difference. The GARCH model is combined with the POT model of extreme value theory to measure the risk. Comparing the failure rates at the three significance levels of 0.05, 0.025 and 0.01, the failure rates are close to the level of significance, which conduct that the GARCH-POT model can measure the risk of Chinese stock market well. Keywords GARCH Model, Extreme Value Theory, VaR 1. Introduction The risk of the stock market has been always widely concerned by researchers and investors. In the past three decades, financial market volatility has become obvious because of economic globalization and financial innovation, which makes financial risk management a necessary tool and capability for business enterprises and financial institutions to operate and manage. Value at Risk (VaR) is a widely used method of measuring financial risk. Domestic and foreign scholars have made fruitful research on the risk measurement of the stock market. Domestic and foreign scholars have done a lot of research on the measurement of risk value. The most traditional VaR model is based on the premise that the hypothesis rate of return follows the normal distribution, while the financial time DOI: 10.4236/ojbm.2019.72065 Apr. 25, 2019 963 Open Journal of Business and Management
S. X. Du et al. DOI: 10.4236/ojbm.2019.72065 series often has volatility clustering and non-normality [1]. In order to better solve the non-normality and volatility aggregation of financial time series, Engle (1982) first proposed an autoregressive conditional heteroskedastic (ARCH) model, which is considered to be a function of past error. The ARCH model is a good indicator of the fluctuating agglomeration of financial markets [2]. On the basis of the ARCH model, Bollerslev (1986) extended it to the generalized autoregressive conditional heteroskedasticity (GARCH) model, and considered that the condi- tional variance is not only a function of the past error, but also a function of the conditional variance of the lag. The GARCH model not only reveals the “fluc- tuating agglomeration” characteristics of financial markets, but also reflects the “thick tail” characteristics. Therefore, the GARCH model can be used to charac- terize financial time series that are thicker than the tail of a normal distribution [3]. Later, in order to reveal the asymmetry and leverage effect of financial data, Nelson (1991) proposed an index GARCH (EGARCH) model, which introduces parameters γ in the conditional variance equation to distinguish between posi- tive and negative external impacts on the price of financial products. Since then, the GARCH family model has continued to grow and develop, such as: TGARCH model, IGARCH model, GJRGARCH model [4]. Gong Rui (2005) compared the GARCH model, EGARCH model and PARCH model of the Shanghai Composite Index, Shanghai Stock Exchange 180 Index and Shenzhen Composite Index under the normal distribution, t distribution and GED distribution, and calculated the corresponding VaR values. It is considered that the VaR value under the t dis- tribution is too conservative, and the GED distribution is more accurate to de- scribe the characteristics of the yield [5]. Dong ran (2016) measured the risks in the national debt market, believed that the Shanghai municipal government bond index had leverage effect, and suggested investors to accurately judge the market information before investing [6]. Extreme events are small probability events, once happen these small probability events can cause huge losses. Extreme Value Theory (EVT) is a theory for mod- eling and analyzing tail data. It is widely used in many fields such as technology engineering, environment, and geological disasters. The two models commonly used in extreme value theory are BMM model and POT model. The BMM model first divides the financial time series into several sub-intervals according to a certain standard, models the maximum (minimum) value of each set of data, and obtains the parameter estimation according to the maximum likelihood estima- tion or other estimation methods. However, the traditional BMM model relies on the choice of the length of the subinterval, and it is easy to ignore some valuable data. At present, the most widely used is the POT model, which can make full use of limited observations. The POT model models all observations that exceed a given threshold, which approximates the Generalized Pareto Distribution. There are many methods for selecting the threshold, such as Hill estimation method, kurtosis method, and over-expectation function graph method. These methods have their own advantages and disadvantages, and there is no unified selection 964 Open Journal of Business and Management
S. X. Du et al. standard. Mcneil (2000) combined the GARCH family model with extreme value theory for the first time and constructed the GARCH-EVT model to predict the VaR of financial products [7]. Tang Yong (2012) used the kernel fitting goodness statistical method and the average excess distribution function graph to select the threshold respectively, and fitted the POT model to the low-frequency data and high-frequency data distribution, and considered that the kernel fitting goodness statistical method was more effective to select the threshold comparison [8]. Zhang Hu (2016) used the extreme value theory to study the risk value of the return rate of stock market [9]. Wang Miao (2017) used extreme value theory to study the fluctuation of the Shanghai stock market and the conditional risk value (CVaR) of the characteristic of thick tail [10]. With the integration of international finance, every stock market is not an in- dependent entity. Due to the increasing economic ties around the world, the stock market in Hong Kong and Taiwan have a longer history than the mainland. At present, there are relatively few literatures on the risk measurement of these four stock markets in China. This paper selects the daily closing price data of Shanghai and Shenzhen 300 Index of Shanghai stock exchange, the Shenzhen 300 Index of Shenzhen stock exchange, the Hang Seng Index of Hong Kong stock exchange market and the Taiwan Weighted Index of Taiwan stock exchange market, and combines the GARCH model with the POT model of extreme value theory to measure risk. 2. Empirical Methodology 2.1. VaR Value at Risk (VaR) is the maximum possible loss of a particular financial asset for a given period of time m at a certain level of confidence 1 p− . Define it as follow: VaR − = p 1 inf { ( ) x F x | m 1 ≥ − } p (1) where the inf indicates the minimum value that satisfies the real number x in the condition. As you can see from the definition, ≤ ) F VaR m p 1 − L m VaR 1 ≥ − t 1 − 1 ≥ − , that is: p (2) Pr p ( ) (   p   ( ) tL m represents the loss random variable of the position of the financial where asset. Therefore, from time t to time t m+ , the probability that the potential loss of the gold financial position holder is less than or equal to is 1 p− . VaR − 1 p 2.2. GARCH Model Common one-dimensional volatility models include ARCH models, GARCH models, and EGARCH models. In general, the return on financial assets has the characteristics of “volatility aggregation”, that is, the volatility is high in a certain period of time, and the volatility is small in other time periods. (p, q) in the GARCH (p, q) model represents the ARCH coefficient p and the GARCH term 965 Open Journal of Business and Management DOI: 10.4236/ojbm.2019.72065
S. X. Du et al. t a t tr , Define r µ= − t ta GARCH(p,q) model if coefficient q. The model consists of two equations: one is the conditional mean equation and the other is the conditional variance equation. For logarithmic rate of return as the disturbance or new interest at time t. We ta satisfies the following formula: call r =  t  = a t   2 a σ α α =  t i − µ + t σε t 2 βσ t − (3) ∑ ∑ a t + + 1 = 1 = p q 0 t t i j j i j where { }tε is independent and identically distributed random variable sequence with mean 0 and variance of 1. jβ ≥ , . For 0 p> i iα = ; for j jβ = . The last constraint for 0 ta is limited. At the same iα β+ tσ is time-varying. This article uses a low-order time, its conditional variance GARCH (1, 1) model: , it must be guarantees that the unconditional variance of iα > , 0 0 q> , it must be ( α β i 0α > , ( ∑ max + < 0 p q , 1 ) 1 = 2 ) i i i a t r µ = +  t t  = a σε t t  2 2 a σ α α βσ− = +  t t 1 1 1 − + 0 1 t t (4) 2.3. POT Model The full name of the POT model is the peak beyond the threshold, independent of the choice of the length of the subinterval, but requires a threshold. For a fixed shape parameter ξ, the mean excess function is a linear function of . De- fine the empirical mean excess function as: u u− 0 ( ) e u T = N 1 u −∑ N = iu 1 x t i ( ) u (5) uN is the number of exceeding the threshold u, named the excess; itx is where the value corresponding to the rate of return. The criterion for selecting a thre- shold by exceeding the excess function is that selecting a threshold that is large enough, and the mean excess plot after the threshold is approximately linear. Suppose the log-return is 1 x i pressed as: = − is the number of excess. The distribution function of the excess is ex- r , the distribution function is ( )F r , r r 2, , n u r i , ( ) F x u = Pr ( r u r | − > u ) = ( F x u + 1 − ( ) ) F u − ( ) F u , x ≥ 0 (6) Transform Equation (6) to get: ( ) F r = ( F x u + ) = − 1   ( ) ( ) F u F x u   + ( ) F u r , > u (7) For a sufficiently large threshold u, uF x can be approximately expressed as Generalized Pareto Distribution (GPD). GPD distribution includes shape para- G x ξβ can be de- meters ξ and scale parameters β. GPD distribution ( ) ; , ( ) DOI: 10.4236/ojbm.2019.72065 966 Open Journal of Business and Management
scribed as: ( ) G x ; ξβ , (  − + 1 1 =  1 exp −  1 − ) x ξ β ( ) − β x S. X. Du et al. ξ ≠ ξ = ξ 0 0 (8) For the distribution function pirical estimate ( )u actual sample as follows: n N − n ( )F r , it can be approximately expressed as em- . Then you can get the distribution function of the ( ) F r 1 = − uN n +  1   ξ ( β r u − − 1 ξ )    , r > u (9) 3. Empirical Study 3.1. Data Description This paper selects the daily closing price data of the Shanghai and Shenzhen 300 Index of the Shanghai stock exchange, the Shenzhen 300 Index of Shenzhen stock exchange, the Hang Seng Index of the Hong Kong stock exchange market and Taiwan Weighted Index of the Taiwan stock market. The time period is from March 7th, 2014 to March 7, 2019. The trading hours of stock markets in mainland China, Hong Kong and Taiwan are not completely consistent. In order to main- tain the consistency of data, the data processing method of Pei Yanhua et al. [11] (2017) is adopted to eliminate data with inconsistent transaction dates to get 1158 data. The logarithm and a difference operation of the daily closing price sequence is performed to obtain a daily logarithmic rate of return sequence, for a total of 1157 data. The article data comes from the CSMAR database, and the analysis of this article is implemented by R software. It can be seen from Table 1 that the daily logarithmic yield series of the Shanghai and Shenzhen 300 Index, the Shenzhen 300 Index, the Hang Seng Index and the Taiwan Weighted Index are all left-biased as the skewness is less than 0, and the kurtosis is significantly greater than 3, both with sharp peaks and thick tails feature. The J-B statistic rejects the null hypothesis of the normal distribution. It can also be seen from Figure 1 that the four indices are not obey to normal distribution as the scatter is away from the line. Both the ADF test and the P-P test passed, indicating that the yields of the four indices were stable. The volatility clustering feature refers to the fact that the volatility is high in a specific period of time while it is low in other periods. It can be seen from the daily logarithmic rate of return timing chart in Figure 2 that all four index yields exhibit volatility ag- gregation characteristics. The ARCH effect test is performed on the four yield series, and the results show that the rate of return series has conditional hete- roscedasticity. 3.2. GARCH Model Establishment According to the above analysis, the Shanghai and Shenzhen 300 Index, the Shenzhen 300 Index, the Hang Seng Index and the Taiwan Weighted Index are 967 Open Journal of Business and Management DOI: 10.4236/ojbm.2019.72065
S. X. Du et al. Table 1. Descriptive statistics of log-returns. Sample size Minimum Maximum Mean Stdev Skewness Kurtosis J-B ADF P-P ARCH-LM HS300 1157 −0.091544 0.077524 0.000489 0.016196 −0.811693 6.594329 2234.2971 −10.095 −33.121 197.25 SZ300 1157 0.000310 0.098733 0.000310 0.017878 −0.785064 5.326195 1494.1584 −10.376 −32.322 241.38 HSI 1157 TWII 1157 −0.060183 −0.065206 0.069870 0.000212 0.011238 −0.248895 3.758039 697.0613 −10.545 −33.393 36.617 0.035175 0.000145 0.008668 −0.919614 6.500682 2211.0504 −10.756 −33.899 89.327 Figure 1. QQ plot of daily log-returns. DOI: 10.4236/ojbm.2019.72065 Figure 2. Time series of daily log-returns. 968 Open Journal of Business and Management
S. X. Du et al. all stationary sequences, and all have an ARCH effect. The GARCH model not only reveals the “fluctuating aggregation” of financial markets, but also reveals the “thick tail” feature. According to previous scholars' research, the low-order GARCH (1, 1) model can well describe the characteristics of financial time series. Therefore, this paper establishes a low-order GARCH (1, 1) model for the Shanghai and Shenzhen 300 Index, the Shenzhen 300 Index, the Hang Seng Index and the Taiwan Weighted Index. Table 2 is the parameter estimation results of the GARCH (1, 1) model, except for the constant term, all the parameters pass the significance test. Among them, the ARCH parameter αand the GARCH parameter β are both greater than zero, which guarantees the positive definiteness of the conditional variance, and also indicates that fluctuations of the stock market are characterized by agglo- meration. The historical fluctuations are positively correlated with the current fluctuations and the speed is gradually slowed down. The large fluctuations often follow with large fluctuations, investors are more speculative in the stock market. α β+ A white noise test is performed on the standardized residual of the model, and the test statistic is the LB statistic. The results are shown in Table 3. Under the 6th and 12th order lags, the P value of the LB statistic is significantly greater than 0.05. It can be considered that the residual of the model belongs to the white noise sequence, that is, the model can well characterize the log yield series. 1 < ensure that the model is a smooth GARCH model. 3.3. POT Model Establishment An estimate of the standard deviation can be obtained from the above GARCH (1, 1) model. Since the GARCH (1, 1) model measures the risks in the normal market, in reality, there are often a few extreme values that can cause huge losses. Therefore, here we establish a GARCH (1, 1)-POT model for the normalized re- sidual of the GARCH (1, 1) model, that is, model the super-threshold data ex- ceeding the threshold u. The most important thing to establish a POT model is the choice of threshold. If the threshold is too large, the sample data exceeding Table 2. Parameter estimation of GARCH (1, 1) model. µ ω α β HS300 0.000683 0.000001 0.073318 0.925682 SZ300 0.000469 0.000002 0.051851 0.943504 HSI 0.000570 0.000002 0.049425 0.935894 TWII 0.000359 0.000006 0.106299 0.816171 Table 3. White noise test of residual. HS300 SZ300 HSI TWII P value Lag = 6 0.6919 Lag = 6 0.7399 Lag = 6 0.9241 Lag = 6 0.7382 Lag = 12 0.3987 Lag = 12 0.9446 Lag = 12 0.8979 Lag = 12 0.7359 969 Open Journal of Business and Management DOI: 10.4236/ojbm.2019.72065
S. X. Du et al. DOI: 10.4236/ojbm.2019.72065 the threshold will be too small, which may increase the variance of the parameter estimation. If the threshold is too small, the excess cannot be guaranteed to obey the generalized pareto distribution, which may result in biased parameter esti- mates. The method of selecting the threshold generally adopts the method of transcending the expectation function, and the selected threshold which value exceeds the threshold part to present a linear characteristic. In addition, there is the 10% principle proposed by Du Mouchel [12]: when the threshold u is allowed, about 10% of the data is selected as the extreme value data to be studied. According to the graph of the transcendental expectation function of Figure 3, based on the 10% principle proposed by Du Mouchel, it can be seen that the mean excess plot shows a linear trend after 1.1 or 1.2. Therefore, the thresholds are se- lected to be 1.2, 1.2, 1.2 and 1.1 respectively. The generalized Pareto distribution (GPD) is fitted to the data exceeding the threshold, and the estimation results of the shape parameters and the scale parameters are shown in Table 4. In order to test the fitting effect of the model, the diagnostic test charts of the GARCH (1, 1)-POT models of the standardized residuals of the four index yield series are shown in Figures 4-7 respectively. The upper left corner is the distri- bution function graph of the overrun, the upper right corner is the distribution tail probability estimation graph, the lower left corner is the residual scatter plot, and the lower right corner is the residual QQ plot. The closer the scatter and the line are, the better the model fits. It can be seen from the four graphs that most of the scatter points are on or near the line, and very few points have slight deviations Figure 3. Mean excess plot. Table 4. Model parameter estimation of GARCH (1,1)-POT model. Threshold Excess ξ β HS300 1.2 117 SZ300 1.2 100 HSI 1.2 110 TWII 1.1 108 −0.009138573 −0.008040062 0.583801381 0.498632659 0.05280263 0.48535986 −0.04282757 0.53193409 970 Open Journal of Business and Management
分享到:
收藏