计量经济学第九章时间序列结构模型课件

1、第九章结构型时间序列模型时间序列回归模型分类：1 .不含外生变量的非结构型模型， 包括单方程模型（如 ARMA 模型） 和多方程模型（如向量自回归模型，VAR ）2 .传统的结构模型，包括含有外生变量的单方程回归模型（如确定性趋势或季节模型、静态模型、分布滞后模型、自回归分布滞后模型等）和联立方程模型3 .协整和误差修正模型等现代时间序列模型第二、三类模型反统称为结构型时间序列模型。本章将对最基本的几种结构型时间序列模型进行简要介绍。第一节 确定性趋势与季节模型确定性趋势与季节模型将经济变量看作是时间的某种函数，用于描述时间序列观测值的长期趋势特征和周期性变动特征。其中的自变量是确定性的时间变

2、量t 或反映季节的虚拟变量。由于自变量是非随机变量，自然是严格外生的，所以不涉及诸如非平稳性、高度持久等问题，一般可以如同横截面数据一样，直接使用经典线性模型的回归分析方法。一、确定性趋势模型（一）种类按照因变量y与时间t的关系不同，常用的确定性趋势模型主要有以下三 类：1. 线性趋势模型yt01tut（ 9.1）当时间序列的逐期增长量（即一阶一次差分ytytyt 1 ）大体相同时，可以考虑拟合直线趋势方程。2. 曲线趋势模型2kyt01t 2t2kt k ut（ 9.2）若逐期增长量的逐期增长量（二阶一次差分2ytytyt 1） 大致相同，可拟合二次曲线yt01t2t2 ut。类似地，如果事

3、物发展趋势有两个拐点，可以拟合三次曲线 yt0it2t23t3 Uto其他更高次的曲线趋势比较少用。3. 指数曲线模型yt0 1 e t（9.3）或ln（ yt）ln 0 （ln 1）tut指数曲线的特点是各期的环比增长速度大体相同（即自然对数的一阶一次差分yt / yt 1 ln ytln yt 1 基本为常数），时间序列的逐期观测值大致按一定的百分比递增或衰减。为了更好地表现事物发展的特征，我们还可以给指数曲线设置发展的上限或下限（渐进线）。常用的带渐进线的指数曲线有以下几种：（1）修正指数曲线模型y k 0 ；e k>0,0w0 , 01 1（9.4）(2)龚伯兹(Gompertz

4、)曲线模型(9.5)(9.6)Vt k o 1teut, k>0, 0V。*1, 0< 1 T(3)逻辑斯蒂(Logistic)曲线模型1，八一一Vt ;TV，k>0，产0，0V 产1k 0 1e(二)模型选择标准如果对同一时间序列有几种趋势线可供选择，在考虑经济意义的基础 上，线性模型可以参照第三章第四节有关多元回归模型的评价准则进行优 选。除此之外，还可以用下列指标作为辅助标准，选择这些指标比较小的模 型。1 .均方根误差(Root Mean Squared Error , RMSE )：RMSE（yt 夕）22 .平均绝对误差（Mean Absolute Error ,

5、 MAE ）: 1-MAE Tyt 又3 .平均绝对百分比误差（Mean Absolute Percent Error , MAPE ）：MAPE -yt yt 100%t I yt I（三）模型估计上述确定性趋势模型要么本身是参数线性的，要么通过变换模型形式可 以使其线性化，所以均属于广义的线性模型范畴。而且自变量（时间变量 t 及其函数）都是非随机变量（肯定符合严格外生条件），所以，可以直接使用OLS进行估计，所有估计与推断方法完全等同于横截面数据。唯一与横截面回归不同的是，时间序列往往存在误差项自相关（将在第 十章讨论）。例9-1改革开放以来我国GDP数据如表9-1所示。试估计我国实际

6、GDP的年平均增长率。表9-1部分年份中国GDP数据（1978年不变价，亿元）年份tGDP年份tGDP197813645.219941716505.9197923922.219951818309.2198034228.719961920141.7198144450.419972022014.2198254853.519982123738.7198365380.319992225547.5198476196.820002327701.5198587031.220012430000.8198697653.320022532725.51987108539.820032636006.419881195

