经济理论与应用

基于GARCH模型的我国创业板市场风险度量*

宋永辉, 陈 晨

(沈阳工业大学 经济学院, 沈阳 110870)

以创业板综合指数的日对数收益率为基本数据,分别在正态分布、T分布与GED分布下总结GARCH参数并计算VAR值,通过实证分析指出序列分布具有尖峰厚尾的特征,利空消息带给创业板市场的波动比等量的利好消息更大,创业板市场存在明显的杠杆效应,波动具有聚集性和持久性;GED分布能够更好地拟合创业板市场综合指数,得出的VAR值精度较高;模型种类并不是VAR值结果的决定性因素。

创业板; 风险度量; GARCH模型; 利空消息; 利好消息; 杠杆效应

对于一个国家的资本市场来说,创业板市场无疑是不可缺少的重要部分,其主要意义在于为成长中的中小企业提供融资渠道和退出机制。国际上成熟的创业板市场已经比较完善,具有丰富的操作经验,而我国的创业板市场还存在一些缺陷与不足。无论对于创业板本身或暴露于创业板的主体而言,对风险的认识都是必不可少的。

一、文献综述

刘建民、张颖(2015)认为,创业板的顺利发展需要借鉴国际上成功创业板的发展经验,加强监管,修正问题,完善制度体系,提高资本市场中各利益相关者的素质[1]。王璇(2015)在总结创业板风险来源的基础上,对创业板的六家公司进行了投资组合的风险计算,并提出在创业板市场风险的控制和防范方面,应挑选高质量的上市公司,增加信息的透明度,建立全面的监管体系,完善退市制度[2]。林博(2013)对GARCH模型和VAR方法的发展历程进行了详细的介绍,在此基础上对上证指数的风险进行了度量分析[3]。薛炜星(2016)探讨了测量VAR的三种基本方法:历史模拟法(historical simulation)、方差协方差法和蒙特卡罗模拟法(Monte Carlo simulation)[4]。沈一凡(2015)介绍了风险度量方法的发展历史,对比了各种方法之间的差异[5]。黄崇珍、曹奇(2017)以华夏沪深300ETF联接为例,对开放式基金的风险度量进行了研究[6]

Neil(2016)的研究采用极值理论的方法估计GARCH模型的残差,度量了存在异方差特性的金融时间序列的尾部风险[7]。Ramazan(2015)对极值理论、GARCH族、历史模拟法等方法进行了比较分析[8]。Samit和Prateek(2017)利用条件极值理论(EVT)构建GARCH模型,并且进行了动态的预测能力模型与独立的GARCH模型的比较[9]

二、我国创业板市场现状

我国的创业板市场成立于2009年10月,自成立以来,在不断暴露出问题的同时也不断地发展与完善。截至2017年7月12日,创业板市场从2009年的28家上市公司扩展到656家上市公司。2009—2017年中国创业板市场基本数据如表1所示。

由表1可以看出,2009年12月—2017年7月,我国创业板的上市公司数量、总发行股本(百万股)、总流通股本(百万股)都在不断地增长。创业板上市公司市价总值已经从2009年创业板刚成立时的1 610亿元连续增加到约49 780亿元,可见其规模不断扩大。但是在规模不断扩大的同时,其各个指标的增速并不一致。2009—2017年总发行股本和总流通股本一直保持着匀速稳定增长的状态,而市价总值和流通市值则在2009—2015年保持增长,但在2015—2017年有下降的趋势。

表120092017年中国创业板市场基本数据

日期上市公司数总发行股本百万股总流通股本百万股上市公司市价总值百万元 上市公司流通市值百万元平均市盈率%2009/12/3136 3460.31 647.79161008.35 29897.03105.382010/12/3115317506.295037.76736521.89200564.1078.532011/12/3128139953.3414222.03743379.22250408.1537.622012/12/3135560089.4724205.22873120.46333528.6332.012013/12/3135576156.3143000.761509197.85821882.7855.212014/12/31406107725.9268768.672185094.931307290.0964.512015/12/31492184045.22116888.785591624.783207867.73109.012016/12/30570263060.93170044.115225450.433053690.2973.212017/07/12656310372.05203926.284977991.912978259.5849.88

