中国股市动态VaR计量模型分析

沪市动态VaR计量模型分析案例
西南财经大学统计学院马丹
一、案例背景
金融风险是指由于经济活动的不确定性而导致资金在筹措和运用中遭受损失的可能性。

在20世纪80年代以前，金融机构所面临的风险还主要是信用风险。

然而，近20多年来金融市场发生了重大变化，全球化的证券市场迅猛发展，资产证券化的趋势越来越强，外汇交易和衍生品交易成了金融市场交易的重要组成部分。

这使得金融机构所面临的主要风险已从信用风险转向了市场风险。

市场风险度量的方法有多种，VaR（Value-at-Risk）方法是目前金融市场风险测量的主流方法。

VaR是一定的概率水平下，证券组合在未来特定一端时间内的最大可能损失。

其优点在于将不同的市场因子、不同市场的风险集成为一个数，较准确测量由不同风险来源及其相互作用而产生的潜在损失，适应了金融市场发展的动态性、复杂性和全球整合性的趋势。

应该注意的是，VaR本身的含义是风险价值，是一个数值，但是它是由一系列统计方法实现的，因此在许多场合它又被称为风险度量方法。

VaR方法是由JP MORGAN公司率先提出来的，并在实践中得到广泛应用。

其在风险测量合管理中的巨大优点已为国际金融监管当局认可和接受，1990年以来的许多金融监管法案和原则都充分强调了基于VaR的风险监管方法。

巴塞尔委员会的巴塞尔银行业有效监管核心原则、欧盟的资本充足法令等都要求金融机构用VaR技术确定内部风险资本要求、风险控制等。

VaR计算方法包括历史模拟法、方差——斜方差法和蒙特卡罗模拟法。

历史模拟法和蒙特卡罗模拟法由于在技术上较为复杂，较难实现，因此，在实践中常用的方法仍然是方差——斜方差法。

方差——斜方差法中有两个需要注意的问题。

一是如何刻画金融数据的尖峰厚尾、波动簇集的时变特征。

二是如何寻找金融数据的分布密度函数。

Risk Mrtrics TM利用指数平滑技术来反映波动时变现象，并假定金融数据服从正态分布。

然而，大量文献资料证明，金融数据有强烈的ARCH效应，其尾部和中间部位集中了大量的概率分布，比正态分布拥有“肥尾”特性。

金融数据的集群性，异方差性等特征显然违背了古典假设，利用传统的基于古典假设的回归模型，OLS建模方法难以刻画出真实的数量规律，所做出的统计推断也是不精确的。

为了对金融数据做出有效的描述，Engle于1982提出了ARCH模型，模拟出数据的集群性特征，他的学生Bollerslev在1986年提出了GARCH （Generalized Auto Regressive Conditional Heteroscedasiticity）模型[1]，将高阶的ARCH模型转化成为简洁的GARCH模型,描绘出金融数据方差项的某种自相关性。

但实际应用中，GARCH模型虽然有助于模拟金融数据分布的宽尾特征，却不能解释金融市场上存在的杠杆效应。

针对这一问题，近年来Nelson等人提出了非对称性（Asymmetric）GARCH模型以刻画出条件方差对正的价格变化反应弱而对负的价格变化反应强这一现象。

这类模型包括TARCH、EGARCH等，收益率序列残差往往假定为正态分布，但实际上正态GARCH模型不能充分描述数据的尖峰厚尾性。

对此，可以假定残差服从t分布、混合正态分布或一般误差分布。

本案例讨论了基于Normol-GARCH、t-GARCH的中国股市动态VaR测度方法，并比较了两种假定下VaR的估计效果。

二、GARCH模型与VaR计算方法简介
1、VaR 计算模型
VaR （Value-at-Risk ）是一定时期内，在一定的置信度下，投资组合可能出现的最大损失，用公式表示为：
VaR=()
(
)
111---+-p F t e
p σμ （1）
其中，()∙F 是收益序列的分布函数，p 为给定的显著性水平，即左尾概率。

如果是静态VaR 只需利用极大似然法估计出参数μ和σ，将之代入计算公式即可。

2、基于GARCH 模型的VaR 计算方法
VaR 估计的条件方差方法属于动态VaR 计算的分析方法,由于实际金融市场中收益率的厚尾性会导致VaR 对风险的低估,因此可以利用GARCH 模型类中的条件方差来度量股票市场VaR 。

这样，VaR t t h z p α1-=，其中t h 是由GARCH 模型估计得到的条件方差， αz 根据收益率分布决定。

