时间序列数据的基本回归分析

1、第十章 时间序列数据的基本回归(hugu)分析共四十页10.1 时间序列(xli)数据的性质我们应该怎样认识时间序列数据的随机性？回答：很明显，经济时间序列满足作为随机变量结果所要求的直观条件，这些变量的结果都无法事先预料到。（例如，我们今天(jntin)不知道道琼斯工业指数在下一个交易日收盘时会是多少，我们也不知道加拿大下一年的年产出增长会是多少。）规范地，一个标有时间脚标的随机变量序列被称为一个随机过程（stochastic process）或时间序列过程（time series process）。共四十页10.2 时间序列回归(hugu)模型的例子1、静态模型我们将有两个变量（例如y和z

2、）的时间序列数据标注(bio zh)相同的时期，将这样的y和z联系起来即为一个静态模型（static model）：“静态模型”的名称来源于我们正在模型化y和z的同期关系的事实。在一个静态回归模型中也可以有几个解释变量。2、有限分布滞后模型在有限分布滞后模型（finite distributed lag model，FDL）中，我们容许一个或多个变量对y的影响有一定时滞。共四十页考察一个二阶FDL：（1）当z发生一个暂时性的提高时， 则表示z在t时期提高一个单位所引起y的即期变化。 通常被称作冲击倾向(qngxing)（impact propensity）或冲击乘数（impact multip

3、lier）。（注意： 分别表示这一暂时变化发生后，下一时期、两个时期、j个时期后y的变化如图10.1）（2）当z从t期开始永久性提高，一期后y提高了 ，两期后y提高了 。 这表明，z的当期和滞后系数之和 ，等于z的永久性提高导致y的长期变化，它被称为长期倾向（long-run propensity, LRP）或长期乘数（long-run multiplier）。共四十页共四十页一个q阶有限分布滞后模型可写成：静态模型是上式的一种(y zhn)特例，当 都为0即可。冲击倾向总是同期z的系数 。长期倾向便是所有变量 的系数之和。共四十页10.3 经典(jngdin)假设下OLS的有限样本性质假定

4、TS.1（线性于参数）假定 TS.2（无完全共线性）：在样本中，没有任何自变量是恒定不变的，或者是其他自变量的一个完全线性组合。假定 TS.3（零条件均值）：假定 TS.4（同方差性）：该假定意味着， 不能依赖于X（只要 和X相互独立就足够了满足TS.3即可），且在所有时期都保持不变。假定 TS.5（无序列相关）：【提问：我们为什么不假定不同横截面观测的误差是无关的呢？答：前述有随机抽样的假定，则以样本中所有解释变量为条件，不同观测的误差是独立的。因此，就我们当前目的而言，序列相关只是时间(shjin)序列和回归中的一个潜在问题。】假定 TS.6（正态性）：误差 独立于X，且具有独立同分布共四

5、十页定理 10.1（OLS的无偏性） 在假定TS.1、TS.2和TS.3下，以X为条件，OLS估计量是无偏的，并因此下式也无条件地成立：定理10.2（OLS的样本方差） 在时间序列高斯-马尔可夫假定TS.1-TS.5下，以X为条件， 的条件方差为： 其中(qzhng)， 是 的总平方和， 为由 对所有其他自变量回归得到的共四十页定理10.3（ 的无偏估计） 在假定TS.1-TS.5下，估计量 是 的一个无偏估计量，其中df=n-k-1定理10.4（高斯-马尔可夫定理） 在假定TS.1-TS.5下，以X为条件，OLS估计量是最优线性无偏估计量。定理10.5（正态抽样分布） 在时间序列的CLM假定

