内容摘要：本文利用GARCH 和VaR模型对上海国债市场的个案进行分析，具体包括GARCH模型的事前估计和分析，模型参数估计与最优选择，GARCH模型计算VaR。结果显示该国债在样本期内的变动性状，并结合实际政策进行了分析和论证。
　　要害词：VaR GARCH 沪市国债
　　
　　本文基于GARCH和VaR模型对沪市国债市场的个案进行了相关分析，以期显示国债在样本期内的变动性状。
　　
　　VaR模型
　　
　　VaR是国际上新近发展起来的一种卓有成效的风险量化技术，它是英文Value-At-Risk的缩写，中文通常译为在险价值， VaR模型通常被称为风险估值模型。对它的一种较为通俗的定义是：未来一定时间内，在给定的条件下，任何一种金融工具和品种的市场价格的潜在最大损失。确切的说，VaR可表示为：
　　Prob△t（△P＞VaR）=1-C
　　其中△P为证券组合在持有期△t内的损失，C为风险测量的置信水平，VaR为此置信水平下处于风险中的价值。根据定义，VaR实际上是要度量估测“正常”情况下资产或资产组合的预期价值与在一定置信区间下的最低价值之差，即
　　VaR=wo[E(r)-r*] （1）
　　其中，wo为持有期初资产组合的价值，r为收益率，E（r）为资产组合的预期收益率，r*为一定置信区间C下最低的收益率。
　　计算VaR最要害也是最困难的问题是确定资产报酬的分布形式，以找到特定分布的一定置信区间C下的最低收益率r*。假如收益率序列{r}服从正态分布，要想求出给定置信水平C下的r*，只要利用正态分布表找到标准正态分布的一个上分位点Z，使得：
　　 ，
　　根据正态分布的性质得出：
　　其中，μ△t为收益率序列的期望，而
　　 为标准差，即可得到与置信度C相对应的r*，即
　　 （2）
　　将公式（2）代入公式（1），可得到一般求解VaR的方差协方差模型，即
　　 （3）
　　由此可知，一般的测算VaR方差协方差模型，通常假定收益率序列服从正态分布。而现实中，国债收益率序列是否满足这一假定，需要进行收益率序列的正态性检验。目前常用的估算VaR模型波动率的方法有移动平均法、GARCH模型和隐含波动率法，移动平均法假设资产收益率的波动服从白噪声过程，这一假设与实际的观察结果经常不一致，隐含波动率则只能运用于有期权产品的资产，因此本文采用GARCH模型方法。
　　
　　GARCH 模型
　　
　　近几年，许多金融实际数据表明：市场在一定时期内有较大波动，为解决实际分析中存在的异方差现象，ENGLE于上世纪80年代初提出了ARCH模型(Autoregressive Conditional Heteroskedasticity)，即自回归条件异方差模型。Bollerslev又于1986年提出广义的ARCH模型，即GARCH模型(Generalized ARCH)，其形式为：
　　 (4)
　　 (5)
　　其中rt为收益率的时间序列，Vart-1(εt)为εt的条件方差。（5）式中系数有关系：
　　 ＜1，α＞0，Gi＞0，Aj＞0
　　我们称这样的模型为GARCH（P，Q）模型，模型中的（4）式描述了模型的条件期望部分，投资者在交易中所得到的信息依靠于过去时刻的收益率以及过去时刻预期收益和实际收益间的误差；（5）式描述了模型的条件方差部分，它不仅是滞后随机扰动项平方的线性函数，还是滞后条件方差的线性函数，表明了过去时刻的波动对未来价格波动有着正向缓解的影响，从而模拟了波动的集群性。
　　
　　实证前的分析
　　
　　由于上证国债指数自2003年1月2日起对外发布，可考察样本期太小。故挑选其中个股考察，经比较，20国债（4）具有较强的代表性，考察其从2001年1月2日到2004年5月11日每日的收盘价，共795个样本。关于收益率公式，本文采用：
　　 （6）
　　统计样本均值为-3.1567×10-5，其标准差为0.0021，其样本偏度为-0.0310，其样本峰度为13.0043，从中可以看出国债收益率序列的经验分布频率与相应的正态分布相比，具有尖峰宽尾的特征，这正是非线性的显著特性。这说明不能用简单的正态分布来模拟国债收益率的变化。考虑到条件正态分布既保留了正态分布的特点又能更好地对收益率进行模拟，在此引入条件正态分布，并应用GARCH模型来计算时变的条件方差ht，代入公式（3），可得出计算我国国债市场的VaR值：
　　VaRt=Pt-1Zα （7）
　　其中，Pt-1为前一日国债价格，Zα是置信度为标准正态分布的临界值，当α=0.05时，Zα=1.645，ht是国债收益率序列的条件方差。 则是标准差。
　　使用自相关函数（ACF）和偏自相关函数(PACF)对收益率序列的自相关性进行分析后，发现在大部分时滞上函数值在置信区间内0四周上下波动，所以收益率序列并不具有自相关性，因此在条件期望模型中不需要引入自相关性的描述部分，于是收益率rt由一个常数项C加上扰动项εt组成，其中Eεt=0，即rt=C εt。
　　尽管收益率序列的ACF值揭示了其弱相关性，但收益率的平方的ACF值却表现出了一定的自相关性，在时滞为1时呈现显著自相关性。最后对相关性进行定量检验，使用ARCH检验和Q检验。
　　在不存在序列相关的假设下，Q统计量和ARCH检验的统计量近似服从x2分布。对收益率序列进行Q检验，发现时滞为10，12，15时的P值分别为0.2747，0.2560，0.2938均大于0.05的显著性水平，说明接收原假设，不存在明显的序列相关性。
　　但是在对收益率平方进行Q检验时，在0.05的显著水平下时滞为10，12，15的P值均为0，收益率方差存在非常明显的序列相关性。结果证实了收益率序列具有GARCH效应，说明国债价格收益率序列具有异方差性，因此使用GARCH建模是合适的。
　　
　　模型参数估计与选择
　　
　　GARCH(1，1)模型参数估计结果如下：rt=-0.00010954 εt；
　　 (-1.7528)
　　σt2 =1.6117×10-6 0.1851σ 2 t-1 0.599εt2
　　 (13.6621) (4.7064) (10.2098)
　　括号中为t统计量，可以看到G1，A1在0.05的显著水平下明显地异于零，说明此国债收益率过去时刻的波动大小对当前波动大小有明显的影响，此外G1 A1<1，但并不是很接近于1，说明条件方差序列具有不是很强的长记忆性，表明了收益率序列的波动的持续性比较高，且投机因素较强，总体风险大，但是并不是非常大，总体来说中等偏上。
　　对模型进行最优选择，使用AIC和BIC信息准则来比较可供选择的GARCH模型。对各种模型的检验如表1，发现随着参数的增加，AIC，BIC值并没有显著增加，直观上说明可选用较为简单的GARCH(1，1)来建模。