创业板的增长态势值得鼓励,但在市盈率走高的背后却是高估值与高业绩增长水平的背离。在创业板市场,投资者们在高收益的引导下脱离基本的理论知识与分析,意识不到高风险的存在,盲目地抱着投机的心态一心想牟取暴利。因此,投资者对创业板市场进行冷静理性的分析、识别投资的利弊、保持冷静与谨慎是非常有必要的[10]

三、GARCH-VAR模型

1. GARCH模型

在金融领域,扰动项μt的条件方差会受到来自之前时刻变量的影响。考虑到的分布滞后模型,其中α0α1,…,αp为方差方程中的p+1个系数,为时刻t-1的残差平方,用一个或两个的滞后值代替许多的滞后值,这就是广义自回归条件异方差模型(GARCH模型)[11]。在GARCH模型中,涉及到两个不同的方程——均值方程和方差方程,分别如式(1)、(2)所示。

一般的GARCH模型可以表示为

(1)

(2)

则称序列服从Garch(p,q)过程。式(1)中:xt为1×(k+1)维外生变量;δ为1×(k+1)维系数向量。式(2)中:k为常数项;可以度量前期的波动信息即ARCH项;为上一期的预测方差即GARCH项;ε和λ为系数。

与波动的聚集效应相同:伴随着较大的波动后可能是更大的波动,较小的波动后伴随着更小的波动。Garch(p,q)模型是ARCH模型的扩展,因此都具备这一特点[12]GARCH模型在描述高阶且计算量不大时优于ARCH模型,因而具有比ARCH更强的适用性。

2.VAR模型

1993年7月G30集团在一个报告中提到用“风险价值系统”来评估金融风险,这份报告是在研究金融衍生工具的基础上提出的。之后Morgan推算出了计算VAR的风险模型用以衡量风险[13]。2004年发布的新巴塞尔协议中进一步主张用VAR模型对商业银行面临的风险进行综合管理[14]。目前,国际金融机构已经广泛采用VAR作为衡量金融法风险的指标。

VAR的基本定义:在正常的市场条件下和给定的持有期内,某一金融资产或证券组合在未来特定一段持有期内的最大可能损失[15]。统计学公式表示为

Prob(ΔPΔt≤VAR)=θ

(3)

式中:ΔP为金融资产在持有期Δt的价值损失额;VAR为给定置信水平θ下的在险价值,即可能的损失上限;θ为给定的置信水平。在确定持有期和置信水平的情况下可以计算出资产或资产组合最大可能的预期损失,也就是VAR值。

四、实证分析

1. 数据的选取

对创业板市场风险进行度量,采用Eviews6.0软件对创业板市场综合指数分析处理。选取了2012年1月7日—2017年6月1日的创业板综指,共1 089个样本数据。在实证分析前,对创业板综指进行日对数收益率的处理,计算出1 088个日对数收益率,表示为其中Pt为创业板综指t日的收盘价;Pt-1为t-1日的收盘价。这种操作既可以消除数据的异方差性,又可以使得收益率具有弹性的概念[16]

2. 数据的基本分析

(1) 数据的统计特征。数据的基本统计特征如表2所示。由表2可知,创业板综合指数日对数收益率均值为0.001 078,偏度为-0.740 38,左偏,峰度为5.134 497。数据在服从正态分布时,通常偏度为0,峰度为3,因此可以看出创业板综指日对数收益率有显著的尖峰厚尾特征,不服从正态分布。

表2数据的基本统计特征

序列观测数量均值中值最大值最小值标准差偏度峰度JB统计量P值R10880.0010780.0021980.067616-0.090130.22016-0.740385.134497305.94230.0000

(2) 单位根检验。遵循AIC和SC最小原则,对创业板日对数收益率进行ADF单位根检验,结果如表3所示。由表3可知,ADF检验值为-30.210 53,P值为0.000 0,小于5%显著性的临界值,故拒绝序列不平稳的原假设,创业板综指的日对数收益率序列是平稳的序列。

表3ADF检验结果