6、TS.1-TS.6下，以X为条件，OLS估计量遵循正态分布。而且，在虚拟假设下，每个t统计(tngj)量服从t分布，F统计量服从F分布，通常构造的置信区间也是确当的。共四十页例10.1 静态菲利普斯曲线研究失业和通货膨胀之间是否存在替代关系。 H0： H1：文件(wnjin)：PHILLIPS.RAW命令：reg inf unem结果： 上述方程并没有表明unem和inf之间存在替代关系（因为 ）分析中可能存在的问题：（1）CLM假定不成立（12章）；（2）静态菲利普斯曲线不是最佳模型（附加预期的菲利普斯曲线）共四十页例10.2 通货膨胀和赤字对利率的影响1948-2003年数据。 i3：三月

7、期国债利率； inf：据消费者价格指数得出的年通货膨胀率 def：联邦赤字占GDP 的百分比文件：INTDEF.RAW命令：reg i3 inf def结果：Inf与def对于i3的影响在统计上十分显著，即通货膨胀上升或赤字相对规模的扩大都会提高短期(dun q)利率。（但前提是CLM假定成立）共四十页10.4 函数形式、虚拟变量(binling)和指数在应用研究中经常出现具有恒定百分比效应的时间序列回归（自然对数形式）将对数函数形式用于分布滞后模型：方程中的冲击倾向 也被称为短期弹性（short-run elasticity）：它度量了GDP增长(zngzhng)1%时货币供给的即期百分比变

8、化；长期倾向 有时也被称为长期弹性（long-run elasticity）：它度量了GDP持久地增长1%，4个月后货币供给的百分比变化。共四十页二值或虚拟自变量在时间序列应用中也相当有用。既然观测单位是时间，所以虚拟变量代表某特定事件在每个时期是否发生。在事件研究（event study）中，二值变量是关键(gunjin)成分。事件研究的目标是为了确定某个特定的事件是否会影响到某项结果。讨论指数（index number）的概念： （1）基期、基值； （2）标准的经济产出都是用真实价值表示的；共四十页例10.3 波多黎各的就业和最低工资研究美国的最低工资对波多黎各就业的影响。 prepopt

9、：波多黎各第t年的就业率（就业人口占总人口的比例）； usgnpt：美国的真实(zhnsh)国民生产总值（以10亿美元计） mincov：度量最低工资相对于平均最低工资的重要性。mincov=(avgmin/avgwage)*avgcov，其中，avgmin是平均最低工资，avgwage是总体平均工资，avgcov是平均工资覆盖率。共四十页文件：PRMINWGE.RAW命令： reg lprepop lmincov lusgnp结果：prepop对mincov的估计弹性(tnxng)是-0.154，而根据t=-2.37，它在统计上是显著的。因此，更高的最低工资降低了就业率，这与古典经济学的预言

10、一样。共四十页例10.4 个人税收豁免对生育率的影响总生育率（gfr）是每个1000个育龄妇女生育孩子的个数。对1913-1984年这段时间，方程 pe：个人税收减免的实际美元金额； ww2：在1941-1945年间为1（第二次世界大战）； pill：从避孕药开始用于控制生育的1963年后一直为1文件(wnjin)：FERTIL3.RAW命令：sum pe reg gfr pe ww2 pill结果：共四十页考虑生育率对pe变化的反应滞后，估计一个包含两期滞后的分布滞后模型命令：reg gfr pe ww2 pill pe_1 pe_2在这个回归(hugu)中，我们只有70次观测，这是因为pe

11、滞后两次减少了2次观测。pe变量的系数估计得很不准确，每一个变量都不是个别显著的。事实上，pet，pet-1和pet-2明显相关，这种多重共线性使得估计每个滞后的影响非常困难。共四十页1、pet，pet-1和pet-2是联合显著的，F统计量的p值为0.012。命令：test pe pe_1 pe_2因此，pe的确对gfr有影响，但我们并没有足够好的估计值判断这种影响是即期的，还是存在一期或者两期的滞后(zh hu)（或都有一些）。2、实际上， pet-1和pet-2不是联合显著的，因而我们使用静态模型还算合理。 命令：test pe_1 pe_2共四十页式（10.19）中估计的（命令：disp

