第七章时间序列分析(计量经济学,潘省初)

1、第七章第七章 时间序列分析时间序列分析(Time Series Analysis)第一节第一节 时间序列分析的基本概念时间序列分析的基本概念 经济分析通常假定所研究的经济理论中涉及的变量之间存在着长期均衡关系。按照这一假定，在估计这些长期关系时，计量经济分析假定所涉及的变量的均值和方差是常数，不随时间而变。 然而，经验研究表明，在大多数情况下，时间序列变量并不满足这一假设，从而产生所谓的“伪伪回归回归”问题(spurious regression problem)。 为解决这类问题，研究人员提出了不少对传统估计方法的改进建议，其中最重要的两项是对变量的非平稳性非平稳性 (non-station

2、arity) 的系统性检验和协整和协整(cointegration)。协整协整 协整分析被认为是上世纪八十年代中期以来计量经济学领域最具革命性的进展。 简单地说，协整分析涉及的是一组变量，它们各自都是不平稳的（含义是随时间的推移而上行或下行），但它们一起漂移。这种变量的共同漂移使得这些变量之间存在长期的线性关系，因而使人们能够研究经济变量间的长期均衡关系。如果这些长时间内的线性关系不成立，则对应的变量被称为是“非协整的” 。误差修正模型误差修正模型 一般说来，协整分析是用于非平稳变量组成的关系式中长期均衡参数估计的技术。它是用于动态模型的设定、估计和检验的一种新技术。 此外，协整分析亦可用于短

3、期或非均衡参数的估计，这是因为短期参数的估计可以通过协整方法使用长期参数估计值，采用的模型是误差修正模型误差修正模型 (error correction model)。 在介绍上述方法之前，下面先介绍所涉及的一些术语和定义。一一 平稳性（平稳性（Stationarity） 严格平稳性严格平稳性 (strict stationarity) 如果一个时间序列Xt的联合概率分布不随时间而变，即对于任何n和k，X1, X2, Xn的联合概率分布与X1+k, X2+k, Xn+k 的联合分布相同，则称该时间序列是严格平稳的。 由于在实践中上述联合概率分布很难确定，我们用随机变量Xt（t=1,2,）的均值

4、、方差和协方差代替之，即所谓的“弱平稳性”。2. 弱平稳性弱平稳性 (weak stationarity)一个时间序列是“弱平稳的”，如果： （1）均值 E(Xt) =，t=1,2, （7.1） （2 ）方差 Var(Xt) = E(Xt -)2 =2，t =1,2,（7.2） （3）协方差 Cov(Xt, Xt+k）= E (Xt -)(Xt+k -) rk， t=1,2,，k0 （7.3）3. 平稳性和非平稳性平稳性和非平稳性 通常情况下，我们所说的平稳性指的就是弱平稳性。一般来说，如果一个时间序列的均值和方差在任何时间保持恒定，并且两个时期t和t+k之间的协方差仅依赖于两时期之间的距离（

5、间隔或滞后）k，而与计算这些协方差的实际时期t无关，则该时间序列是平稳的。 只要这三个条件不全满足，则该时间序列是非平稳的。事实上，大多数经济时间序列是非平稳的。例如，在图7.1中，某国的私人消费（CP）和个人可支配收入（PDI）这两个时间序列都有一种向上的趋势，几乎可以断定它们不满足平稳性条件（7.1），因而是非平稳的。图7 .1 某 国私人消费和个人可支配收入，1 9 6 0 1 9 9 5 年度数据单位：百万美元（1 9 7 0 年不变价）10000020000030000040000050000060000019601965197019751980198519901995CPPDI二二

6、 几种有用的时间序列模型几种有用的时间序列模型1、白噪声（、白噪声（ White noise） 白噪声通常用t表示，是一个纯粹的随机过程，满足:（1）E(t) = 0 , 对所有t成立；（2）V ar(t) = 2，对所有t成立；（3）Cov (t, t+k) = 0，对所有t和k0成立。 白噪声可用符号表示为： tIID(0, 2) （7.4）注：这里IID为Independently Identically Distributed（独立同分布)的缩写。2、随机漫步（、随机漫步（Random walk） 随机漫步是一个简单随机过程，由下式确定： Xt = Xt1+t （7.5） 其中t为白噪