ADF检验值1%临界值5%临界值10%临界值P值结论-30.21053-3.436149-2.863989-2.5681250.0000平稳

(3) 自相关检验。对创业板综合指数进行自相关检验,结果如表4所示。通过Eviews系统分析总结出创业板综合指数36阶滞后的AC值和PAC值。P值接近零,高阶后等于零,可见数据之间不存在明显的自相关性,也就是说创业板日对数收益率之间基本不存在自相关。

表4自相关检验结果

阶数AC值PAC值Q统计量P值1 0.087 0.0878.17450.00410-0.089-0.08828.64600.001200.0250.02343.38100.002360.0160.02572.49200.000

(4) 正态分布检验。Q-Q图是最常用的视觉检测是否服从正态分布的常规方法,在数据满足正态分布时,它的Q-Q散点图会是一条直线。根据自相关检验结果绘制数据的Q-Q散点图,得出Q-Q散点图两边有凸起,中间弯曲凹陷,因此明显看出创业板综指日对数收益率是不服从正态分布的。

(5) LM异方差检验。创业板综指日对数收益率的波动集中在(-0.08,0.08),但在这个区间内价格波动较大,时高时低,有时激烈波动,有时平缓,可以看出其更容易在较大的波动后出现更大的波动,较小的波动幅度后随之而来的是更小的波动幅度,因此创业板综指日对数收益率存在显著的集聚性效应。

对创业板综指的残差平方进行自相关检验,自相关与偏相关系数AC和PAC显著不为零,Q统计量显著,因此拒绝不存在异方差的原假设,我国创业板市场综指的日对数收益率有明显的异方差性。

对数据进行ARCH效应检验,结果如表5所示。由表4可知,P值为零,也就是χ2分布的伴随概率P为0,创业板综指存在ARCH效应,拒绝不存在ARCH效应的原假设。

表5ARCH效应检验结果

阶数qF统计量观测值个数与R2乘积P值147.7747345.84250.0000350.79126134.02850.0000537.18601159.41950.0000726.71753160.42960.0000922.14887169.55860.0000

3. GARCH-VAR模型的建立

对数据进行ADF单位根检验、正态分布检验、自相关检验和异方差性检验后,可以判断出创业板综合指数日对数收益率是平稳的,在不满足正态分布的同时不存在自相关性,存在异方差并具有ARCH效应。因此,设收益方程为一般的均值回归方程,即rt=μ+σt,其中,μ为均值,σt为均值方程的残差。由于数列存在明显的异方差性,遵循AIC最小原则,经过反复的测试得出滞后阶数(pq)为(1,1)时AIC值取到最小,所以本文模型均为GARCH(1,1)类模型。

(1) 基于正态分布的GARCH族模型参数及VAR检验值。每个模型对应两行内容,第一行为模型参数值,第二行为Z检验值(见表6)。由表6可知,模型参数在5%的显著性水平下都是显著的。因为常数项不影响效果的估计,故没有残留明显的异方差现象。这说明模型可以较好地拟合创业板综合指数日对数收益率的异方差现象。EGARCH模型的非对称项系数γ1=-0.004 779显著小于零,根据指数模型的对数性质,可以得出利空消息带给创业板市场的波动比等量的利好消息更大,创业板市场存在明显杠杆效应的结论。

表6基于正态分布的GARCH类模型参数值及Z

模型ω1α1β1γ1ρ1ψ1GARCH(1,1)TARCH(1,1)EGARCH(1,1)CARCH(1,1)2.24×10-60.0418240.952642 0.0001000.0000000.0000002.04×10-60.0508830.956113-0.022221 0.0006000.0000000.0000000.087900-0.1100970.0854830.994440-0.004779 0.0000000.0000000.0000000.596500 0.0004030.0168710.6314680.9946570.040881 0.0126000.5545000.3966000.0000000.000000

注:ωc值,α为RESID(-1)2值,β为GARCH(-1)值,γ为RESID(-1)2与RESID(-1)<0值的乘积,ρC(3)值,ψC(4)值,下同。

