计量经济学--时间序列数据分析.doc

时间序列数据的计量分析方法1.时间序列平稳性问题及处理方案1.1序列平稳性的定义从平稳时间序列中任取一个随机变量集，并把这个序列向前移动h个时期，那么其联合概率分布仍然保持不变。平稳时间序列要求所有序列间任何相邻两项之间的相关关系有相同的性质。1.2不平稳序列的后果可能两个变量本身不存在关系而仅仅因为有相似的时间趋势而得出它有关系，也就是出现伪回归；破坏回归分析的假设条件，使得回归结果和各种检验结果不可信。1.3平稳性检验方法：ADF检验1.3.1ADF检验的假设：辅助回归方程：（是否有截距和时间趋势项在做检验时要做选择）原假设：H0：p=0,存在单位根备择假设：H1：P0，不存在单位根结果识别方法：ADF Test Statistic 值小于显著性水平的临界值，或者P值小于显著性水平则拒绝原假设并得出结论：所检测序列不存在单位根，即序列是平稳序列。1.3.2实例对1978年2008年的中国GDP数据进行ADF检验，结果如表一。表一ADF检验结果Augmented Dickey-Fuller test statistict-Statistic Prob.*3.0636211Test critical values:1% level-3.6998715% level-2.97626310% level-2.62742从结果可以看出，ADF的t统计量值大于10%显著性水平上的临界值，P值为，接受原假设，说明所检测的GDP数据是不平稳序列。1.4不平稳序列的处理方法1.4.1方法如果所要分析的数据是不平稳序列，可以对序列进行差分使其变成平稳序列，但是这样做的后果是使新得出的数据丧失了许多原序列的特征，我们能从数据中得到的信息会变少，通常差分的次数不能超过两次。经验表明，存量数据是二阶单整，做二次差分可以使其平稳，流量数据是一阶单整，做一次差分可以使其平稳，增量数据通常就是平稳序列。1.4.2实例将上述GDP数据作二阶差分，然后进行ADF检验，结果见表二。表二二阶差分数据ADF检验结果Augmented Dickey-Fuller test statistict-Statistic Prob.*-2.7347750.0809Test critical values:1% level-3.6891945% level-2.97185310% level-2.625121从结果中可以看出，ADF的t统计量小于10%显著性水平的临界值，P值为0.08，如果选择10%的显著性水平则可以认为序列是平稳序列。说明：在单位根检验的EVIEWS软件界面上有多种检验平稳性的方法，如果选择ADF检验，还需要选择要检验的序列的形式，原序列、一阶差分序列或是二阶差分序列；还可以选择辅助回归函数的形式：只有截距项；有截距项和时间趋势项；没有截距项也没有时间趋势项。如果对回归结果的残差进行ADF检验（不包含时间趋势和截距）可以确定回归方程的变量是否存在协整关系，或者确定模型的设定是否合理。单位根检验是为数不多的对时间序列数据和残差都适用的检验方法，其余的各种检验都是基于回归结果的残差进行的，虽然自相关性检验对变量可以检验，但通常没有什么意义。2.异方差性问题及其解决方案2.1 异方差的定义对于不同的样本点，方差不再是一个固定的常数，而有些是互不相同的，则称模型存在异方差2.2 异方差产生的原因模型中缺少某些解释变量；测量误差。2.3 异方差产生的后果参数诂计量非有效；统计量失真；模型的预测失效。2.4 异方差的检验，2.4.1检验方法图示检验法；Goldfeld-Quandt(戈德菲尔德-匡特)检验；Glejser（戈里瑟）检验；Park(帕克)检验；White（怀特）检验；Remark：说明：怀特检验和戈德菲尔德-匡特检验要求样本数量比较大；图示检验和怀特检验只能检验出是否存在异方差，但不能给出异方差与解释变量的函数关系；帕克检验和戈里瑟检验可以给出异方差如何依赖具体的解释变量和具体的函数形式。2.4.2 G-Q检验a.检验假设H0:1 =2,H1:1 2；结果识别方法：F远大于为递增方差，F接近于为同方差；若计算得到的F值小于所设显著性水平的临界值则认为总体不存在异方差；F值大于所设显著性水平则认为存在递增异方差。b.检验统计量FRSS2/RSS1C.步骤将解释变量按升序或降序排列，去除中间的四分之一，对小数值的四分之一点五回归得SSR，对大数量的四分之一点五回归得SSR2。获得F统计量临界值用命令scalar f=qfdist(0.95,10,10)（5%的显著性水平，两样本的容量都为12）。说明：G-Q检验要求样本数量较大，并且不能知道异方差与解释变量的具体函数形式。只对解释变量排序，排序功能在工作文件的Procs里。2.4.3White检验的假设回归方程：辅助回归方程：（是否包括交叉项在做检验时要做选择）原假设：H0：不存在异方差（辅助回归方程除截距项外其余系数均为零）备择假设：H1：存在异方差结果识别方法：White检验提供了三个统计量，三个统计量对应的P值大于所设定的显著性水平则说明不存在异方差，否则存在异方差。但是同时也要考虑辅助回归方程的有效性，需要通过回归结果的F统计量结果来帮助判断检验结果是否可信。White检验的各个变量之间往往呈现多重共线性，因此我们关注的是方程的总体上是否显著（F统计量），而不关心某个单个变量是否显著（t统计量）。White检验能回答是否存在异方差，但不能给出异方差的函数形式White检验可以粗略认为通过所有解释变量的二阶泰勒级数展开逼近任何非线性函数，如果在White检验情况下没有发现异方差，方程存在异方差的可能性一般会很小，而Park检验和戈里瑟检验则不一定2.4.2实例用模型，进行OLS回归，对残差进行White检验，检验结果如表三。表三White检验结果Heteroskedasticity Test: WhiteF-statistic33.72139 Prob. F(1,29)0.0000 Obs*R-squared16.66677 Prob. Chi-Square(1)0.0000 Scaled explained SS9.277265 Prob. Chi-Square(1)0.0023 VariableCoefficientStd. Errort-StatisticProb. C37350022132087342.8276760.0084 XF20.0232610.0040065.8070120.0000 表三的上部分是White检验结果，共列出了三种统计量，用F统计量表示回归方程系数同时为零的概率。F统讲量的P值很小，说明系数不太可能同时为0，则拒绝原假设，接受存在异方差的观点。表三的下部分是辅助回归方程的各种统计统数，其中各系数的t统计量很小，说明各系数不太可能为零，这与检验结果是一致的。所以检验结果可信。2.4.3 Glejser戈里瑟检验a.辅助回归方程b.方法第一步，产生残差的绝对值序列第二步，做最小二乘回归，尝试不同的函数形式，这依赖于散点图第三步，用t统计量检验变量的显著性来判断是否存在异方差结果识别方法：通过回归结果相关系数的t统计量来判断是否存在异方差，并可以知道方差存在的具体形式。2.4.4 Park检验a.辅助回归方程b.检验假设H 0: =0结果识别方法：回归结果的相关系数显著为0，则说明不存在异方差。说明：Park检验不存在异方差并不一定真正不存在异方差。C方法第一步，获取回归方程的残差第二步，以残差平方的对数为因变量，以比例因子为唯一的解释变量，建立辅助回归方程第三步，采用t检验看是否存在异方差说明：选择比例因子是Park检验的关键，也就是说可能有多个解释变量但只选其中的一个。怀特检验是广义的Park检验可以构造权序列对原始方程进行加权，以为权数，但这种方法不一定能完全消除异方差，原因在于异方差可能还有其他解释变量有关2.5异方差的解决方案2.5.1 对数变换法 对各时间序列数据取对数，缩小解释变量的尺度从而缩小异方差的倍数。2.5.2加权最小二乘法对模型进行加权，使其成为一个不存在异方差性的模型，然后用普通最小二乘法进行估计。通常采用残差绝对值的倒数作为权数。说明：进行加权最小二乘回归时首先进行普通最小二乘回归生成残差序列，然后新建残差变量，可命名为e1，然后最进加权最小二乘回归。在EVIEWS软件的残差检验界面中提供了多种检验方法，如果选择了WHITE检验则需要选择辅助方程中是否包括交叉项，选择交叉项会由于解释变量太多而减少自由度，因此在多元回归中建议不选择包括交叉项。3 序列相关问题及其解决方案3.1 序列相关的定义模型中的随机于扰项不相互独立，也就是模型回归后的残差序列存在着自相关。3.2序列相关产生的原因及后果原因：经济现象固有的惯性；模型设定偏误（漏掉了重要的解释变量或是模型数形式误设，如二次关系误设成只是一次关系）；数据处理的影响。后果：参数估计无偏但估计量无效；预测功能失效。3.3序列相关的检验3.3.1D-W检验这种检验的结果在一般的回归结果中都直接给出，识别方法是，如果其值为2或接近于2则说明不存在序列相关。D-W检验存在的问题：只能检难残差的一阶自相关；模型中要有截距项；模型右边不能含滞后项；无数据缺失。3.3.2 相关图、Q统计量检验 主要用于确定ARMA模型的具体形式及滞后阶数。3.3.3BG(L-M)检验3.3.3.1L-M检验的假设辅助回归方程：原假设：H0：（直到阶滞后不存在序列相关，为预先定义好的整数。）备择假设：H1：存在（存在自相关）识别方法：L-M检验结果的F统计量的P值大于设定的显著性小平则接受原假设，接受不存在序列相关；从一阶开始不断调整检验阶数，并通过辅助回归函数系数的t统计量确定存在几阶自相关。3.3.3.2 实例用模型，进行OLS回归，对残差进行L-M检验，检验结果如表四。表四序列相关检验结果Breusch-Godfrey Serial Correlation LM Test:F-statistic119.3309 Prob. F(3,26)0.0000 Obs*R-squared28.901 Prob. Chi-Square(3)0.0000 VariableCoefficientStd. Errort-StatisticProb. C-78.72442757.602-0.1039130.9180 XF0.0089150.0244260.3649630.7181 RESID(-1)1.4971250.1952767.6667010.0000 RESID(-2)-0.4397630.349507-1.2582380.2195 RESID(-3)-0.1517190.249972-0.6069430.5492 表四上部所示，F统计量的值近似为零，说明辅助回函数的系数在1%的显著平拒绝原假设，接受存在序列相关的假设。表下部是辅助回归函数的相关统计量，其中相关系数的t统计值可帮助确定存在几阶滞后相关。表四中确定不存在三阶滞后相关，存在一阶滞后相关，对于是否存在二阶滞后相关要将LM检验设定为二阶来做判断。实际上自相关性检验是从一阶开始，而不是直接检难三阶自相关。3.4序列相关的处理方案：最常用的方法是广义最小二乘法（GLS:Generalized least squares）一阶差分法（First-Order Difference)广义差分法(Generalized Difference)( ARMA模型)修正标准误方法（HAC）(可以在OLS回归中直接实现)3.4.1用自相关和偏自相关系数来识别AR或者是MA模型MA(q)的识别：MA(q)的自相关系数的q步以后是截尾的；MA(q)的偏自相关系数呈现出某种衰减的形式是拖尾的；可以通过识别一个序列的偏自相关系数的拖特征，大致确定模型是否服从MA(q)过程。AR(p)的识别：AR(p)模型的自相关系数会呈现出某种形式的衰减；可以通过识别模型的偏自相关系数的个数，来确定AR(p)模型的阶数p。（也可以通过模型的LM检验来确定滞后的阶数）3.4.2实例用模型，进行OLS回归，对残差进行自相关性检验。结果显示，残差的自相关性呈现拖尾，偏自相关系数呈现一阶截尾。我们可以初步判断模型应为AR(1)模型。但是通过残差的LM检难显示为二阶序列相关，因此我们将模型定义为AR(2)。修正后的回归结果如表五。表五修正序列相关的AR模型的回归结果Method: Least SquaresVariableCoefficientStd. Errort-StatisticProb. C-2633.2890 6911.2710 -0.3810 0.7064 XF2.4631 0.2272 10.8403 0.0000 AR(1)1.9346 0.2224 8.6971 0.0000 AR(2)-1.0231 0.2535 -4.0357 0.0005 R-squared0.999511F-statistic17033.39Adjusted R-squared0.999452Prob(F-statistic)0.0000 Durbin-Watson stat2.153095上表中的相关系数都通过了t检验，并且DW值为2.15，可以认为模型不存在一阶自相关。实际上，此时对残差进行LM检验应就不存在序列相关。说明：标准误经过HAC修正后，DW统计量失效；与HC修正后，White检验失效具有异曲同工之处。4. 多重共线性问题及解决方案4.1 多重共线性的定义如果两个或多个解释变量之间出现了较强的近拟相关性并且是线性相关性则称为多重共线性。4.2 多重共线性出现的原因经济变量之间内在着内在联系；经济变量在时间上有相关的共同趋势;解释变量中含有滞后项。4.3 多重共线性的影响参数估计值不稳定，方差很大；回归系数的符号有误，经济含义不合理；变量的显著性检验失去意义。4.4 多重共线性的诊断R方值很大，t值较小（P很大）；理论性强，检验值弱；新引入变量后方差增大；4.5 多重共线性的解决方案先验信息法：根据经济理论或其它研究提前确定模型参数之的关系；用相对数代替决定数；删去模型中次要的或可替代的解释变量；差分法；增大样本容量；逐步回归法。4时间序列的计量经济分析4.1多元回归分析4.1.1数据来源本文从中国统计年鉴2009 及2001中获得了从1978年至2008年的中国宏观经济数据，包括消费、投资、政府购买、进口、出口，并以此为根据计算出了国民生产总值。4.1.2分析模型a.为了分析消费、投资与政府购买对国内生产总值的影响，构构建了如下模型：b.将上述模型中的数据进行单位根检验，结果显示所有数据均为二阶单整。此刻应该对数据进行处理使其变为平稳序列，通常的方法是做差分，做二阶差分可以做数据变为平稳序列，但是这样数据会失去对于经济现象的解释能力。在这里暂不做处理。c.模型回归后结果进行自相关性检验，发现残差存在着二阶滞后相关，进一步检验显示残差的自相关性呈现一阶截尾，偏自相关系数呈现二阶尾。因此为修正模型存在的自相关性，用ARMA(2,1)模型。d.进行线性回归。回归结果见表六。e.对上述回归结果进行异方差性检验、序列相关性检验并对残差进行单位根检验。结果见表七、表八和表九。表六ARMA(2,1)模型回归结果Method: Least SquaresVariableCoefficientStd. Errort-StatisticProb. C1.1900580.08298914.339870.0000 LOGXF0.3170940.0433497.3149910.0000 LOGI0.4142090.03570911.599550.0000 LOGG0.2671790.0658034.060270.0005 AR(1)1.3855920.11075812.510040.0000 AR(2)-0.932420.117263-7.9515310.0000 MA(1)-0.9999130.119961-8.3353150.0000 R-squared0.999926Adjusted R-squared0.999906F-statistic49709.31 Durbin-Watson stat2.337958Prob(F-statistic)0.0000 结果显示，模型中各变量的系数通过了t检验，但是拟合较高，说明模型中变量存在线性相关。表七ARMA(2,1)模型回归的异方差检验结果Heteroskedasticity Test: WhiteF-statistic1.44882 Prob. F(7,21)0.2389Obs*R-squared9.444254 Prob. Chi-Square(7)0.2223Scaled explained SS3.676272 Prob. Chi-Square(7)0.8162VariableCoefficientStd. Errort-StatisticProb. C0.0244830.0358030.6838230.5016GRADF_0124.10E-054.19E-050.9766260.3399GRADF_022-2.44E-062.43E-06-1.000310.3285GRADF_032-2.85E-063.02E-06-0.946170.3548GRADF_0425.97E-063.95E-061.5105420.1458GRADF_0520.0067830.0261550.2593210.7979GRADF_0620.0163320.0147811.1049390.2817GRADF_072-3.84E-055.61E-05-0.683480.5018R-squared0.325664Durbin-Watson stat2.000288Adjusted R-squared0.100885Prob(F-statistic)0.238901表七分为三个部分，上部分为异方差检验的三种统计量，结果大于10%的显著性水平，因此接受原假设不存在异方差。表的中部为辅助回归方程的系数，t统计量的p值大于10%的显著性水平，接受回归系数均为零的原假设，这些统计值可以帮助判断White检验统计量的可信度。表的下部为辅助回归方程回归结果的统计量，可以判断辅助回归结果的有效性，并帮助判断White检验统计量的可信度。通常辅助加归方程有效性高，White检验统计量的可信度才高。表八ARMA(2,1)模型回归的自相关性检验结果Breusch-Godfrey Serial Correlation LM TestF-statistic1.664913 Prob. F(1,21)0.211Obs*R-squared2.050869 Prob. Chi-Square(1)0.1521VariableCoefficientStd. Errort-StatisticProb. C-0.039990.08959-0.4463560.6599LOGXF0.0239110.047460.5038150.6196LOGI-0.005310.039672-0.1337520.8949LOGG-0.016750.070903-0.2362320.8155AR(1)0.029480.1180020.249830.8051AR(2)-0.03390.121122-0.2798710.7823MA(1)-0.000280.000379-0.7408920.467RESID(-1)-0.259050.259587-0.9979220.3297R-squared0.07072Durbin-Watson stat2.110225Adjusted R-squared-0.23904Prob(F-statistic)0.973851如上表所示，由于经过进行了自相关性的修正，因此只检验了回归结果的一阶滞后相关。上表分为三部分，上部分为LM检验的三个统计值，结果的p值大于10%的显著性水平，因此接受在滞后p阶内不存在自相关的原假设，表的中部和下部分别为辅助回归方程的系数的相关统计量，可以判断辅助回归方程回归结果的有效性，进而帮助判断LM检验结果的可信度。表九ARMA(2,1)模型回归残差的单位根检验结果Augmented Dickey-Fuller test statistict-Statistic Prob.*Augmented Dickey-Fuller test statistic-6.2662030.0000 Test critical values:1% level-2.6501455% level-1.95338110% level-1.609798VariableCoefficientStd. Errort-StatisticProb. E11(-1)-1.1770090.187834-6.2662030.0000 R-squared0.592447Durbin-Watson stat2.061427Adjusted R-squared0.592447上表第一部分是ADF检验结果，结果P值为0，说明原假设存在单位根不能接受，所检验的残差序列不存在单位根。表的中部和下部是辅助回归方程的相关统计量，协助判断ADF统计结果的可信度。综上所述，所建立的ARMA(2,1)模型通过了异方差、序列相关检验，并且残差序列平稳，虽然模型中各变量是I(2)序列，但是回归结果依然是有效的。因为如果经济变