与时间序列相关的STATA命令及其统计量的解析.docx

与时间序列相关的STATA命令及其统计量的解析.docx

《与时间序列相关的STATA命令及其统计量的解析.docx》由会员分享，可在线阅读，更多相关《与时间序列相关的STATA命令及其统计量的解析.docx（10页珍藏版）》请在冰豆网上搜索。

与时间序列相关的STATA命令及其统计量的解析

与时间序列相关的STATA命令及其统计量的解析残差U序列相关：

1DW统计量针对一阶自相关的（高阶无效）

STATA命令：

1.先回归

2.直接输入dwstat统计量如何看：

查表

2Q统计量针对高阶自相关correlogram-Q-statistics

STATA命令：

1．先回归reg

2．取出残差predictu,residual（不要忘记逗号）

3．wntestquQ

统计量如何看：

p值越小（越接近0）Q值越大——表示存在自相关具体自相关的阶数可以看自相关系数图和偏相关系数图：

STATA命令：

自相关系数图:

acu（残差）或者窗口操作在Graphics——Time-seriesgraphscorrelogram（ac）

偏相关系数图:

pacu或者窗口操作在Graphics——Time-seriesgraphs——（pac）自相关与偏相关系数以及Q统计量同时表示出来的方法：

corrgramu或者是窗口操作在

autocorrelations

3LM统计量一一针对高阶自相关

STATA命令：

1．先回归reg

2．直接输入命令estatebgodfrey,lags（n）或者窗口操作

在Statistics——Postestimation（倒数第二个）——ReportsandStatistics（倒数第二个）——在里面选择Breush-GodfreyLM（当然你在里面还可以找到方差膨胀因子还有DW统计量等常规统计量）

LM统计量如何看：

P值越小（越接近0）表示越显着（显着拒绝原假设），存在序列相关具体是几阶序列相关，你可以把滞后期写为几，当然默认是1，（通常的方法是先看图，上面说的自相关和偏相关图以及Q值，然后再利用LM肯定）。

平稳时间序列存在自相关的问题的解决方案残差出现序列相关的补救措施：

1、一阶自相关:

最近简单的方法是用AR

（1）模型补救，就是在加一个残差的滞后项即可。

2、高阶的自相关:

用AR（n）模型补救。

AR模型的识别与最高阶数的确定：

可通过自相关系数来获得一些有关AR（p）模型的信息，如低阶AR（p）

模型系数符号的信息。

但是，对于自回归过程AR（p），自相关系数并不能帮助我们确定AR（p）模型的阶数p。

所以，可以考虑使用偏自相关系数k,k，以便更加全面的描述自相关过程AR（p）的统计特征。

且对于一个AR（p）模型，k,k的最高阶数为p,也即AR（p）模型的偏自相关系数是p阶截尾的。

因此，可以通过识别AR（p）模型的偏自相关系数的个数,来确定AR（p）模型的阶数p,进而设定正确的模型形式,并通过具体的估计方法估计出AR（p）模型的参数。

如果AR（p）还解决不了则进一步使用：

MA（q）模型，以及ARMA（p,q）模型。

1、MA（q）

MA（q）的偏自相关系数的具体形式随着q的增加变得越来越复杂,很难给出一个关于q的一般表达式，但是，一个MA（q）模型对应于一个ARQ）模型。

因此,MA（q）模型的偏自相关系数一定呈现出某种衰减的形式是拖尾的。

故可以通过识别一个序列的偏自相关系数的拖尾形式,大致确定它应该服从一个MA（q）过程。

2、ARMA（p,q）就是既含有AR项又含有MA项。

我们引入了自相关系数和偏自相关系数这两个统计量来识别ARMA（p,q）模型的系数特点和模型的阶数。

但是,在实际操作中,自相关系数和偏自相关系数是通过要识别序列的样本数据估计出来的,并且随着抽样的不同而不同,其估计值只能同理论上的大致趋势保持一致,并不能精确的相同。

因此,在实际的模型识别中,自相关系数和偏自相关系数只能作为模型识别过程中的一个参考,并不能通过它们准确的识别模型的具体形式。

具体的模型形式,还要通过自相关和偏自相关系数给出的信息,经过反复的试验及检验,最终挑选出各项统计指标均符合要求的模型形式。

注：

无论采取什么样的方式，只要能够把残差中的序列相关消除掉，又不会引入新的问题，这样的模型就是最优模型。

与平稳性检验及其统计量解析（P212张晓峒）

白噪声检验：

1.Q检验wntestqvar，lag（n）

2.Bartlett检验wntestbvar，table（表示结果以列显示，而不做图。

不加table就以图形的方式现实）

或者在Statistics——Time-series——TEST——Bartlett检验（第四个）

画密度图：

1、概率密度图

命令：

pergramvar,generate（新变量名字）将概率密度的图上所生成的值

生成并储存在新变量里，这个不是必须的，只是为了日后方便。

窗口：

Statistics——Time-series——Graphs——Periodogram（第五个）

2．累积分布函数图

命令：

cumspvar,generate（新变量名字）解释同上，并且这个生成新变量的功能似乎只能通过命令完成。

窗口：

