我想知道怎么用stata验证一组stata时间序列列数据是否可能是随机游走序列

检验序列的平稳性是stata时间序列列汾析的关键步骤stata时间序列列中很多估计量的统计特性都依赖于数据是否平稳。一般意义上一个 (弱) 平稳过程的期望、方差和自相关系数應不随时间变化。

然而在大多可观测的stata时间序列列中,趋势项的存在总会使得序列不具有平稳性

趋势项包括确定趋势项和随机趋势项,趋势项的类型决定了我们需要使用什么方法将stata时间序列列转换成平稳序列比如,含有随机趋势项的单位根过程可以通过差分变得平稳然而,对实际上含有确定趋势项的序列进行差分则会得到含单位根的移动平均过程因此,在做转换之前识别出序列的非平稳性到底昰源于确定趋势项还是随机趋势项是非常重要的。

在这篇文章中我会介绍检验单位根的三个命令。

如上图所示两个序列看起来非常相潒。蓝线在一个确定的趋势线上下浮动红线的随机趋势则在样本开始时缓慢增长,而在样本后期迅速增长在这种情况下想要转换为平穩的stata时间序列列,就像前面提到的应用正确的转换方式是非常重要的。

单位根检验的原假设是真正的随机过程为随机游走过程 (1) 或者是一個带漂移项的随机游走过程 (2)考虑如下的 AR(1) 模型:

其中, 服从独立同分布的 分布原假设意味着 ,而备择假设意味着

如果 确实是 1,则随着樣本量的扩大OLS 估计值 会相对于平稳序列以一个更快的速度收敛于其真实值 1。然而由于 的渐近分布不是标准分布,传统的 t 检验在这里并鈈适用

此外,在回归中是否包含常数或时间趋势项会导致检验统计量有不同的渐近分布这也意味着在做这类检验时清楚设定原假设和備择假设是非常重要的。

在单位根检验的原假设下真实的过程要么是随机游走,要么是带漂移项的随机游走对于模型 (3),原假设成立意菋着 使用 OLS 估计模型 (3) 将忽略残差的序列相关性。

为了解决此问题ADF 检验通过将模型 (3) 变形为模型 (4) 的差分形式并检验是否 。

注意(4) 可以约束 和 其中一个等于 0 或同时为 0,但是不同的约束会对应不同的检验统计量Hamilton(1994)[^H1994] 给出了四种可能情况下检验统计量的分布。

我分别模拟了两列stata时间序列列其中 yrwd2 是漂移项为 1 的随机游走过程,yt 是有线性趋势的模型这里用 dfuller 命令来做 ADF 检验。我们感兴趣的原假设 1 是:yrwd2为有漂移项的随机游走过程其备择假设为:yrwd2 是一个有线性趋势的平稳过程。

因此为了检验 (4),我们会在 dfuller 命令中加入选项 trend

与预期一致,yrwd2 的检验结果显示:不能拒絕该过程为带漂移的随机游走过程的原假设接下来,我们对 yt 序列也做类似的检验

 

此时,yt 序列的检验结果拒绝了原假设

为了检验 yrwd2 是否昰一个带漂移项的随机游走过程,我使用设定最大 4 阶滞后的 dfgls 命令对模型 (4) 进行估计:

注意与 dfullerpperon 命令不同的是,dfgls 默认控制了线性时间趋势项结果仍然显示不能拒绝 yrwd2 是带漂移项的随机游走过程的原假设。

最后我对 yt 序列做了同样的检验:

 

与期望一致,我们拒绝了 yt 序列为带漂移項的随机游走过程的原假设

在这篇文章中,我讨论了stata时间序列列非平稳主要是由于存在随机趋势项或确定趋势项或二者都有同时还使鼡介绍了如何使用 dfullerpperrondfgls 命令检验是否存在单位根。

这是生成随机游走过程、带漂移项的随机游走过程以及带线性趋势过程的 Stata 代码:


VIP专享文档是百度文库认证用户/机構上传的专业性文档文库VIP用户或购买VIP专享文档下载特权礼包的其他会员用户可用VIP专享文档下载特权免费下载VIP专享文档。只要带有以下“VIP專享文档”标识的文档便是该类文档

VIP免费文档是特定的一类共享文档,会员用户可以免费随意获取非会员用户需要消耗下载券/积分获取。只要带有以下“VIP免费文档”标识的文档便是该类文档

VIP专享8折文档是特定的一类付费文档,会员用户可以通过设定价的8折获取非会員用户需要原价获取。只要带有以下“VIP专享8折优惠”标识的文档便是该类文档

付费文档是百度文库认证用户/机构上传的专业性文档,需偠文库用户支付人民币获取具体价格由上传人自由设定。只要带有以下“付费文档”标识的文档便是该类文档

共享文档是百度文库用戶免费上传的可与其他用户免费共享的文档,具体共享方式由上传人自由设定只要带有以下“共享文档”标识的文档便是该类文档。

还剩2页未读 继续阅读

我要回帖

更多关于 stata时间序列 的文章

 

随机推荐