第六章 分布滞后模型与自回归模型

1、1 第六章 自回归模型和 分布滞后模型 2 第一节 分布滞后模型和自回归模型的概念 第二节 分布滞后模型的估计 第三节 部分调整模型和适应预期模型 第四节 自回归模型的估计 第五节 阿尔蒙多项式分布滞后 第六节 格兰杰因果关系检验 3 第一节 分布滞后模型和自回归模型的概念 很多经济过程的实现需要若干周期的时间，因此 需要在我们的计量经济模型中引入一个时间维，通 常的作法是将滞后经济变量引入模型中。让我们用 两个简单的例子说明之。 4 例6.1 yt = +xt-1 + ut, t = 1,2,n 本例中y的现期值与x的一期滞后值相联系，比较一 般的情况是： yt = +0xt +1xt-1

2、+sxt-s + ut, t = 1,2,n 即y的现期值不仅依赖于x的现期值，而且依赖于x的 若干期滞后值。这类模型称为分布滞后模型分布滞后模型，因为x 变量的影响分布于若干周期。 5 例6.2 yt = +yt-1 + ut, t = 1,2,n 本例中y的现期值与它自身的一期滞后值相联系， 即依赖于它的过去值。一般情况可能是： yt = f (yt-1, yt-2, , x2t, x3t, ) 即y的现期值依赖于它自身若干期的滞后值，还依赖 于其它解释变量。 在本例中，滞后的因变量（内生变量）作为解释变 量出现在方程的右端。这种包含了内生变量滞后项 的模型称为自回归模型自回归模型。 6

3、动态经济模型 我们上面列举了模型中包含滞后经济变量的两种 情况。第一种是仅包含滞后外生变量的模型，第二 种是包含滞后内生变量的模型。在两种情况下，都 通过一种滞后结构将时间维引入了模型，即实现了 动态过程的构模。 7 第二节 分布滞后模型的估计 我们在上一节引入了分布滞后模型： yt =+0xt +1xt-1 +sxt-s + ut (1) 在这类模型中，由于在x和它的若干期滞后之间 往往存在数据的高度相关，从而导致严重多重共线 性问题。因此，分布滞后模型极少按（1）式这样的 一般形式被估计。通常采用对模型各系数j施加某 种先验的约束条件的方法来减少待估计的独立参数 的数目，从而解决多重共线性

4、问题。这方面最著名 的两种方法是科克(koyck)方法和阿尔蒙(almon)方 法。 8 一、科克分布滞后模型 科克方法简单地假定解释变量的各滞后值的系数 （有时称为权数）按几何级数递减，即： 其中 01 这实际上是假设无限滞后分布，由于0，即长期影响大 于短期影响。 ,x ,y (1)(8)yxy (9) 1 yx 15 从实践的观点来看，科克变换模型很有吸引力， 一个ols回归就可得到、和的估计值。这显 然比前面介绍的格点搜索法要省时很多，大大简化 了计算。 可是，科克变换后模型的扰动项为 ut-ut-1 ， 这带来了自相关问题（这种扰动项称为一阶移动平 均扰动项）。并且，解释变量中包含了

5、yt-1，它是 一个随机变量，从而使得高斯马尔柯夫定理的解 释变量非随机的条件不成立。此问题的存在使得 ols估计量是一个有偏和不一致估计量。这可以说 是按下葫芦起了瓢。我们将在第四节中讨论科克模 型的估计问题。 16 第三节 部分调整模型和适应预期模型 有两个著名的动态经济模型，它们最终可化成与 上一节（2）式相同的几何分布滞后形式，因此都 是科克类型的模型。它们是： 部分调整模型 （partial adjustment model, pdm） 适应预期模型 （adaptive expectations model, aem） 17 一、部分调整模型 在部分调整模型中，假设行为方程决定的是因

