复旦大学-经济学院-谢识予-计量经济学第七章-误差序列相关

1、1 第七章第七章 误差序列相关误差序列相关 2 第一节 误差序列相关的性质和原因 第二节 误差序列相关的发现和判断 第三节 误差序列相关的克服和处理 本章结构 3 第一节 误差序列相关的性质和原因 n两变量和多元线性回归模型都要求模型 的误差项不存在序列相关性，即： 对任意 都成立。 n这条假设的含义是误差项是纯粹的微小 外来扰动因素，不同期之间相互独立， 不包含任何有规律性、趋势性的因素。 0)()()( jijjii EEEE ji 4 n这条假设对线性回归分析也十分重要， 最小二乘估计的最小方差性和一致估计， 得到残差方差的无偏估计，以及进行各 种统计推断等，也都以这条假设为基础。 n但

2、误差项无序列相关的假设也不是总能 成立。由于误差项包含的因素常常有时 间趋势，数据处理也会导致不同期数据 产生内在联系，因此误差序列往往是有 自相关性的。 5 n这种问题称为线性回归模型的“误差序 列相关”，表现为： 对至少部分 成立。 n当线性回归模型存在误差序列相关性时， 参数估计的有效性和一致性都不成立或 无法证明，残差方差和参数估计量方差 的估计无法得到，从而各种统计推断和 预测分析也同样会遇到困难。 0)()()( jijjii EEEE ji 6 n误差序列相关可以有多种不同的情况， 其中相邻两期误差项之间的相关性，也 就是误差项 受前一期误差项 的影响， 称为误差项的“一阶自回归

3、”。 n一阶自回归可以表示为 ， 其中 满足 ，称“一阶自回归系 数”， 是均值为0的独立同分布随机变 量。 i 1i iii 1 10 i 7 第二节 误差序列相关的发现和判断 n处理和克服误差序列相关性的基础是判 断该问题的存在和类型。 n对于一阶自回归性，就是要判断一阶自 回归系数 的大小和符号的正负性。 n回归残差序列分析也是发现和检验误差 序列相关性的基本方法，残差序列分析 包括残差序列图分析和杜宾瓦森检验 。 8 一、残差序列图分析 n以i为横轴，以残差e或e/S为纵轴，画出 残差序列的分布图。 n不存在误差序列相关问题，同时也不存 在系统偏差和异方差性的模型，回归残 差序列的分布

4、形态应该如图7.1（a）那 样，无规律而均匀地分布在横轴上下的 一定区域内。 9 S e i S e i 10 n如果模型存在误差序列相关问题，那么残差序 列的分布会呈现相应的规律性。 n例如当误差项有一阶正自相关问题时，残差序 列的分布形态会出现类似图7.1（b）的情况， 因为相邻残差之间的正相关关系会使它们出现 相同趋势和符号的机会增大。 S e i 11 n如果误差项有一阶负自相关性，那么相 邻误差项异号的机会增大，因此残差分 布会出现类似图7.1（c）的形态。 n还可以用另一种残差序列图发现误差序 列相关问题。以 为横轴，以 为纵轴， 用相邻残差项构成坐标 ，然后观察这些 坐标的分布情

5、况。 1i e i e 12 图7.2 分析误差序列相关残差分布图 000 c a b i e i e 1i e 1i e 1i e 13 n如果这些坐标如图7.2（a）那样均匀地 分布在四个象限内，应认为不存在误差 序列相关问题。 n如果坐标分布如图7.2（b）和（c），那 么应分别判断有一阶正自相关性和一阶 负自相关性，因为（b）图坐标分布落在 一、三象限多意味着相邻残差同号的居 多，而（c）图坐标分布落在二、四象限 多则意味着相邻残差异号居多。 14 二、杜宾瓦森检验（DW检验） nDW检验也是一种残差序列分析，因为用于检 验的DW统计量是根据回归残差序列计算的。 nDW检验的方法是，首

6、先假设线性回归模型 有一阶自回归问题，即 。 n然后检验一阶自回归系数 的显著性。如果检 验结果是显著的，那么认为误差项有一阶自回 归性，否则认为误差项没有一阶自回归性。 KK XXY 110 iii 1 15 n要检验一阶自回归系数是否有显著性， 首先必须对它的值进行估计。 n为此我们考察相邻误差项之间的协方差 公式。根据 和 的性质有 111 222 111 () ()()()() iiiii iiiii EE EEEE i i 16 n因此： n由于模型误差项的数值无法得到，因此 的真实值是无法得到的。 n但可以根据误差项与回归残差的关系， 用残差序列构造下列统计量： )( )( 2 1

7、 i ii E E i i n i ii e ee 2 1 1 17 n作为误差序列一阶自回归系数的估计。 n更进一步，杜宾和瓦森考虑用已知分布 且与 有密切关系的DW统计量 来代替 。 i i n i ii e ee DW 2 2 2 1 18 n两者之间的关系是： n因此 =0（无一阶自回归性）对应 DW=2， 1（误差项有强正自相关） 对应DW0， -1（误差项有强负自相 关）对应DW4。 1222 2 2 2 1 2 2 1 2 2 1 2 2 i i n i ii i i n i ii n i i n i i e ee e eeee DW 19 nDW的精确分布实际上也不清楚，而且分

8、 布情况与解释变量的取值有关。但杜宾 和瓦森证明对于解释变量的任意情况， DW统计量有一个上限和一个下限，在一 定条件下它们服从 分布。 n杜宾和瓦森计算了对应显著性水平 和0.01，样本容量 ，解释变量 个数 时，判断误差序列存在一阶正 自相关性的上下限 和 的临界值表， 作为检验的基本工具。 05. 0 10015 n 5K L d U d 20 n由于正自相关性和负自相关性之间的对 称性，实际上只需要有检验正自相关性 的方法即可，检验负自相关性时可先作 变换 ，然后用检验正自相关 性的临界值进行判断。 n检验误差序列正自相关的前提是DW ，接受 ，认为误差项不存 在一阶自相关性；若DW

9、认为存在 一阶正自相关；若 DW2，则根据检验负自相关性。 n上述临界值和判断方法也可以总结为如 图7.3中的5个判断区域。 DW 0 24 U d4 L d4 U d L d 正自相关 无法判断 无自相关 无法判断 负自相关 23 DW检验的缺点 n首先是它只适用于一阶自回归性的检验， 而且在样本数较小或解释变量数较大时 不适用； n其次是当解释变量有随机性（分布滞后 模型或联立方程组模型中普遍）时不适 用； n此外是DW检验存在无法判断的DW值区 间。 24 n不过，当样本容量n较大而解释变量数K 较小时，无法判断的区间也比较小，因 此可以通过增加样本容量，换新的样本 或修改模型等克服，也

10、可以当作存在误 差序列相关加以处理克服无法判断的困 难。 n例71。详见Eviews演示。 25 第三节 误差序列相关的克服和处理 一、一阶差分法 二、广义差分法 三、柯奥迭代法 四、杜宾两步法 26 一、一阶差分法 n设线性回归模型为 ，而 且已经经过检验知 有很强的一阶正自 相关性，即 。其中 是均值为0的独立同分布、无自相关性 的随机变量。 n把滞后一期的观测值代入变量关系，得 方程 iii XY 10 i iii 1 i 1 11101 iii XY 27 n（7-1）式减去（7-2）式得到： n由于 ，因此有 n令 为Y的一阶差分， 为X 的一阶差分，则有 1111 iiiiii X

11、XYY iii 1 11111 iiiiiii XXYY 1 iii YYY 1 iii XXX 111 iiiii XY 28 n因为 ，因此上式近似为 n用该Y和X的一阶差分模型进行回归分析， 可以避免原模型误差序列相关性的困难， 得到参数 估计值 。 n原模型的另一个参数 的估计值，则可 以根据 计算。 n这种克服误差序列相关性的方法称“一 阶差分法”。 1 iii XY 1 1 1 b 0 XbYb 10 29 n变量的增长率与一阶差分有密切关系，用变量 的增长率进行回归也能起到消除误差项强正自 相关性的作用。 n一阶差分模型在克服误差序列相关性方面的作 用有较大的局限，只适用 接近1

12、的一阶正自 相关性，而且如果模型没有误差序列相关性、 有负自相关性或只轻微正自相关性，运用一阶 差分模型反而会导致更强的误差序列相关性。 30 二、广义差分法 n设模型为 ，且已知模型 的误差项有一阶自相关性，即 其中 n同样把滞后一期的观测值代入模型，可 得 （7-6） n把（7-5）式减去（7-6）式与 的乘积， 可得 iii XY 10 iii 1 10 11101 iii XY 1111001 iiiiii XXYY 31 n我们称 为Y 的“广义差分”， 为X 的“广义差分”。 n再根据 ，有 n如果记 ，则有 n由于 不存在误差序列相关问题，因此 用广义差分模型进行回归分析，可以避

13、 免原模型的误差序列相关问题。 1 * iii YYY 1 * iii XXX 1ii * 10 * 1 ii XY 1 0 A iii XAY * 1 * i 32 n可以利用模型直接得到 的参数估计 值 ，并得到A的估计值 。 n再利用该估计值得到原模型参数 的估 计值 。 n这种克服误差序列相关问题的参数估计， 称“广义差分法”。 1 1 bA 0 1 0 Ab 33 广义差分法的缺点 n首先是差分变换会减少一个样本容量。 这通常可以通过对第一组数据作变换后 保留的方法解决。 n另一个问题是假设已知的一阶自回归系 数 实际上是无法知道的，只能根据原 模型的回归残差序列求其估计值。既然 原

14、模型存在误差序列相关，那么回归残 差就会受到影响，从而估计值也有偏差。 34 三、柯奥迭代法 n“柯-奥迭代法”是克服广义差分法无法得 到一阶自回归系数 较准确估计困难的 方法。 n首先对原模型运用普通最小二乘估计， 估计参数并计算出回归残差序列，再根 据回归残差序列计算 的第一个估计值 i i n i ii e ee 2 1 1 35 n这个 的准确程度可能有问题，但它至 少可以反映 的大部分情况，用这个 作为估计值进行广义差分法处理，可以 消除掉模型的大部分误差序列相关性。 n作广义差分变换： n再用 和 进行线性回归分析得到估计 值 和 ，并计算相应的回归残差序列。 1 * iii YY

15、Y 1 * iii XXX * i Y * i X A 1 b 36 n用 和 的回归残差进行DW检验，如果 不存在误差序列相关性问题，说明广义 差分已经消除了原模型误差序列相关的 影响，把 和 作为原模型两 个参数的估计值。 n如果仍然有误差序列相关性，则可以用 新的回归残差序列重新计算 的估计值， 再进行广义差分变换和回归，计算相应 的回归残差序列 * i Y * i X 1 0 Ab 1 b 37 n这样的过程可以反复进行，直到检验结 果不存在误差序列相关性。 n通常迭代一至两次就能实现目标，或 就会收敛。我们把最后得到的一组估计 量，作为原模型两个参数的估计。 38 四、杜宾两步法 n我们直接从两变量模型的广义差分式 出发。如果我们对该式稍作调整，可得 这是一个带滞后变量的多元线性回归模 型，其中 为均值为0的独立同分布随机 变量。 1111001 iiiiii XXYY iiiii XXYY 11