《《计量经济分析方法与建模》第二版课件-第06章条件异方差模型讲义资料》由会员分享，可在线阅读，更多相关《《计量经济分析方法与建模》第二版课件-第06章条件异方差模型讲义资料（99页珍藏版）》请在金锄头文库上搜索。
1、1,第六章 条件异方差模型,EViews中的大多数统计工具都是用来建立随机变量的条件均值模型。本章讨论的重要工具具有与以往不同的目的建立变量的条件方差或变量波动性模型。 我们想要建模并预测其变动性通常有如下几个原因: 首先，我们可能要分析持有某项资产的风险；其次，预测置信区间可能是时变性的，所以可以通过建立残差方差模型得到更精确的区间；第三，如果误差的异方差是能适当控制的，我们就能得到更有效的估计。,2,6.1 自回归条件异方差模型 自回归条件异方差(Autoregressive Conditional Heteroscedasticity Model， ARCH)模型是特别用来建立条件方差模
2、型并对其进行预测的。 ARCH模型是1982年由恩格尔(Engle, R.)提出，并由博勒斯莱文(Bollerslev, T., 1986)发展成为GARCH (Generalized ARCH)广义自回归条件异方差。这些模型被广泛的应用于经济学的各个领域。尤其在金融时间序列分析中。 按照通常的想法，自相关的问题是时间序列数据所特有，而异方差性是横截面数据的特点。但在时间序列数据中，会不会出现异方差呢？会是怎样出现的？,3,恩格尔和克拉格（Kraft, D., 1983）在分析宏观数据时，发现这样一些现象：时间序列模型中的扰动方差稳定性比通常假设的要差。恩格尔的结论说明在分析通货膨胀模型时，大
3、的及小的预测误差会大量出现，表明存在一种异方差，其中预测误差的方差取决于后续扰动项的大小。,5,6.1.1 ARCH模型 为了说得更具体，让我们回到k -变量回归模型： (6.1.1) 如果 ut 的均值为零，对 yt 取基于(t-1)时刻的信息的期望，即Et-1(yt)，有如下的关系： (6.1.2) 由于 yt 的均值近似等于式（6.1.1）的估计值，所以式（6.1.1）也称为均值方程。,6,假设在时刻 ( t 1 ) 所有信息已知的条件下，扰动项 ut 的条件分布是： (6.1.7) 也就是，ut 遵循以0为均值，(0+1u2t-1 )为方差的正态分布。,7,由于(6.1.7)中 ut
4、的方差依赖于前期的平方扰动项，我们称它为ARCH(1)过程： 通常用极大似然估计得到参数0, 1, 2, , k, 0, 1的有效估计。 容易加以推广，ARCH (p)过程可以写为： (6.1.8) 这时方差方程中的(p+1)个参数0, 1, 2, , p也要和回归模型中的参数0, 1, 2, , k一样，利用极大似然估计法进行估计。,8,如果扰动项方差中没有自相关，就会有 H0 ： 这时 从而得到扰动项方差的同方差性情形。 恩格尔曾表明，容易通过以下的回归去检验上述虚拟假设： 其中，t 表示从原始回归模型（6.1.1）估计得到的OLS残差。,9,在 ARCH(p) 过程中，由于 ut 是随机
5、的，ut2 不可能为负，所以对于 ut 的所有实现值，只有是正的，才是合理的。为使 ut2 协方差平稳，所以进一步要求相应的特征方程 （6.1.9） 的根全部位于单位圆外。如果 i（i = 1, 2, , p）都非负，式（6.1.9）等价于 1 + 2 + + p 1。,10,6.1.2 ARCH的检验,下面介绍检验一个模型的残差是否含有ARCH效应的两种方法：ARCH LM检验和残差平方相关图检验。 1. ARCH LM检验 Engle在1982年提出检验残差序列中是否存在ARCH效应的拉格朗日乘数检验（Lagrange multiplier test），即ARCH LM检验。自回归条件异方
6、差性的这个特殊的设定，是由于人们发现在许多金融时间序列中，残差的大小与最近的残差值有关。ARCH本身不能使标准的OLS估计无效，但是，忽略ARCH影响可能导致有效性降低。,11,ARCH LM检验统计量由一个辅助检验回归计算。为检验原假设：残差中直到q阶都没有ARCH，运行如下回归： 式中 t 是残差。这是一个对常数和直到 q 阶的滞后平方残差所作的回归。这个检验回归有两个统计量： （1）F 统计量是对所有残差平方的滞后的联合显著性所作的一个省略变量检验； （2）TR2 统计量是Engles LM检验统计量，它是观测值个数 T 乘以回归检验的 R2 ；,12,普通回归方程的ARCH检验都是在残
7、差检验下拉列表中进行的，需要注意的是，只有使用最小二乘法、二阶段最小二乘法和非线性最小二乘法估计的方程才有此项检验。,Breusch-Pagan-Godfrey Harvey Glejser ARCH White Custom Test Wizard,图6.4 普通方程的ARCH检验列表,13,2. 残差平方相关图 显示直到所定义的滞后阶数的残差平方t2的自相关性和偏自相关性，计算出相应滞后阶数的Ljung-Box统计量。残差平方相关图可以用来检查残差自回归条件异方差性（ARCH）。如果残差中不存在ARCH，在各阶滞后自相关和偏自相关应为0，且Q统计量应不显著。可适用于使用LS，TSLS，非线
8、性LS估计方程。在图6.4中选择Residuals Tests/ Correlogram Squared Residuals项，它是对方程进行残差平方相关图的检验。单击该命令，会弹出一个输入计算自相关和偏自相关系数的滞后阶数设定的对话框，默认的设定为36，单击OK按钮，得到检验结果。,14,例6.1 沪市股票价格指数波动的ARCH检验 为了检验股票价格指数的波动是否具有条件异方差性，本例选择了沪市股票的收盘价格指数的日数据作为样本序列，这是因为上海股票市场不仅开市早，市值高，对于各种冲击的反应较为敏感，因此，本例所分析的沪市股票价格波动具有一定代表性。在这个例子中，我们选择的样本序列sp是19
9、96年1月1日至2006年12月31日的上海证券交易所每日股票价格收盘指数，为了减少舍入误差，在估计时，对sp进行自然对数处理，即将序列ln(sp)作为因变量进行估计。,15,由于股票价格指数序列常常用一种特殊的单位根过程随机游动（Random Walk）模型描述，所以本例进行估计的基本形式为： (6.1.12) 首先利用最小二乘法，估计了一个普通的回归方程，结果如下： (6.1.13) (2.35) (951) R2= 0.997,16,可以看出，这个方程的统计量很显著，而且，拟合 的程度也很好。但是需要检验这个方程的误差项是否存在条件异方差性，。,17,图6.1 股票价格指数方程回归残差,
10、观察上图，该回归方程的残差，我们可以注意到波动的“成群”现象：波动在一些较长的时间内非常小，在其他一些较长的时间内非常大，这说明残差序列存在高阶ARCH效应。,18,因此，对式(6.1.26)进行条件异方差的ARCH LM检验，得到了在滞后阶数p = 3时的ARCH LM检验结果如下。此处的P值为0，拒绝原假设，说明式（6.1.26）的残差序列存在ARCH效应。,可以计算式（6.1.26）的残差平方t2的自相关（AC）和偏自相关（PAC）系数，结果说明式（6.1.26）的残差序列存在ARCH效应。,19,例6.2 中国CPI模型的ARCH检验 本例建立CPI模型，因变量为中国的消费价格指数（上
11、年同月=100）减去100，记为cpit；解释变量选择货币政策变量：狭义货币供应量M1的增长率，记为m1rt；3年期贷款利率，记为Rt，样本期间是1994年1月2007年12月。由于是月度数据，利用X-12季节调整方法对 cpit 和 m1rt 进行了调整，结果如下： t = (19.5) (-5.17) (2.88) (-2.74) R2=0.99 对数似然值 = -167.79 AIC = 2.045 SC =2.12,20,这个方程的统计量很显著，拟合的程度也很好。但是观察该回归方程的残差图，也可以注意到波动的“成群”现象：波动在一些时期内较小，在其他一些时期内较大，这说明误差项可能具有
12、条件异方差性。,21,从自相关系数和偏自相关系数可以看出：残差序列存在着一阶ARCH效应。再进行条件异方差的ARCH LM检验，得到了在滞后阶数p = 1时的ARCH LM检验结果：,因此计算残差平方t2的自相关（AC）和偏自相关（PAC）系数，结果如下：,22,从自相关系数和偏自相关系数可以看出：残差序列存在着一阶ARCH效应。因此利用ARCH(1)模型重新估计模型（6.1.14），结果如下： 均值方程： z = (12.53) (-1.53) (4.72) (-3.85) 方差方程： z = (5.03) (3.214) R2=0.99 对数似然值 = -151.13 AIC = 1.87
13、 SC = 1.98 方差方程中的ARCH项的系数是统计显著的，并且对数似然值有所增加，同时AIC和SC值都变小了，这说明ARCH(1)模型能够更好的拟合数据。,23,再对这个方程进行条件异方差的ARCH LM检验，得到了残差序列在滞后阶数p=1时的统计结果： 此时的相伴概率为0.69，接受原假设，认为该残差序列不存在ARCH效应，说明利用ARCH(1)模型消除了式（6.1.14）的残差序列的条件异方差性。式（6.1.15）的残差平方相关图的检验结果为： 自相关系数和偏自相关系数近似为0。这个结果也说明了残差序列不再存在ARCH效应。,24,6.1.3 GARCH模型 扰动项 ut 的方差常常
14、依赖于很多时刻之前的变化量（特别是在金融领域，采用日数据或周数据的应用更是如此）。因此 必须估计很多参数，而这一点很难精确的做到。但是如果我们能够意识到方程(6.1.8)不过是 t2 的分布滞后模型， 我们就能够用一个或两个 t2 的滞后值代替许多 ut2的滞后值，这就是广义自回归条件异方差模型(generalized autoregressive conditional heterosce-dasticity model，简记为GARCH模型)。在GARCH模型中，要考虑两个不同的设定：一个是条件均值，另一个是条件方差。,25,在标准化的GARCH(1,1)模型中： 均值方程：(6.1.17
15、) 方差方程：(6.1.18) 其中：xt 是 (k+1)1维外生变量向量， 是(k+1)1维系数向量。 (6.1.17)中给出的均值方程是一个带有扰动项的外生变量函数。由于t2是以前面信息为基础的一期向前预测方差 ，所以它被称作条件方差，式（6.1.18）也被称作条件方差方程 。,26,(6.1.18)中给出的条件方差方程是下面三项的函数： 1常数项（均值）： 2用均值方程(6.1.11)的扰动项平方的滞后来度量从前期得到的波动性的信息： ut2-1（ARCH项）。 3上一期的预测方差： t2-1 （GARCH项）。 GARCH(1,1)模型中的(1,1)是指阶数为1的GARCH项（括号中的第一项）和阶数为1的ARCH项（括号中的第二项）。一个普通的ARCH模型是GARCH模型的一个特例，GARCH(0,1)，即在条件方差方程中不存在滞后预测方差t2-1的说明。,27,在EViews中ARCH模型是在扰动项是条件正态分布的假定下，通过极大似然函数方法估计的。例如，对于GARCH(1,1)，t 时期的对数似然函数为： (6.1.19) 其中 (6.1.20) 这个说明通常可以在金融领域得到解释，因为代理商或贸易商可以通过建立长期均值的加权平均（常数），上期的预期方差（GARCH项）和在以前各
