- 1 -
酿酒食品板指数日收益率波动性特征及其
分析
戴屹,雷钦礼**
(暨南大学经济学院,广州 510632) 5
摘要:本文采用 GARCH模型,对我国酿酒食品板指数的日收益率进行研究。研究表明:投
资酿酒食品板指数的平均收益率大于投资大盘指数的;酿酒食品板指日收益率序列是平稳序
列,具有尖峰厚尾特征,波动呈现聚集性和异方差性;非对称的 TARCH(1,1)模型模拟
效果最好,酿酒食品板指存在负的杠杆效应。投资者投机目的强,但投机成分在减少,日益
变得理性。 10
关键词:酿酒食品板指;波动性;GARCH模型
中图分类号:
The Characteristics and Analysis of the Daily Yield Rate of
Chinese Brewing Food Plate Index 15
DAI Yi, LEI Qinli
(Economy School, Jinan University, GuangZhou 510632)
Abstract: In this paper, We studied the daily yield rate of Chinese brewing food plate Index by
using the GARCH model. The results show that the average rate of return on investing brewing
food plate index outweighs the market index. Daily yield rate series of brewing food plate Index is 20
stationary, characterized by High Kurtosis and Fat Tails, and its volatility is also clustered and
heteroscedastic. The effect of simulation is the best by using asymmetric TARCH (1, 1) model,
Brewing food plate exists negative leverage effect. Investors’ speculation purposes are strong, but
speculation component are reducing, and becoming more rational.
Keywords: Brewing food plate Index; Volatility; GARCH Model 25
0 引言
随着居民可支配收入的增加,人们对酿酒食品类消费支出逐渐加大。近年来,酿酒食品
行业迅速发展,使其在资本市场屡受机构投资者青睐,更被广大投资者追捧。而投资就面临
市场波动的风险,因此研究酿酒食品板块的波动性特征,对认识酿酒食品行业的整体风险,30
做出正确的投资决策,以及为市场监管者防范和控制风险提供依据,具有重要的意义。
波动性作为金融资产风险大小的度量方法之一,一直以来是金融领域研究的重点。
Engle(1982) [1] 提出方差随时间变化的自回归条件异方差模型(ARCH),为解决传统回归模型
碰到的异方差问题提供了新思路。Bollerslev(1986) [2]在 Engle 的基础上,借助 ARMA 模型思
想,对异方差的表现形式进行线性扩展,形成应用更为广泛的 GARCH 模型。为捕捉金融市35
场波动的非对称性,研究正、负资产收益率的非对称效应,Nelson(1991) [3] 提出了指数
GARCH(EGARCH)模型。随后,许多经济学家对上述模型进一步拓展完善,形成了以
EGARCH 模型、TARCH 模型、GARCH-M 模型等为主要模型的 GARCH 族模型。在实证研
究方面,我国学者岳朝龙(2006)[4]利用 GARCH 模型族实证分析了上海股市收益率的波动特
征,得出了上海股市收益率不仅具有条件异方差性,而且还存在杠杆效应。王明照,郭冰40
(2004)[5]借助 GARCH 族模型研究了沪深股市收益率波动性特征,结果表明沪深股市收益率
- 2 -
序列不服从正态分布,波动确实存在集群性,尖峰厚尾性显著,我国股市的波动性与收益率
显著正相关。聂富强,宋国军(2007)[6]运用 GARCH 模型的在引入两个虚拟变量的情况下对
上证指数股票收益率序列建模,得出股市不具有杠杠效应。孟卫东等(2010) [7]采用 GARCH
类模型对沪深地产指数收益率的波动进行了估计和预测,发现沪深地产指数收益率的波动不45
存在杠杆效应,投资者投机目的较强,GARCH(1,1)-M 模型的样本外预测效果最好。
目前,国内学者对于 GARCH 族模型的应用,大多以沪深两市大盘指数收益率作为研究
对象,对行业板块指数的研究则显得不足。在具体 GARCH 模型的运用上,也大多假定均值
方程的残差分布为高斯正态分布,其他分布形式考虑较少。为刻画酿酒食品板市场波动的风
险,分析其独特的波动性特征,本文力图应用 GARCH 族模型对我国酿酒食品板块日收益率50
序列波动性进行相应的实证检验,得出研究结论,希望对股市行业板块波动实证研究有所贡
献。
1 理论模型简介
ARCH 模型
自回归条件异方差模型,即 ARCH 模型由 Engle(1982)开创性的提出,它的基本思想为55
干扰项 tε 的方差 2tσ 取决于其前期值的大小, q 阶 ARCH 模型 ARCH(q)可表示如下:
2
0
22
22
2
110
2 )(... tqtqttttt
ttt
ttt
B
xy
εααεαεαεαασ
νσε
ελ
+=++++=
=
+=
−−−−−
其 中 为 ty 为 被 解 释 变 量 , tx 为 解 释 变 量 , tε 为 残 差 项 ,
)1,0(,0,0 ,,0 N
dii
tt ⎯⎯→⎯≥> ναα 。
GARCH 模型和 GARCH-M 模型 60
GARCH 模型就是广义的条件异方差,它是 ARCH 模型的拓展,GARCH 通过用一个或
者两个 2tσ 的滞后值代替许多 2tε 的滞后值,更好的描述了波动的持续性,在 GARCH 模型中,
考虑两个不同的设定:一个是均值条件,另一个是条件方差,GARCH(p, q)模型表示如下:
22
11
22
22
2
110
2 ...... ptpttqtqttttt
ttt
ttt xy
−−−−−−−− +++++++=
=
+=
σβσβεαεαεαασ
νσε
ελ
= ∑ ∑
= =
−− ++
q
i
p
j
jtjiti
1 1
22
0 σβεαα 65
其中 q 是 ARCH 项的阶数,p 是自回归 GARCH 项的阶数, 00 >α 且 0,0 >> ij αβ ,
a(B)和 b(B)为后移算子多项式,在大多数实证研究中,一般取 2,2 ≤≤ qp 既满足模型的灵
活性和简约性,也可以满足对收益率的解释。
在资本市场中,证券的收益率与它的波动率成正比,波动率越大,期望的收益也越大。
在这种认识的基础上,Engle 等人把条件方差项 2tσ 加入均值方程中,建立了 GARCH-M 模70
型:
- 3 -
中国科技论文在线
2
0
22
22
2
110
2
2
)(...
)(
tqtqttttt
ttt
tttt
B
fxy
εααεαεαεαασ
νσε
εσρλ
+=++++=
=
++=
−−−−−
其中,参数ρ 叫做风险溢价参数,代表了风险和收益之间的一种权衡。当 1>ρ 时,表
明投资者是风险厌恶型, 1<ρ 表明投资者是风险偏好型, 1=ρ 则属于风险中性者。
TARCH 模型 75
对金融资产的价格研究发现,在相等的价格上涨和下跌幅度上,收益率的下降往往出现
更大的波动,为解释这一波动率的非对称性,Zakoian(1990)提出了门限 ARCH(Threshold
ARCH)模型,在条件方差方程中引入非对称项 kL ,模型表示如下:
∑∑ ∑
=
−
−
−−
= =
−− +++=
r
k
ktktktk
q
i
p
j
jtjitit IL
1
2
1 1
22
0
2 εεσβεαασ
其中, −−ktI 是一个虚拟变量,出现利空消息时, 0<−ktε , 1=−−ktI ;出现利好消息80
时 0>−ktε , 0=−−ktI 。当 0>kL 时,利空消息会比利好消息带来更大的波动,就存在杠
杠效应。
EGARCH 模型
上述模型的条件分布可以是正态分布,t 分布,进一步的改进允许具有比前面假设的二次
方程映射更加灵活的关系。Nelson(1991)在条件方差中采用自然对数形式,提出另外一个非85
对称模型—EGARCH 模型:
)ln()ln( 2
11
0
2
jt
p
j
j
n
k kt
kt
k
it
it
it
it
t LE −
== −
−
−
−
−
− ∑∑ ++⎟⎟⎠
⎞
⎜⎜⎝
⎛+= σβσ
ε
σ
ε
σ
εασ
对条件方差取对数,这意味着 2tσ 非负,杠杠效应呈指数形式。从方程表达式中可以看
出,只要 0≠kL ,信息冲击的影响就存在着非对称性, 0<kL 说明杠杠效应显著。同样,
可以看出 EGARCH 模型不要求条件方差方程右端系数都为正,这样在求解过程中更为简洁、90
灵活。
2 实证分析
数据来源及基本统计特征分析
(1)数据来源。本文以酿酒食品板块指数为研究对象,选取 2006 年 8 月 30 日至 2012
年 7 月 2 日共 1418 个交易日收盘指数的数据(数据由大智慧软件导出)。指数收益率由相95
邻两天收盘价自然对数的一阶差分来表示,收益率 tR 计算公式为 1loglog −−= ttt PPR ,其中
tP 为日指数的收盘价,本文采用 Eviews[8]作为数据分析软件。
(2)数据的描述分析。对收益率序列做基本统计分析,发现样本期内酿酒食品板块的
收益率序列的均值为 %,标准差为 %,偏度系数为,小于零,呈负偏,
意味着投资者在此期间投资酿酒食品的收益看涨的可能性小于看跌的可能性,即收益率存在100
非对称性。另外序列峰度系数为 ,大于 3,说明酿酒食品收益率序列存在“尖峰厚
- 4 -
中国科技论文在线
尾”特性,同时 J-B 统计量的值非常大,为 ,其尾概率为 0,这些充分拒绝收益率
序列的正态分布假设。对同时期的上证综指日收益率序列做统计分析,知大盘收益率均值为
%,标准差为 %。
表 1 酿酒食品板指与上证综指日收益率基本统计量 105
The Basic statistics of the Daily Yield Rate of Chinese Brewing Food Plate Index and Shanghai Composite
Index
统计量 酿酒食品板指日收益率 上证综指日收益率
均值
标准差
偏度
峰度
J-B 值
从上述两组数据的对比中可以发现,该时期股市投资,整体收益率偏低,在标准差几乎
同等变动条件下,投资酿酒食品行业的收益率是大盘收益率的近 5 倍。究其原因在于,本文110
的研究时段正好经历了酿酒食品行业飞速发展,市场对酒类尤其是高端白酒消费需求旺盛,
上市公司利润大幅增长。此后大盘受美国金融危机和欧债危机的影响,股市下跌,危机对整
个股市的负面冲击比较大,主要表现为平均收益低和波动性较大。然而酿酒食品行业是传统
的保守型行业,在众多行业板块中,上市公司相对较多,市场需求比较刚性,受外部影响也
要小得多,所以在危机发生之后,经历短暂的低落,表现远好过大盘。 115
平稳性检验及均值方程
(1)平稳性检验。对酿酒食品快板指数收益率做单位根检验,通过 ADF 检验来测试其
稳定性。ADF 检验假定收益率在 0 上下波动,因此在计算 ADF 统计量过程中假定回归方程
不含常数项和时间趋势项。计算得到 ADF 统计量为,显著的小于临界值值。故
拒绝收益率序列存在单位根的原假设,可见酿酒食品板块指数收益率是平稳序列。 120
(2)均值方程。由于收益率序列是平稳的,因此可以对其建立时间序列模型。通过对
收益率做自相关和偏相关分析,发现 tR 与其滞后 38 阶存在明显序列相关关系。应用最小二
乘法对酿酒食品指数日收益率进行回归,并剔除不显著项,建立均值方程如下:
squared-R )-(F Prob
() )(
38-t
===
++= ttr ε
方程下方括号内的值为参数估计对应的t值,由此可知各参数都通过显著性检验。从 2R125
的数值来看,方程的OLS拟合的效果不是很好,说明误差项可能存在条件异方差。通过估计
方程的残差分析,发现残差项存在着波动性聚集现象,即大的波动集中在几个时段,小的变
化也集中在几个时段,大波动时常伴随着小波动,导致这种变化很可能原因在于股票市场的
波动性容易受国外危机、谣言、政局变动的影响,而且有时候突然的政府货币政策与财政政
策变化也容易引起市场波动,从而有理由相信误差项的条件方差不是某个自变量的函数,而是130
随时间变化并且依赖于过去误差的大小。
ARCH 检验及误差分布估计
(1)ARCH检验。运用GARCH 模型的一个前提是所分析的时间序列具有ARCH效应,
ARCH效应检验最常用的是LM 检验。通过对残差做ARCH-LM检验,发现F和LM统计量
对应的P值均为0,故拒绝原假设,认为存在ARCH效应,该序列可用GARCH模型进行分析。 135
- 5 -
中国科技论文在线
(2)误差分布估计。利用GARCH族模型对条件异方差进行检验和估计时,通常的一个
假设是残差分布服从高斯正态分布,然而在分析金融时间序列的过程中,正态分布不再是常
态,往往会出现多种分布形式,比如t分布,广义误差分布(GED)等。因此在模型建立之
前,有必要对残差分布进行检验。通过对残差做基本统计分析,发现J-B统计量值为,
对应的P值为0。由此可知,均值方程的残差不服从正态分布,对残差序列作Q-Q散点图,比140
较可知,t分布下的拟合情况更好。
-6
-4
-2
0
2
4
6
.00 .05 .10
N
or
m
al
Q
ua
nt
ile
Quantile-Quantile of RESID
图1 正态分布下残差Q-Q散点图
Residual Q - Q Scatter diagram Under Normal Distribution
145
-6
-4
-2
0
2
4
6
.00 .05 .10
St
ud
en
t's
t
Q
ua
nt
ile
Quantile-Quantile of RESID
图2 t分布下残差Q-Q散点图
Residual Q - Q Scatter Diagram Under Student’s t Distribution
GARCH 族模型估计分析 150
由上述分析可知,假定残差分布为 t 分布更为合适,因此可构建残差分布基于 t 分布的
GARCH 族模型。本文采用在实证研究中最有代表性的 GARCH、GARCH-M、TGARCH、
EGARCH 模型对酿酒食品板块指数日收益率的波动性进行分析。在 2,2 ≤≤ qp 条件下,根
据 AIC 和 SC 越小以及 ARCH-LM 检验的统计量相伴概率越大,模型解释变量滞后阶数越合
适的原则,确定模型阶数,结合方程中估计参数的显著性,得出模型估计结果见表 1。 155
- 6 -
中国科技论文在线
表 2 t 分布下 GARCH 族模型的参数估计
The Parameter Estimation of GARCH under T Distribution
模型
参数
GARCH
(1,1)
GARCH—M
(1,1)
TARCH
(1,1)
EGARCH
(1,1)
均值方程
λ
[]
[]
[]
[] ρ
[]
方差方程
0α 510−×
[]
510−×
[]
510−×
[]
510−×
[]
1α [] [] [] []
1β [] [] [] []
kL []
[]
AIC 值
对数似然函数值
ARCH( 2)检验
Obs*R-squared
[]
[]
[]
[]
注:[ ]内为各参数对应的p值。
160
从表 1 中可以看出,除 EGARCH(1,1)条件均值方程中的系数λ不满足 5%检验水平
下显著外,各模型的其他参数估计结果都是显著的。对残差序列进行 ARCH 检验,结果显示,
残差序列不存在 ARCH 效应,说明利用上述 GARCH 族模型消除了酿酒食品板指日收益率
序列的条件异方差性。
对比各模型的估计结果发现,在 GARCH(1,1)模型中,方差方程 ARCH 项系数 1α 与165
GARCH 项系数 1β 之和为 ,满足参数约束条件,模型是平稳的,可认为该方程对数据
有较好拟合效果,系数非常接近 1,说明酿酒食品板收益率对外部冲击的反应以一个相对较
慢的速度递减,一次冲击对股价造成的影响将有无限延伸下去的趋势,波动具有较长时间的
记忆性。另外,ARCH 项和 GARCH 项的系数都很显著,说明酿酒食品板块指数收益率有很
明显的异方差,且波动呈“聚集性”特征。在 GARCH-M(1,1)模型中,风险溢价系数ρ 等170
于 ,在 0 与 1 之间,说明期望收益与风险成正相关关系,投资者偏好风险,投机目
的较强。在 TARCH 模型中,杠杆效应系数 kL =>0,这表示当出现利好消息时,会
给板块价格指数带来一个 倍的正面冲击,而当出现利空消息时,则会带来一个
+ 倍的负面冲击。 在 EGARCH(1,1)模型中,杠杆效应系数 kL =
<0,符号与 TARCH(1,1)模型的相反,但根据对数的性质,得出的结论是一致的,175
当出现利好消息时,会给板块价格指数带来一个 倍的正面冲击,而当出
现利空消息时,则会带来一个 + 倍的负面冲击,进一步证明存酿酒食品板
块股票价格波动存在杠杆效应。
通过对酿酒食品板指数对数收益率序列建立的 GARCH 族模型分析可以发现,TARCH
(1,1)模型不仅各个参数的系数都很显著,而且有更小的 AIC 值和更大的对数似然值,180
能很好的拟合我国酿酒食品板指日收益率的时间序列,故本文将 TARCH(1,1)模型方程
作为对我国酿酒食品板块指数日收益率序列的波动方程。
- 7 -
中国科技论文在线
3 结论
本文通过 GARCH 族模型对我国酿酒食品板块指数日收益率的波动性进行实证分析,得
出的主要结论和相应的解析如下: 185
第一、我国股市酿酒食品板块收益率序列不服从高斯正态分布,存在尖峰厚尾特征,波
动呈现聚集性和异方差性。对比同时期的大盘收益率序列发现,投资酿酒食品板块股票的风
险与投资大盘股风险相近,但收益率却是大盘的近 5 倍。
第二、GARCH 族模型的持续性参数接近于 1,说明酿酒食品板块股票价格的波动受过
去冲击影响较强,收益率的波动具有长记忆性。 190
第三、GARCH-M 模型的均值方程中风险溢价系数为 ,收益率与市场波动存在
正相关性。较大的风险溢价系数意味着投资酿酒食品板块股票的投资者风险意识不断加强,
投资者正由过去的盲目跟从日趋走向理性投资,对承担市场风险要求得到更多的风险补偿。
第四、酿酒食品板指数收益率具有明显的负杠杠效应,利空消息引起的波动比同等大小
的利多消息引起的波动要大。究其原因,一是消息的传导机制并不是每次都有效的,二是我195
国证券市场尚未达到弱式有效。
第五、TARCH(1,1)模型是所有 GARCH 族模型中对酿酒食品板指波动性拟合效果
最好的。
[参考文献] (References) 200
[1] Engle. Auto-regressive Conditional Heteroskedasticity with Estimates of the Variance of Inflation[J].
Econometrica, 1982, 5(50):987-1008.
[2] Bollerslev,T. Generalized Autoregressive Conditional Heteroskedasticity[J].Journal of
Econometrics,1986,12(31): 307-327.
[3] Nelson,. Conditional heteroskedasticity in Asset Returns: A New Approach[J],Econometrica, 205
1991,20(59):347-370.
[4] 岳朝龙.上海股市收益率 GARCH 模型族的实证研究 [J]. 2006,06(10):126-129.
[5] 王明照,郭冰.沪深两市波动性特征的实证研究[J].统计与决策. 2004,12(15):98-99.
[6] 聂富强,宋国军. 沪、深股市波动不对称性的实证分析[J]. 数理统计与管理;2007,12(5):178-183.
[7] 孟卫东 , 宋丽伟 , 阳军 . 基于 GARCH 模型的沪深地产股波动性分析及预测 [J]. 统计与决210
策.2010,12(3):68-71.
[8] 高铁梅.计量经济分析方法与建模:Eviews 应用及实例[M].北京:清华大学出版社,2009.