763
个编辑
更改
创建页面，内容为“{{#seo: |keywords=格兰杰因果关系,时间序列分析 |description=统计检验,时间序列 }} File:GrangerCausalityIllustration.png|360px|thumbnail|当时…”
{{#seo:
|keywords=格兰杰因果关系,时间序列分析
|description=统计检验,时间序列
}}
[[File:GrangerCausalityIllustration.png|360px|thumbnail|当时间序列X可以作为格兰杰原因引起时间序列Y时，X中的模式在经过一段时间迟滞后与Y大致重复（两个示例用箭头表示）。因此，X的过去值可用于预测Y的未来值。]]
格兰杰因果关系检测是一种[https://en.wikipedia.org/wiki/Statistical_hypothesis_testing 统计假设检验方法]，用于确定一个时间序列是否可以用来预测另一个时间序列，于1969年首次提出。<ref>{{cite journal |jstor=1912791 |pages=424–438 |last1=Granger |first1=C. W. J. |title=Investigating Causal Relations by Econometric Models and Cross-spectral Methods |volume=37 |issue=3 |journal=Econometrica |year=1969 |doi=10.2307/1912791}}</ref> 一般来说，回归反映的是“纯粹的”相关性，但'''克莱夫 · 格兰杰 Clive Granger'''认为，可以通过衡量使用一个时间序列的先验值来预测另一个时间序列未来值的能力来检验经济学中的因果关系。由于“真正因果关系”的问题是一个深刻的哲学问题，并且存在着后此谬误（后发者因此而发）的问题，即：因为A先于B，所以A引起B这种假设本身存在着谬误，因此计量经济学者断言格兰杰因果检验只能发现“预测因果关系”。<ref>{{cite book |last=Diebold |first=Francis X. |title=Elements of Forecasting |location=Cincinnati |publisher=South Western |year=2001 |edition=2nd |page=254 |isbn=978-0-324-02393-0 }}</ref>
时间序列X可以作为'''格兰杰原因 Granger-cause'''引起时间序列Y，通常当且仅当通过对<math>X</math>的滞后值（包括<math>Y</math>的滞后值）进行一系列t检验和F检验之后，这些<math>X</math>值可以提供关于<math>Y</math>的未来值的显著统计信息。
格兰杰本人还强调，在经济学以外的领域，有一些研究使用“格兰杰因果关系”检验得出了“荒谬”的结论。 他在诺贝尔奖演讲中说道：“当然，出现了许多荒谬的论文。”。 <ref>{{cite journal|last1=Granger|first1=Clive W. J|title=Time Series Analysis, Cointegration, and Applications|journal=American Economic Review|date=2004|volume=94|issue=3|pages=421–425|doi=10.1257/0002828041464669|url=https://www.nobelprize.org/uploads/2018/06/granger-lecture.pdf|accessdate=12 June 2019|citeseerx=10.1.1.370.6488}}</ref> 然而，由于计算简单，格兰杰因果关系检验仍然是时间序列因果关系分析的一种常用方法。格兰杰因果关系的原始定义没有考虑到潜在的混杂效应，也没有捕捉到瞬时的和非线性的因果关系，不过现在已经提出了几种扩展方法来解决这些问题。<ref name=Eichler>{{cite book|last1=Eichler|first1=Michael|editor1-last=Berzuini|editor1-first=Carlo|title=Causality : statistical perspectives and applications|date=2012|publisher=Wiley|location=Hoboken, N.J.|isbn=978-0470665565|pages=327–352|edition=3rd|chapter-url=http://researchers-sbe.unimaas.nl/michaeleichler/wp-content/uploads/sites/31/2014/02/causalstatistics.pdf|chapter=Causal Inference in Time Series Analysis}}</ref>
==问题思路==
我们认为，在预测Y值时，如果综合使用Y和X的过去值优于仅使用Y的过去值，则随时间演化的变量X可以导致另一个演化变量Y。
==基本原则==
格兰杰基于两个原则定义了因果关系：<ref name ="Eichler"/><ref name=Granger80>{{cite journal|last1=Granger|first1=C.W.J.|title=Testing for causality: A personal viewpoint|journal=Journal of Economic Dynamics and Control|date=1980|volume=2|pages=329–352|doi=10.1016/0165-1889(80)90069-X}}</ref>
# 自变量发生在因变量之前。
# 这一自变量对于被其影响的因变量的未来值有着独特的信息。
鉴于这两个关于因果关系的假设，格兰杰提出了如下假说，用来衡量<math>𝑋</math>对<math>𝑌</math>的因果影响：
:<math>
\mathbb{P}[Y(t+1) \in A\mid \mathcal{I}(t)] \neq \mathbb{P}[Y(t+1) \in A\mid \mathcal{I}_{-X}(t)],
</math>
其中<math>\mathbb{P}</math>是概率，<math>𝐴</math>是任意的非空集合，<math>\mathcal{I}(t)</math>和<math>\mathcal{I}_x(t)</math>分别表示整个宇宙中截止到时间<math>𝑡</math>的可用信息和排除了<math>𝑋</math>的修正宇宙中的可用信息。如果上述假设成立，我们可以说<math>X</math>是<math>Y</math>的格兰杰原因（<math>X</math>可以格兰杰引起<math>Y</math>）。
==格兰杰因果分析方法==
如果[[时间序列]]是平稳过程，则使用两个（或更多）变量的水平值执行测试。如果变量是非平稳的，则使用一阶差分（或高阶差分）进行测试。通常使用信息准则来选择要包含的滞后数，例如'''赤池信息量准则 the Akaike information criterion'''或'''[https://en.wikipedia.org/wiki/Schwarz_information_criterion 施瓦兹信息准则]the Schwarz information criterion'''。当满足以下条件时，变量之一的任意特定滞后值在回归中被保留，即（1）根据'''t检验 t-test'''它是显著的，并且（2）根据'''F检验 F-test '''，它和变量的其他滞后值共同增加了模型的[https://en.wikipedia.org/wiki/Explanatory_power 解释力]。只有当且仅当回归中没有保留解释变量的滞后值时，才不否定无格兰杰因果关系的零假设。
实践中可能会出现两个变量均不互为格兰杰原因，或两个变量均互为格兰杰原因的情况。
===数学表述===
设<math>y</math>和<math>x</math>均为平稳时间序列。为了检验<math>x</math>不是<math>y</math>的格兰杰原因这一零假设，首先要<math>y</math>的适当滞后值，以包含在<math>y</math>的单变量[https://en.wikipedia.org/wiki/Autoregression 自回归]中：
:<math>y_t = a_0 + a_1y_{t-1} + a_2y_{t-2} + \cdots + a_my_{t-m} + \text{error}_t.</math>
接下来，使用包含<math>x</math>的滞后值来扩展这一自回归:
:<math>y_t = a_0 + a_1y_{t-1} + a_2y_{t-2} + \cdots + a_my_{t-m} + b_px_{t-p} + \cdots + b_qx_{t-q} + \text{error}_t.</math>
在这个回归中，保留了根据t统计量所有单个显著的<math>x</math>的滞后值，前提是根据F检验，它们共同增加了回归的解释力(F检验的零假设不是x联合增加的解释力)。 在上述增广回归的式子中，对于显著的<math>x</math>的滞后值，<math>p</math>是最短的滞后期，<math>q</math>是最长的滞后期。
当且仅当回归中没有保留<math>x</math>的滞后值时，<math>x</math>不是<math>y</math>的格兰杰原因的零假设被接受。
===多元分析===
多元格兰杰因果关系分析通常是通过对时间序列拟合一个'''[https://en.wikipedia.org/wiki/Vector_autoregression 向量自回归模型] VAR'''来完成的。特别地，对 <math>t=1, \ldots, T</math> 是一个 <math>d</math> 维多元时间序列，设 <math>X(t) \in \mathbb{R}^{d\times 1}</math> 。格兰杰因果关系可以通过拟合一个带<math> L</math> 时滞的 VAR 模型来执行:
:<math>
X(t) = \sum_{\tau=1}^L A_{\tau}X(t-\tau) + \varepsilon(t),
</math>
其中，<math> 𝜀(t)</math> 是一个'''白高斯随机向量 white Gaussian random vector'''， <math>A_\tau</math> 是每个 <math>\tau</math> 的矩阵。如果其中至少一个元素 <math>A_{\tau}(j, i)</math> 对于 <math>\tau = 1, \ldots, L</math> （绝对值）显著大于零，则时间序列 <math>X_i</math> 被称为另一个时间序列 <math>X_j</math> 的格兰杰原因。<ref>{{cite book|last1=Lütkepohl|first1=Helmut|author-link = | title=New introduction to multiple time series analysis|date=2005|publisher=Springer|location=Berlin|isbn=978-3540262398|pages=41–51|edition=3}}</ref>
===非参数检验===
上述线性的方法适用于检验格兰杰因果关系的均值。但是，这种方法无法在高阶矩中检测格兰杰因果关系，例如方差。格兰杰因果关系的非参数检验就是为了解决这个问题而设计的。<ref>{{cite journal |last=Diks |first=Cees |last2=Panchenko |first2=Valentyn |year=2006 |title=A new statistic and practical guidelines for nonparametric Granger causality testing |journal=Journal of Economic Dynamics and Control |volume=30 |issue=9 |pages=1647–1669 |url=http://research.economics.unsw.edu.au/vpanchenko/papers/2006_GC_JEDC.pdf |doi=10.1016/j.jedc.2005.08.008}}</ref>这些检验中的格兰杰因果关系为一般定义，不涉及任何建模假设，如线性自回归模型。格兰杰因果关系的非参数检验可以作为诊断工具，用来建立包括高阶矩和/或非线性的更优参数模型。<ref>{{cite journal |last=Francis |first=Bill B.|last2=Mougoue|first2=Mbodja |last3=Panchenko |first3=Valentyn |year=2010 |title=Is there a Symmetric Nonlinear Causal Relationship between Large and Small Firms? |journal=Journal of Empirical Finance |volume=17 |issue=1 |pages=23–28 |url=http://research.economics.unsw.edu.au/vpanchenko/papers/2009_SLFIRMS_JEF.pdf |doi=10.1016/j.jempfin.2009.08.003}}</ref>
==局限性==
顾名思义，格兰杰因果关系并不一定是真正的因果关系。实际上，格兰杰因果关系检验只能满足[https://en.wikipedia.org/wiki/Humean_definition_of_causality 休谟因果关系定义]，而休谟定义用经常性联结来定义因果关系。<ref name=":0">{{Cite journal|last=Mariusz|first=Maziarz|date=2015-05-20|title=A review of the Granger-causality fallacy|url=http://hrcak.srce.hr/index.php?show=clanak&id_clanak_jezik=229894|journal=The Journal of Philosophical Economics: Reflections on Economic and Social Issues|language=en|volume=VIII.|issue=2|issn=1843-2298}}</ref>如果,<math>x</math>和<math>y</math>都是由具有不同滞后的共同第三过程驱动的，人们可能仍然无法拒绝格兰杰因果关系的备择假设。然而，操纵其中一个变量不会改变另一个变量。实际上，格兰杰因果检验是为了处理成对变量而设计的，当真实关系涉及三个或更多变量时，还可能会产生误导性的结果。尽管如此，有人认为，考虑原因的概率视角的话，格兰杰因果关系仍然可以被认为是真正的因果关系，特别是当考虑到赖欣巴哈的概率因果关系中的“屏蔽”概念时。<ref>{{cite journal |doi=10.1016/j.plrev.2015.09.002 |pmid=26429630 |title=Foundational perspectives on causality in large-scale brain networks |journal=Physics of Life Reviews |volume=15 |pages=107–23 |year=2015 |last1=Mannino |first1=Michael |last2=Bressler |first2=Steven L |bibcode=2015PhLRv..15..107M }}</ref> 其他可能导致产生误导检验结果的原因有:
# 采样不够频繁或过于频繁;
# 非线性因果关系
# 时间序列的非平稳性和非线性
# 理性预期的存在。
一个包含更多变量的类似检验可以应用[https://en.wikipedia.org/wiki/Vector_autoregression 向量自回归]。
== 扩展==
从误差项是正态分布的假设中，发展出了一种对偏差不敏感的格兰杰因果关系分析方法。<ref>{{cite journal |pages=1489–1500 |doi=10.1080/00036840500405763|title=Tests for causality between integrated variables using asymptotic and bootstrap distributions: Theory and application|journal=Applied Economics|volume=38|issue=13|year=2006|last1=Hacker|first1=R. Scott|last2=Hatemi-j|first2=A.|url=https://ideas.repec.org/a/taf/applec/v38y2006i13p1489-1500.html}}</ref> 由于许多金融变量是非正态分布的，因此这种方法在金融经济学中有着很大的用处。<ref>{{cite journal |pages=394–419 |doi=10.1086/294632 |title=The Variation of Certain Speculative Prices |journal=The Journal of Business |volume=36 |issue=4 |year=1963 |last1=Mandelbrot |first1=Benoit }}</ref> 最近，文献中提出了非对称因果关系检验，以区分积极变化和消极变化的因果影响。<ref>{{cite journal |pages=447–456 |doi=10.1007/s00181-011-0484-x |title=Asymmetric causality tests with an application |journal=Empirical Economics |volume=43 |year=2012 |last1=Hatemi-j |first1=A. |url=https://ideas.repec.org/a/spr/empeco/v43y2012i1p447-456.html}}</ref> 同时，格兰杰(非)因果关系检验在面板数据上的扩展也是可行的。<ref>{{cite journal |title=Testing for Granger non-causality in heterogeneous panels |year=2012 |last1=Dumistrescu |first1=E.-I. |last2=Hurlin |first2=C. |journal=Economic Modelling |volume=29 |issue=4 |pages=1450–1460 |doi=10.1016/j.econmod.2012.02.014|citeseerx=10.1.1.395.568 }}</ref>
==在神经科学领域==
关于神经功能的一个长期观念认为，大脑的不同区域是用来执行特定任务的；与特定区域局部的结构性连结以某种方式决定了该部分的功能。通过收集多年来完成的工作，已经有了一种不同的、以网络为中心的方法来描述大脑中的信息流。对功能的解释开始包括在大脑不同层次和不同位置存在的网络这一概念。<ref name=pmid17569852>{{cite journal |doi=10.1126/science.1144677 |pmid=17569852 |title=NEUROSCIENCE: Neural Networks Debunk Phrenology |journal=Science |volume=316 |issue=5831 |pages=1578–9 |year=2007 |last1=Knight |first1=R. T }}</ref> 这些网络的行为可以用随时间演化的非确定性过程来描述。也就是说，给定相同的输入刺激，网络的输出会是不同的。这些网络的动态是由概率控制的，所以我们把它们当作随机过程，这样我们就可以捕捉到大脑不同区域之间的这种动态。
过去人们曾探索过从神经元及其周围集合的放电活动中获得某种信息流度量的不同方法，但是这些方法可以得出的结论种类有限，对信息的方向性流动、其影响大小以及它如何随时间变化提供的见解很少。<ref name="Kim">{{cite journal |doi=10.1371/journal.pcbi.1001110 |pmid=21455283 |pmc=3063721 |title=A Granger Causality Measure for Point Process Models of Ensemble Neural Spiking Activity |journal=PLoS Computational Biology |volume=7 |issue=3 |pages=e1001110 |year=2011 |last1=Kim |first1=Sanggyun |last2=Putrino |first2=David |last3=Ghosh |first3=Soumya |last4=Brown |first4=Emery N |bibcode=2011PLSCB...7E1110K }}</ref> 最近，格兰杰因果关系已经应用于解决其中一些问题，并取得了巨大成功。<ref name=pmid20202481>{{cite journal |doi=10.1016/j.neuroimage.2010.02.059 |pmid=20202481 |title=Wiener–Granger Causality: A well established methodology |journal=NeuroImage |volume=58 |issue=2 |pages=323–9 |year=2011 |last1=Bressler |first1=Steven L |last2=Seth |first2=Anil K }}</ref> 简单地说，就是研究如何最好地预测一个神经元的未来：要么使用整个神经元集合，要么使用除了某个目标神经元之外的整个集合。如果把目标神经元排除在外使预测变得更糟，那么我们说它与当前神经元存在“g-因果”关系。
===点过程模型的扩展===
以往的格兰杰因果关系分析方法只能对连续值数据进行分析，因此对'''神经[https://en.wikipedia.org/wiki/Action_potential 尖峰序列] neural spike train'''记录的分析涉及到了转换，这些转换最终改变了数据的随机属性，从而间接地改变了由此得出的结论的有效性。 然而，在2011年，一个新的通用格兰杰因果关系框架被提出，使用这种框架可以直接操作任何模态，包括神经尖峰序列。<ref name="Kim" />
神经尖峰序列数据可以被建模为一个'''[https://en.wikipedia.org/wiki/Point_process 点过程] Point process'''。时间点过程是一个连续时间内发生的二元事件的随机时间序列。它只能在每个时间点采用两个值，表示事件是否实际发生。由于单个神经元的动作电位具有典型的波形，这种信息的二值表示形式适合于[https://en.wikipedia.org/wiki/Neuronal_ensemble 神经元集群]的活动。通过这种方式，携带从神经元输出的实际信息是“尖峰”的出现，以及连续尖峰之间的时间。利用这种方法，可以将神经网络中的信息流简单概括为每个神经元在一个观察周期内的峰值时间。点进程可以用尖峰本身的计时，尖峰之间的等待时间以及使用计数过程来表示，或者，如果时间离散到足以确保每个空当只有一个事件有可能发生，即，一个时间段只能包含一个事件，作为1和0的集合，也可以被表示为非常类似于二进制的方式。
最简单的神经尖峰模型之一是'''[https://en.wikipedia.org/wiki/Poisson_point_process 泊松过程] Poisson process'''。然而，这是有限的，因为它是无记忆性的。在计算当前的放电概率时，它不考虑任何峰值历史。然而，神经元通过其相对和绝对不应期表现出基本的(生物物理)历史依赖性。为了解决这个问题，一个条件强度函数被用来表示一个神经元在其历史条件下出现尖峰的概率。条件强度函数表示瞬时放电概率，并隐式定义点过程的完整概率模型。它定义了每个单位时间的概率。因此，如果这个单位时间足够小，以确保在那个时间窗口中只能发生一个尖峰，那么我们的条件强度函数完全指定了一个给定神经元在一定时间内放电的概率。
==参阅==
* [https://baike.baidu.com/item/%E6%A0%BC%E5%85%B0%E6%9D%B0%E5%9B%A0%E6%9E%9C%E5%85%B3%E7%B3%BB%E6%A3%80%E9%AA%8C/2485970?fr=aladdin 百度百科]
* [https://zhuanlan.zhihu.com/p/79936875 知乎：是前因还是后果?——在真格量化中进行格兰杰因果检验（含Python代码）]
* [https://blog.csdn.net/luciazxx/article/details/44224145 CSDN：格兰杰因果关系检验]
==参考文献==
{{reflist|30em}}
==扩展阅读==
* {{cite book |last=Enders |first=Walter |title=Applied Econometric Time Series 应用计量经济学时间序列|location=New York |publisher=Wiley |edition=Second |year=2004 |isbn=978-0-471-23065-6 |pages=[https://book.douban.com/subject/1480917/] |url-access=registration |url=https://book.douban.com/subject/1480917/}}
* {{cite book |last=Gujarati |first=Damodar N. |last2=Porter |first2=Dawn C. |title=Basic Econometrics |location=New York |publisher=McGraw-Hill |year=2009 |edition=Fifth international |isbn=978-007-127625-2 |chapter=Causality in Economics: The Granger Causality Test 经济学中的因果关系：格兰杰因果检验|pages=652–658 }}
* {{cite book |last=Hoover |first=Kevin D. |authorlink= |chapter=Granger-causality |pages=[https://archive.org/details/newclassicalmacr0000hoov/page/168 168–176] |title=The New Classical Macroeconomics 新古典宏观经济学|location=Oxford |publisher=Basil Blackwell |year=1988 |isbn=978-0-631-14605-6 |url=https://archive.org/details/newclassicalmacr0000hoov/page/168 }}
* {{cite encyclopedia |year= 2008 |title = Granger–Sims causality |last= Kuersteiner|first=Guido|encyclopedia=The New Palgrave Dictionary of Economics 新帕尔格雷夫经济学大词典 |url=https://book.douban.com/subject/1481407/}}
*Kleinberg, S. and Hripcsak, G. (2011) [http://people.dbmi.columbia.edu/samantha/papers/kleinberg_jbi11_preprint.pdf "生物医学信息学因果推理综述"] 于2012年4月30日[https://web.archive.org/web/20120430094823/http://people.dbmi.columbia.edu/samantha/papers/kleinberg_jbi11_preprint.pdf 存档]于[https://en.wikipedia.org/wiki/Wayback_Machine Wayback Machine] J. Biomed Informatics
==编者推荐==
[[File:GLJ1.jpg|200px|缩略图|右|格兰杰计量经济学文集封面]]
*[https://book.douban.com/subject/2369017/ 格兰杰计量经济学文集]
::格兰杰是最具创造力的现代计量经济学家之一，他对帮助我们理解因果关系、建模方法、非平稳性、季节性和预测做出了重要的贡献。一些被广泛应用的时间序列计量经济学概念就来自于他的研究成果，这是因为格兰杰善于将观测到的经济结果系统地用模型来表述，而且这些模型很好地反映了经济体复杂和进化的本质，这本书收录了他发表过的作品，而且把主要主题下相关联的研究集合到一起，这样便更容易找到相关的论文，里面的很多文章值得我们认真反复研读。
----
本中文词条由[[用户:Yillia Jing|Yillia Jing|]]参与编译，[[用户:苏格兰|苏格兰]]审校，[[用户:乐多多|乐多多]]编辑，欢迎在讨论页面留言。
'''本词条内容源自wikipedia及公开资料，遵守 CC3.0协议。'''
[[Category:多元时间序列]]
[[Category:时间序列统计检验]]
|keywords=格兰杰因果关系,时间序列分析
|description=统计检验,时间序列
}}
[[File:GrangerCausalityIllustration.png|360px|thumbnail|当时间序列X可以作为格兰杰原因引起时间序列Y时，X中的模式在经过一段时间迟滞后与Y大致重复（两个示例用箭头表示）。因此，X的过去值可用于预测Y的未来值。]]
格兰杰因果关系检测是一种[https://en.wikipedia.org/wiki/Statistical_hypothesis_testing 统计假设检验方法]，用于确定一个时间序列是否可以用来预测另一个时间序列，于1969年首次提出。<ref>{{cite journal |jstor=1912791 |pages=424–438 |last1=Granger |first1=C. W. J. |title=Investigating Causal Relations by Econometric Models and Cross-spectral Methods |volume=37 |issue=3 |journal=Econometrica |year=1969 |doi=10.2307/1912791}}</ref> 一般来说，回归反映的是“纯粹的”相关性，但'''克莱夫 · 格兰杰 Clive Granger'''认为，可以通过衡量使用一个时间序列的先验值来预测另一个时间序列未来值的能力来检验经济学中的因果关系。由于“真正因果关系”的问题是一个深刻的哲学问题，并且存在着后此谬误（后发者因此而发）的问题，即：因为A先于B，所以A引起B这种假设本身存在着谬误，因此计量经济学者断言格兰杰因果检验只能发现“预测因果关系”。<ref>{{cite book |last=Diebold |first=Francis X. |title=Elements of Forecasting |location=Cincinnati |publisher=South Western |year=2001 |edition=2nd |page=254 |isbn=978-0-324-02393-0 }}</ref>
时间序列X可以作为'''格兰杰原因 Granger-cause'''引起时间序列Y，通常当且仅当通过对<math>X</math>的滞后值（包括<math>Y</math>的滞后值）进行一系列t检验和F检验之后，这些<math>X</math>值可以提供关于<math>Y</math>的未来值的显著统计信息。
格兰杰本人还强调，在经济学以外的领域，有一些研究使用“格兰杰因果关系”检验得出了“荒谬”的结论。 他在诺贝尔奖演讲中说道：“当然，出现了许多荒谬的论文。”。 <ref>{{cite journal|last1=Granger|first1=Clive W. J|title=Time Series Analysis, Cointegration, and Applications|journal=American Economic Review|date=2004|volume=94|issue=3|pages=421–425|doi=10.1257/0002828041464669|url=https://www.nobelprize.org/uploads/2018/06/granger-lecture.pdf|accessdate=12 June 2019|citeseerx=10.1.1.370.6488}}</ref> 然而，由于计算简单，格兰杰因果关系检验仍然是时间序列因果关系分析的一种常用方法。格兰杰因果关系的原始定义没有考虑到潜在的混杂效应，也没有捕捉到瞬时的和非线性的因果关系，不过现在已经提出了几种扩展方法来解决这些问题。<ref name=Eichler>{{cite book|last1=Eichler|first1=Michael|editor1-last=Berzuini|editor1-first=Carlo|title=Causality : statistical perspectives and applications|date=2012|publisher=Wiley|location=Hoboken, N.J.|isbn=978-0470665565|pages=327–352|edition=3rd|chapter-url=http://researchers-sbe.unimaas.nl/michaeleichler/wp-content/uploads/sites/31/2014/02/causalstatistics.pdf|chapter=Causal Inference in Time Series Analysis}}</ref>
==问题思路==
我们认为，在预测Y值时，如果综合使用Y和X的过去值优于仅使用Y的过去值，则随时间演化的变量X可以导致另一个演化变量Y。
==基本原则==
格兰杰基于两个原则定义了因果关系：<ref name ="Eichler"/><ref name=Granger80>{{cite journal|last1=Granger|first1=C.W.J.|title=Testing for causality: A personal viewpoint|journal=Journal of Economic Dynamics and Control|date=1980|volume=2|pages=329–352|doi=10.1016/0165-1889(80)90069-X}}</ref>
# 自变量发生在因变量之前。
# 这一自变量对于被其影响的因变量的未来值有着独特的信息。
鉴于这两个关于因果关系的假设，格兰杰提出了如下假说，用来衡量<math>𝑋</math>对<math>𝑌</math>的因果影响：
:<math>
\mathbb{P}[Y(t+1) \in A\mid \mathcal{I}(t)] \neq \mathbb{P}[Y(t+1) \in A\mid \mathcal{I}_{-X}(t)],
</math>
其中<math>\mathbb{P}</math>是概率，<math>𝐴</math>是任意的非空集合，<math>\mathcal{I}(t)</math>和<math>\mathcal{I}_x(t)</math>分别表示整个宇宙中截止到时间<math>𝑡</math>的可用信息和排除了<math>𝑋</math>的修正宇宙中的可用信息。如果上述假设成立，我们可以说<math>X</math>是<math>Y</math>的格兰杰原因（<math>X</math>可以格兰杰引起<math>Y</math>）。
==格兰杰因果分析方法==
如果[[时间序列]]是平稳过程，则使用两个（或更多）变量的水平值执行测试。如果变量是非平稳的，则使用一阶差分（或高阶差分）进行测试。通常使用信息准则来选择要包含的滞后数，例如'''赤池信息量准则 the Akaike information criterion'''或'''[https://en.wikipedia.org/wiki/Schwarz_information_criterion 施瓦兹信息准则]the Schwarz information criterion'''。当满足以下条件时，变量之一的任意特定滞后值在回归中被保留，即（1）根据'''t检验 t-test'''它是显著的，并且（2）根据'''F检验 F-test '''，它和变量的其他滞后值共同增加了模型的[https://en.wikipedia.org/wiki/Explanatory_power 解释力]。只有当且仅当回归中没有保留解释变量的滞后值时，才不否定无格兰杰因果关系的零假设。
实践中可能会出现两个变量均不互为格兰杰原因，或两个变量均互为格兰杰原因的情况。
===数学表述===
设<math>y</math>和<math>x</math>均为平稳时间序列。为了检验<math>x</math>不是<math>y</math>的格兰杰原因这一零假设，首先要<math>y</math>的适当滞后值，以包含在<math>y</math>的单变量[https://en.wikipedia.org/wiki/Autoregression 自回归]中：
:<math>y_t = a_0 + a_1y_{t-1} + a_2y_{t-2} + \cdots + a_my_{t-m} + \text{error}_t.</math>
接下来，使用包含<math>x</math>的滞后值来扩展这一自回归:
:<math>y_t = a_0 + a_1y_{t-1} + a_2y_{t-2} + \cdots + a_my_{t-m} + b_px_{t-p} + \cdots + b_qx_{t-q} + \text{error}_t.</math>
在这个回归中，保留了根据t统计量所有单个显著的<math>x</math>的滞后值，前提是根据F检验，它们共同增加了回归的解释力(F检验的零假设不是x联合增加的解释力)。 在上述增广回归的式子中，对于显著的<math>x</math>的滞后值，<math>p</math>是最短的滞后期，<math>q</math>是最长的滞后期。
当且仅当回归中没有保留<math>x</math>的滞后值时，<math>x</math>不是<math>y</math>的格兰杰原因的零假设被接受。
===多元分析===
多元格兰杰因果关系分析通常是通过对时间序列拟合一个'''[https://en.wikipedia.org/wiki/Vector_autoregression 向量自回归模型] VAR'''来完成的。特别地，对 <math>t=1, \ldots, T</math> 是一个 <math>d</math> 维多元时间序列，设 <math>X(t) \in \mathbb{R}^{d\times 1}</math> 。格兰杰因果关系可以通过拟合一个带<math> L</math> 时滞的 VAR 模型来执行:
:<math>
X(t) = \sum_{\tau=1}^L A_{\tau}X(t-\tau) + \varepsilon(t),
</math>
其中，<math> 𝜀(t)</math> 是一个'''白高斯随机向量 white Gaussian random vector'''， <math>A_\tau</math> 是每个 <math>\tau</math> 的矩阵。如果其中至少一个元素 <math>A_{\tau}(j, i)</math> 对于 <math>\tau = 1, \ldots, L</math> （绝对值）显著大于零，则时间序列 <math>X_i</math> 被称为另一个时间序列 <math>X_j</math> 的格兰杰原因。<ref>{{cite book|last1=Lütkepohl|first1=Helmut|author-link = | title=New introduction to multiple time series analysis|date=2005|publisher=Springer|location=Berlin|isbn=978-3540262398|pages=41–51|edition=3}}</ref>
===非参数检验===
上述线性的方法适用于检验格兰杰因果关系的均值。但是，这种方法无法在高阶矩中检测格兰杰因果关系，例如方差。格兰杰因果关系的非参数检验就是为了解决这个问题而设计的。<ref>{{cite journal |last=Diks |first=Cees |last2=Panchenko |first2=Valentyn |year=2006 |title=A new statistic and practical guidelines for nonparametric Granger causality testing |journal=Journal of Economic Dynamics and Control |volume=30 |issue=9 |pages=1647–1669 |url=http://research.economics.unsw.edu.au/vpanchenko/papers/2006_GC_JEDC.pdf |doi=10.1016/j.jedc.2005.08.008}}</ref>这些检验中的格兰杰因果关系为一般定义，不涉及任何建模假设，如线性自回归模型。格兰杰因果关系的非参数检验可以作为诊断工具，用来建立包括高阶矩和/或非线性的更优参数模型。<ref>{{cite journal |last=Francis |first=Bill B.|last2=Mougoue|first2=Mbodja |last3=Panchenko |first3=Valentyn |year=2010 |title=Is there a Symmetric Nonlinear Causal Relationship between Large and Small Firms? |journal=Journal of Empirical Finance |volume=17 |issue=1 |pages=23–28 |url=http://research.economics.unsw.edu.au/vpanchenko/papers/2009_SLFIRMS_JEF.pdf |doi=10.1016/j.jempfin.2009.08.003}}</ref>
==局限性==
顾名思义，格兰杰因果关系并不一定是真正的因果关系。实际上，格兰杰因果关系检验只能满足[https://en.wikipedia.org/wiki/Humean_definition_of_causality 休谟因果关系定义]，而休谟定义用经常性联结来定义因果关系。<ref name=":0">{{Cite journal|last=Mariusz|first=Maziarz|date=2015-05-20|title=A review of the Granger-causality fallacy|url=http://hrcak.srce.hr/index.php?show=clanak&id_clanak_jezik=229894|journal=The Journal of Philosophical Economics: Reflections on Economic and Social Issues|language=en|volume=VIII.|issue=2|issn=1843-2298}}</ref>如果,<math>x</math>和<math>y</math>都是由具有不同滞后的共同第三过程驱动的，人们可能仍然无法拒绝格兰杰因果关系的备择假设。然而，操纵其中一个变量不会改变另一个变量。实际上，格兰杰因果检验是为了处理成对变量而设计的，当真实关系涉及三个或更多变量时，还可能会产生误导性的结果。尽管如此，有人认为，考虑原因的概率视角的话，格兰杰因果关系仍然可以被认为是真正的因果关系，特别是当考虑到赖欣巴哈的概率因果关系中的“屏蔽”概念时。<ref>{{cite journal |doi=10.1016/j.plrev.2015.09.002 |pmid=26429630 |title=Foundational perspectives on causality in large-scale brain networks |journal=Physics of Life Reviews |volume=15 |pages=107–23 |year=2015 |last1=Mannino |first1=Michael |last2=Bressler |first2=Steven L |bibcode=2015PhLRv..15..107M }}</ref> 其他可能导致产生误导检验结果的原因有:
# 采样不够频繁或过于频繁;
# 非线性因果关系
# 时间序列的非平稳性和非线性
# 理性预期的存在。
一个包含更多变量的类似检验可以应用[https://en.wikipedia.org/wiki/Vector_autoregression 向量自回归]。
== 扩展==
从误差项是正态分布的假设中，发展出了一种对偏差不敏感的格兰杰因果关系分析方法。<ref>{{cite journal |pages=1489–1500 |doi=10.1080/00036840500405763|title=Tests for causality between integrated variables using asymptotic and bootstrap distributions: Theory and application|journal=Applied Economics|volume=38|issue=13|year=2006|last1=Hacker|first1=R. Scott|last2=Hatemi-j|first2=A.|url=https://ideas.repec.org/a/taf/applec/v38y2006i13p1489-1500.html}}</ref> 由于许多金融变量是非正态分布的，因此这种方法在金融经济学中有着很大的用处。<ref>{{cite journal |pages=394–419 |doi=10.1086/294632 |title=The Variation of Certain Speculative Prices |journal=The Journal of Business |volume=36 |issue=4 |year=1963 |last1=Mandelbrot |first1=Benoit }}</ref> 最近，文献中提出了非对称因果关系检验，以区分积极变化和消极变化的因果影响。<ref>{{cite journal |pages=447–456 |doi=10.1007/s00181-011-0484-x |title=Asymmetric causality tests with an application |journal=Empirical Economics |volume=43 |year=2012 |last1=Hatemi-j |first1=A. |url=https://ideas.repec.org/a/spr/empeco/v43y2012i1p447-456.html}}</ref> 同时，格兰杰(非)因果关系检验在面板数据上的扩展也是可行的。<ref>{{cite journal |title=Testing for Granger non-causality in heterogeneous panels |year=2012 |last1=Dumistrescu |first1=E.-I. |last2=Hurlin |first2=C. |journal=Economic Modelling |volume=29 |issue=4 |pages=1450–1460 |doi=10.1016/j.econmod.2012.02.014|citeseerx=10.1.1.395.568 }}</ref>
==在神经科学领域==
关于神经功能的一个长期观念认为，大脑的不同区域是用来执行特定任务的；与特定区域局部的结构性连结以某种方式决定了该部分的功能。通过收集多年来完成的工作，已经有了一种不同的、以网络为中心的方法来描述大脑中的信息流。对功能的解释开始包括在大脑不同层次和不同位置存在的网络这一概念。<ref name=pmid17569852>{{cite journal |doi=10.1126/science.1144677 |pmid=17569852 |title=NEUROSCIENCE: Neural Networks Debunk Phrenology |journal=Science |volume=316 |issue=5831 |pages=1578–9 |year=2007 |last1=Knight |first1=R. T }}</ref> 这些网络的行为可以用随时间演化的非确定性过程来描述。也就是说，给定相同的输入刺激，网络的输出会是不同的。这些网络的动态是由概率控制的，所以我们把它们当作随机过程，这样我们就可以捕捉到大脑不同区域之间的这种动态。
过去人们曾探索过从神经元及其周围集合的放电活动中获得某种信息流度量的不同方法，但是这些方法可以得出的结论种类有限，对信息的方向性流动、其影响大小以及它如何随时间变化提供的见解很少。<ref name="Kim">{{cite journal |doi=10.1371/journal.pcbi.1001110 |pmid=21455283 |pmc=3063721 |title=A Granger Causality Measure for Point Process Models of Ensemble Neural Spiking Activity |journal=PLoS Computational Biology |volume=7 |issue=3 |pages=e1001110 |year=2011 |last1=Kim |first1=Sanggyun |last2=Putrino |first2=David |last3=Ghosh |first3=Soumya |last4=Brown |first4=Emery N |bibcode=2011PLSCB...7E1110K }}</ref> 最近，格兰杰因果关系已经应用于解决其中一些问题，并取得了巨大成功。<ref name=pmid20202481>{{cite journal |doi=10.1016/j.neuroimage.2010.02.059 |pmid=20202481 |title=Wiener–Granger Causality: A well established methodology |journal=NeuroImage |volume=58 |issue=2 |pages=323–9 |year=2011 |last1=Bressler |first1=Steven L |last2=Seth |first2=Anil K }}</ref> 简单地说，就是研究如何最好地预测一个神经元的未来：要么使用整个神经元集合，要么使用除了某个目标神经元之外的整个集合。如果把目标神经元排除在外使预测变得更糟，那么我们说它与当前神经元存在“g-因果”关系。
===点过程模型的扩展===
以往的格兰杰因果关系分析方法只能对连续值数据进行分析，因此对'''神经[https://en.wikipedia.org/wiki/Action_potential 尖峰序列] neural spike train'''记录的分析涉及到了转换，这些转换最终改变了数据的随机属性，从而间接地改变了由此得出的结论的有效性。 然而，在2011年，一个新的通用格兰杰因果关系框架被提出，使用这种框架可以直接操作任何模态，包括神经尖峰序列。<ref name="Kim" />
神经尖峰序列数据可以被建模为一个'''[https://en.wikipedia.org/wiki/Point_process 点过程] Point process'''。时间点过程是一个连续时间内发生的二元事件的随机时间序列。它只能在每个时间点采用两个值，表示事件是否实际发生。由于单个神经元的动作电位具有典型的波形，这种信息的二值表示形式适合于[https://en.wikipedia.org/wiki/Neuronal_ensemble 神经元集群]的活动。通过这种方式，携带从神经元输出的实际信息是“尖峰”的出现，以及连续尖峰之间的时间。利用这种方法，可以将神经网络中的信息流简单概括为每个神经元在一个观察周期内的峰值时间。点进程可以用尖峰本身的计时，尖峰之间的等待时间以及使用计数过程来表示，或者，如果时间离散到足以确保每个空当只有一个事件有可能发生，即，一个时间段只能包含一个事件，作为1和0的集合，也可以被表示为非常类似于二进制的方式。
最简单的神经尖峰模型之一是'''[https://en.wikipedia.org/wiki/Poisson_point_process 泊松过程] Poisson process'''。然而，这是有限的，因为它是无记忆性的。在计算当前的放电概率时，它不考虑任何峰值历史。然而，神经元通过其相对和绝对不应期表现出基本的(生物物理)历史依赖性。为了解决这个问题，一个条件强度函数被用来表示一个神经元在其历史条件下出现尖峰的概率。条件强度函数表示瞬时放电概率，并隐式定义点过程的完整概率模型。它定义了每个单位时间的概率。因此，如果这个单位时间足够小，以确保在那个时间窗口中只能发生一个尖峰，那么我们的条件强度函数完全指定了一个给定神经元在一定时间内放电的概率。
==参阅==
* [https://baike.baidu.com/item/%E6%A0%BC%E5%85%B0%E6%9D%B0%E5%9B%A0%E6%9E%9C%E5%85%B3%E7%B3%BB%E6%A3%80%E9%AA%8C/2485970?fr=aladdin 百度百科]
* [https://zhuanlan.zhihu.com/p/79936875 知乎：是前因还是后果?——在真格量化中进行格兰杰因果检验（含Python代码）]
* [https://blog.csdn.net/luciazxx/article/details/44224145 CSDN：格兰杰因果关系检验]
==参考文献==
{{reflist|30em}}
==扩展阅读==
* {{cite book |last=Enders |first=Walter |title=Applied Econometric Time Series 应用计量经济学时间序列|location=New York |publisher=Wiley |edition=Second |year=2004 |isbn=978-0-471-23065-6 |pages=[https://book.douban.com/subject/1480917/] |url-access=registration |url=https://book.douban.com/subject/1480917/}}
* {{cite book |last=Gujarati |first=Damodar N. |last2=Porter |first2=Dawn C. |title=Basic Econometrics |location=New York |publisher=McGraw-Hill |year=2009 |edition=Fifth international |isbn=978-007-127625-2 |chapter=Causality in Economics: The Granger Causality Test 经济学中的因果关系：格兰杰因果检验|pages=652–658 }}
* {{cite book |last=Hoover |first=Kevin D. |authorlink= |chapter=Granger-causality |pages=[https://archive.org/details/newclassicalmacr0000hoov/page/168 168–176] |title=The New Classical Macroeconomics 新古典宏观经济学|location=Oxford |publisher=Basil Blackwell |year=1988 |isbn=978-0-631-14605-6 |url=https://archive.org/details/newclassicalmacr0000hoov/page/168 }}
* {{cite encyclopedia |year= 2008 |title = Granger–Sims causality |last= Kuersteiner|first=Guido|encyclopedia=The New Palgrave Dictionary of Economics 新帕尔格雷夫经济学大词典 |url=https://book.douban.com/subject/1481407/}}
*Kleinberg, S. and Hripcsak, G. (2011) [http://people.dbmi.columbia.edu/samantha/papers/kleinberg_jbi11_preprint.pdf "生物医学信息学因果推理综述"] 于2012年4月30日[https://web.archive.org/web/20120430094823/http://people.dbmi.columbia.edu/samantha/papers/kleinberg_jbi11_preprint.pdf 存档]于[https://en.wikipedia.org/wiki/Wayback_Machine Wayback Machine] J. Biomed Informatics
==编者推荐==
[[File:GLJ1.jpg|200px|缩略图|右|格兰杰计量经济学文集封面]]
*[https://book.douban.com/subject/2369017/ 格兰杰计量经济学文集]
::格兰杰是最具创造力的现代计量经济学家之一，他对帮助我们理解因果关系、建模方法、非平稳性、季节性和预测做出了重要的贡献。一些被广泛应用的时间序列计量经济学概念就来自于他的研究成果，这是因为格兰杰善于将观测到的经济结果系统地用模型来表述，而且这些模型很好地反映了经济体复杂和进化的本质，这本书收录了他发表过的作品，而且把主要主题下相关联的研究集合到一起，这样便更容易找到相关的论文，里面的很多文章值得我们认真反复研读。
----
本中文词条由[[用户:Yillia Jing|Yillia Jing|]]参与编译，[[用户:苏格兰|苏格兰]]审校，[[用户:乐多多|乐多多]]编辑，欢迎在讨论页面留言。
'''本词条内容源自wikipedia及公开资料，遵守 CC3.0协议。'''
[[Category:多元时间序列]]
[[Category:时间序列统计检验]]