7、03.120042739637.71989129889.220052843773.019901310268.920062948871.219911411211.420073055243.019921512808.020083160189.519931614596.620093265426.0资料来源：中国统计年鉴2010首先，描出GDP及其对数的时间序列图:70000 *60000 -50000 -.PDG. 40000 -.30000 r/*,20000 -* *.*10000 -011111105101520253035T3511.210.810.410.0PDGNL305 2 20 5

8、1 O 1 5可见，GDP与时间t之间不存在线性关系，而其对数与时间存在线性关 系。即GDP呈指数趋势发展:GDPt0 1teut线性化的趋势模型形式为ln GDPt ln 0 （ ln 1） t ut ，即 ln GDPt01t ut利用 OLS 估计，并消除误差项自相关（自相关问题见第十章），输出结果如下：表9-2 EViews输出结果Dependent Variable: LNGDPVariableCoefficientStd. Error t-StatisticProb.C8.0704290.014908541.35080.0000T0.0942070.000778121.06280.

9、0000AR(1)1.2815660.1311449.7721790.0000AR(2)-0.733140.128869 -5.6890800.0000R-squared0.999647Mean dependent var9.718339Adjusted R-squared0.999606S.D. dependent var0.833467-5.24254S.E. of regression0.016540Akaike info criterion6-5.05572Sum squared resid0.007113Schwarz criterion0Log likelihood82.63820

10、F-statistic24538.49Durbin-Watson stat2.013367Prob(F-statistic)0.000000Inverted AR Roots .64-.57i.64+.57i可见，消除误差项自相关后，回归方程为in GDP 8.0700.094 t回归结果说明，样本内我国实际 GDP年均增长率为0.094,即9.4%用水平法计算的GDP年均增长率为3165422/3645.2 1 0.098 9.8% 二、季节变动模型与横截面数据相比，时间序列往往受到季节变动的影响。所谓季节变动 是指经济变量因受自然因素或社会经济因素影响，从而形成的有规律的周期 性变动。这里

11、的 季节”一词是广义的，泛指任何一种有规律的、按一定周期（如 季、月、旬、周、日）重复出现的变化如果忽视其季节因素，回归中就会犯遗漏变量的错误，影响参数估计的 无偏性和一致性。根据回归分析的目的不同，对含有季节变动的数据有不同 的处理方法：1 .季节差分。如季度数据和月度数据分别采用 4阶和12阶差分。2 .季节修匀。即通过移动平均、指数平滑等方法，去除时间序列中的周 期性变化。3 .引入季节虚拟变量。例9-2某企业最近6年的商品销售额的季度数据表 9-3。试用虚拟变量 方法建立季节波动模型。表9-3某企业商品销售额的季度数据年份季度yttD2tD3tD4t200611132.63310002

12、1298.058210031207.310301041116.035400111580.5915000200721797.940610031709.866701041650.171800111856.4849000200822062.6241010032022.3101101042090.4561200112261.06813000200922518.5681410032440.0671501042465.8851600112673.70417000201022887.9031810032896.7401901042913.3852000113163.66521000201123486.161

13、2210033205.5292301043348.30724001数据的变化轨迹见图9-3。图9-3销售额变化轨迹图可见，商品销售额既有明显的线性增长趋势，又含有明显的季节波动。 故应构造含有季节虚拟变量的线性趋势模型。为避免虚拟变量陷阱”，只对第二、三、四季度设置虚拟变量，记为D2、D3和D4。EViews回归结果如下： 表9-4EViews输出结果Dependent Variable: Y Method: Least Squares Sample: 2006Q1 2011Q4 Included observations: 24Coefficie Std. Error t-Statistic

14、 Prob.C975.974134.3654928.399830.0000T103.21671.95921652.682630.0000D2127.3019D3-70.8203D4-156.96737.9063238.0579138.309233.358330-1.860857-4.0973760.00330.07830.0006R-squared0.993293Adjusted R-squared0.991881S.E. of regression65.56791Sum squared resid81683.87-131.645Mean dependent var2241.061S.D. d

15、ependent var727.6855Akaike info criterion11.38710Schwarz criterion11.63253Hannan-QuinnLog likelihood2 criter.11.45221F-statistic703.4771 Durbin-Watson stat 1.893665Prob(F-statistic) 0.000000故样本回归方程如下:% 975.974 103.217t 127.302D2t 70.820D3t 156.967D4t(9.7) 在使用确定性趋势和季节模型时，我们实际上是用时间变量和季节虚拟变量作为因变量的所有影响因

16、素的代理变量，简化了回归分析的假定，在经济预测方面应用比较广泛。但有时失之于过于简单化，不能对影响时间序列变动的具体社会经济因素进行分析，不利于验证经济理论、分析经济结构和政策评价。另外，这类模型的一个潜在问题是往往存在误差项自相关，这也是需要在应用时加以注意的。第二节静态模型静态模型(Static Model )是时间序列回归分析的基础模型之一。实际上，上一节介绍的确定性趋势模型和季节模型都属于静态模型。在静态模型中，被解释变量的值只取决于各解释变量和随机误差项的当期值。静态模型一般形式如下：ytkxktut9.8)这实际上是将多元线性回归模型应用于时间序列数据。一、静态模型对数据的要求1

17、.如果模型中的自变量xj (j 1,2,., k) 是非随机变量，或尽管是随机变量，但严格外生于随机误差项ut ，即E(ut |X )=E(ut |X1 ,X 2,.Xk )=0 ， 其中 X j ( , xj(t 2) ,xj(t 1) , xjt , xj(t 1) , xj (t 2) , )则静态模型( 9.8) 可以完全按照横截面数据的经典线性回归模型方法进行参数估计和统计推断，如本章第一节介绍的确定性趋势(季节)回归模型。2 . 如果模型中的自变量xj(j 1,2,., k) 是随机变量，但不严格外生于随机误差项，而是(宽)外生于随机误差项，即仅有E(ut|Xt)=E(ut|x1t

18、,x2t,.xkt)=0 ，这时时间序列的特殊性将给回归分析带来一些问题。但如果xjt , yt 都是产生自平稳、遍历过程，应用OLS 在大样本下仍然可以获得参数的一致和渐进正态的估计量，因此也可以按照横截面数据回归分析的方法进行参数估计和统计推断。 对于非平稳和高度持久的时间序列，由于OLSE 的一致性和渐进正态性受到破坏，容易产生伪回归问题。一般而言，如果生成时间序列的随机过程是非平稳的（或高度持久的），必须将数据平稳化（或弱相依化），然后再进行 OLS 回归。由于在一般情况下，大多数非平稳时间序列都是高度持久的，所以时间序列平稳化的过程同时也是弱相依化的过程。二、趋势平稳序列的回归什么是

19、趋势平稳序列？时间序列中确定性趋势导致的伪回归称为第一种类型的伪回归。如果作 为时间序列中包含了确定性的时间趋势，为了避免第一种类型的伪回归，有 等价的两种处理方法。（一）在模型中引入时间变量如果x、y中含有确定性趋势，若直接将y对x进行回归，相当于遗漏了 重要变量t,所以，应该在模型中将其作为独立的自变量引入进来。这样 x 的回归系数才反映在时间t固定（保持不变）条件下，x对y的 纯净”影响， 避免伪回归。例9-3表9-5中，x是某地人均寿命（岁）,y是该地区稻米产量（千 克/亩），t是时间变量（以1991年作为起始年份）。表9-5 某地人均寿命与稻米产量历史数据年份 1tl x 1yli

20、年份 1tl x | y1991177.05517.1020011177.50630.781992273.00525.8320021277.72610.611993369.55534.2020031378.98634.811994476.78550.9220041478.48651.921995568.60548.7920051579.23661.271996676.05561.5620061676.00642.991997774.72571.4620071780.25613.981998878.55582.4020081874.79619.941999975.52600.3220091976

21、.00678.3320001076.51603.1820102081.28627.13如果直接将y对x进行回归，结果如下:夕 29.034 8.2199xt(9.9)2t (16.13)(5.23)R2 0.305回归系数高度显著，但在经济上无法解释。考虑到 x和y都是t的函数,x?t 70.145 0.302t(9.10)2t (31.46) (2.881)R2 0.316和yt 450.311 7.223t(9.11)2t (26.47) (9.043)R2 0.819一个自然而然的猜测是，正是由于两个变量存在的长期趋势形成了( 9.9)的结果。将t引入(9.9),观察x对y的 纯净”效应

22、(偏效应)，得回归模型 如下：y?t 382.37 0.969xt 6.931t(9.12)t (2.943) (0.527) (7.033)R2 0.822可见，消除了t 的因素， x 对 y 的偏效应在统计上不显著。说明，( 9.9)确实是一个伪回归。(二)去除原序列中的确定性趋势剔除”掉时间序列中的确定性趋势，使时间序列平稳化，然后再应用OLS 进行回归分析。例9-4在例9-3中，我们已经验证出x和y序列中都含有显著的线性 趋势，我们将其从原序列 剔除”出去，即作变换x' xt 5? xt (70.145 0.302t)y' yt yt yt (450.311 7.223

23、t)变换后的数据见表9-6。表9-6变换后的时间序列数据年份x'y'年份x'y'19916.63359.56920014.063101.01319922.28161.06920023.98173.6211993-1.47162.22120034.94090.60219945.45771.72020044.138100.4911995-3.02562.36920054.586102.61819964.12367.90720061.05377.10719972.49170.59020075.00140.88119986.01974.3022008-0.76139.6

24、1919992.68785.00520090.14790.77720003.37680.64020105.12532.354对变换后的序列进行 ADF检验，可以发现，xI(0), y'1(0)。使用OLS,回归结果如下：y' 69.275 0.968x'(9.13t(0.543)R2 0.016可见，二者的线性关系不显著，说明人均寿命与稻米亩产量之间没有实质性的系统联系。去势平稳化变换对回归系数的影响与在原序列的回归模型中加入时间趋势 t 的效果完全一样。二者反映的都是剔除了t 的影响后，x 对 y 的 “ 纯净效应。三、差分平稳序列回归什么是差分平稳序列？两个差分平稳

25、序列的回归往往产生第二种类型的伪回归。我们可以通过差分处理“剔除 ”其中的随机趋势使时间序列平稳化，然后再进行回归。例9-5纽约证券交易所（NYSE）证券市场综合指数y和斯里兰卡人口数x的数据如表9-7所示。试分析两个序列之间的经济关系。表9-7 NYSE证券市场综合指数和斯里兰卡人口数年份y (%)x（千人）年份y (%)x（千人）196643.711440197657.913717196753.811702197752.513942196858.911992197853.614184196951.512252197961.914471197050.212516198077.91473819

26、7154.612608198171.114988197264.512816198281.115189197351.813091198395.215417197436.113284198496.715599197547.6134961985121.515837节选自多米尼克 萨尔瓦多、德里克 瑞杰统计学与计量经济学（第 版）（复旦大学出版社，2008）第291页，表11.10.如果直接将y对x进行OLS回归，结果如下：(9.14)y?t313.01 0.03xt2 t(8.72)R2 0.75将两个序列进行一阶差分变换，分别记为xt 和yt ，重新进行OLS 回归，回归方程如下：?yt 291.80 0.004xt(9.15)2t( 0.93)R2 0.03可见，消除了随机趋势，x 对 y 的效应在统计上不显著。说明，( 9.14)确实是一个伪回归。在单整序列之间进行回归分析时，我们常常遇到不同序列的单整阶数不同，比如我们研究y 与x、z 之间的关系，但yt I(1)，xt I(2)，zt I (3) 。这就需要对不同序列分别进行不同次数的差分使之平稳化，再利用差分变量yt、2Xt和3Zt进行回归。例9-6讨论EViews6.0软件自带的示例文件 demo.wfl中GDP与RS 的关系。打开 EV