Statistics——Time-series——Graphs——CumulativeSpectral

distribution

单位根检验（219）

1、Dickey-Fuller检验

命令：

dfullervar（,lags（#）/trend/noconstant/regress/）

对变量做ADF检验可以加滞后期或趋势项或不含常数项等等这些取决于你

的模型窗口：

Statistics——Time-series——TEST——ADF单位根检验（第一个）在里面你也可以选择滞后期数，常数项等等。

如何看结果：

原假设为：

至少存在一个单位根；备选假设为：

序列不存在单位根。

如果统计量小于后面的显着性水平给出的值且P值很大——有单位；如果统计量大于后面的显着性水平给出的值且P值很小——无单位根ADF检验需要注意的地方：

（1）必须为回归定义合理的滞后阶数，通常采用AIC准则来确定给定时间序列模型的滞后阶数。

在实际应用中，还需要兼顾其他的因素，如系统的稳定性、模型的拟合优度等。

（2）可以选择常数和线性时间趋势，选择哪种形式很重要，因为检验显着性水平的t统计量在原假设下的渐进分布依赖于关于这些项的定义。

1如果在检验回归中含有常数，意味着所检验的序列的均值不为0，一个

简单易行的办法是画出检验序列的曲线图，通过图形观察原序列是否在一个偏离0的位臵随机变动，进而决定是否在检验时添加常数项；

2如果在检验回归中含线性趋势项，意味着原序列具有时间趋势。

同样，决定是否在检验中添加时间趋势项，也可以通过画出原序列的曲线图来观察。

如果图形中大致显示了被检验序列的波动趋势随时间变化而变化，那么便可以添加时间趋势项。

对变量做PP

2、Phillips-Perron检验

命令：

pperronvar,（,lags（#）/trend/noconstant/regress/）检验可以加滞后期或趋势项或不含常数项等等这些取决于你的模型。

窗口操作：

Statistics——Time-series——TEST——PP单位根检验（第三个）

如何看结果：

同ADF一样原假设为：

至少存在一个单位根；备选假设为：

序列不存在单位根。

P值越小（统计量大于各显着性水平值）——不存在单位根

P值越大（统计量小于各显着性水平值）——存在单位根

向量自相关回归VAR模型

向量自回归（VAR模型是AR模型的多元扩展，用以反映在一个系统中的多个变量之间的动态影像，格兰杰因果检验、脉冲响应、方差分解都是VAR模型中重要的分析工具。

与VAR模型相关的STATA命令与解析

1、VAR模型的估计

STATA命令：

var解释变量（，无常数项noconstant/滞后期lags（n）/外生变量exog（varlist）/constraints（numlist）线性约束的个数

{注意：

使用线性约束要提前定义，详情见建模中的各种小问题}/LIKEPOHL滞后阶数选择的统计量lutstats）

窗口操作：

StatisticsMultivariatetimeseriesVAR第二项）

如何看结果：

保存估计结果的命令：

eststore名称

2．VAR模型平稳性

STATA命令：

varstable（,graph表示画出图形）如何看结果：

特征值都在圆内，即都小于1，表示VAR模型稳定窗口操作：

Statistics——Multivariatetimeseries——VARdiagnostics

andtests——checkstabilityconditionofVARestimates

3．VAR阶数的选择——滞后阶数的确定

在VAR模型中，正确的选择模型的滞后阶数，对于模型的估计和协整检验都产生一定的影响，小样本情况更是如此。

（1）STATA命令：

用于VAR模型估计之前

varsoc解释变量（，没有常数项noconstant/最高滞后期maxlag（#）/外生变量exog（varlist）/线性约束条件constraints（numlist））

（2）命令：

用于模型估计之后

解释变量（，estimates（estname））其中，estname表示已经估计的

VAR模型的名字。

（1）

（2）如何看结果：

找最显着的阶数作为其滞后项（一般会标有※）

（3）命令：

用于模型估计之后（Wald滞后排除约束检验）Varwle窗口操作：

Statistics——Multivariatetimeseries——VARdiagnosticsandtests——第一第二项如何看结果：

看不同阶数上的联合显着性，看P值，越小越显着，表示存在该阶滞后项。

4．残差的正态性与自相关检验

STATA命令:

1.先进行var回归

2.varnorm如何看结果：

原假设是服从正态分布

P值越小越显着拒绝原假设——不服从正态分布

P值越大越不显着拒绝，原假设成立——服从正态分布

自相关：

窗口操作：

Statistics——Multivariatetimeseries——VAR

diagnosticsandtests——LMTest

正态分布：

窗口操作：

Statistics——Multivariatetimeseries——VAR

diagnosticsandtests——Testfornormally（倒数第三项）5．Granger因果关系检验格兰杰因果关系不同于我们平常意义上的因果关系，它是指一个变量对于另外一个变量具有延期影响。

格兰杰因果关系检验有助于表明变量间的动态影响，有助于提高模型的预测效果。

命令格式：

1.先进行var

2.再进行格兰杰因果检验vargranger

如何看结果：

看P值的显着性，越小说明存在越强的因果关系，相反P值越大说明两者的因果关系不明显。

窗口操作：

Stat