6、变 量的理想值（desired value）或均衡值yt*，而不 是其实际值yt： yt* =+xt+ut （1） 由于yt*不能直接观测，因而采用 “部分调整假 说” 确定之，即假定因变量的实际变动（ytyt- 1）,与其理想值和前期值之间的差异（yt* yt-1） 成正比： yt yt-1=（yt* - yt-1） (2） 01, 称为调整系数。 18 从（3）式可看出，yt是现期理想值和前期实际值 的加权平均。的值越高，调整过程越快。如果 =1，则yt=yt*,在一期内实现全调整。若=0，则 根本不作调整。 （2）式 yt yt-1=(yt* - yt-1) (2） 可改写为： yt =

7、yt* +(1-) yt-1 (3） 19 (1）式 yt* =+xt+ut 代入（3）式 yt =yt* +(1-) yt-1 ，得到 yt=+xt+(1-)yt-1+ut （4） 用此模型可估计出、和的值。 与科克模型类似，这里也存在解释变量为随机变 量的问题（yt-1）.区别是科克模型中,yt-1与扰动项 （ut-ut-1）同期相关，而部分调整模型不存在同期 相关。在这种情况下，用ols法估计，得到的参数估 计量是一个一致的估计量。 20 不难看出，（4）式 yt=+xt+(1-)yt-1+ut （4） 与变换后的科克模型的形式相似，我们也不难通过对 （4）式中yt-1进行一系列的置换化

8、为几何分布滞后的 形式: （4）式两端取一期滞后，得 将此式代入（4）式，得到（为简单起见，省略扰动 项）： 1121 (1)(5) tttt yxyu 2 12 1 (1)(1)(1) tttt yxxy 21 我们可以用同样的方法置换yt-2，以及随后的yt-3, yt-4,直至无穷，结果是将yt表示为x的当前值和滞 后值的一个滞后结构，系数为科克形式的几何递减权 数，具体形式为： ttttt xxxy .)1 ()1 ( 2 2 1 .)1 ()1 ( 2 2 1 tttt uuu 与上节（2）式形式完全一样。 令=1-，=，则得 其中 2 12 .(6) ttttt yxxx 22 例

9、6.3 林特纳（lintner）的股息调整模型 jlintner建立的股息调整模型是应用部分调整模 型的一个著名实例。 在对公司股息行为的研究中，lintner发现，所有 股份公司都将其税后利润的一部分以股息的形式分配 给股东，其余部分则用作投资。 当利润增加时，股息一般也增加，但通常不会将增 加的利润都用作股息分配，这是因为： （1）利润的增加可能是暂时的； （2）可能有很好的投资机会。 23 为了建立一个描述这种行为的模型，lintner假设各 公司有一个长期的目标派息率，理想的股息dt*与现 期利润t有关，其关系为 dt*=t tttt uddd )( 1 * ttttt udddd )

10、( 1 * 1 ttt ud 1 其中ut为扰动项。因此 而实际股息服从部分调整机制 24 使用美国公司部门19181941年数据，得到如下回 归结果： 1 70. 015. 03 .352 ttt dd 各系数在1%显著水平下都显著异于0。 从回归结果可知，（1-）的估计值为0.70，因而 调整系数的估计值为0.30，即调整速度为0.30。由 于t的系数是的估计值，除以0.30，则得到长 期派息率（）的估计值为0.50。 tttt udd 1 )1 (即 25 1、在模型中考虑预期的重要性 预期（expectation）的构模往往是应用经济学家 最重要和最困难的任务，在宏观经济学中更是如此。

11、 投资、储蓄等都是对有关未来的预期很敏感的。如 果政府实施一项扩张政策，这将影响工商界人士有 关未来经济总状况的预期，特别是关于盈利能力的 预期，因而影响他们的投资计划。 例如，如果存在很可观的失业，则政府支出增加 被认为是有益的，并将刺激投资。另一方面，如果 经济正接近充分就业，则政府的扩张政策被认为将 导致通货膨胀，结果是工商界的信心受挫，投资下 降。 二、适应预期模型 26 2、适应预期模型 由上所述，可知在模型中考虑预期的重要性。不 幸的是，在宏观经济领域，不存在令人满意的直接 计量预期的方法。作为一种权宜之计，某些模型使 用一种称为适应预期过程适应预期过程的间接方法。 11 ()01

12、(8) eee tttt xxxx 适应预期过程是一种简单的学习过程，其机制是， 在每一时期中，将所涉及变量的当前观测值与以前 所预期的值相比较，如果实际观测值大，则将预期 值向上调整，如果实际观测值小，则预期值向下调 整。调整的幅度是其预测误差的一个分数，即： 27 （8）式可写成 1 (1)01(9) ee ttt xxx 适应预期和部分调整之间当然有很多明显的类 似之处，可是从适应预期模型的最初形式导出仅 包含可观测变量的模型（可操作模型）不象在部 分调整模型的情况那么简单。 上式表明，x的预期值是其当前实际值和先前预期 值的加权平均。的值越大，预期值向x的实际发 生值调整的速度越快。

13、28 假设你认为因变量yt与某个变量x的预期值xte有关， 则可写出模型 )10( t e tt uxy )11()1 ( 211 e tt e t xxx 我们用“降阶”法来解决这个问题。如果（9）式 成立，则对于t-1期，它也成立，即： 若假定xte 用适应预期机制确定，这就是一个适应 预期模型，其中解释变量xte是不可观测的，必须用可 观测变量取代之。 29 将（11）式 代入（9）式 ，得 )12()1 ()1 ( 2 2 1 e ttt e t xxxx , 2 e t x )13(.)1 ()1 ( 2 2 1 ttt e t xxxx 我们可以用类似的方法，消掉（12）式中的 这

14、一过程可无限重复下去，最后得到： e tt e t xxx 211 )1 ( e tt e t xxx 1 )1 ( 30 将（13）式 代入（10）式 ，得 )1 ()1 ( 11 ttttt uuyxy .)1 ()1 ( 2 2 1 ttt e t xxxx t e tt uxy 不难看出，此式与上节中科克分布（2）的形式相 同。该模型的参数可用上一节介绍的非线性方法估 计。对（14）式施加科克变换，将简化模型的数学 形式，但由于与科克模型同样的理由，不宜直接用 ols法估计。施加科克变换的适应预期模型为： 2 12 (1)(1).(14) ttttt yxxxu 31 例6.4 fri

15、edman的持久收入假说 1957年，弗里德曼对传统消费函数提出批评， 提出了自己的消费模型。在他的模型中，第i个消 费者在第t期的消费与持久性收入（permanent income）yitp有关，而不是与当期的收入yit有关。 持久性收入是一种长期收入概念，它表示在考 虑了各种可能的波动的情况下，某人大体上可以依 靠的收入。 32 持久收入是根据最近的经验和有关未来的预期而 主观决定的，由于是主观的，因而无法直接计量。任 何一年中的实际收入可能高于或低于持久收入，取决 于该年中的特别因素。实际收入和持久收入之差称为 暂时性收入 (transitory income)，记为yitt ， 我们有

16、： )15( t it p itit yyy 33 他以同样方式区分了持久性消费、实际消费和暂 时性消费的概念。持久性消费是与持久性收入的水 平相对应的消费水平。实际消费可能与持久消费有 差异，原因是出现了某些特殊的未预料到的情况 （如未预料到的医疗费用），或者是冲动性购买的 结果。二者之差称为暂时性消费，记为citt: )16( t it p itit ccc yitt和citt被假定为具有0均值和常数方差的随机 变量，它们相互独立，且与yitp和citp无关。 弗里德曼进一步假定持久消费与持久收入成正比： 34 上式中持久收入yitp不可观测，为解决这一问题， 弗里德曼假设持久收入遵从适应

17、预期过程，也就是说 ，如果某人的现期收入高于（或低于）其先前的持久 收入概念，则他将增加（或减少）后者，增加（或减 少）的幅度是二者之差乘以： )17( p it p it yc )18()( 1 p itit p it yyy 一般位于0和1之间。因此人们在实际收入增加时 将调整他们的持久收入概念，但不会做全额调整， 这是因为认识到实际收入的变动或许有一部分是由 于收入的暂时分量变动的结果。 35 （18）式可改写为： )19()( 11 p itit p it p it yyyy )20()1 ( 1 p itit p it yyy 此式表明，在第t年，消费者将持久收入估计为 实际收入和以

18、前的持久性收入概念的加权平均。 如果接近于1，则该消费者将绝大部分权重给了 实际收入，yp迅速向y调整，若接近0，则很小 部分权重给了实际收入，调整过程将很缓慢。 即 36 至此，我们得到了实际消费和持久收入之间的关系 式，即消费函数的弗里德曼模型。式中citt起着扰动 项的作用。 )21( t it p itit cyc 将（17）式 代入（16）式 我们有： )17( p it p it yc )16( t it p itit ccc 37 为了估计这个模型，弗里德曼用（20）式（适应 预期机制）将持久收入表示成实际收入的现期值和 各期滞后值： p itit p it yyy 1 )1 (

19、 p ititit yyy 2 2 1 )1 ()1 ( )22()1 ()1 ()1 ( 3 3 2 2 1 itititit yyyy 若01，这就是一个合理的假设，现期收入的权数 最大，上一年次之，随着时间往回推，影响逐年衰减 。最后，权数变得非常之小，使得无需考虑该年之前 那些过去值。 38 弗里德曼采用的估计方法是我们前面介绍过的非 线性方法，即首先试位于0和1区间内的大量值， 为每个值计算相应的持久收入时间序列，然后用 消费对每个持久收入数据集回归，根据r2选出最佳 值。 为了与传统消费函数相比较，弗里德曼用美国 19051951（战争期间除外）的人均实际消费和人 均可支配收入数据

20、进行了回归。在格点搜索计算中， 他将持久收入计算为现期收入和16个滞后收入项的 加权平均值，的最优值为0.37，得到消费函数中 的估计值为0.88。 39 第四节 自回归模型的估计 上两节中，我们讨论了下列三个模型： 科克模型 部分调整模型 适应预期模型 ) 1 ()()1 ( 11 ttttt uuyxy )2()1 ( 1tttt uyxy )3()1 ()1 ( 11 ttttt uuyxy 40 这种解释变量中包括因变量的滞后值的模型称为 自回归模型。由于在解释变量中包含了因变量的滞 后值,我们就可以动态地考察该变量在若干周期中的 变动,因此称为动态模型。 在自回归模型（4）中，由于随

21、机解释变量的存 在和序列相关的可能性这双重原因，ols法不能直 接应用，因此我们必须研究这类模型的估计问题。 )4( 1210tttt vyxy 这三个模型具有一种共同的形式，即： 41 一、自回归模型的估计问题 ols法的应用，要求解释变量xt为非随机的。在自 回归模型中，由于yt-1作为解释变量，这一条件已无 法满足，这是因为，由于 因此： 这表明，yt-1是随着随机扰动项vt-1的变动而变 动的，即yt-1部分地由vt-1决定，因而yt-1是随机变量。 tttt vyxy 1210 1221101 tttt vyxy 42 1. 解释变量为随机变量时ols估计量的统计性质 由第三章，当x

22、为非随机量这一条不满足时 （1）若每一个xt都独立于所有的扰动项ut，即 cov(xs,ut)=0, s=1,2,n t=1,2,n 则ols估计量仍为无偏估计量。 （2）若解释变量xt独立于相应的扰动因素ut，即随 机解释变量与扰动项同期无关 ： cov(xt,ut)=0, t=1,2,n 则ols估计量为一致估计量。 （3）若上述两条均不满足，则ols估计量既是有偏 的，又是不一致的。 43 2、自回归模型的估计问题 在自回归模型的情况下，第（1）条已无法满足，因 为yt-1显然可以表示为vt-1,vt-2,v1等的函数，因而依 赖于vt-1和所有早期的扰动因子。 现在让我们来看是否有可能

23、满足解释变量与扰动项 同期无关的条件，从而得到一个一致的估计量。 44 在自回归模型（4）的情况下： 也就是要求yt-1独立于vt,或 cov(yt-1,vt)=0 不难看出，只要扰动项vt是序列独立的（即自回 归模型（4）的各期扰动项相互独立），我们就可以 假定yt-1独立于所有未来的扰动因子（包括vt），在 这种假定下，yt-1与vt无关，我们对（4）式应用ols 得到的参数估计量是一致估计量。 0121tttt yxyv 45 让我们回到本节开始时列出的三个模型，看看我们 关于yt-1独立于所有未来的扰动因子，特别是yt-1与vt 无关的假定是否能成立。 在科克模型和适应预期模型中，扰动

24、因子序列独立 的条件不成立，以科克模型为例，扰动项 vt = ut-ut-1 假定ut满足标准假设条件，则容易证明 21 22 121 ttttttt uuuuuuue 2 2 )()( 2111 tttttt uuuuevve 该式非0，即vt序列相关。 46 我们还不难证明 即yt-1与vt相关。适应预期模型的情况与此类似。 )( ,),( 111 ttttt uuycovvycov 2 因此，对于科克模型和适应预期模型，应用ols 法不仅得不到无偏估计量，而且也得不到一致估计 量。也就是说，即使样本容量无限增大，参数估计 量也不趋向于其总体值。因此，不宜采用ols法估计 上述两种模型。

25、47 但是，部分调整模型不同，在该模型中， vt=ut，若ut满足标准假设条件，则vt也满足。因 此，可用ols法直接估计部分调整模型，将产生一致 估计值。 综上所述，ols法可用于部分调整模型的估计，并 提供一致的估计值。而科克模型和适应预期模型，则 由于其扰动项存在序列相关，用ols进行估计得到的 估计量既是有偏的，也是不一致的。 48 ols法不能应用于科克模型和适应预期模型的原 因是解释变量yt-1与扰动项vt相关，如果这种相关 能够被消除的话，我们就可以用ols得到一致估计 值。如何实现这一点呢？利维顿（liviatan）提出 的工具变量法是一种解决方法。 二、自回归模型的估计 49

26、 工具变量法的基本思路是当扰动项u与解释变量x 相关时，设法找到另一个变量z，z与x高度相关，而 与扰动项u不相关，在模型中，用z替换x，然后用 ols法估计，变量z称为工具变量工具变量。 只要工具变量的选取能够保证z与x高度相关，而与 u不相关，则我们得到的将是一致估计量。z与x的相 关程度越高，这种替代的效果就越好。我们下面回到 科克模型和适应预期模型，研究工具变量的选取。 50 我们的模型为 这里x是唯一的外生变量，而y的行为部分地依赖 于x的行为，yt-1的取值部分地取决于xt-1的数值。因 此，这里xt-1就是一个比较理想的工具变量，即用滞 后外生变量作为滞后内生变量的工具： zt=

27、xt-1 , t=1,2,n 来估计 tttt vyxy 1210 ntvzxy tttt ,.,2 , 1 210 51 在实践中，自回归模型还可以用极大似然法估计， 得到的估计量是一致估计量。 当然，对于本节所涉及的三种模型，由于它们都 是几何滞后模型，因而都可以用前面介绍的非线性 方法进行估计，该方法尽管费时，但没有估计问题。 52 第五节 阿尔蒙多项式分布滞后 （almon polynomial distributed lags） 科克分布假定滞后解释变量的系数按几何级数 递减。对于很多应用问题来说，这是一种令人满 意的近似，但对于另一些应用问题，这种假设就 未必符合现实情况。例如，在

28、某些情况下较现实 的假设是，因变量对解释变量变动的响应是，开 始小，然后随时间变大，尔后再次衰减，如下图 所示。 53 t 滞后的权数 54 阿尔蒙滞后分布的基本假设 阿尔蒙滞后分布为这类行为的构模提供了灵活 的选择，同时使待估计的参数数目大大减少。 基本假设是，如果y依赖于x的现期值和若干期 滞后值，则权数由一个多项式分布给出。由于这个 原因，阿尔蒙滞后也称为多项式分布滞后。最简单 的例子是二次和三次多项式的情况，如下图所示： 55 图2 二次函数 56 图3 三次函数 57 一般情况下，在分布滞后模型 中，假定： 其中p为多项式的阶数，如图2中p=2，图3中p=3。 也就是用一个p阶多项式

29、来拟合分布滞后，该多项 式曲线通过滞后分布的所有点。 由用户选择最大滞后周期m和多项式阶数p。 p pi iaiaiaa 2 210 tmtmttt uxxxy 110 58 例6.5 若你根据一实际问题设定下面的模型： 2 210 iaiaa i 00 a 2101 aaa 2102 42aaa 2103 93aaa 2104 164aaa ) 1 ( 443322110ttttttt uxxxxxy 我们有： 这表明，你所选择的最大滞后周期m=4，模型中共 有6个参数。若决定用二次式进行拟合，即p=2，则 59 代入原模型，得 t i iti ux 4 0 4 0 2 210 )( i t

30、it uxiaiaa t i it i it i it uxiaixaxa 4 0 2 2 4 0 1 4 0 0 令： z0t=xt-i, z1t=ixt-i, z2t=i2xt-i ttttttt uxxxxxy 443322110 60 （2）式中有4个参数，比（1）式的6个少了两个， 估计出，a0,a1,a2 的值之后，我们可以转换为i的 估计值,公式为： 2 210 iaiaa i ) 2( 221100ttttt uzazazay 显然，z0t,z1t和z2t可以从现有观测数据中得出， 使得我们可用ols法估计下式： 61 在实践中，人们期望m尽量小一些，如果有10年的 数据，通常

31、滞后取二至三期。 应用阿尔蒙滞后的关键在于如何选择最大滞后周期 m和多项式的阶数p。 一般说来，采用高阶多项式，拟合效果要好一些， 但出现多重共线性问题的可能性要比二阶、三阶多项 式大。一般情况下，三次多项式是一个不错的选择。 对于p，我们可直接由（2）式用t检验法检验 h0: ap = 0,如果接受原假设，我们就可以去掉ap，然后用 （p-1）阶来估计（2）式，如果h0: ap = 0被拒绝，我 们可以试（p+1）阶，并检验 h0: ap+1= 0，等等。 62 *第六节第六节 格兰杰因果关系检验格兰杰因果关系检验 (granger causality test) 相关并不意味着必定存在因果

32、关系。然而，许多经济相关并不意味着必定存在因果关系。然而，许多经济 变量之间存在单向或双向的影响关系。例如，宏观经变量之间存在单向或双向的影响关系。例如，宏观经 济学经常会问到的一个问题是：是济学经常会问到的一个问题是：是gdp“引起引起”货币供货币供 给给m，还是货币供给，还是货币供给m “引起引起” gdp？ 格兰杰因果关系检验法要检验的就是这类的因果关格兰杰因果关系检验法要检验的就是这类的因果关 系是否存在，比如说，系是否存在，比如说，x是否为是否为y的因？检验思路是，的因？检验思路是， 用当前的用当前的y对对y的若干期滞后及的若干期滞后及x的若干期滞后回归，的若干期滞后回归， 然后检验

33、然后检验x的这些滞后变量作为一个整体是否改善了回的这些滞后变量作为一个整体是否改善了回 归结果，如果回答是肯定的，则归结果，如果回答是肯定的，则x被称为被称为y的格兰杰原的格兰杰原 因因 （x granger causes y）。）。 63 值得注意得是，这种因果关系往往是双向的，值得注意得是，这种因果关系往往是双向的， x 是是y的格兰杰原因，的格兰杰原因，y又是又是x的格兰杰原因。的格兰杰原因。 格兰杰因果关系检验要求进行以下两个回归格兰杰因果关系检验要求进行以下两个回归： 1 11 2 11 (1) (2) mm tit iit it ii mm tit iit it ii yxyu xyxu 下面分四种