基于GARCH模型参数,计算得到VAR的统计特征值及Kupiec检验结果(见表7)。由表7可知,在95%和99%置信水平下,通过正态分布GARCH类模型得出的VAR均值差别不大,而最大值和最小值的差别却比较大。由表7可知,LR检验的结果:在显著水平为5%的情况下,TARCH(1,1)和CARCH(1,1)模型的LR检验值显著小于3.41,也就是TARCH(1,1)、CARCH(1,1)模型通过了LR检验,相对于没有通过检验的模型结果更加准确。在1%的显著水平下,GARCH(1,1)模型和EGARCH(1,1)模型LR检验的结果小于6.635,所以这两类模型都通过了检验,即GARCH(1,1)和EGARCH(1,1)模型下计算VAR检验值得到的结果较准确。

表7基于正态分布的GARCH类模型VAR返回检验值

模型置信度/%最小值最大值均值标准差失败天数失败率LR检验值GARCH(1,1)TARCH(1,1)EGARCH(1,1)CARCH(1,1)95-0.0574-0.0340-0.02790.0310870.07994.7262099-0.0743-0.0156-0.20860.0132250.02306.4836095-0.0617-0.0087-0.02760.01071030.09460.02783990.0819-0.0127-0.03460.0149470.04329.4680095-0.0691-0.0089-0.06930.0102740.068011.3945099-0.0832-0.0147-0.04370.0156170.01565.0767095-0.0607-0.0298-0.02780.01521040.09551.4734099-0.1480-0.0295-0.04810.0195430.039511.03940

(2) 基于T分布的GARCH族模型参数及VAR检验值。模型参数在5%的显著性水平下结果均较为显著,检验残差的LM效果估计和残差检验异方差性见表8、9,结果表明没有残留明显的异方差现象,因此模型可以较好地拟合创业板综合指数日对数收益率的异方差现象,并且从参数估值结果可以看出创业板指数存在明显的非对称效应。

表8基于T分布的GARCH类模型参数值及Z

模型ω2α2β2γ2ρ2ψ2DF值GARCH(1,1)TARCH(1,1)EGARCH(1,1)CARCH(1,1)1.44×10-6 0.058894 0.9397497.346026 0.2378000.0000000.0000000.0000001.40×10-60.0733780.943506-0.0322457.313138 0.2333000.0025000.0000000.1844000.000000-0.1195010.1173670.9963850.0001017.314141 0.0009000.0000000.0000000.9946000.000000 0.001986-0.055709-0.4519270.9993390.0588587.129285 0.9057000.0409000.3216000.0000000.0000000.000000

表9基于T分布的GARCH类模型VAR返回检验值

模型置信度/%最小值最大值均值标准差失败天数失败率LR检验值GARCH(1,1)TARCH(1,1)EGARCH(1,1)CARCH(1,1)95-0.0640-0.0131-0.02290.0107570.052321.476699-0.1900-0.0188-0.04670.0226210.01937.537695-0.0469-0.0106-0.04270.0124490.045023.629199-0.1299-0.0263-0.07200.0211150.01389.215695-0.0386-0.0153-0.03260.0136470.043227.076599-0.1310-0.0179-0.06140.0199190.01746.742695-0.0600-0.0235-0.02950.0240560.051424.748699-0.1201-0.0354-0.05100.0153100.009211.0947

由表9可知,在95%和99%置信水平下,各种计算出的VAR均值之间差距较最大值与最小值之间差距大。在5%的显著水平下GARCH(1,1)、EGARCH(1,1)、TARCH(1,1)、CARCH(1,1)模型的LR检验值均大于3.841;在1%的显著水平下,模型的LR检验值均大于6.635,所以,在显著性水平为1%和5%时都要拒绝原假设。可见,基于T分布的创业板综合指数计算出的VAR值过于保守,存在着风险被缩小的可能。综合而言,在T分布假设下,根据创业板综合指数日收益率各种GARCH类模型计算的参数值和Z检验结果,得出结论:创业板综合指数日收益率在服从T分布时是不合理的。

