GARCH模型案例文档良心出品

1、GARCH模型案例1 .数据选取与时段选择本案例以上证指数为例,通过ARCH/GARCH 模型研究我国证券市场的波动性规律.尽管上证指数从1990年12月9日开始公布,但由于在开始阶段,进入流通的样本股票数量少,而且交易制度不完善,股票 投机性强,所以股市异常波动性太大.1996、1997年以后,这种异常波动趋于平稳,上证指数方差变化指小于0.03.考虑到我国股市制度变化对收益变化有很大影响,因此在时段选择上还要考虑股市交易制度的 变化.为了保证股市稳定,预防过度投机行为,中国股票市场交易1996年12月6日开始实行T+1交易制度,以及实施涨跌停板限制.综合以上因素,把数据分析时段选择为 19

2、98年1月1日至2007年9月28日,共2350个数据.2 .波动率及其特征金融资产收益率的波动性在证券、期权交易中是一个重要因素,它是标的资产的条件方差.波动率在 风险治理中也是重要的,它为计算资产的在险价值 VaR提供了一个简单的方法.一般来说,波动率不能 被直接观测到,但它也具有一些特征值得研究.这些特征包括：1波动率存在聚类性,也就是波动率可能在一些时间段上较高,而在另一些时间段上较低；2波动率以连续方式变化,波动率的跳跃现象是少见的；3波动率是平稳的,不会发散到无穷,而是在一定范围内随时间连续变化；4波动率对利好消息和利空消息的反响是不同的,即存在杠杆效应.3 .数据根本分析本案例研

3、究的收益率形式为日对数收益率,即 k=她.酶斗对收益率数据进行初步分析得当的结果如下表所示:均值最大值最小值标准差偏度峰度J-B检验值0.0006540.09404-0.0925540.01490.0455338.0032451.654从表中数据可以看出,股指日对数收益率的均值很小,可以认为是 0.收益率的分布具有正的偏度,所 以分布的尾部略向右拖,说明盈利的概率要大于亏损的概率.峰度值大于正态分布的峰度正态分布的峰度为3),这反映了收益率分布具有尖峰厚尾的特征.下面再进行上证指数时序特征分析.观察上证指数时序图,收益率确实存在明显的聚类效益(即一次大的波动后往往伴随着另一次大的波动).对其进

4、行不包含截距项和趋势项的 ADF单位根检验,得到的结果如表所示,说明日对数收益率序列在各个显著性水平下都是平稳的,图()是序列的自相关图和偏自相关图,Q统计量显示序列在 3阶以后根本不存在自相关现象,因此可以用AR(3)模型拟合收益率序列,使用软件估计模型得到：ft = 0.058452rt_a + ut(2.835)R: = o.om D.W.=1.96ADF 统#<1%临界值5%临界值10%临界值-47.48441-2.565947-1.940959-1.616609Autocorrelation Partial Correlation AC PAC Q-Stat Prob|111