7、声。 Xt的均值： E(Xt）= E(Xt-1+t）= E(Xt1) + E(t) = E(Xt1) 这表明Xt的均值不随时间而变。 为求Xt的方差，对（7.5）式进行一系列置换： Xt = Xt1+t = Xt2+t-1+t = Xt3+t-2+t-1+t = = X0+1+2+t = X0+t 其中X0是Xt的初始值，可假定为任何常数或取初值为0，则 2110)()()(tVarXVarXVarttttttt 这表明Xt的方差随时间而增大，平稳性的第二个条件（7.2）不满足，因此，随机漫步时间序列是非平稳时间序列。可是，若将（7.5）式 Xt = Xt1+t写成一阶差分形式： Xt=t （

8、7.6) 这个一阶差分新变量Xt是平稳的，因为它就等于白燥声t，而后者是平稳时间序列。3、带漂移项的随机漫步、带漂移项的随机漫步 (Random walk with drift) Xt=+Xt1+t （7.7） 其中是一非0常数，t为白燥声。 之所以被称为“漂移项”，是因为（7.7）式的一阶差分为 Xt = XtXt-1 =+t 这表明时间序列Xt向上或向下漂移，取决于的符号是正还是负。显然，带漂移项的随机漫步时间序列也是非平稳时间序列。4、自回归过程、自回归过程 随机漫步过程（7.5）（ Xt = Xt1+t）是最简单的非平稳过程。它是 Xt=Xt1+t （7.8）的特例，（7.8）称为一阶

9、自回归过程 (AR(1)，该过程在11时是平稳的，其他情况下，则为非平稳过程。 更一般地，（7.8）式又是 Xt=1Xt1+2Xt2+qXt-q+t （7.9）的特例，（7.9）称为q阶自回归过程 (AR(q)。可以证明，如果特征方程 11L2L23L3qLq = 0 （7.10）的所有根的绝对值均大于1，则此过程（7.9）是平稳的，否则为非平稳过程。三三 单整的时间序列（单整的时间序列（Integrated series） 从（7.6）可知，随机漫步序列的一阶差分序列Xt = XtXt-1是平稳序列。在这种情况下，我们说原非平稳序列Xt是“一阶单整的”，表示为I(1)。 与 此 类 似 ，

10、若 非 平 稳 序 列 必 须 取 二 阶 差 分(2Xt=XtXt-1)才变为平稳序列，则原序列是“二阶单整的”，表示为I(2)。 一般地，若一个非平稳序列必须取d阶差分才变为平稳序列，则原序列是“d阶单整的”(Integrated of order d)，表示为I(d)。 由定义不难看出，I(0)表示的是平稳序列，意味着该序列无需差分即是平稳的。另一方面，如果一个序列不管差分多少次，也不能变为平稳序列，则称为“非单整的”。 第二节第二节 平稳性的检验平稳性的检验 平稳性检验的方法可分为两类：传统方法和现代方法。前者使用自相关函数(Autocorrelation function)，后者使用

11、单位根(Unit roots)。单位根方法是目前最常用的方法，因此本节中，我们仅介绍单位根方法。一一 单位根单位根 考察（7.8）式的一阶自回归过程，即 Xt=Xt1+t （7.11） 其中t为白噪声，此过程可写成 XtXt1=t 或（1L）Xt = t （7.12） 其中L为滞后运算符，其作用是取时间序列的滞后，如Xt 的一期滞后可表示为L(Xt），即 L(Xt）= Xt1 由上节所知，自回归过程Xt平稳的条件是其特征方程的所有根的绝对值大于1。由于这里特征方程为1L=0，该方程 仅有一个根L=1/ ，因而平稳性要求11。 因此，检验Xt的平稳性的原假设和备择假设为： H0：1 Ha：1 接

12、受原假设H0表明Xt是非平稳序列，而拒绝原假设（即接受备择假设Ha）则表明Xt是平稳序列。单位根检验方法的由来 在=1的情况下，即若原假设为真，则（7.11）就是随机漫步过程（7.5），从上节得知，它是非平稳的。因此，检验非平稳性就是检验=1是否成立，或者说，就是检验单位根是否存在。换句话说，单位根是表示非平稳性的另一方式。这样一来，就将对非平稳性的检验转化为对单位根的检验，这就是单位根检验方法的由来。（7.11）式 Xt=Xt1+t 两端各减去Xt-1，我们得到 XtXt1= Xt1Xt1+t即 Xt= Xt1+t （7.13） 其中是差分运算符，=1。 假设为正（绝大多数经济时间序列确实如

13、此），前面的假设 H0：1 Ha：1 可写成如下等价形式： H0：0 Ha：0 在=0的情况下，即若原假设为真，则相应的过程是非平稳的。 换句话说，非平稳性或单位根问题，可表示为=1或=0。从而我们可以将检验时间序列Xt的非平稳性的问题简化成在方程（7.11）的回归中，检验参数=1 是否成立或者在方程（7.13）的回归中，检验参数=0是否成立。这类检验可用t检验进行，检验统计量为： 或 （7.14）其中， 和 分别为参数估计值 和 的标准误差，即 这里的问题是，（7.14）式计算的t值不服从t分布，而是服从一个非标准的甚至是非对称的分布。因而不能使用t分布表，需要用另外的分布表。1tStSSS

14、()SSe( )SSe二二 Dickey-Fuller检验（检验（DF检验）检验） 迪奇（Dickey） 和福勒（Fuller）以蒙特卡罗模拟为基础，编制了（7.14）中t统计量的临界值表，表中所列已非传统的t统计值，他们称之为统计值。这些临界值如表7.1所示。后来该表由麦金农（Mackinnon）通过蒙特卡罗模拟法加以扩充。表7.1 Dickey-Fuller 统计量临界值表取更小值的概率样本容量0.010.0250.050.100.900.950.9750.99无常数项无时间项(统计量)25-2.66-2.26-1.95-1.600.921.331.712.1650-2.62-2.25-1

15、.95-1.610.911.311.662.08100-2.60-2.24-1.95-1.610.901.291.642.03250-2.58-2.23-1.95-1.620.891.291.632.01500-2.58-2.23-1.95-1.620.891.281.622.00-2.58-2.23-1.95-1.620.891.281.622.00有常数项无时间项(统计量)25-3.75-3.33-3.00-2.62-0.370.000.340.7250-3.58-3.22-2.93-2.60-0.40-0.030.290.66100-3.51-3.17-2.89-2.58-0.42-0.

16、050.260.63250-3.46-3.14-2.88-2.57-0.42-0.060.240.62500-3.44-3.13-2.87-2.57-0.43-0.070.240.61-3.43-3.12-2.86-2.57-0.44-0.070.230.60有常数项有时间项(统计量)25-4.38-3.95-3.60-3.24-1.14-0.80-0.50-0.1550-4.15-3.80-3.50-3.18-1.19-0.87-0.58-0.24100-4.04-3.73-3.45-3.15-1.22-0.90-0.62-0.28250-3.99-3.69-3.43-3.13-1.23-0

17、.92-0.64-0.31500-3.98-3.68-3.42-3.13-1.24-0.93-0.65-0.32-3.96-3.66-3.41-3.12-1.25-0.94-0.66-0.33 有了表，我们就可以进行DF检验了，DF检验按以下两步进行： 第一步：对（7.13）式执行OLS回归，即估计 Xt=Xt-1+t （7.15） 得到常规t值。 第二步：检验假设 H0：= 0 Ha：0 用上一步得到的t值与表7.1中查到的临界值比较，判别准则是： 若 t， 则接受原假设H0，即Xt非平稳。 若t，则拒绝原假设H0，Xt为平稳序列。 Dickey和Fuller注意到临界值依赖于回归方程的类型

18、。因此他们同时还编制了与另外两种类型方程中相对应的统计表，这两类方程是： Xt=+Xt-1+t （7.16）和 Xt=+t+Xt-1+t （7.17） 二者的临界值分别记为和T。这些临界值亦列在表7.1中。尽管三种方程的临界值有所不同，但有关时间序列平稳性的检验依赖的是Xt-1的系数，而与、无关。 在实践中，经济数据一般不用（7.15）式那样的无常数项的形式。带漂移项的时间序列通常采用（7.17）式，而不带漂移项的时间序列采用（7.16）式。 例7.1 检验某国私人消费时间序列的平稳性。表7 . 2 某 国 私 人 消 费 和 个 人 可 支 配 收 入 ( 1 9 7 0 百 万 美 元 )

19、年 份 私 人 消 费 个 人 可 支 配 收 入 消 费 价 格 指 数1 9 6 01 0 7 8 0 8 . 01 1 7 1 7 9 . 20 . 7 8 3 1 4 21 9 6 11 1 5 1 4 7 . 01 2 7 5 9 8 . 90 . 7 9 1 6 8 41 9 6 21 2 0 0 5 0 . 01 3 5 0 0 7 . 10 . 8 0 1 7 5 81 9 6 31 2 6 1 1 5 . 01 4 2 1 2 8 . 30 . 8 2 8 6 8 81 9 6 41 3 7 1 9 2 . 01 5 9 6 4 8 . 70 . 8 4 7 1 8 51 9

20、 6 51 4 7 7 0 7 . 01 7 2 7 5 5 . 90 . 8 8 5 8 2 81 9 6 61 5 7 6 8 7 . 01 8 2 3 6 5 . 50 . 9 1 6 5 0 51 9 6 71 6 7 5 2 8 . 01 9 5 6 1 1 . 00 . 9 3 4 2 3 21 9 6 81 7 9 0 2 5 . 02 0 4 4 7 0 . 40 . 9 4 1 1 9 31 9 6 91 9 0 0 8 9 . 02 2 2 6 3 7 . 50 . 9 6 9 6 3 01 9 7 02 0 6 8 1 3 . 02 4 6 8 1 9 . 01 . 0

21、0 0 0 0 01 9 7 12 1 7 2 1 2 . 02 6 9 2 4 8 . 91 . 0 3 3 7 2 71 9 7 22 3 2 3 1 2 . 02 9 7 2 6 6 . 01 . 0 6 8 0 6 41 9 7 32 5 0 0 5 7 . 03 3 5 5 2 1 . 71 . 2 2 8 1 5 61 9 7 42 5 1 6 5 0 . 03 1 0 2 3 1 . 11 . 5 1 7 7 9 51 9 7 52 6 6 8 8 4 . 03 2 7 5 2 1 . 31 . 7 0 1 1 4 71 9 7 62 8 1 0 6 6 . 03 5 0 4 2

22、 7 . 41 . 9 2 9 9 0 61 9 7 72 9 3 9 2 8 . 03 6 6 7 3 0 . 02 . 1 5 9 8 7 21 9 7 83 1 0 6 4 0 . 03 9 0 1 8 8 . 52 . 4 3 6 3 6 41 9 7 93 1 8 8 1 7 . 04 0 6 8 5 7 . 22 . 8 3 8 4 5 31 9 8 03 1 9 3 4 1 . 04 0 1 9 4 2 . 83 . 4 5 9 0 3 01 9 8 13 2 5 8 5 1 . 04 1 9 6 6 9 . 14 . 0 8 1 8 4 41 9 8 23 3 8 5 0 7

23、. 04 2 1 7 1 5 . 65 . 1 1 4 1 6 91 9 8 33 3 9 4 2 5 . 04 1 7 9 3 0 . 36 . 0 6 7 8 3 51 9 8 43 4 5 1 9 4 . 04 3 4 6 9 5 . 77 . 1 3 0 6 0 21 9 8 53 5 8 6 7 1 . 04 5 6 5 7 6 . 28 . 4 3 5 2 8 51 9 8 63 6 1 0 2 6 . 04 3 9 6 5 4 . 11 0 . 3 0 0 8 1 01 9 8 73 6 5 4 7 3 . 04 3 8 4 5 3 . 51 1 . 9 1 9 5 0 01 9

24、 8 83 7 8 4 8 8 . 04 7 6 3 4 4 . 71 3 . 6 1 4 4 8 01 9 8 93 9 4 9 4 2 . 04 9 2 3 3 4 . 41 5 . 5 9 2 8 5 01 9 9 04 0 3 1 9 4 . 04 9 5 9 3 9 . 21 8 . 5 9 5 3 9 01 9 9 14 1 2 4 5 8 . 05 1 3 1 7 3 . 02 2 . 0 9 1 1 6 01 9 9 24 2 0 0 2 8 . 05 0 2 5 2 0 . 12 5 . 4 0 1 2 2 01 9 9 34 2 0 5 8 5 . 05 2 3 0 6 6

25、 . 12 8 . 8 8 3 4 6 01 9 9 44 2 6 8 9 3 . 05 2 0 7 2 7 . 53 2 . 0 0 3 8 5 01 9 9 54 3 3 7 2 3 . 05 1 8 4 0 6 . 93 4 . 9 8 0 8 5 0 用表7.2中的私人消费（Ct）时间序列数据，估计与（7.16）和（7.17）相对应的方程，分别得到如下估计结果：(1) (1) =12330.48-0.01091 C =12330.48-0.01091 Ct-1 t-1 R R2 2=0.052=0.052 (t:) (5.138) (-1.339) DW=1.765 (t:) (5.1

26、38) (-1.339) DW=1.765(2) (2) =15630.83+346.4522t-0.04536C =15630.83+346.4522t-0.04536Ct-1t-1 R R2 2=0.057=0.057 (t:) (1.966) (0.436) (-0.5717) DW=1.716 (t:) (1.966) (0.436) (-0.5717) DW=1.716 两种情况下，t值分别为 -1.339和 -0.571，二者分别大于表7.1中从0.01到0.10的各种显著性水平下的值和值。因此，两种情况下都不能拒绝原假设，即私人消费时间序列有一个单位根，或换句话说，它是非平稳序列

27、。tCtC 下面看一下该序列的一阶差分（Ct）的平稳性。做类似于上面的回归，得到如下结果：(3) (3) 2 2 = 7972.671-0.85112= 7972.671-0.85112C Ct-1 t-1 R R2 2=0.425=0.425 (t:) (4.301) (-4.862) DW=1.967 (t:) (4.301) (-4.862) DW=1.967(4) (4) 2 2 =10524.35-114.461t-0.89738=10524.35-114.461t-0.89738C Ct-1 t-1 R R2 2=0.454=0.454 (t:) (3.908) (-1.294)

28、(-5.073) DW=1.988 (t:) (3.908) (-1.294) (-5.073) DW=1.988其中其中2 2C Ct t= =C Ct t- -C Ct-1t-1。tCtC 两种情况下，t值分别为 -4.862和-5.073，二者分别小于表7.1中从0.01到0.10的各种显著性水平下的值和T值。因此，都拒绝原假设，即私人消费一阶差分时间序列没有单位根，或者说该序列是平稳序列。 综合以上结果，我们的结论是： Ct是平稳序列，CtI(0）。 而Ct是非平稳序列，由于CtI(0），因而 CtI(1）。第三节第三节 协整协整 让我们考察弗里德曼的持久收入假设：私人总消费（让我们考

29、察弗里德曼的持久收入假设：私人总消费（Ct）是）是持久私人消费和暂时性私人消费（持久私人消费和暂时性私人消费（t）之和，持久私人消费与）之和，持久私人消费与持久个人可支配收入（持久个人可支配收入（Yt）成正比。则消费函数为：）成正比。则消费函数为： (7.18) 其中其中011。 用表用表7.2中数据对此消费函数进行中数据对此消费函数进行OLS估计，假定持久个人估计，假定持久个人收入等于个人可支配收入，我们得到：收入等于个人可支配收入，我们得到： = 0.80969Yt R2=0.9924 (t:) (75.5662) DW=0.8667tttPttYcC1tC 除除DW值低以外，估计结果很好

30、。值低以外，估计结果很好。t值很高表明回归值很高表明回归系数显著，系数显著，R2也很高，表明拟合很好。可是，由于方也很高，表明拟合很好。可是，由于方程中的两个时间序列是趋势时间序列或非平稳时间序程中的两个时间序列是趋势时间序列或非平稳时间序列，因此这一估计结果有可能形成误导。结果是，列，因此这一估计结果有可能形成误导。结果是，OLS估计量不是一致估计量，相应的常规推断程序不估计量不是一致估计量，相应的常规推断程序不正确。正确。 这种结果看上去非常好但涉及的变量是趋势时间序这种结果看上去非常好但涉及的变量是趋势时间序列的回归被列的回归被Granger 和和 Newbold 称为称为“伪回归伪回归

31、” (Spurious regression)。 当回归方程中涉及的时间序列是非平稳时间序列时，当回归方程中涉及的时间序列是非平稳时间序列时，OLS估计量不再是一致估计量，相应的常规推断程序估计量不再是一致估计量，相应的常规推断程序会产生误导。这就是所谓的会产生误导。这就是所谓的“伪回归伪回归”问题。问题。 他们指出，如果在时间序列的回归中他们指出，如果在时间序列的回归中DW值低于值低于R2，则应怀疑有伪回归的可能。我们上面的结果正，则应怀疑有伪回归的可能。我们上面的结果正是如此（是如此（R2 = 0.9924 DW = 0.8667）。）。 考虑到经济学中大多数时间序列是非平稳序列，考虑到经

32、济学中大多数时间序列是非平稳序列，则我们得到伪回归结果是常见的事。避免非平稳性则我们得到伪回归结果是常见的事。避免非平稳性问题的常用方法是在回归中使用时间序列的一阶差问题的常用方法是在回归中使用时间序列的一阶差分。可是，使用变量为差分形式的关系式更适合描分。可是，使用变量为差分形式的关系式更适合描述所研究的经济现象的短期状态或非均衡状态，而述所研究的经济现象的短期状态或非均衡状态，而不是其长期或均衡状态，描述所研究经济现象的长不是其长期或均衡状态，描述所研究经济现象的长期或均衡状态应采用变量本身。期或均衡状态应采用变量本身。 由上面的讨论，自然引出了一个明显的问题：我们由上面的讨论，自然引出了

33、一个明显的问题：我们使用非均衡时间序列时是否必定会造成伪回归？使用非均衡时间序列时是否必定会造成伪回归？ 对此问题的回答是，如果在一个回归中涉及的趋势对此问题的回答是，如果在一个回归中涉及的趋势时间序列时间序列“一起漂移一起漂移”，或者说，或者说“同步同步”，则可能没，则可能没有伪回归的问题，因而取决于有伪回归的问题，因而取决于t检验和检验和F检验的推断也检验的推断也没有问题。这种非均衡时间序列的没有问题。这种非均衡时间序列的“同步同步”，引出了，引出了我们下面要介绍的我们下面要介绍的“协整协整”概念。概念。一协整的概念一协整的概念 在方程（在方程（7.18） 中，持久收入假设中，持久收入假设

34、要求两时间序列要求两时间序列Ct和和Yt的线性组合，即时间序列的线性组合，即时间序列Ct1Yt必须是平稳的，这是因为此序列等于必须是平稳的，这是因为此序列等于t，而暂时，而暂时性私人消费（性私人消费（t）按定义是平稳时间序列。）按定义是平稳时间序列。 可是，可是，Ct和和Yt都是非平稳时间序列，事实上，不难都是非平稳时间序列，事实上，不难验证：验证：CtI(1），），YtI(1）。）。 也就是说，尽管也就是说，尽管CtI(1），），YtI(1），但持久收），但持久收入假设要求它们的线性组合入假设要求它们的线性组合t=Ct1Yt是平稳的，即是平稳的，即t=Ct1YtI (0）。在这种情况下，我们

35、说时间序）。在这种情况下，我们说时间序列列Ct和和Yt是协整的是协整的(Cointegrated)。下面给出协整。下面给出协整(Cointegration)的正式定义。的正式定义。1tttCY协整的定义协整的定义 如果两时间序列如果两时间序列YtI(d)，XtI(d)，并且这，并且这两个时间序列的线性组合两个时间序列的线性组合a1Yt+a2Xt 是是(d-b)阶单阶单整的，即整的，即a1Yt+a2XtI(d-b)（db0），则），则Yt 和和Xt被称为是（被称为是（d, b）阶协整的。记为）阶协整的。记为 Yt, XtCI(d , b)这里这里CI是协整的符号。构成两变量线性组合的是协整的符号

36、。构成两变量线性组合的系数向量（系数向量（a1，a2）称为）称为“协整向量协整向量”。 下面给出本节中要研究的两个特例。下面给出本节中要研究的两个特例。 1、Yt, XtCI(d, d） 在这种情况下，在这种情况下，d=b，使得，使得a1Yt+a2XtI(0），即两时），即两时间序列的线性组合是平稳的，因而间序列的线性组合是平稳的，因而 Yt, XtCI(d, d）。）。 2、Yt, XtCI(1, 1) 在这种情况下，在这种情况下，d=b=1，同样有，同样有a1Yt+a2XtI(0)，即两，即两时间序列的线性组合是平稳的，因而时间序列的线性组合是平稳的，因而 Yt, XtCI(1, 1）。）

37、。 让我们考虑下面的关系让我们考虑下面的关系 Yt = 0+1Xt （7.19） 其中，其中，YtI(1），），XtI(1）。）。 当当0= Yt01Xt时，该关系处于长期均衡状态。时，该关系处于长期均衡状态。 对长期均衡的偏离，称为对长期均衡的偏离，称为“均衡误差均衡误差”，记为，记为t： t = Yt01Xt 若长期均衡存在，则均衡误差应当围绕均衡值若长期均衡存在，则均衡误差应当围绕均衡值0波动。波动。也就是说，均衡误差也就是说，均衡误差t应当是一个平稳时间序列，即应当是一个平稳时间序列，即应有应有 tI(0），），E(t）= 0。 按照协整的定义，由于按照协整的定义，由于 YtI(1），

38、），XtI(1），且线性组合），且线性组合 t=Yt01XtI(0） 因此，因此，Yt 和和Xt是（是（1，1）阶协整的，即）阶协整的，即 Yt，XtCI(1, 1）协整向量是（协整向量是（1, 0, 1） 综合以上结果，我们可以说，两时间序列之间的协综合以上结果，我们可以说，两时间序列之间的协整是表示它们之间存在长期均衡关系的另一种方式。整是表示它们之间存在长期均衡关系的另一种方式。因此，若因此，若Yt 和和Xt是协整的是协整的, 并且均衡误差是平稳的且并且均衡误差是平稳的且具有零均值，我们就可以确信，方程具有零均值，我们就可以确信，方程 Yt =0+1Xt+t （7.20）将不会产生伪回归

39、结果。将不会产生伪回归结果。 由上可知，如果我们想避免伪回归问题，就应该在由上可知，如果我们想避免伪回归问题，就应该在进行回归之前检验一下所涉及的变量是否协整。进行回归之前检验一下所涉及的变量是否协整。二协整的检验二协整的检验 我们下面介绍用于检验两变量之间协整的两种简我们下面介绍用于检验两变量之间协整的两种简单方法。单方法。1、Engle-Granger法法 步骤步骤1.用上一节介绍的单位根方法求出两变量的单整的阶，用上一节介绍的单位根方法求出两变量的单整的阶，然后分情况处理然后分情况处理, 共有三种情况：共有三种情况：（1）若两变量的单整的阶相同，进入下一步；若两变量的单整的阶相同，进入下

40、一步；（2）若两变量的单整的阶不同，则两变量不是协整若两变量的单整的阶不同，则两变量不是协整的；的；（3）若两变量是平稳的，则整个检验过程停止，因若两变量是平稳的，则整个检验过程停止，因为你可以采用标准回归技术处理。为你可以采用标准回归技术处理。 步骤步骤2. 若两变量是同阶单整的，如若两变量是同阶单整的，如I(1），则用），则用OLS法估法估计长期均衡方程（称为协整回归）：计长期均衡方程（称为协整回归）： Yt=0+1Xt+t并保存残差并保存残差et，作为均衡误差，作为均衡误差t的估计值。的估计值。 应注意的是，虽然估计出的协整向量（应注意的是，虽然估计出的协整向量（1， ， ）是）是真实协

41、整向量（真实协整向量（1，0，1）的一致估计值，这些系数的标）的一致估计值，这些系数的标准误差估计值则不是一致估计值。由于这一原因，标准误差估准误差估计值则不是一致估计值。由于这一原因，标准误差估计值通常不在协整回归的结果中提供。计值通常不在协整回归的结果中提供。01步骤步骤3. 对于两个协整变量来说，均衡误差必须是平稳的。对于两个协整变量来说，均衡误差必须是平稳的。为检验其平稳性，对上一步保存的均衡误差估计值为检验其平稳性，对上一步保存的均衡误差估计值（即协整回归的残差（即协整回归的残差et）应用单位根方法。具体作法）应用单位根方法。具体作法是将是将DickeyFuller检验法用于时间序列

42、检验法用于时间序列et，也就是，也就是用用OLS法估计形如下式的方程：法估计形如下式的方程： et =et-1 + t （7.21） 有两点须提请注意：有两点须提请注意：（1）（）（7.21）式不包含常数项，这是因为）式不包含常数项，这是因为OLS残差残差et应以应以0为为中心波动。中心波动。（2）DickeyFuller统计量不适于此检验，表统计量不适于此检验，表7.3提供了用于提供了用于协整检验的临界值表。协整检验的临界值表。表7.3 协 整 检 验EG或A EG的 临 界 值变 量 个 数m=2m=3m=4显 著 性 水 平样 本 容 量0.010.050.100.010.050.100

43、.010.050.1025-4.37-3.59-3.22-4.92-4.10-3.71-5.43-4.56-4.1550-4.12-3.46-3.13-4.59-3.92-3.58-5.02-4.32-3.98100-4.01-3.39-3.09-4.44-3.83-3.51-4.83-4.21-3.89 -3.90-3.33-3.05-4.30-3.74-3.45-4.65-4.10-3.81例7.2 某 国 私 人 消 费 和 个 人 可 支 配 收 入 的 协 整 。 第 一 步 ： 求 出 两 变 量 的 单 整 的 阶私 人 消 费 变 量 ：tC=12330.48-0.01091C

44、t-1 (7.22) (t:)(5.138) (-1.339) R2=0.052 DW=1.7652tC=7972.671-0.85112 Ct-1 (7.23) (t:) (4.301) (-4.862) R2=0.425 DW=1.967个人可支配收入变量：tY=19903.93-0.02479Yt-1 (7.24)(t:) (3.054) (-1.387) R2=0.055 DW=2.2702tY=12889.39-1.11754Yt-1 (7.25) (t:) (3.983) (-6.270) R2=0.551 DW=2.014 由表由表7-37-3中可见，中可见，C Ct t和和Y

45、Yt t都是非平稳的，而都是非平稳的，而CCt t和和YYt t都是平稳的。这就是说，都是平稳的。这就是说， C CtI(1），），YtI(1）因而我们可以进入下一步。因而我们可以进入下一步。第 二 步 ， 进 行 协 整 回 归 ， 结 果 如 下 ：tC = 11907.23+0.779585Yt (7.26)( t:) (3.123) (75.566)R2=0.994 DW =1.021同 时 我 们 计 算 并 保 存 残 差 （ 均 衡 误 差 估 计 值 ） et。 第 三 步 ， 检 验et的 平 稳 性 。te = -0.51739et-1 (7.27)(t:) (-3.150

46、)R2=0.224 DW =1.948 第四步，得出有关两变量是否协整的结论。第四步，得出有关两变量是否协整的结论。 用用t3.150与表与表73中的临界值相比较（中的临界值相比较（m=2），），采用显著性水平采用显著性水平=0.05，t大于临界值大于临界值，因而接受，因而接受et非平稳的原假设，意味着两变量不是协整的，我们非平稳的原假设，意味着两变量不是协整的，我们不能说在私人消费和个人可支配收入之间存在着长不能说在私人消费和个人可支配收入之间存在着长期均衡关系。期均衡关系。 可是，如果采用显著性水平可是，如果采用显著性水平=0.10，则，则3.150与表与表73 中的临界值大致相当，因而可

47、以预期，若中的临界值大致相当，因而可以预期，若=0.11，t将小于临界值将小于临界值，我们接受，我们接受et为平稳的备择为平稳的备择假设，即两变量是协整的，或者说两变量之间存在假设，即两变量是协整的，或者说两变量之间存在着长期均衡关系。着长期均衡关系。2、Durbin-Watson法法 此方法非常简单，步骤如下：此方法非常简单，步骤如下：步骤步骤1. 估计协整回归方程估计协整回归方程 Yt=0+1Xt+t 保存残差保存残差et，计算，计算DW统计值（现称为统计值（现称为“协整回协整回归归”DurbinWatson统计值（统计值（CRDW），）， 即即 CRDW= 其中其中 为残差的算术平均值。

48、为残差的算术平均值。221)()(eeeettte步骤步骤2. 根据下述原假设和备择假设得出有关两变量协整根据下述原假设和备择假设得出有关两变量协整的结论：的结论： H0：et非平稳，即非协整非平稳，即非协整 H1： et平稳，平稳， 即协整即协整 若若CRDWd，则接受原假设，则接受原假设H0； 若若CRDWd，则拒绝原假设，则拒绝原假设H0。 这里原假设成立的临界这里原假设成立的临界d值为值为d = 0，对应于显著性，对应于显著性水平为水平为0.01，0.05和和0.10的临界值分别为的临界值分别为 0.511，0.386和和0.322。例例7.3 某国私人消费和个人可支配收入的协整某国私

49、人消费和个人可支配收入的协整 将将CRDW应用于上例。应用于上例。 第一步：由上例中（第一步：由上例中（7.26）式知）式知CRDW=1.021 第二步：因为第二步：因为CRDW=1.021大于上面提到的临界值大于上面提到的临界值, 故拒绝原假设，接受备择假设，因此得出结论：故拒绝原假设，接受备择假设，因此得出结论： 私人消费和个人可支配收入可以协整。私人消费和个人可支配收入可以协整。*三误差修正模型（三误差修正模型（ECM）的估计）的估计 协整分析中最重要的结果可能是所谓的协整分析中最重要的结果可能是所谓的“格兰杰代格兰杰代表定理表定理”（Granger representation the

50、orem）。按照）。按照此定理，如果两变量此定理，如果两变量Yt和和Xt是协整的，则它们之间存是协整的，则它们之间存在长期均衡关系。在长期均衡关系。 当然，在短期内，这些变量可以是不均衡的，扰动当然，在短期内，这些变量可以是不均衡的，扰动项是均衡误差项是均衡误差t。两变量间这种短期不均衡关系的动。两变量间这种短期不均衡关系的动态结构可以由误差修正模型（态结构可以由误差修正模型（error correction model）来描述，来描述，ECM模型是由模型是由Sargan提出的。这一联系两提出的。这一联系两变量的短期和长期行为的误差修正模型由下式给出：变量的短期和长期行为的误差修正模型由下式给

51、出： Yt = 滞后的（滞后的（Yt，Xt）+t-1 + vt （7.28） 10 其中其中 YtI(1），），XtI(1） Yt ，XtCI (1，1） t= Yt01XtI（0） vt=白噪声，白噪声，为短期调整系数。为短期调整系数。（7.28）式是）式是ECM模型的一般形式，实践中可根据模型的一般形式，实践中可根据情况建立具体的情况建立具体的ECM模型。最简单的是一阶模型。最简单的是一阶ECM模型，形式如下：模型，形式如下：011ttttYXv不难看出，在（不难看出，在（7.28）中，所有变量都是平稳的，）中，所有变量都是平稳的，因为因为 YtI (1), XtI (1)YtI (0),

52、 XtI (0) Yt, XtCI (1, 1) tI (0)因此，有人或许会说，该式可用因此，有人或许会说，该式可用OLS法估计。但事法估计。但事实上不行，因为均衡误差实上不行，因为均衡误差t不是可观测变量。因而不是可观测变量。因而在估计该式之前，要先得到这一误差的值。在估计该式之前，要先得到这一误差的值。Engle 和和 Granger建议采用下述两步方法估计方程建议采用下述两步方法估计方程(7.28):第一步：估计协整回归方程第一步：估计协整回归方程 Yt=0+1Xt+t得到协整向量的一致估计值（得到协整向量的一致估计值（1， ， ），用），用它得出均衡误差它得出均衡误差t的估计值的估计

53、值 et= Yt Xt第二步：用第二步：用OLS法估计下面的方程法估计下面的方程 Yt = 滞后的（滞后的（Yt, Xt）+et-1+vt （7.29）0110例例7.4 估计某国私人消费和个人可支配收入之间的估计某国私人消费和个人可支配收入之间的误差修正模型。误差修正模型。第一步第一步 ：由例：由例7.2 中中7.26式协整回归的结果：式协整回归的结果： = 11907.23 + 0.779585Y= 11907.23 + 0.779585Yt t (7.30)(7.30) (t:) (3.123) (75.566) (t:) (3.123) (75.566) R R2 2=0.994 DW

54、=1.021=0.994 DW=1.021 我们得到残差我们得到残差et。tC 第二步：估计误差修正模型，结果如下：第二步：估计误差修正模型，结果如下： = 5951.557+0.28432Y= 5951.557+0.28432Yt t 0.19996e0.19996et-1 t-1 (7.31)(7.31) (t:) (7.822) (6.538) ( (t:) (7.822) (6.538) (2.486)2.486) R R2 2=0.572 DW=1.941=0.572 DW=1.941 （7.31）中的结果表明个人可支配收入）中的结果表明个人可支配收入Yt的短期变的短期变动对私人消费

55、存在正向影响。此外，由于短期调整系动对私人消费存在正向影响。此外，由于短期调整系数是显著的，表明每年实际发生的私人消费与其长期数是显著的，表明每年实际发生的私人消费与其长期均衡值的偏差中的均衡值的偏差中的20%（0.19996）被修正。）被修正。tC 以韩德瑞(D.F.Hendry)为代表的动态建模方法（也称为伦敦经济学院(LSE)方法）指出自回归分布滞后模型（ADL）是最通用的线性模型形式。当变量为非平稳时间序列时ADL模型尤为适用，因为只要模型包括了足够多的滞后项，就一定能摆脱单位根的困扰。当变量间存在协整关系时，ECM模型便成为ADL模型的一个特例。以一阶ADL模型（7.32)为例，He

56、ndry对模型变量进行了等价变换，得到（7.33）所示的ECM模型。012131(7.32)tttttYXYXtttttttttttttttttXYXXYXXYXXYXYYY123120211231121013112101312101t)11)(1()1)(1()() 1() 1(即013t12122(1)()(7.33)11tttYXYX 式（7.33）将 依次分解为三个具有不同含义的部分：短期扰动、非均衡项和白噪声。 称为负反馈系数。 当YtI(1)，XtI(1）时，式（7.33）方程左边YtI (0），方程右边XtI (0），tI (0）。如果非均衡项 I (0)，则 Yt 与 Xt 存

57、在（1，1）阶协整关系。tY0) 1(2123120)11(tXYHendry论证了 对应经济理论模型中的长期均衡解，它自身不含任何变动的趋势。当外生变量的波动引起 时，该相对于长期均衡解的非均衡项在负反馈系数的作用下引起的延迟波动，促使重新回到其长期均衡解，因此称式（7.33）为“均衡修正模型”或“误差修正模型”。ttXY23120110)11(123120tXY 实际建模中，Hendry的动态建模方法主张从“一般到特殊”的原则，从包含被解释变量的最广泛影响因素的ADL模型开始，逐级约化，每一步约化都需要满足各项检验标准，力求在数据信息损失最小的情况下得到包含被解释变量长期均衡关系的最简洁的

58、ECM模型，有效避免了“伪回归”问题。这一动态建模方法已成为当今主流经济计量建模方法之一。例7.5 运用动态建模方法估计某国私人消费和个人可支配收入之间的误差修正模型。 第一步：确定私人消费和个人可支配收入的单整阶数，由例7.1知: CtI(1）和YtI(1）。 第二步：建立ADL模型。取ADL模型滞后阶数为2时，运用OLS法，方程估计通过自相关、异方差、正态分布等各项检验。表明可以从滞后阶数为2的ADL模型开始对方程进行约化。 1212278400.81520.07690.272( :)(2.84) (4.31)(0.461)(5.91)0.09420.0957(7.34)( 1.25)(

59、1.40)0.9991.953387.68ttttttCCCYtYYRDWSigma1129255 0.79940.27890.125(7.35)( :)(4.72) (9.83)(6.19)( 1.96)0.9991.953355.41ttttCCYYtRDWSigma第三步：逐级约化ADL模型为最简化模型，原则是在通过各项检验标准的条件下，运用OLS法逐步略去方程中t检验值最不显著的变量。具体将式（7.33）中的Ct-2和Yt-2分别依次略去，最终得到各变量均显著的最简化模型（7.35）。可以看出，拟合方程的标准差由式（7.34）的3387.68下降为式（7.35）的3355.41，方程得

60、到了优化。第四步：将最简化模型因变量改写为一阶差分形式，并设定ECM项。1129255 0.27890.20060.1538(7.36)( :)(4.72) (6.19)( 2.47)(2.42)0.5771.953355.41ttttCYCYtRDWSigma依据式（7.33）的推导，令ECMt = Ct - 0.767*Yt，其中-0.767 =0.1538/(-0.2006)。由于式（7.36）只是对式（7.35）进行了变量的等价变换，因此方程的标准差没有发生变化，均是3355.41。可以验证ECMtI(0)，即Ct与Yt存在协整关系。 第五步：运用OLS法得到包含ECM项的误差修正模型