(3) 基于GED分布的GARCH族模型参数及VAR检验值。检验残差的LM效果估计和残差检验异方差性见表10、11,同样没有残留明显的异方差现象,因此模型可以较好地拟合创业板综合指数日对数收益率的异方差现象,创业板指数存在明显的非对称效应。GED分布下各模型估值结果的尾部参数即V值在1.45上下波动,因此,相较于其他分布,GED分布可以更好地刻画创业板指数日对数收益率数据的厚尾特征。

表10基于GED分布的GARCH类模型参数值及Z

模型ω3α3β3γ3ρ3ψ3VGARCH(1,1)TARCH(1,1)EGARCH(1,1)CARCH(1,1)1.98×10-6 0.050531 0.9453911.452988 0.0448000.0000000.0000000.0000001.91×10-60.0631200.948965-0.0284821.454149 0.0534000.0017000.0000000.1806000.000000-0.1191420.1029190.995025-0.0022421.457718 0.0001000.0000000.0000000.8685000.000000 0.000500-0.036820-0.4870200.9961710.0503351.446115 0.2970000.2371000.4288000.0000000.0000000.000000

注:V为GED PARAMETER值。

表11基于GED分布的GARCH类模型VAR返回检验值

模型置信度/%最小值最大值均值标准差失败天数失败率LR检验值GARCH(1,1)TARCH(1,1)EGARCH(1,1)CARCH(1,1)95-0.0701-0.0191-0.02760.0124680.06243.619999-0.1190-0.0106-0.03340.0163170.01560.571895-0.0480-0.0102-0.02480.0119840.07711.434299-0.0816-0.0132-0.03910.0104290.02660.681695-0.0517-0.0087-0.02510.0109510.04681.434299-0.0725-0.0281-0.04290.0127260.02390.571895-0.0628-0.0304-0.02500.0129830.07623.304099-0.0930-0.0810-0.04960.0159190.01740.5427

由表11可知,基于GED分布与基于正态分布假设的估计的结果比较相似,检验残差的LM效应,结果显示残差不再有显著的异方差现象,GED分布下的模型较好地拟合了数据的异方差现象。在5%的显著水平下,GARCH(1,1)模型、EGARCH(1,1)模型、TARCH(1,1)模型、CARCH(1,1)的LR检验值均小于3.841,可见4个模型全部通过了LR检验,意味着在GED分布下GARCH类模型得出的VAR值结果都比较精确。而在显著性为l%的情况下,LR检验值远小于6.635。所以,在两种置信水平下,GARCH(1,1)模型、EGARCH(1,1)模型、TARCH(1,1)模型、CARCH(1,1)模型都能够比较好地拟合创业板市场风险。

五、结 论

选取中国创业板市场综合指数,分别在正态分布、T分布、GED分布假设下,运用GARCH类模型计算VAR值并进行Kupiec检验,对创业板市场风险进行度量。对比分析不同显著性水平下的计算结果,得出以下结论:

第一,分析创业板块日收盘价对数收益率序列的统计特征,得出创业板综指日对数收益率有显著的尖峰厚尾特征,不服从正态分布,且其波动呈现出聚集性和持久性。

第二,正态分布下采用TARCH(1,1)模型和CARCH(1,1)模型计算VAR值时得到的结果较准确,而基于T分布的创业板综合指数计算出的VAR值过于保守,存在着风险被缩小的可能。只有GED分布在两种显著性水平下都能够很好地拟合数据,在GED分布下计算得出的VAR值的结果相对科学、精度高。四类GARCH模型计算得到的VAR值、失败率、失败天数、VAR标准差的差距都不大,从而可以看出模型种类并不是VAR值结果的决定性因素。

第三,创业板市场存在明显的杠杆效应即波动的非对称性,从模型参数可以看出利空消息带给创业板市场的波动比等量的利好消息更大,即负的冲击带来的波动大于正的冲击。

第四,在实际的投资决策过程中,投资者应选择适当的GARCH类模型才能更好地规避风险。投资者在创业板市场进行投资选择时,可以利用VAR进行风险评估和预测,根据理论基础选择投资观念,根据分析来实时更新投资策略。研究表明,VAR值随着置信水平的不同而有所不同,设定置信水平越高,风险价值VAR值也就越准确,也就意味着投资者可以在更准确地意识到风险的基础上有更大可能避免风险带来的损失。同时,创业板市场的杠杆效应也提醒投资者理性投资,不要盲目跟风。所以,投资者要根据实际情况为自己量身定制投资方案,及时更新投资方法和策略,灵活投资,这样才能保证风险和收益的最佳分配。