GARCH 模型的一般表达式可写成：
t
t t i t q
i i j t p
j j t t
t h h h x r νεεαβωεψδβ=++=++=-=-=∑∑21
1
/ （2）
其中，t h 为条件方差,t v 为独立同分布的随机变量,t h 与t v 互相独立。

一般常假定t v 为标准正态分布。

但正态性往往不足以反映金融数据的厚尾性，Nelson 和Hamilton 提出用广义误差分布（GED ）和t 分布来刻画尖峰厚尾特性。

t v 服从t 分布时，相应的对数似然函数为：
()()()()∑=⎥⎦
⎤⎢⎣⎡⎪⎪⎭⎫ ⎝⎛-++---⎪⎭⎫ ⎝⎛⎪⎭⎫ ⎝⎛Γ-⎪⎭⎫ ⎝⎛⎪⎭⎫ ⎝⎛+Γ=N
t t t v h v h v v v L 1221log 1212log 21
2log 21log ε
这里的()∙Γ是GAMMA 函数，v 是自由度。

似然函数可通过对偶牛顿算法或信赖域算
法极大化得到。

由于股市收益率的波动回随信息的变化而出现非对称性的特点，利空消息引
起的波动一般更大，为了刻画波动的非对称性，Nelson 等人提出了非对称性（Asymmetric ）GARCH 模型，包括EGARCH 及TGARCH 等模型。

在EGARCH 模型中，条件方差t h 为延迟扰乱i t -ε的反对称函数：
()()()∑∑=--=++=p
j j t i t q i i t h z g h 1
1
ln ln βαω，()[]t i t t i t z E z z z g -+=--γθ，t
t
t h z ε=
在TGARCH 模型中，条件方差t h 表示为：
∑∑=---=-+++=q
i i t J t t p j j t t d h h 1
212
11
εαγεβω 其中 0,0;0,11111>=<=----t t t t d d εε
另外，以上模型中
∑∑==+q
i i p
j j
1
1
αβ
<1，蕴涵GARCH 过程为宽平稳。

如果将此限制放宽
即是积分型GARCH 模型，即 IGARCH 模型。

对于IGARCH 模型而言，无条件方差并不存在。

同时，考虑到股票市场的风险溢价，在均值方程中加入条件方差的函数，得到GARCH-M 模型。

一般来说，有效的股票市场上股价与其风险是正向相关的。

GARCH 模型建模程序包括：
（1） 用Augunmented-Dikey Fuller 方法检验序列平稳性。

（2） 用Q 统计量识别序列自相关或偏自相关的阶数。

（3） 用Lagrange multiplier 方法检验序列残差项的异方差性以确定条件方差方程数。

（4） 用Most —Likelihood 方法得到方程中参数的最大似然估计。

（5） 评价模型。

3、VaR 模型检验方法
对VaR 模型可采用Kupiec 检验。

设N 为检验样本中损失高于VaR 的次数，T 为检验样本总数，P=1-C ，C 是即定的置信水平。

则检验的假设为：
似然比统计量为：
在原假设下，lr 服从于自由度为1的2
χ分布。

在大样本条件下，也可以用正态分布来逼近，同样有较好的检验效果。

三、沪市动态VaR 模型的实证分析
1、数据选取和统计描述
本案例以上海证券市场综合指数为例，说明VaR 的计算。

考虑到中国股市是从1996年12月16日起实行涨跌停板限价交易制度的，采用此后的数据可以减少异常值的干扰，提高模型的精度。

因此收集的数据为1996年12月16日至2002年10月22日的上证综合指数收盘价格,共计1404个。

数据来源于中国股票市场交易数据库和证券之星（，即设第t 日的收益率为t r ，则1ln ln --=t t t p p r ，t p 是t 日上证综合指数的收盘价格。

以下模型中的参数估计均是用SAS6.12统计软件实现的。

如图1所示，上证综合指数收益率集中在（-0.05，0.05）之间，在波动较大的几个时段上收益率绝对值的峰值达到了0.1,在到达顶峰后又迅速回落。

在样本期间，平均收益率为0.00284，偏度为-0.256529，峰度为9.401051，数据分布向左偏移，有一个沉重的尾巴。

这说明从整体上看，收益率低于其均值的时候较多，波动性强，尾部和中间包含了大量的统计信息。

J-B 正态性检验得到的统计量值为2558.392，也说明收益率序列显著异于正态分布。

2、 GARCH 模型族的建立 利用Q 检验和LM 检验识别序列的异方差性，结果见表1 。