5、0 018 0.018 0.7296 0.393(112 -0.035 QD35 3.5316 0.1711103 0,057 0,058 11.163 0.01111114 0.040 0.037 14.993 0.00511115 0.000 0.D03 14.993 0.0100116 -0.053 -0.054 21 622 0.001111j7 0.029 0.027 23.692 0.00111118 0.006 -0.001 23.670 0.D0311119 0.004 0.012 23.703 0.005111110 0.000 0.001 23.704 0.00811111

6、 0.027 0.D26 25472 0.D0811J12 0.056 0.054 33.551 0.001在此根底上进行条件异方差检验,主要使用Engle提出的ARCH-LM检验方法.ARCH-LM检验的原假设是残差序列中直到 p阶都不存在 ARCH效应,能够得到两个检验统计量：(1) F统计量是对所有残差平方的滞后的联合显著性所作的一个省略变量检验；(2)统计量是Engle's LM检验统计量,一般情况下渐近服从分布.在本案例中,检验结果如卜表所不,说明日对数收益率序列确实存在条件异方差现象.ARCH(q)LM统计量概率结论128.971340.000000存在异方差249.648

7、320.000000存在异方差3100.50290.000000存在异方差4195.84110.000000存在异方差4.建立GARCH模型通过对对数收益率数据的简单分析,收益率序列存在明显的异方差性,应该考虑建立ARCH/GARCH模型.由于GARCH模型具有比ARCH更优秀的性质,例如它比 ARCH模型需要更小的滞后阶数,并有与 ARMA模型相似的结构,因此本案例直接建立GARCH模型.建立GARCH模型通常包括三个步骤：(1)建立均值方程(例如 ARMA ),以别离出数据中任何线性相关成分；确定GARCH (p, q)模型的结构.与ARCH模型不同,GARCH模型的阶数较难确定,在实际应

8、用时,一般只会用到低阶的GARCH模型,如GARCH(1,1)、GARCH(1,2)、GARCH(2,1)等.在针对金融时间序列分析的研究中,有两类模型经常用到,包才GARCH(1,1)模型和GARCH(1,1)t模型.其中后者是指随机扰动项,服从t分布的情况,能更好地刻画一些金融时间序列的波动性和顶峰厚尾现象；(3)检验拟合的 GARCH模型,对其进行必要的修正.本案例的建模过程也遵循这一步骤,直接使用上面的 AR(3)作为均值方程,在此根底上建立GARCH(1,1) 模型,模型形式为：%+与IL情二% + %+ 网熄T阳ffQ,肉,由于待估参数,使用软件估计得到：( ft = 0.0425

9、43_a +stj(2.085)R2= 0.001345 D.W.= 1.96(睹=瓦虱 X 10-6 + 0.12364eL + 0.85544g,(6.55)(12.897)(82.64)需要特别指出的是,均值方程伴有GARCH方程以后,方程中的某些项通常常失去显著性.为了进行比拟,表给出了分别在正态分布、t分布和广义误差分布假设下的估计结果：分布类型均值方程GARCH方程正态分布1=B.042543Jf_j + 层=6,41X_10-6 + 0,12364e_i+ 0r85544ufi1(2.085)(6.55)(12.897)(82.64)t分布n = 10728rE_a - %醴=7

10、.21 x IO'6 + 0.10783f= + lfi654U=(3.576)(3.82)(6.266)(44.84)广义误差分布.；-；-.二,二 ：一 ：二二-(3.530)(3.756)(12.897)(82.64)不难看出,在随机干扰项 0服从t分布或者广义误差分布的假设下,均值方程的参数显著性都比上服从正态分布假设条件下要高,进一步验证了金融时间序列具有顶峰厚尾的特点.5 .模型结果分析在本案例针对上证指数日对数收益率拟合的一组GARCH(1,1)模型都具有以下几个共同特征：(1) 模型中的 良系数都较大,并且通过了显著性检验,说明指数波动具有“长期记忆性,即过去价格的波动

11、与其无限长期价格波动的大小都有关系.(2) GARCH方程中见卡氏接近于1,说明条件方差函数具有单位根和单整性,也就是说条件方差波动具有持续记忆性,说明证券市场对外部冲击的反响以一个相对较慢的速度递减,股市一旦出现大的波动 在短时期内很难消除.(3) GARCH方程中 附+用1, 说明收益率条件方差序列是平稳的,模型具有可预测性.6 . GARCH模型改良在进行金融时间序列分析时还经常用到GARCH-M模型和EGARCH模型.前者将条件方差项引入均值方程,说明了风险与预期收益存在紧密联系；后者那么为研究条件方差的非对称性提供了便利的分析工具.此外,EGARCH模型放宽了传统 GARCH模型中对

12、参数的非负限制,而且由于EGARCH模型的条件方差方程是4的一次函数,可以比拟好地判断波动源的持续性.首先建立GARCH(1,1)-M 模型,模型形式为：(1)正态分布假设下：1 = 0 刖1865%+0,04328械 + 球(2.054)(2.312)户二 0刖5641 D.M=L967W = U7x MT6+0122021 + 0.85822%(6.(460) (12.382)(81.854)(2)在t分布假设下：ft = 007105%t + 0,M4 优 + 品(2.054) (2.233) = 0.00549 5,=1967砰二710xl0-6 + 0.10945七 + 0.8649

13、36!(3.(781) (6.27)(44.82)(3)在广义误差分布假设下：ft = 0,0651 + 0.0488+(3.406)(2.632) R: = 0.0059 DM= 1968砂=655 x 1O7 + o,iu2% + 0.8631(3.708)(6.583)(45.33)在各种分布假设下,均值方程中的.都具有显著性,说明了证券市场的预期收益与风险成比例,即风险 越大,期望收益越大.下面再建立EGARCH(1,1)模型,模型形式为：h = %+免使用软件估计得到:A 0.03344镐 + %(1.987)R2 = 0,0012 DM= L96X = -042634-0.0378 |邑+ 0.2345 IgT+ 0S7