参考文献

[1]刘建民,张颖.创业板发展现状及完善对策 [J].人民论坛,2015(2):77-79.

[2]王璇.中国创业板市场风险分析 [J].北方经贸,2015(6):175-176.

[3]林博.基于GARCH模型的VAR方法在上证市场中的应用 [J].时代金融,2013(7):250-253.

[4]薛炜星.VAR在我国金融风险管理中的应用研究 [J].现代经济信息,2016(5):276-278.

[5]沈一凡.基于GARCH-VAR模型对我国股票市场的风险度量研究 [J].时代金融,2015(3):130-132.

[6]黄崇珍,曹奇.基于GARCH-VAR模型的开放式基金风险度量 [J].统计与决策,2017(1):152-155.

[7]Neil M.A comparison of GARCH models for VAR estimation [J].Expert Systems with Application,2016,39(3):47-58.

[8]Ramazan G.The use of GARCH models in VAR esti-mation [J].Statistical methodology,2015,8(15):4-12.

[9]Samit P,Prateek S.Improved VAR forecasts using extreme value theory with the realized GARCH model [J].Studies in Economics and Finance,2017,11(8):238-259.

[10]田岚,刘澄皇,甫玉婷.投资者情绪与新股长期表现研究:来自我国创业板市场的证据 [J].河北经贸大学学报,2016(4):84-91.

[11]王艳秋,张青,王甲山.基于VAR模型的中国石油消费影响因素分析 [J].沈阳工业大学学报(社会科学版),2017(1):64-70.

[12]洪晶晶,吴鹏跃.基于GARCH模型的金融市场波动性分析与预测 [J].价值工程,2016(1):46-47.

[13]任继勤,单晓彤,梁策.中国主板与创业板市场风险比较分析:基于GARCH-VAR方法 [J].财贸研究,2015(3):128-132.

[14]Carol A,Emese L,Silvia S.Forecasting VAR using analytic higher moments for GARCH processes [J].International Review of Financial Analysis,2013,5(6):36-45.

[15]任继勤,单晓彤,梁策.中国主板与创业板市场风险比较分析:基于GARCH-VAR方法 [J].财贸研究,2015(3):128-132.

[16]奉颖香.VAR测度股市风险的实证研究 [J].现代经济信息,2017(4):304-305.

RiskmeasurementofChineseGrowthEnterpriseMarketbasedonGARCHmodel

SONG Yong-hui, CHEN Chen

(School of Economics, Shenyang University of Technology, Shenyang 110870, China)

AbstractUsing the logarithm daily yields of the Chinese Growth Enterprise Market (GEM)comprehensive index as the basic data, the GARCH parameters are summed up and the VAR values are calculated under the normal distribution, theTdistribution and the GED distribution. Through the empirical analysis, it is pointed out that the sequence distribution has the characteristics of spikes and thick tails, the bad news can bring greater volatility to the growth enterprise market than the same amount of good news, the GEM market has obvious leverage effect, and the volatility has aggregation and persistence. The GED distribution can better fit the comprehensive index of the GEM market, and the accuracy of the VAR value is higher; the type of the model is not the decisive factor of the result of the VAR value.

KeywordsGrowth Enterprise Market (GEM); risk measurement; GARCH model; bad news; good news; leverage effect

网络出版地址:http:∥kns.cnki.net/kcms/detail/21.1558.C.20180320.1042.026.html

*本文已于2018-03-20 13∶06在中国知网优先数字出版。

作者简介宋永辉(1964-),男,辽宁沈阳人,教授,博士,主要从事风险投资、国际投融资等方面的研究。

基金项目辽宁省教育厅智库项目(201654025)。

收稿日期2017-10-20

doi:10.7688/j.issn.1674-0823.2018.05.08

文章编号:1674-0823(2018)05-0428-07

文献标志码:A

中图分类号F 425

(责任编辑:张 璐)