格兰杰检验的基本步骤_【Stata教程】格兰杰因果检验-程序员宅基地

原标题：【Stata教程】格兰杰因果检验

“社会科学中的数据可视化”第432篇推送

引言

在实证分析中，我们经常需要确定因果关系是x导致y，还是y导致x。对此，Granger提出了一种解决方法:如果x是y的原因，且不存在反向因果，则x过去值可以预测y未来值，反之则不然。具体来说，我们建立时间序列模型如下，并提出假设H0： βm=0，m=1,2…p。如果接受该假设，则意味着x过去值不能够预测y未来值；如果拒绝该假设，则可以，即x是y的格兰杰因(Granger cause)。

本次推送中，我们介绍如何在Stata中进行格兰杰因果检验。

基本命令

Stata中格兰杰因果检验有三种方法:前两种方法基本思路大致相同，它们均是先确定最佳滞后期，随后再进行格兰杰因果检验；第三种方法则是先拟合VAR模型，之后再进行检验。具体来说，第一种方法的基本命令为：

reg y l.y l.x

对滞后一期的变量进行回归，也就是说此时

式中p、q的值都为1。

estat ic

显示AIC 与BIC 取值，以便选择最佳滞后期。

reg y l.y l.x l2.y l2.xestat ic

再次显示AIC与BIC取值。

……

根据信息准则确定p, q后，检验所用的命令为test。值得注意的是，这种方法中p和q的取值可以是不同的，以便获得最有说服力的结果。

第二种方法的基本命令为：

ssc install gcause

下载格兰杰因果检验程序gcause

gcause y x, lags( 1)

对滞后一期的变量进行回归。

estat ic

显示AIC与BIC 取值，以便选择最佳滞后期

gcause y x, lags( 2)

滞后2期

estat ic

再次显示AIC与BIC取值，以便选择最佳滞后期。

在选定滞后期后，我们可以使用F检验或者卡方检验进行因果检验。

第三种方法的主要命令是：

vary x

进行向量自回归。

vargranger 操作实例

在这一部分中，我们以第二种方法为例，演示如何在stata中进行因果检验。首先我们在联网状态下导入数据，输入命令如下：

use http: //www.stata-press.com/data/imeus/ukrates, clear

安装gcause格兰杰因果检验程序

ssc install gcause

输出结果如下：

gcause r20 rs, lags( 1)estat icgcause r20 rs, lags( 2)estat ic

依次对滞后一期、滞后两期等变量进行回归，根据AIC及BIC的取值确定最佳的滞后期。在本例中，我们发现p=q=3时AIC及BIC的值最小，因此我们将p和q都赋值为3。滞后三期回归结果如图所示：

我们发现F检验和卡方检验得出一致结论，接受原假设，即rs不是r20的格兰杰因。

注：操作平台为Stata14.0

文章来源：格兰杰因果检验stata操作指南，部分文字进行删改。更多请点击“阅读原文” 返回搜狐，查看更多

责任编辑：

本文链接：https://blog.csdn.net/weixin_39675215/article/details/111730121

原作者删帖不实内容删帖广告或垃圾文章投诉

智能推荐

艾美捷Epigentek DNA样品的超声能量处理方案-程序员宅基地

文章浏览阅读15次。空化气泡的大小和相应的空化能量可以通过调整完全标度的振幅水平来操纵和数字控制。通过强调超声技术中的更高通量处理和防止样品污染，Epigentek EpiSonic超声仪可以轻松集成到现有的实验室工作流程中，并且特别适合与表观遗传学和下一代应用的兼容性。Epigentek的EpiSonic已成为一种有效的剪切设备，用于在染色质免疫沉淀技术中制备染色质样品，以及用于下一代测序平台的DNA文库制备。该装置的经济性及其多重样品的能力使其成为每个实验室拥有的经济高效的工具，而不仅仅是核心设施。

11、合宙Air模块Luat开发：通过http协议获取天气信息_合宙获取天气-程序员宅基地

文章浏览阅读4.2k次，点赞3次，收藏14次。目录点击这里查看所有博文本系列博客，理论上适用于合宙的Air202、Air268、Air720x、Air720S以及最近发布的Air720U（我还没拿到样机，应该也能支持）。先不管支不支持，如果你用的是合宙的模块，那都不妨一试，也许会有意外收获。我使用的是Air720SL模块，如果在其他模块上不能用，那就是底层core固件暂时还没有支持，这里的代码是没有问题的。例程仅供参考！..._合宙获取天气

EasyMesh和802.11s对比-程序员宅基地

文章浏览阅读7.7k次，点赞2次，收藏41次。1 关于meshMesh的意思是网状物，以前读书的时候，在自动化领域有传感器自组网，zigbee、蓝牙等无线方式实现各个网络节点消息通信，通过各种算法，保证整个网络中所有节点信息能经过多跳最终传递到目的地，用于数据采集。十多年过去了，在无线路由器领域又把这个mesh概念翻炒了一下，各大品牌都推出了mesh路由器，大多数是3个为一组，实现在面积较大的住宅里，增强wifi覆盖范围，智能在多热点之间切换，提升上网体验。因为节点基本上在3个以内，所以mesh的算法不必太复杂，组网形式比较简单。各厂家都自定义了组_802.11s

线程的几种状态_线程状态-程序员宅基地

文章浏览阅读5.2k次，点赞8次，收藏21次。线程的几种状态_线程状态

stack的常见用法详解_stack函数用法-程序员宅基地

文章浏览阅读4.2w次，点赞124次，收藏688次。stack翻译为栈，是STL中实现的一个后进先出的容器。要使用 stack，应先添加头文件include<stack>，并在头文件下面加上“ using namespacestd;"1. stack的定义其定义的写法和其他STL容器相同, typename可以任意基本数据类型或容器：stack<typename> name;2. stack容器内元素的访问..._stack函数用法

2018.11.16javascript课上随笔（DOM）-程序员宅基地

文章浏览阅读71次。<li> <a href = "“#”>-</a></li><li>子节点：文本节点（回车），元素节点，文本节点。不同节点树：　　节点（各种类型节点）childNodes:返回子节点的所有子节点的集合，包含任何类型、元素节点（元素类型节点）：child。node.getAttribute(at...

随便推点

layui.extend的一点知识第三方模块base 路径_layui extend-程序员宅基地

文章浏览阅读3.4k次。//config的设置是全局的layui.config({ base: '/res/js/' //假设这是你存放拓展模块的根目录}).extend({ //设定模块别名 mymod: 'mymod' //如果 mymod.js 是在根目录，也可以不用设定别名 ,mod1: 'admin/mod1' //相对于上述 base 目录的子目录}); //你也可以忽略 base 设定的根目录，直接在 extend 指定路径（主要：该功能为 layui 2.2.0 新增）layui.exten_layui extend

5G云计算：5G网络的分层思想_5g分层结构-程序员宅基地

文章浏览阅读3.2k次，点赞6次，收藏13次。分层思想分层思想分层思想-1分层思想-2分层思想-2OSI七层参考模型物理层和数据链路层物理层数据链路层网络层传输层会话层表示层应用层OSI七层模型的分层结构TCP/IP协议族的组成数据封装过程数据解封装过程PDU设备与层的对应关系各层通信分层思想分层思想-1在现实生活种，我们在喝牛奶时，未必了解他的生产过程，我们所接触的或许只是从超时购买牛奶。分层思想-2平时我们在网络时也未必知道数据的传输过程我们的所考虑的就是可以传就可以，不用管他时怎么传输的分层思想-2将复杂的流程分解为几个功能_5g分层结构