1、0学生姓名：何雨芹 学号：41108065学校：西南财经大学课程：金融统计分析 1我国银行间同业拆借市场利率风险度量基于 VaR 模型的实证研究摘 要本文利用 VaR 模型通过 2013 年 1 月 4 日至 2014 年 10 月 30 日我国银行间同业拆借市场每日加权平均利率进行实证研究，建立了基于 GARCH 模型的我国银行间同业拆借市场利率风险测度 GARCH 族模型（ GARCH(1，1)/TARCH(1，1)/EGARCH(1，1)）,得出以下结论: t 分布不适合描述我国银行间同业拆借利率序列的分布状况,广义误差分布能较好刻画我国银行间同业拆借利率序列的分布；根据样本数据，现阶段
2、我国银行间同业拆借利率风险也较低。 关键词：VaR 模型 同业市场拆借利率 GARCH 族模型 2一、 文献综述同业拆借市场是金融机构之间进行短期、临时性头寸调剂的市场，是货币市场的重要组成部分。1996 年以来, 我国银行同业拆借市场在中国人民银行的监督管理下稳步发展, 目前已形成全国统一的银行同业拆借市场格局，生成了全国统一的中国银行间同业拆借利率(China inter -bank offered rate,CHIBOR)，它是我国货币市场最早市场化的利率，也是目前唯一直接的市场利率，能够十分灵敏的反应市场上货币资金的供求状况，因此可视为我国货币市场上的基准利率。随着利率市场化的深人，利
3、率结构体系的完善和合理，我国金融市场的成熟，金融衍生工具的丰富，各商业银行在利率风险管理方面经验的逐渐成熟。商业银行可以逐步向先进的利率风险度量模型演进。从而跟国际先进风险管理水平接轨，增强自身的市场竞争力和抵御风险的能力。而银行间同业拆借是我国利率市场化改革的前沿阵地，银行间同业拆借市场利率市场化改革始于 1996 年，同业拆借利率市场化程度已经较高，具备了运用VaR 模型的客观条件。再加上 VaR 模型已经是一个比较成熟的模型，用它来研究我国银行间同业拆借市场应当是有一定研究价值的。VaR 模型源自马科维兹于 1952 年创立的基本均值方差模型，蒂尔.古尔迪曼被视为“风险价值”这一术语的创
4、立者，该理论一经提出就迅速得到学者的关注。国外已有很成熟的关于 VaR 的理论研究和实证研究， Jeremy Berkowitz（ 1999）提出了新的评价 VaR 的方法，Tean-Philippe Bouchaud 和 Marc Poters（1999）提出了如何利用金融资产波动的正态性去简化计算复杂的非线性组合 VaR；大部分学者在计算风险价值 VaR 值时，都以金融时间序列数据服从正态分布和无条件方差为假设前提，但是大量的实证研究表明，金融时间序列数据并不严格符合这一假定，为了解决这一问题，随着研究的不断深入，又有学者提出了半参数模型和广义条件异方差模型（GARCH 模型）等模型，大大
5、丰富了 VaR 的计算方法。Kees Koedijk（2001）将 VaR 风险管理模型应用于资产组合选择和资本资产定价，并指出由于资产组合收益率呈现出尖峰肥尾的特征，这会导致传统的均值-方差模型存在低估风险资产组合所面临的风险，可能会导致投资风险。在国内，近几年关于 VaR 的实证研究已经越来越丰富和深入，早在 2000 年初，国内就有学者王春峰、万海辉和李刚指出用蒙特卡洛模拟法计算 VaR 值所存在的缺陷，并提出基于马尔科夫链蒙特卡洛的计算方法。之后也有一批学者相继提出了对 VaR 计算方法的改进，同时，VaR 方法的应用研究开始受到重视。杜海涛（2000）在VaR 方法在证券风险管理中的
6、应用一文中在市场指数风险度量、单个证券的风险度量、基金管理人员绩效评价及确定配股价格等方面运用了 VaR 方法。他认为沪深两市的指数。单个证券、投资基金的收益都服从正态分布，在这一前提下去计算资产的 VaR 值并进行模型检验，得出了较好的结果。以他的研究为代表，早期的关于 VaR 的实证研究多集中在证券市场。迄今为止，将 VaR 方法运用于银行间同业拆借市场的研究还不太丰富，但就我刚才所说，在我国面临重大金融市场改革的前提下，银行间的同业市场越来越重要，对它进行风险度量分析时非常必要的。现有研究下，郑尧天和杜子平（2007）选择隔夜拆借利率为研究对象并分别用组合正态 VaR 方法和蒙特卡罗模拟
7、法对其进行建模，经后验区间检验发现蒙特卡罗模拟法的估计结果更为理想；冯科和王德全（2009）以 2002 年 6 月 4 日至 2009 年 3 月 31 日期间我国银行同业拆借利率为研究对象分别建立了隔夜拆借和 7 天拆借品种的预测模型，并度量了其利率风险。得出通过选择适当滞后阶数的 ARMA-GARCH 类模型，可以有效地刻画同业拆借利率的动态特性：t-分布和 g-分布下的模型能更好地捕捉同业拆借利率序列的尖峰厚尾3性，同业拆借利率存在显著的自相关性、风险溢价效应和波动的反杠杆效应，即利率上升时的波动更大，VaR 方法可以有效地预测同业拆借利率风险。多数研究都是围绕 GARCH模型展开的拓
8、展性研究，并且得出了不同的结论。二、 我国银行间同业拆借市场利率现状的分析与模型建立1996 年 6 月 1 日，中国人民银行取消了对同业拆借利率的上限管理，同业拆借利率由交易双方根据市场资金的供求状况自行确定，开始了我国同业拆借市场市场化改革。1996年拆借市场的交易量仅为 5871.58 亿元，至 2003 年交易量已达 24113 亿元，增幅惊人。银行间市场发展极为迅速。随着市场交易量的激增，同业拆借利率的不确定性增加。参与的商业银行面临巨大的利率风险。虽然拆借资金只是短期使用，但现实中商业银行的拆借资金已经拆除了调剂头寸的需求，同业拆借资金处于利率风险之中。因此，用风险价值 VaR模型
9、测度我国商业银行在同业拆借市场中的利率风险是很有价值的。（一）数据及其来源由于选取的是中国银行同业拆借市场中的隔夜拆借为研究对象，因此我选择了最新的数据，是2013年1月4日至2014年10月30日的同业拆借市场每日加权利率这457个观测值为样本数据，其中有些天数的数据是缺失的（数据来源为东方财富网数据中心http:/ 。（二）描述性统计分析在运用 VaR 模型测度利率风险之前，须首先对数据进行描述性统计分析和检验同业拆借利率序列的正态性、平稳性、自相关性和条件异方差性。1、 描述性统计首先对数据做描述性统计分析，运用软件 SPSS17.0，得出结果如下图所示：可以看出，同业拆借利率中最小的利
10、率为1.6922% ，最大利率为13.8284% ，波动幅度非常的大，这也印证前面所说同业拆借市场的风险也很大的结论。同业拆借利率的均值为3.109752%。另一方面，可以从表中看到同业拆借利率的偏度为3.669，为右偏分布；峰度为24.529，数据为尖峰分布。2、正态性检验对样本数据做正态 QQ 图和无趋势正态概率图如下，通过正态概率图可以看出 , 在正态直线以外散布着大量的点，数据点组成的线呈曲线状，且两端有摆动, 说明 CHIBOR 的实际分布两侧具有厚尾现象。通过无趋势正态概率图可知, 大多数散点并不是随机分布在通过零点的水平直线周围，而是呈现明显的抛物线形状，所以，正态 QQ 图和无
11、趋势正态概率图分析，可初步得出数据服从正态分布假设不成立。4对 CHIBOR 做单样本 K-S 检验，结果如下图，从结果可以看出，样本数据均值为3.109752，标准差为1.1130473，K-S 的 Z 统计量为3.304，对应的相伴概率为0.000，小于0.05的显著性水平，因此，正态性假设被拒绝，及同业拆借利率的实际分布与正态分布存在显著性差异。5通过上述检验, 我国银行间同业拆借利率分布是非正态的，如果仍然假设其服从正态分布，虽然可以大大简化运算，但误差会比较大。因而，我们要用其它分布形式代替正态分布以提高模型的拟合度。3、 平稳性检验对CHIBOR做单位根检验，检验结果如下，可以看出
12、ADF统计量值为-3.978563，P值为0.0017，因此拒绝原假设，这说明CHIBOR数据是平稳的。Null Hypothesis: CHIBOR has a unit rootExogenous: ConstantLag Length: 7 (Automatic - based on SIC, maxlag=17)t-Statistic Prob.*Augmented Dickey-Fuller test statistic -3.978563 0.0017Test critical values: 1% level -3.4446915% level -2.86775710% leve
13、l -2.5701454、 自相关检验自相关检验的目的是检验各期的拆借利率是否具有相关性。对于自相关检验的方法, 目前主要有以 Dickey-Fuller 检验为代表的序列相关法和以 Ljung-Box 统计量为代表的序列相关检验法。运用 EViews 求 CHIBOR 序列的自相关系数和偏相关系数及 Ljung-Box 统计量 Q-t 结果如下:6可以看出，CHIBOR序列具有较强的自相关性，因此求CHIBOR一阶差分序列的自相关系数、偏相关系数和Ljung-Box统计量Q-t 如下图：由图可知，CHIBOR一阶差分序列存在较弱的自相关现象，结合上一条结论，我们可以知道我国同业拆借市场利率存
14、在序列相关，即各期的同业拆借利率是相互影响的。5、 条件异方差检验金融时间序列变量的方差是计算 VaR 过程中最重要的参数，大多数学者计算 VaR 时，假设方差是常数，即无条件分布假设，然而，金融变量往往与无条件分布假设相异，方差具有时变性。由于统计学中常常以方差和标准差描述金融变量的波动性，为了准确测定金融时间序列的波动性，经济学家们努力构建各种模型，其中以 ARCH 模型和 GARCH 模型的应用最为广泛。大量前人的实证证明 GRACH 模型适合预测金融时间序列的波动性。因此，这里我也采用 GRACH 模型来估计我国银行间同业拆借利率间的波动性。7对我国银行间同业拆借利率求一阶差分后做时序
15、图（如上图） ，由图中可以看出波动具有明显的时变性，且不同时期波动性的大小也不相同，同时波动出现聚集现象，因而银行间同业拆借利率序列存在条件异方差。下面我们分析各期方差之间的关系，时序图所显示的波动聚集性说明相邻的波动性之间存在一定的相关性。根据近似估计，我们可以用以 代替方差 检验各期方差之E(chibor2) 2间是否存在相关性。下图是 CHIBOR 的平方序列的自相关检验结果：可以看出，CHIBOR 的平方序列依然与其滞后阶数存在比较强的自相关，存在自相关现象。8（三）模型建立通过检验同业拆借利率序列的正态性、平稳性、自相关性和条件异方差。我们可以初步断定银行间同业拆借利率序列存在自相关
16、、异方差等特性。GARCH 族模型包括很多模型，其中比较为人所知有 GARCH 模型，EGARCH 模型和 TGARCH 模型。这些模型都是有 ARCH 模型发展而来的，都适用于估计存在自回归条件异方差的时间序列。要进行模型的建立，首先就要确定滞后的阶数，分别将这三个模型用于数据，用赤池信息准则 AIC 来确定最有滞后阶数。同业拆借利率波动模型试算结果滞后阶数 GARCH AIC TARCH AIC EGARCH AIC（1,1） -2.6411 -2.5075 -2.6498（2,1） -2.6309 -2.4880 -2.6106（1,2） -2.6125 -2.4034 -2.5845（
17、2,2） -2.5797 -2.4422 -2.5981由于 AIC 准则要求仅当所增加的解释变量能够减少 AIC 值时才在原模型中增加该解释变量。在上述结果中，AIC 值都是随滞后阶数的增加增大（TRACH 和 EARCH 在（2,2）出现了 AIC 值的减小，但幅度很小） ，考虑到自由度，选取 GARCH 族群滞后（1,1）的模型。三、实证分析及模型检验在上文确定了模型滞后阶数的情况下，我们要进行实证检验，因为我们采用了GARCH 族模型中的三种模型，这意味着在这里我们也就有三个 VaR 的估计结果。我们在风险价值定义的基础上，将利率波动的风险价值 VaR 用公式表示为：VaRt+1=t(
18、E(r)-r*)=tct+1其中， 是 t 时期预测下一天的风险价值； 是利率的期望值；VaRt+1 E(r)是最低利率； 是在给定置信水平 c 时对应的分位数； 是 t 时t是 t时 期的 资产 价 值 ；r* c t+1期预测的下一天的同业拆借利率的波动性。分位数的取值依赖于利率的概率分布函数，在不同概率分布函数下，相同的置信水平c 对应的分位数存在差异。标准差口分为条件标准差和无条件标准差，无条件标准差假设利率的波动始终不变，条件标准差是对利率随时间波动的不一致性的反映，它是对利率风险的动态度量，能准确的把握风险的变化。本文对 VaR 的计算主要是通过建立估计条件标准差的动态模型而实现的
19、。关于 GARCH 族模型中的残差分布，一般会有 3 个假设: 正态分布、学生 t 分布和广义误差分布( GED )。选择哪种分布要根据时间序列的特征来定。正态分布是一种重要的概率分布，正态分布所具有的对称性、易加性等良好的统计特性决定了它在金融分析中的重要作用。本章在采用其它分布形式分析我国银行间同业拆借利率序列的风险价值 VaR 之前, 首先在正态分布假设下对拆借利率序列的风险价值 VaR 进行了计算，作为比较的依据。大量的金融实证研究表明，金融市场的市场因子（如价格、利率）序列的分布具有尖峰厚尾的特征，这与正态分布存在较大差异。通过对我国银行间同业拆借利率序列的正态性检验可知，拆借利率序
20、列并不服从正态分布且具有厚尾性。t 分布是比正态分布的概率密9度函数具有更厚尾部的一种分布，适合于分析具有厚尾特性的分布，本章在计算正态分布下的 VaR 之后，以 t 分布为假设，计算银行间同业拆借利率序列的风险价值 VaR。广义误差分布(GED) 是一种更为灵活的分布形式，通过对参数的调整可以适应不同的分布情形。它是由 J.P.Morgan 提出的，其概率密度函数既可以分析尾部比正态分布和 t 分布厚的情形，也可以分析尾部比正态分布薄的情形。在这里用参数法计算 VaR 时，分别假设了模型的残差服从正态分布、t 分布和滚滚一误差分布。下面就是计算结果。运用EView 软件分析GARCH族模型在
21、不同的分布情况下同业拆借利率的波动率（标准差），再将估计值带入VaR的计算公式，在 95%的置信水平下计算同业拆借利率的日均VaR，得到下表。模型估计结果模型 参数 正态分布 t 分布 广义误差分布0.000378 0.000195 0.000197(常数)（9.4768） （2.0794） （2.0097）0.768224 0.479672 0.438659（利好的影响系数）（17.9872） （6.1157） （6.136954）0.608172 0.730903 0.809832（系数）（25.7843） （24.8821） （19.8425）分位数（95%） 2.4679 2.7458
22、 2.3964GARCH(1，1)模型日均 VaR（95%）0.0067943Pt 0.0069801Pt 0.0055031Pt0.000301 0.000209 0.000193(常数)(4.6793) （1.8693） （1.6493）0.841312 0.559931 0.463320（利好的影响系数）(12.9783) （5.0032） （5.3248）-0.637952 -0.569327 -0.173856（刻画杠杆效 应 ）（-6.5548） （-1.7698） （-1.5099）0.647079 0.854671 -0.173856（系数）（35.9941） （24.3644
23、） （18.3291）分位数（95%） 2.4591 2.8643 2.5068TARCH(1，1)模型日均 VaR（95%）0.0061921Pt 0.7092304Pt 0.0064698Pt-0.983054 -0.973054 -0.649545(常数)（-12.0180） (-12.9843) (-6.4740)0.672970 0.684521 0.499100（利好的影响系数）（14.8752） (15.0081) (5.9481)0.297472 0.279451 0.178364（刻画杠杆效 应 ）（8.0949） (8.8495) (3.4215)1.003430 1.03
24、5698 0.954678（系数）（70.9844） (73.9151) (60.5125)分位数（95%） 2.6691 2.8161 2.4633EGARCH(1，1)模型日均 VaR（95%）0.0063571Pt 0.0074891Pt 0.0062871Pt10四、结论通过建立基于 GARCH(1，1)族模型的 VaR 模型，实证分析我国银行间同业拆借市场的利率风险，在分析前，我们已经可以通过基本的描述性统计看出我国银行间同业拆借利率波动是非常剧烈，这一方的原因是市场化改革的不断深入，市场化加深；另一方面却也反映了市场的不稳定。而我们通过确切的模型则可以得出更为数据化的结论：1、通过
25、比较 GARCH(1，1)模型、TARCH(1，1)模型和 EGARCH(1，1)模型的估计结果，可以看出，实际上三类模型得出的结论是没有太大的差别的。但 t 分布下三类模型计算的 VaR 值都大于正态分布和广义误差分布，且同时结合之前的正态性检验可知，t 分布不适合描述我国银行间同业拆借利率的分布状况，而广义误差分布能很好的刻画我国银行间同业拆借利率的分布。2、在分析第二个 GARCH 模型 TARCH(1，1) 时，我们可以从表中看到，由于获得的并不显著，说明我国银行间同业拆借利率不存在明显的杠杆效应；而在 EGARCH(1，1)中，系数 ，说明我国银行间同业拆借利率存在杠杆效应且利坏消息
26、影响大于利好消0且 0息。而两个模型的结论不同则说明是否存在杠杆效应并不明确。3、三个 GARCH 模型在计算 VaR 值时，几乎都得出了相近的 VaR 值，大约都在0.006 左右，这说明在数据选取区间我国银行间同业拆借利率风险较低，大概在 0.6%左右。Pt11参考文献1. 肖春来, 宋然. VaR 理论及其应用研究 J. 数理统计与管理, 2003, (2): 23-452. 沈悦. VAR 宏观计量经济模型的演变与最新发展基于 2011 年诺贝尔经济学奖得主 Smis 研究成果的拓展脉络N. 数量经济技术经济研究, 2012.10(3).3. 张娜, 黄新飞. 我国同业拆借市场利率的波
27、动性J. 统计与决策 , 2006, (8): 45-474. 谭晖. 基于 VaR 模型的中国商业银行利率风险测度的实证研究 以同业拆借市场为例D. 重庆 :重庆师范大学, 2012.5. 辛俚. 我国商业银行间同业拆借市场利率风险 VaR 度量实证研究 D. 成都:西南财经大学, 2012.6. 王德全. ARMA-GARCH 模型及 VaR 方法在我国银行间同业拆借市场中的应用研究J. 系统工程 , 2009, (5): 13-167. 李成, 马国校. VaR 模型在我国银行同业拆借市场中的应用研究J. 金融研究, 2007, (5): 46-488. 李志辉, 李志辉. 我国商业银行利率风险的度量研究以同业拆借市场为例EB/OL. http:/ 徐炜. GARCH 模型与 VaR 的度量研究EB/OL. http:/ 李良松. 上海银行间同业拆放利率 VaR 的有效性研究J. 金融研究, 2009, (9): 56-5811. 郑尧天, 杜子平 . 基于 VaR 模型的银行同业拆借利率风险估计 J. 工业技术经济, 2007, (12): 105-10712. 房小定, 吕鹏. 基于 GARCH 模型的上海同业拆借利率风险度量N. 西安电子科技大学学报(社会科学版), 2013(4).
