与时间序列相关的STATA

1、与时间序列相关的STATA命令及其统计量的解析残差U序列相关：DW统计量一一针对一阶自相关的（高阶无效）STATA命令：先回归直接输入dwstat统计量如何看：查表Q统计量针对高阶自相关correlogram-Q-statisticsSTATA命令：先回归reg取出残差predict u,residual（不要忘记逗号）wntestq u Q统计量如何看：p值越小（越接近0） Q值越大一一表示存在自相关具体自相关的阶数可以看自相关系数图和偏相关系数图：STATA命令：自相关系数图：ac u（ 残差）或者窗口操作在 GraphicsTime-series graphscorrelogram（ac

2、）偏相关系数图：pac u 或者窗 口操作在 GraphicsTime-series graphs（pac）自相关与偏相关系数以及Q统计量同时表示出来的方法：corrgram u或者是窗口操作在StatisticsTime-seriesGraphsAutocorrelations&Partial autocorrelationsLM统计量一一针对高阶自相关STATA命令：先回归reg直接输入命令estate bgodfrey,lags（n）或者窗口操作在 StatisticsPostestimation（倒数第二个） Reports and Statistics（倒数第二个） 在里面选择Bre

3、ush-Godfrey LM （当然你在里面还可以找到方差膨胀因子还有DW统计量等常规 统计量）LM 统计量如何看：P值越小（越接近0）表示越显著（显著拒绝原假设），存在序列相关具体是几阶序列相关，你可以把滞后期写为几，当然默认是1,（通常的方法是先看图，上 面说的自相关和偏相关图以及Q值，然后再利用LM肯定）。平稳时间序列存在自相关的问题的解决方案残差出现序列相关的补救措施：1、一阶自相关：最近简单的方法是用AR（1）模型补救，就是在加一个残差的滞后项即可。2、高阶的自相关：用AR（n）模型补救。AR模型的识别与最高阶数的确定：可通过自相关系数来获得一些有关AR(p)模型的信息，如低阶AR(

4、p)模型系数符号的 信息。但是，对于自回归过程AR(p)，自相关系数并不能帮助我们确定AR(p)模型的阶数p。 所以，可以考虑使用偏自相关系数pk,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的增加变得

5、越来越复杂，很难给出一个关于 q的一般表达式，但是，一个MA(q)模型对应于一个AR(-)模型。因此，MA(q)模型的偏 自相关系数一定呈现出某种衰减的形式是拖尾的。故可以通过识别一个序列的偏自相关系数 的拖尾形式，大致确定它应该服从一个MA(q)过程。2、ARMA(p,q)就是既含有AR项又含有MA项。我们引入了自相关系数和偏自相关系数这两个统计量来识别ARMA(p,q)模型的系数特 点和模型的阶数。但是，在实际操作中，自相关系数和偏自相关系数是通过要识别序列的样 本数据估计出来的，并且随着抽样的不同而不同，其估计值只能同理论上的大致趋势保持一 致，并不能精确的相同。因此，在实际的模型识别中

6、，自相关系数和偏自相关系数只能作为 模型识别过程中的一个参考，并不能通过它们准确的识别模型的具体形式。具体的模型形式， 还要通过自相关和偏自相关系数给出的信息，经过反复的试验及检验，最终挑选出各项统计 指标均符合要求的模型形式。注：无论采取什么样的方式，只要能够把残差中的序列相关消除掉，又不会引入新的问 题，这样的模型就是最优模型。与平稳性检验及其统计量解析(P212张晓炯)白噪声检验：Q 检验 wntestq var， lag(n)Bartlett检验wntestb var ，table (表示结果以列显示，而不做图。不加table就以图形 的方式现实)或者在 StatisticsTime-

7、seriesTESTBartlett 检验(第四个)画密度图：1、概率密度图命令：pergram var ,generate(新变量名字)将概率密度的图上所生成的值生成并储存在新变 量里，这个不是必须的，只是为了日后方便。窗口： StatisticsTime-seriesGraphsPeriodogram(第五个)累积分布函数图命令：cumsp var ,generate(新变量名字)解释同上，并且这个生成新变量的功能似乎只能通 过命令完成。窗口： StatisticsTime-seriesGraphsCumulative Spectral distribution单位根检验(219)1、Di

8、ckey-Fuller 检验命令：dfuller var （,lags（#）/trend/noconstant/regress/）对变量做ADF检验可以加滞后期或趋势项或不含常数项等等这些取决于你的模型。窗口： StatisticsTime-seriesTESTADF单位根检验（第一个）在里面你也可以选择 滞后期数，常数项等等。如何看结果：原假设为：至少存在一个单位根；备选假设为：序列不存在单位根。如果统计量小于后面的显著性水平给出的值且P值很大一一有单位；如果统计量大于后面的显著性水平给出的值且P值很小一一无单位根ADF检验需要注意的地方：（1）必须为回归定义合理的滞后阶数，通常采用AIC准

