《恒生指数收益率建模》由会员分享,可在线阅读,更多相关《恒生指数收益率建模(21页珍藏版)》请在金锄头文库上搜索。
1、基于GARCH模型的金融资产收益波动率的研究以恒生指数为例【摘要】在金融领域中,金融资产收益率的条件方差通常作为对资产风险的一种度量,它是金融资产定价的数学理论以及在险价值计算的重要组成部分。本文以恒生指数的月收益率作为样本数据,利用GARCH模型对金融资产收益的波动率建立模型,并对2013年1月到11月的波动率进行预测,与样本外数据作对比,从而发现用GARCH模型拟合金融资产收益波动率的恰当性。【关键词】波动率 聚集性 GARCH模型 尖峰厚尾1研究目的与意义波动率研究是资产定价方面的重要研究方向。波动率估计准确与否直接关系到模型运用是否得当,投资策略是否成立。关于波动率的研究相对较长的历史
2、,在1982年,Engle提出了ARCH模型的思路,为波动率研究开启了新的篇章。在1986年,Bollerslev在ARCH的基础上将历史波动率加入模型,提出了GARCH模型,更好地抓住了股票收益波动率的长记忆性。之后,关于波动率的量化分析都以GARCH模型为基础。本文利用GARCH模型对恒生指数(HSI)1989年12月1日到2013年11月30日月收盘价观察值进行实证分析,比较GARCH各类模型对恒生指数对数收益波动率的分析结果。本文主要分为四个章节。第一章介绍本文的研究背景、意义、方法及目的,介绍本文研究思路和框架。第二章介绍了GARCH族模型的起源与发展,介绍了ARCH模型、GARCH
3、模型提出的背景以及它们的表示方式,和模型识别、诊断。第三章对恒生指数月收盘价样本做ARCH效应分析和检验,包括自相关性检验、平稳性检验、异方差性检验,得出上证指数收益率可以使用ARCH模型的结论,使用GARCH模型进行建模。第四章对全文结论作出总结,并提出了本文的不足及改进建议。结论发现,尽管GARCH中均有不少模型能对样本数据进行分析,但分析结果梅有较大差异,对所使用的数据,GARCH模型能较好地捕捉到上证股指收益率整体样本序列的尖峰厚尾特征和波动的聚集性。对于监管者、从业者、研究者、投资者来说,波动率的实证研究至关重要,希望本文能对各市场相关者提供有意义的参考,为共同建设一个功能更加健全、
4、更加灵活、更加富有弹性的金融市场做出贡献。2理论依据2.1 ARCH模型对金融资产收益序列,可以同时考虑其前二阶矩的时变特征。对于给定的信息集,它的条件均值和条件方差分别为: (2.1.1) (2.1.2)其中用于描述金融资产期望收益的时变性,揭示金融资产收益能力;用于描述金融资产方差的时变性,揭示金融资产的风险特征。Engle(1982)提出的ARCH模型可表示为: (2.1.3)其中为解释变量组成的向量,可以由的滞后项或者其他外生变量组成;是t期的扰动项,它为独立同分布的白噪声过程,表示偶发因素的作用;,保证条件方差严格为正,且保证ARCH过程平稳。由于,往往呈现出波动聚集性,即“大幅波动
5、往往集中在某些时段上,而小幅波动集中在另外一些时段上”。因此ARCH模型可以很好地揭示金融市场的“波动聚集”这一特征。经典的ARCH效应检验方法为拉格朗日乘子检验(Lagrange multiplier,LM)。对ARCH模型中的均值方程,若随机变量为独立白噪声过程,且,这时在式(6.3)的方差方程中有,而为一个常数。如果随机变量服从ARCH过程,那么中至少有一个不为零。因此ARCH模型的原假设和备择假设为:不存在ARCH效应 不全为零取,利用LM检验的ARCH效应,其检验统计量为 (2.1.4)在成立时,统计量LM的极限分布为。利用对数似然函数对所讨论的参数向量求一阶和二阶偏微分,在样本数充
6、分大时给出的LM算式为: (2.1.5)式中,T为样本数,取残差,是对()进行回归所得到的拟合优度,。ARCH模型最常用的估计方法是极大似然估计,其对数似然函数为 (2.1.6)式中,利用优化方法即可得到参数向量的一致估计量。2.2 GARCH模型一般的,GARCH(p,q)模型可以表为 (2.2.1)GARCH模型不仅能揭示金融市场的“波动聚集性”特征,而且可以揭示“厚尾”特征。例如一个GARCH(1,1)模型的波动方程: (2.2.2)式中保证条件方差的正定性;保证该模型为平稳的GARCH模型。可以证明,若,容易得到模型表示的无条件峰度满足: (2.2.3)从而可以刻画那些比正态分布更厚尾
7、部的金融时间序列。一个GARCH(p,q)模型可以通过两步法来确定其滞后阶数。对于残差序列,由于是的无偏估计,可以使用的偏自相关函数来确定ARCH部分的最优滞后阶数。在给定的前提下,可以使用AIC或者BIC准则确定GARCH部分的最优滞后阶数。GARCH模型最常用的估计方法是极大似然估计,定义,则,在随机变量服从正态分布下,其条件密度函数为: (2.2.4)对于T个观测值下的对数似然函数为: (2.2.5)参数向量的最大似然估计为方程组的解。若GARCH模型建立恰当,它的标准化残差序列应该服从独立同分布,且具有零均值和单位方差。可以使用Ljung-Box统计量分别检验和来估计GARCH模型中均
8、值方程与方差方程设定的正确性;通过的Q-Q图识别分布假设的正确性。在GARCH模型预测中,最核心的为条件方差(波动率)预测。对于式(2.2.1)定义的GARCH(p,q)模型,为得到其向前步预测,可以将方差方程向前递推步:因此,对方程两边同时使用条件期望,可得:这样可以通过递归方式对条件方差进行求解,从而得到波动率的向前步预测结果。3实证分析实证分析部分以恒生指数1989年12月1日到2012年12月31日数据作为内部数据建立模型,以2013年1月1日到2013年11月30日数据作为实验外部数据。恒生指数是香港股市价格的重要指标,指数由若干只成份股(即蓝筹股)市值计算出来的,代表了香港交易所所
9、有上市公司的12个月平均市值涵盖率的63%。图表1展示了HSI数据的时间序列图。数据呈现出总体上升的趋势,并且具有随着股价增加波动幅度更大的迹象。令表示时间序列,第t天的收益率定义为:分别计算样本的月收益和日收益,并画出收益的时间序列图。如图表2和图表3,该图显示出收益率在一些时期波动更大,特别是在1997年到1998年间,这一现象是由亚洲金融危机引发的。图表 1HSI价格的时间序列图图表 2 HSI月收益率图表 3HSI日收益率利用Ljung-Box对其月收益率自相关性进行检验,结显示月收益不存在序列自相关。X-squareddfp-value15.348120.2229随后用同样方法检验平
10、方月收益,存在序列自相关。X-squareddfp-value35.311120.利用LM检验月收益序列的ARCH效应,p值小于0.05,存在ARCH效应。Chi-squareddfp-value24.685120.01639计算均值调整对数收益,如图表4,除了在根据平方的自相关和偏自相关函数的拖尾现象,可以使用ARMA(1,1)模型作为均值方程。一般地,使用低阶的GARCH模型就能很好的捕捉到金融序列的波动聚集性,所以这里我们使用GARCH(1,1)作为条件异方差的方程。图表 4 为了比较不同模型拟合的效果,我建立了六种GARCH模型,分别为:GARCH(1,1)-N模型,GARCH(1,2
11、)-N模型,GARCH(1,1)-t模型,GARCH(1,1)-st模型,GARCH(1,1)-GED模型,GARCH(1,1)-SGED模型。将运算结果整理到如下表格1中,具体运算结果见附录。上述所有的模型都能消除异方差;对残差相关性检验,发现残差都可以消除自相关特性,说明残差项不具有自相关的特征。根据AIC,BIC准则,选择AIC和BIC较小的模型,在GARCH(1,1)-N模型和GARCH(1,1)-GED模型中选择。由于只有GARCH(1,1)-N模型的参数估计均通过了显著性检验,而其他模型均有部分参数估计没有通过显著性检验,所以选择GARCH(1,1)-N模型拟合数据。表格 1 估计
12、结果表GARCH(1,1)-NGARCH(1,2)-NGARCH(1,1)-tGARCH(1,1)-stGARCH(1,1)-GEDGARCH(1,1)-SGEDmu1.109591.175871.253551.128051.276281.15146omega3.442854.8452.961922.952593.094713.09745alpha10.232730.124850.19020.184410.207720.19861alpha20.18343beta10.724620.630010.769880.773330.750890.75667Log Likelihood-959.9309
13、-958.4313-955.989-955.5339-955.7724-955.2982AIC6.6.6.6.6.6.BIC6.6.6.6.777096.6.J-B Test13.198857.480818.2936418.823715.7490216.52389Ljung-Box Test9.8.9.035269.9.9.LM Arch Test13.856889.14.1496214.2162913.9493514.07678提取GARCH(1,1)-N模型得到的波动率估计和标准化残差,并画出下图:Q-Q显示出直线的模式Jarque-Bera检验的p值为0.00136,Shapiro-Wi
14、lk检验的p值为0.0276,因此正态假设不能被拒绝。如果GARCH的模型被正确识别,那么标准残差近似独立同分布,信息的独立同分布可以通过他们样本ACF来检验。如上图,分别为波动率、波动率平方、标准残差和标准残差平方样本的ACF图。在数据独立同分布假设下,图形中临界限是基于倍的名义方差。即使模型被正确识别,该名义值与残差平方的自相关函数的真实方差也可能存在很大区别。不过从图形大体可以看出波动率、波动率平方、标准残差和标准残差平方序列不相关。GARCH(1,1)模型可以较好地拟合HSI数据,用其预测未来11个月的条件方差。这里用到了样本外数据与预测数据进行对比,结果如表2和表3所示:表格 2 恒生指数月对数收益的波动率预测预测步长1234567891011GARCH(1,1)-N4.589 4.858 5.103 5.326 5.532 5.722 5.898 6.062 6.215 6.357 6.491 GARCH(1,2)-N4.578 4.604 4.969 5.231 5.480 5.704 5.909 6.097 6.269 6.42