12、lay _bpe+_bpe_1+_bpe_2） LRP=0.073-0.0058+0.034=0.101但我们从式（10.19）中无法得到(d do)这个估计值的标准误。为得到LRP估计值标准误的技巧：令 表示LRP，并将 代入模型便得到共四十页基于上式，可通过将gfrt对pet，（pet-1-pe），（pet-2-pet），ww2t和pillt进行(jnxng)回归而得到 及其标准差。命令： gen dif1=pe_1-pe gen dif2=pe_2-pe reg gfr pe dif1 dif2 ww2 pill说明 在较小的显著性水平上异于0。本例说明：即使 都不是个别显著的，但LRP

13、非常显著。共四十页例10.5 反倾销调查和化学产品进口在美国开展反倾销调查，而后制定反倾销生产条例的过程中，一些有意思的问题：（1）在反倾销调查前的一段时期进口量异常吗？（2）反倾销调查后进口有明显的变化(binhu)吗？（3）有利于美国产业的决策执行后，进口究竟减少了多少？定义的3个虚拟变量：befile6：在开始调查前的六个月为1；affile6：表示开始调查后的六个月；afdec6：代表调查结束并确认构成倾销行为后的六个月；因变量chnimp：从中国进口的数量（取对数形式）；解释变量包括：（1）化工产量指标chempi；（2）石油产量gas；（3）汇率指标rtwex；（均使用对数形式）共

14、四十页文件：BARIUM.RAW命令： reg lchnimp lchempi lgas lrtwex befile6 affile6 afdec6结果：计算出准确的百分比变化（决策执行后，进口减少(jinsho)的比率）命令：display 100*(exp(_bafdec6)-1)共四十页例10.6 选举结果和经济形势(xngsh)费尔利用1916-1992年（每4年一次）的数据 得到的20次观测，解释了两党选举中民主党候选人获得选票的比例。估计费尔模型的一个简化形式：其中： demvote：两党选举中民主党候选人获得选票的比例； partyWH：虚拟变量，民主党在白宫执政时取值为1， 共

15、和党执政时取值为-1； incum：民主党在任总统参加竞选时定义为1，共和党在任总统参加竞选时定义为-1，其他情况为0； gnews：现任政府执政的前15个季度中，人均真实产出增长率超过2.9%的季度数； inf：本届政府前15个季度的年均通货膨胀率。当partyWH为1时， 度量了好的经济消息对执政党的影响(yngxing)； 度量了通货膨胀对执政党的影响(yngxing)。关注交互项哦！共四十页文件：FAIR.RAW命令： reg demvote partyWH incum pWHgnews pWHinf if year1996结果：解释：除了partyWH外，所有变量都在5%的水平上显著

16、。处于执政党的位置可以带来相当于所得(su d)选票份额5.4%的选票。另外，好的经济消息会产生正的影响：每个季度的好消息相当于1.1个百分点。通货膨胀会产生负的影响：如果平均年通货膨胀率上升2个百分点，执政党在选举中会失去1.5个百分点的选票。共四十页利用上述方程(fngchng)的预测：display _b_cons+_bpartyWH+_bincum+_bpWHgnews*3+_bpWHinf*3.019预测结果：共四十页10.5 趋势(qsh)和季节性描述有趋势的时间序列很多经济(jngj)时间序列都有随着时间而上升的共同趋势。忽略两个序列按相同或相反趋势延伸的事实，会导致如下错误结论

17、：认为一个变量的变化由另一个变量的变化所致。在很多情况下，两个时间序列过程表现出相关性仅仅是因为，由于某些无法观测因素的作用，二者具有共同的时间趋势而已。共四十页线性时间趋势（linear time trend）：各期变化(binhu)值相同共四十页指数趋势（exponential trend）： 各期具有相同(xin tn)的平均增长率共四十页在回归分析中使用趋势变量仅因为每个变量都随着时间的推移而增长，便得到两个或多个趋势变量之关系的现象，便是谬误回归（spurious regression problem）考虑一个yt受两个可观测因素xt1和xt2影响的模型。除了这两个变量以外，还有一些