9、则来确定给定时间序列模型的滞后 阶数。在实际应用中，还需要兼顾其他的因素，如系统的稳定性、模型的拟合优度等。（2）可以选择常数和线性时间趋势，选择哪种形式很重要，因为检验显著性水平的t统计 量在原假设下的渐进分布依赖于关于这些项的定义。如果在检验回归中含有常数，意味着所检验的序列的均值不为0，一个简单易行的办法 是画出检验序列的曲线图，通过图形观察原序列是否在一个偏离0的位谿随机变动，进而 决定是否在检验时添加常数项；如果在检验回归中含线性趋势项，意味着原序列具有时间趋势。同样，决定是否在检验 中添加时间趋势项，也可以通过画出原序列的曲线图来观察。如果图形中大致显示了被检验 序列的波动趋势随时

10、间变化而变化，那么便可以添加时间趋势项。2、Phillips-Perron 检验命令：pperron var , （,lags（#）/trend/noconstant/regress/）对变量做 PP 检验可以加滞后期或趋 势项或不含常数项等等这些取决于你的模型。窗口操作：StatisticsTime-seriesTESTPP 单位根检验（第三个）如何看结果：同ADF 一样原假设为：至少存在一个单位根；备选假设为：序列不存在单位根。P值越小（统计量大于各显著性水平值）一一不存在单位根P值越大（统计量小于各显著性水平值）一一存在单位根向量自相关回归VAR模型向量自回归（VAR）模型是AR模型的多

11、元扩展，用以反映在一个系统中的多个变量之间的 动态影像，格兰杰因果检验、脉冲响应、方差分解都是VAR模型中重要的分析工具。与VAR模型相关的STATA命令与解析1、VAR模型的估计STATA命令：var解释变量 （，无常数项noconstant/滞后期lags（n）/外生变量 exog（varlist）/constraints（numlist 线 性约束的个数注意：使用线性约束要提前定义，详情见建模中的各种小问题/LIKEPOHL滞后阶数选 择的统计量lutstats）窗口操作：StatisticsMultivariate time seriesVAR（第二项）如何看结果：保存估计结果的命令：

