《【2017年整理】计量经济学经典eviews时间序列回归》由会员分享，可在线阅读，更多相关《【2017年整理】计量经济学经典eviews时间序列回归（5页珍藏版）》请在金锄头文库上搜索。
1、22时间序列回归本章讨论含有 ARMA 项的单方程回归方法，这种方法对于分析时间序列数据（检验序列相关性，估计 ARMA 模型，使用分布多重滞后，非平稳时间序列的单位根检验）是很重要的。13.1 序列相关理论时间序列回归中的一个普遍现象是：残差和它自己的滞后值有关。这种相关性违背了回归理论的标准假设：干扰项互不相关。与序列相关相联系的主要问题有：一、一阶自回归模型最简单且最常用的序列相关模型是一阶自回归 AR(1)模型定义如下： tttuxy1参数 是一阶序列相关系数，实际上，AR(1)模型是将以前观测值的残差包含到现观测值的回归模型中。二、高阶自回归模型：更为一般，带有 p 阶自回归的回归，
2、AR(p)误差由下式给出：tttuxytptu21AR(p)的自回归将渐渐衰减至零，同时高于 p 阶的偏自相关也是零。13.2 检验序列相关在使用估计方程进行统计推断（如假设检验和预测）之前，一般应检验残差（序列相关的证据） ，Eviews 提供了几种方法来检验当前序列相关。1Dubin-Waston 统计量 D-W 统计量用于检验一阶序列相关。2相关图和 Q-统计量 计算相关图和 Q-统计量的细节见第七章3序列相关 LM 检验 检验的原假设是：至给定阶数，残差不具有序列相关。13.3 估计含 AR 项的模型随机误差项存在序列相关说明模型定义存在严重问题。特别的，应注意使用 OLS 得出的过分
3、限制的定义。有时，在回归方程中添加不应被排除的变量会消除序列相关。1一阶序列相关在 EViews 中估计一 AR(1)模型，选择 Quick/Estimate Equation 打开一个方程，用列表法输入方程后，最后将 AR(1)项加到列表中。例如：估计一个带有 AR(1)误差的简单消费函数ttt uGDPcCS21tu应定义方程为： cs c gdp ar(1)2高阶序列相关估计高阶 AR 模型稍稍复杂些，为估计 AR(k)，应输入模型的定义和所包括的各阶 AR 值。如果想估计一个有 1-5 阶自回归的模型ttt uGDPcCS21ttu5应输入： cs c gdp ar(1) ar(2)
4、ar(3) ar(4) ar(5)3存在序列相关的非线性模型EViews 可以估计带有 AR 误差项的非线性回归模型。例如：估计如下的带有附加 AR(2)误差的非线性方程tctt uGDPCS2123tttt ucu2413使用 EViews 表达式定义模型，在后面的方括号内描述 AR 修正项，对每一阶 AR 滞后项都应包括一个系数，每项之间用逗号隔开。cs=c(1)+gdp c(2)+ar(1)=c(3),ar(2)=c(4)EViews 通过 差分来转换这种非线性模型且使用 Gauss-Newton 迭代法来估计转换后的非线性模型。4存在序列相关的两阶段回归模型通过把二阶段最小二乘法或二阶
5、段非线性最小二乘法和 AR 项结合起来，对于在回归因子和扰动项存在相关性的情况和残差存在序列相关一样估计模型。5AR 估计输出含有 AR 项的模型有两种残差：第一种是无条件残差 ，bxyuttt通过原始变量以及估计参数 算出。在用同期信息对 yt 值进行预测时，这些残差是可以观测出的误差，但要忽略滞后残差中包含的信息。通常，除非有特别的原因来检验这些残差，Eviews 不能自动计算下面的估计。第二种残差是估计的一期向前预测误差 。如名所示，这种残差代表预测误差。一般 AR(p)平稳条件是：滞后算子多项式的根的倒数在单位圆内。 EViews 在回归输出的底部给出这些根：Inverted AR R
6、oots。如果存在虚根，根的模应该小于 1。6EViews 如何估计 AR 模型EViews 估计 AR 模型采用非线性回归方法。这种方法的优点在于：易被理解，应用广泛，易被扩展为非线性定义的模型。注意：非线性最小二乘估计渐进等于极大似然估计且渐进有效。13.4 ARIMA 理论ARIMA（自回归单整动平均）模型是 AR 模型的一般化，EViews 使用三种工具来为干扰项的序列相关建模：自回归 AR、单整 I、动平均 MA。13.5 估计 ARIMA 模型为建立 ARIMA 模型，需要： 差分因变量，确定差分阶数； 描述结构回归模型（因变量和回归因子） ，加入 AR 或 MA 项。一、ARMA
7、 项 模型中 AR 和 MA 部分应使用关键词 ar 和 ma 定义。二、季节 ARMA项 对于带有季节移动的季度数据，Box and Jenkins(1976)建议使用季节自回归 SAR和季节动平均 SMA。三、ARIMA 估计输出 存在 AR 或 MA 定义的估计输出和 OLS 是一样的，只是增加了一个 AR，MA多项式的倒根的下部程序块。四、ARMA 估计选择 带有 AR 或 MA 的模型用非线性最小二乘法估计。非线性估计方法对所有系数估计都要求初值。作为缺省 Eviews 决定初值。用户可设置初值，EViews 使用 C 系数向量。也可使用命令安排 C 向量值定义，例如下面方程的系数Y
8、 c X ma(2) ma(1) sma(4) ar(1)可定义为 param c(1) 50 c(2 ) 0.8 c(3) 0.2 c(4) 0.6 c(5) 0.1 c(6) 0.5初值：常数是 50， X 系数的初值是 0.8， ar(1)、ma(2)、ma(1)、sma(4) 系数的初值分别是 0.2 ， 0.6，0.1，0.5。13.6 诊断检验如果 ARMA 模型定义正确，模型残差将为白噪声。这意味着残差中应不存在序列相关。D-W 统计量24是当方程右边没有滞后变量时对一阶序列相关的检验。如上所述，对残差中序列相关更多的检验可以如：View/Residual Tests/Corre
9、logram-Q-Statistic 和 View/Residual Tests/Serial correlation LM Test。13.7 多项分布滞后（PDLs）一个分布滞后算子如下（13.37）tkttttt xxy 10系数 描述 x 对 y 作用的滞后。在模型中解释变量与随机误差项不相关的情况下，可以直接使用OLS 估计参数。在其它情形下，x 的当前和滞后值具有高共线性时，直接估计失败。可以使用多项式分布滞后（PDLS）来减少要估计的参数个数，以此来平滑滞后系数。平滑就是要求系数服从一个相对低阶的多项式。P 阶 PDLS 模型限制 系数服从如下形式的 p 阶多项式j = 0 ,
10、1 , 2 , , k （13.38）ppj cjcjj )()()( 12321 是事先定义常数：c是 偶 数是 奇 数（ kc/)(PDLS 有时被称为 Almon 分布滞后模型。常数 仅用来避免共线性引起的数值问题，不影响 的估c 计。这种定义允许仅使用参数 p 来估计一个 x 的 k 阶滞后的模型（如果 p k，将显示“近似奇异“错误信息） 。如果定义一个 PDL 模型，EViews 用(13.38) 式代入到(13.37)式，将产生如下形式方程(13.40)tptt zzy 121其中(13.41)ktptptpp ktttk xcxcxczz )()1()(1121 一旦从（13.
11、40）式估计 ，利用（13.38）式就可得到 的各系数。这一过程很明了，因为 是 的线性变换。定义一个 PDLs 有三个元素：滞后长度 k，多项式阶数（多项式最高次幂数）p 和附加的约束。13.8 非平稳时间序列上述 ARMA 估计理论都是基于平稳时间序列。如果一个序列的均值和自协方差不依赖于时间，就说它是平稳的。非平稳序列的典型例子是随机游动 ， 是平稳随机扰动项。序列 y 有一个ttty1t常数预测值，方差随时间增长。随机游动是差分平稳序列，因为 y 一阶差分后平稳。，差分平稳序列称为单整，记为 I(d)，d 为单整阶数。单整阶数是序列中单位tttt yLy)1(根数，或者是使序列平稳而差
12、分的阶数。对于上面的随机游动，有一个单位根，所以是 I(1)，同样，平稳序列是 I(0)。13.9 单位根检验EViews 提供两种单位根检验：Dickey-Fuller(DF)、增广 DF(ADF)检验和 Phillips-Perron（PP）检验。一、ADF 检验为说明 ADF 检验的使用，先考虑一个 AR(1)过程（13.46）ttty1是参数， 假设为白噪声。如果-1 1，y 平稳序列。如果 =1，y 是非平稳序列（带漂移的,t25随机游动） 。如果这一过程在一些点开始，y 的方差随时间增长趋于无穷。如果 的绝对值大于 1，序列发散。因此，一个序列是否平稳，可以检验 是否严格小于 1。
13、DF 和 PP 都用单位根作为原假设。因为发散序列没有经济学含义，所以备选假设为单边假设 。1:0H 1:H从方程两边同时减去 1tytty1其中 （13.47）所以原假设和备选假设可改为 （13.48）0:1H单位根检验可以看作对 进行 t 检验。EViews 将 DF，ADF 检验都看成为 ADF 检验。ADF 检验考虑如下三种回归形式： tpiititt yy1tiittttpiittt yay1210即通过在模型中增加 的滞后项，以消除残差的序列相关性。在检验回归中包括常数，常数和线性t趋势，或二者都不包含。二、Phillips-Perron(PP)检验Phillips 和 Perro
14、n（1988）提出一种非参数方法来控制序列中高阶序列相关。对 AR(1)的 PP 检验为：（13.51）ttty1ADF 检验通过在方程右边添加滞后差分项来修正高阶序列相关。PP 检验 参数的 t 统计量来修正 AR(1)的序列相关。这种修正方法是非参数的，因为我们使用 在零频率的谱估计。零频率对未知形式的异方 差性和自相关性较稳健。EViews 使用 Newey-West 异方差自相关一致估计（13.52）qj j102)(（13.53） Tjtjtj1q 是截断滞后值。PP 统计量由下式计算：（13.54）stt bbp2)(0210是 t 统计量； 是 的标准差； s 是检验回归标准差。
15、PP 统计量渐进分布同 ADF 的 t 统计量一样。bbsEViews 显示 Mackinnon 临界值。对 PP 检验，必须为 Newey-West 纠正定义截断滞后因子 q，即要包括的序列相关期数。对话框开始包括 N-W 自动截断滞后选择（ floor 函数返回的是不超过括号中数的最大整数）)10/(492Tflorq这仅基于检验回归中使用的观测值数，也可定义为任何整数。13.10 命 令命令 equation eq_gdp.ls gdp c ar(1) ar (2) ma(1) ma(2)用来用一个 arma(2,2)模型拟和序列 GDP 并把结果储存26在方程 EQ_GDP 中。命令 eq1.auto(4) 用来检验方程 EQ!残差序列直到四阶的相关系数。命令 eq1.correlogram(12)用来显示方程直到 12 阶的残差相关图。命令 equation eq2.ls gdp c pdl(m1,12,3) 使