18、无法观测的因素也随着时间的推移而系统地增长或缩减。满足(mnz)以上特征的模型为：它可以理解成xt3=t是的多元线性回归。如果上式省略掉t而只做yt对xt1和xt2的回归，一般会得到 和 的偏误估计值。以下例说明时间趋势如何导致谬误回归。共四十页例10.7 住房投资与价格对美国1947-1988年住房投资和住房价格指数的年度观测。文件：HSEINV.RAW变量含义(hny)： invpc：真实人均住房投资（以千美元计）； price：住房价格指数（1982=1）。命令1：reg linvpc lprice结果1： 人均投资对价格的弹性非常(fichng)大，且统计上显著；但我们要小心此处inv

19、pc和price都有上升的趋势。共四十页命令(mng lng)2：reg linvpc t reg lprice t结果2：命令3：reg linvpc lprice t结果3：趋势系数和标准误（虽然不一定(ydng)可靠）揭示了上升趋势现在结论大不相同：估计出的价格弹性是负的，而且在统计上也非显著异于0。因而前一回归方程为invpc和price之间的谬误关系。共四十页在有些情形中，若自变量和因变量有不同类型的趋势（比如一个向上而另一个向下），增加一个时间趋势可使关键解释变量更显著(xinzh)，但自变量围绕其趋势线的变动会导致因变量偏离其趋势线的变动。以例10.8来解释。共四十页例10.8

20、生育方程（基于例10.4的FERTIL3.RAW）（1）在生育方程中添加一个线性时间(shjin)趋势：命令：reg gfr pe ww2 pill t（2）采用二次趋势（观察总生育率在1913-1984年间表现出先上升后下降的趋势。）命令：reg gfr pe ww2 pill t tsq（1）、（2）中pe系数的估计值不断增大，并且更加(gnji)显著共四十页对包含时间趋势回归的除趋势解释在回归模型中引进时间趋势，相当于在回归分析中，在使用原始数据之前(zhqin)，便将它们除趋势（detrending）对于拟合方程 中的 和 可通过如下除趋势化的步骤得到：（1）将yt、xt1和xt2分别

21、对常数项和时间趋势t回归，并记录残差 、 和 ，t=1，2，n。例如，对于 和 的解释与此类似。（2）做 对 和 的回归。这个回归得出的 和 与前式中相同。这意味着，我们最感兴趣的系数估计值（ 和 ）来自于一个没有时间趋势的回归，在这个回归中，我们首先除去了因变量和所有自变量的趋势。共四十页例10.9 波多黎各的就业（基于例10.3的PRMINWGE.WAGE）加入一个线性趋势，估计结果为：log(usgnp)的系数发生了显著变化：从不显著的-0.012提高到非常显著的1.06。最低工资的系数只发生微小变化，然而(rn r)标准误明显变小了，而标准误的变小使log(mincov)比以前更加显著

22、。共四十页因变量有趋势时R2的计算当因变量含有趋势时，时间序列回归(hugu)中的普通或调整R2可能会认为地变大。为解决这一问题，我们首先做yt对t的回归，得到残差 。然后， 将 对xt1，xt2，t回归这一回归中的R2能够更好地反应出xt1，xt2能在多大程度上解释yt，因为它过滤掉了时间趋势的影响。共四十页例10.10 住房投资（基于例10.7的HSEINV.RAW）除去log(invpc)中的趋势，并将得到的变量对log(price)和t做回归(hugu)：predict linvpch, resreg linvpch lprice t对比除趋势之前的回归结果：reg linvpc lprice t除趋势后的R2表明，log(price)围绕其趋势的变动(bindng)，对log(