12、est store名称VAR模型平稳性STATA命令：varstable(, graph表示画出图形)如何看结果：特征值都在圆内，即都小于1，表示VAR模型稳定窗口操作： StatisticsMultivariate time seriesVAR diagnostics and testscheckstability condition of VAR estimatesVAR阶数的选择一一滞后阶数的确定在VAR模型中，正确的选择模型的滞后阶数，对于模型的估计和协整检验都产生一定的影 响，小样本情况更是如此。STATA命令：用于VAR模型估计之前varsoc解释变量(，没有常数项noconsta

13、nt/最高滞后期maxlag(#)/外生变量exog(varlist)/ 线性约束条件 constraints(numlist)命令：用于模型估计之后解释变量(，estimates(estname)其中，estname表示已经估计的VAR模型的名字。(1)(2)如何看结果：找最显著的阶数作为其滞后项(一般会标有浓) 命令：用于模型估计之后(Wald滞后排除约束检验)Varwle窗口操作：StatisticsMultivariate time seriesVAR diagnostics and tests第一第二项如何看结果：看不同阶数上的联合显著性，看P值，越小越显著，表示存在该阶滞后项。残差

14、的正态性与自相关检验STATA命令：先进行var回归varnorm如何看结果：原假设是服从正态分布P值越小越显著拒绝原假设一一不服从正态分布P值越大越不显著拒绝，原假设成立一一服从正态分布自相关：窗口操作：StatisticsMultivariate time seriesVAR diagnostics and testsLMTest正态分布：窗口操作：StatisticsMultivariate time seriesVAR diagnostics andtestsTest for normally 倒数第三项)Granger因果关系检验格兰杰因果关系不同于我们平常意义上的因果关系，它是指一

15、个变量对于另外一个变量具有 延期影响。格兰杰因果关系检验有助于表明变量间的动态影响，有助于提高模型的预测效果。命令格式：先进行var再进行格兰杰因果检验vargranger如何看结果：看P值的显著性，越小说明存在越强的因果关系，相反P值越大说明两者的 因果关系不明显。窗口操作：StatisticsMultivariate time seriesGranger causality test脉冲响应与方差分解（223）脉冲响应与方差分解是一个问题的两个方面。脉冲响应是衡量模型中的内生变量对一个变量的脉冲（冲击）做出的响应一对多，一个变 量 向下所引起的其他变量的变动，而方差分解则是如何将一个变量的

16、响应分解到模型中的内 生变量多对一，一个变量的变动向上追溯引起该变动的若干原因。STATA 的 irf命令用于计算VAR、SVAR、VEC模型的脉冲响应、动态乘子和方差分解。注意：该方法的操作使用于var、svar、vec估计之后。（1）创建irf文件STATA命令：irf create irfname，set（名字）（先进行var,然后使用这条命令就可以 直接把刚刚var的结果保存到该irf文件里，并且只有这条命令是最好用的，其他 命令即使可以建立irf文件但是不能把var的结果保存进去，那也是没用的。）激活irf文件显示当前处于活动状态的irf文件：STATA 命令：irf set激活（或

17、创建）irf文件：STATA命令：irf set文件名称创建新的irf文件并替换正在活动的irf文件：STATA 命令：irf set 文件名称，replace清除所有活动的irf文件：STATA 命令：irf set，clear窗口操作：StatisticsMultivariate time seriesManage IRF results and files（2）用irf文件作图（223）对于VAR、SVAR、VEC模型，脉冲响应函数（IRF）的类型包括简单脉冲响应、正交脉冲响 应、动态乘子三种，方差分解包括Cholesky分解和结构分解两种。没种模型可以采用不同 的分析工具。窗口操作：S

18、tatisticsMultivariate time seriesIRF and FEVD analysis简单的 IRF：（VAR/SVAR/VEC 之后）命令：irf graph irf （，使用哪个文件set（文件名）/脉冲变量impulse（变量名）/响应变量 response （内生变量名）如果不加约束就是默认当前打开的文件动态乘子：（VAR之后）命令：irf graph dm （，使用哪个文件set（文件名）/脉冲变量 impulse（变量名）/响应变量response （内生变量名）方差分解：（VAR/SVAR/VEC之后）命令：irf graph fevd （，使用哪个文件se

19、t（文件名）/脉冲变量impulse（变量名）/响应变量 response （内生变量名）联合图表：将多个脉冲响应图或方差分解图结合起来）命令：irf cgraph （irfname脉冲变量响应变量方差分解的方法fevd/IRF的方法irf） （irfname脉冲变量响应变量IRF方法irf/方差分解的方法fevd）叠加图表：（将多个脉冲响应图或方差分解图叠加起来）命令：irf ograph（irfname脉冲变量响应变量方差分解的方法fevd/IRF的方法irf） （irfname脉冲变量响应变量IRF方法irf/方差分解的方法fevd)Irf列表STATA命令：irf table IRF方

20、法irf/方差分解方法fevd联合列表：(将多个脉冲响应图或方差分解列表结合起来)命令：irf ctable (irfname脉冲变量响应变量方差分解的方法fevd/IRF的方法irf) (irfname脉冲变量响应变量IRF方法irf/方差分解的方法fevd)Irf其他命令 命令：irf describeirf describe，detail7. VAR模型的预测227窗口操作：StatisticsMultivariate time seriesDynamic forecast命令格式1(对于VAR、SVAR模型)：fcast compute prefix命令格式2 (对于VECM模型)：f

21、cast compute prefix对预测进行作图命令：fcast graph prefixvar(prefix 变量名)小结大概流程：估计VAR模型var y x zest store VAR1根据信息准则确定VAR模型的最优滞后结束，根据结果重新估计varsoc x z , maxlag(#)var *(全部变量，或者ln*所有的对数变量)，lags(1/3)(比如最优的滞后期为3,滞后 期 123)est store VAR2考察VAR模型的平稳性varstable, estimates(VAR2) graph dlabel (画图并标出具体数值)检验VAR模型残差的正态分布特征和自相

22、关特征varnorm, jbera estimates(VAR2)对各变量进行Granger因果关系检验vargranger (，estimates (VAR2)绘制脉冲响应图以及预测误差方差分解var y x z，lags(1/3)irf create irfname， set (名称)irf graph irf (, estimates (名称)irf table fevd(, estimates (名称)/预测区间n8step(n)根据VAR模型的估计结果进行预测预测n期(n8)fcast compute prefix(, step(n)fcast compute f_(, step(n)将VAR模型与IRF相结合的窗口操作：StatisticsMultivariate time seriesBasic VAR约翰逊协整检验协整检验是对非平稳变量进行回归的必要前提。只有存在协整关系，协整回归才有意义。在各种协整检验方法中，Johansen(1998)在 VAR框架下的特征值检验和迹检验应用最为普通。命令格式为：vecrank var1 var2 (,lag(n),trend