查看“︁格蘭傑因果關係”︁的源代码

{{noteTA
|G1=Math
|1=zh:軟件;zh-cn:软件;zh-tw:軟體;
|4=zh:滯後值;zh-tw:落後期;zh-hant:落後期;zh-hans:落後期;
|5=zh:追蹤資料;zh-cn:面板資料;
}}

'''格蘭傑因果關係檢驗'''（{{lang-en|Granger causality test}}）是一種[[假設檢定]]的統計方法，檢驗一組[[時間序列]]<math>x</math>是否為另一組時間序列<math>y</math>的原因。它的基礎是[[迴歸分析]]當中的[[自迴歸模型]]。迴歸分析通常只能得出''不同'' [[變量]]間的''同期'' [[相關性]]；自迴歸模型只能得出''同一'' 變量''前後期'' 的相關性；但[[諾貝爾經濟學獎]]得主[[克萊夫·格蘭傑]]（Clive Granger）於1969年論證<ref name="Granger">{{cite journal | title=Investigating Causal Relations by Econometric Models and Cross-spectral Methods | first=Clive W. J. | last=Granger | year=1969 | journal={{tsl|en|Econometrica|計量經濟學 (期刊)|Econometrica}} | volume=37 | issue=3 | pages=424-438 | publisher=[[經濟計量學會|The Econometric Society]] | url=http://www.jstor.org/stable/1912791 | access-date=2015-04-20 | archive-date=2016-03-06 | archive-url=https://web.archive.org/web/20160306202011/http://www.jstor.org/stable/1912791 | dead-url=no }}</ref> ，在自迴歸模型中透過一系列的檢定進而揭示不同變量之間的時間落差相關性是可行的。

格蘭傑本人在其2003年獲獎演說中強調了其引用的局限性，以及“很多荒謬論文的出現”（Of course, many ridiculous papers appeared）。格蘭傑因果關係檢驗的結論只是一種統計估計，不是真正意義上的[[因果關係]]，不能作爲肯定或否定因果關係的根據。同時，格蘭傑因果關係檢驗也有一些不足之處，如並未考慮[[干擾因素]]的影響，也未考慮時間序列間非線性的相互關係。一些基於格蘭傑因果關係檢驗的方法一定程度上解決了這些問題<ref name=LiuBah>{{cite web | url= http://www-bcf.usc.edu/~liu32/sdm_theory.pdf | title= An Examination of Practical Granger Causality Inference (Manuscript) | last1= Bahadori | first1= Mohammad Tahalast | last2= Liu | first2= Yan | year= 2014 | website= University of Southern California | accessdate= 2015-04-20 | archive-date= 2015-04-27 | archive-url= https://web.archive.org/web/20150427143132/http://www-bcf.usc.edu/~liu32/sdm_theory.pdf | dead-url= no }}</ref>。

==核心概念==
'''過去值'''（lag value，或稱滯後值）：同一變項比當期'''時間上'''更早的值。例如：當期為<math>y_{10}</math>，它的滯後值為<math>y_{i<10}</math>。

格蘭傑因果關係檢驗的基本觀念在於：'''未來的事件不會對目前與過去產生因果影響，而過去的事件才可能對現在及未來產生影響'''。<ref name="Granger"></ref><ref>{{cite journal|url=http://www.scholarpedia.org/article/Granger_causality|author=Anil Seth|year=2007|title=Granger causality|work=Scholarpedia|volume=2|issue=7|page=1667|doi=10.4249/scholarpedia.1667|access-date=2011-03-07|archive-date=2011-08-10|archive-url=https://web.archive.org/web/20110810082322/http://www.scholarpedia.org/article/Granger_causality|dead-url=no}}</ref>
也就是說，如果我們試圖探討變數<math>x</math>是否對變數<math>y</math>有因果影響，那麼只需要估計<math>x</math>的滯後值是否會影響<math>y</math>的現在值，因為<math>x</math>的未來值不可能影響<math>y</math>的現在值。'''假如在控制了<math>y</math>變數的過去值以後，<math>x</math> 變數的過去值仍能對Y 變數有顯著的解釋能力'''，我們就可以稱<math>x</math>能「Granger 影響」（Granger-cause）<math>y</math>。<ref name="police">{{cite journal | title=警力增加能導致竊盜犯罪率降低嗎？ - 台灣縣市1998-2007 動態追蹤資料之Granger 因果分析 | author=劉孟奇，張其祿，盧敬植 | journal=公共行政學報 | issue=34 | year=2010 | pages=1–27 | url=http://pa.nccu.edu.tw/download.php?filename=6_4ef12231.pdf&dir=publish&title=%E5%85%A8%E6%96%87%E4%B8%8B%E8%BC%89 | format=[[PDF]] | accessdate=2015-04-20 | archive-date=2016-03-04 | archive-url=https://web.archive.org/web/20160304070909/http://pa.nccu.edu.tw/download.php?filename=6_4ef12231.pdf&dir=publish&title=%E5%85%A8%E6%96%87%E4%B8%8B%E8%BC%89 | dead-url=no }}</ref>

===侷限性和改進===
最初版的格蘭傑因果測試，有時候無法發現真正的因果關係。因為雖然對於認定因果關係而言，理論上還必須控制其他可能的[[偽關係|干擾因素]]，但在 Granger 最初提出這套因果測試的版本中，並未納入干擾變數的分析，而是假設其他可能解釋變數的資訊包含在<math>y</math>的落後值中。如果事實上帶來因果關係的是第三變數（干擾變數），亦即若事實上操控 <math>x</math> 並無法改變 <math>y</math> ，格蘭傑因果關係的[[零假設]]仍然可能被拒絕。因此標準版的格蘭傑因果測試結果可能會產生誤導性。

1980年代由其他的計量經濟學家對Granger測試加以修改、擴充，將可能的第三（以上）變數納入測試，成為使用{{link-en|追蹤資料|panel data}}的[[向量自迴歸模型]]（{{Lang-en|panel data VAR model}}）。相較於最初版的 Granger 測試，擴充版可以產生更有效的估計結果。<ref name="police"></ref>

==步驟描述==
#準備工作：一開始要用幾個落後期來建立模型，需要研究者的評估，通常使用 [[赤池信息量準則]]（{{Lang-en|Akaike information criterion、簡稱AIC}}） 或 {{link-en|貝氏訊息量準則|Bayesian information criterion}}（{{Lang-en|Bayesian information criterion、簡稱BIC}}）來判斷。
#格蘭傑因果關係檢驗的第一步是建立用<math>y</math>的滯後值來預測<math>y</math>的自迴歸模型。此際，如果時間序列<math>y</math>是[[平穩過程|廣義平穩]]的，則可以直接使用滯後值。如果不平穩，就必須對不平穩的時間序列先做（一階或更多階）[[差分]]，直到得出平穩時間數列。
#如果發現 <math>y</math> 的某期滯後值 (1) 在迴歸分析中具有顯著性（根據 [[學生t檢驗|t檢定]]的[[p值]]來判斷），且 (2) 這期滯後值加入模型後可提高迴歸模型的解釋力（根據迴歸分析的[[迴歸分析|F檢定]]），這個滯後值便被留在模型中。
#然後進一步加入<math>x</math>（或Δ<math>x</math>）的滯後期來擴充迴歸模型。關於平穩時間序列的要求、某期滯後值留在模型中的條件，同上述<math>y</math>的處理。
#若且唯若（充分必要）沒有任何解釋變項<math>x</math>（或Δ<math>x</math>）的滯後值被留在模型中，便無法拒絕無格蘭傑因果關係的[[零假設]]。

研究人員希望發現明顯的證據，比如<math>x</math>是<math>y</math>的格蘭傑原因但反之不成立，便能做出因果關係的推論。然而在實際操作中也可能會發現沒有變量是對方的格蘭傑原因，或者<math>x</math>和<math>y</math>兩個變量互為格蘭傑原因。

===數學定義===
1. 令<math>x</math>和<math>y</math>為[[平穩過程#廣義平穩（弱平穩）|廣義平穩]]序列。如要檢測<math>x</math>非<math>y</math>的格蘭傑原因之[[零假設]]，首先引入<math>y</math>的滯後值建立<math>y</math>的[[自迴歸模型]]（AR model on <math>y</math>）：
:<math>y_t = a_0 + a_1y_{t-1} + a_2y_{t-2} + ... + a_my_{t-m} + residual_t.</math>

:所有的<math>y</math>滯後值中：(1) 在迴歸分析中具有顯著性（根據[[學生t檢驗|t-統計值]]的[[p值]]來判斷）的，且 (2) 這期滯後值加入模型後可提高迴歸模型的解釋力（根據迴歸分析的[[迴歸分析|F檢定]]）的－－將被留在模型中。<math>m</math>表示的是<math>y</math>變量滯後期中檢定為顯著的時間上最早一個。

2. 接著，引入<math>x</math>的滯後值建立增廣迴歸模型：
:<math>y_t = a_0 + a_1y_{t-1} + a_2y_{t-2} + ... a_my_{t-m} + b_px_{t-p} + ... + b_qx_{t-q} + residual_t.</math>

:所有的<math>x</math>滯後值中：(1) 在迴歸分析中具有顯著性（根據[[t檢驗]]的[[p值]]來判斷）的，且 (2) 這期滯後值加入模型後可提高迴歸模型的解釋力（根據迴歸分析的[[迴歸分析|F檢定]]）的－－將被留在模型中。在以上增廣迴歸模型中，<math>p</math>代表<math>x</math>變量滯後值中檢定為顯著的時間上最早一個，<math>q</math>則是<math>x</math>變量滯後值中檢定為顯著的時間上最近一個。 

3. 如果'''沒有任何<math>x</math>的滯後值'''被留在模型中，無格蘭傑因果關係的零假設就成立。

===軟件運行===
一些統計軟體可以執行Granger causality test。例如：[[Stata]]、[[SPSS]]<ref>{{cite web|url=http://answers.google.com/answers/threadview/id/27621.html|title=Q: Granger causality test|work=Google Answers|accessdate=2013-03-08|archive-date=2013-05-27|archive-url=https://web.archive.org/web/20130527110100/http://answers.google.com/answers/threadview/id/27621.html|dead-url=no}}</ref>、[[EViews]]<ref>{{cite web|url=http://www.youtube.com/watch?v=YOS4nOmAteE|title=Eviews. Granger. CAUSALITY BETWEEN MONEY AND INTEREST RATE IN CANADA|work=YouTube|accessdate=2013-03-08|archive-date=2016-09-24|archive-url=https://web.archive.org/web/20160924144555/https://www.youtube.com/watch?v=YOS4nOmAteE|dead-url=no}}</ref>、[[R語言]]。

這裡舉個R語言中lmtest程序庫裡grangertest()指令的例子：

<pre>
Granger causality test

Model 1: fii ~ Lags(fii, 1:5) + Lags(rM, 1:5)
Model 2: fii ~ Lags(fii, 1:5)
  Res.Df  Df      F  Pr(>F)
1    629
2    634   5 2.5115 0.02896 *
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Granger causality test

Model 1: rM ~ Lags(rM, 1:5) + Lags(fii, 1:5)
Model 2: rM ~ Lags(rM, 1:5)
  Res.Df  Df      F Pr(>F)
1    629
2    634   5 1.1804 0.3172
</pre>

模型1檢驗將滯後的rM從解釋FII的回歸模型中移除是否可行，答案是不可行的（因為p值 = 0.02896）。但由模型1和模型2的組合可發現從解釋rM的模型中移除FII的滯後值是可能的。我們可以由此斷定rM是FII的格蘭傑原因，反之則不成立。

==延伸==
承繼著回歸模型的基本性質，格蘭傑因果關係分析也假設實際值與預測值之間的誤差呈常態分佈，若實際現象不呈常態分佈將嚴重影響推論的有效性。

Hacker & Hatemi-J (2006)發展出一種不必在乎誤差項是否呈正態分佈的格蘭傑因果關係研究方法<ref>Hacker R.S. and Hatemi-J A. (2006) [http://ideas.repec.org/a/taf/applec/v38y2006i13p1489-1500.html "Tests for causality between integrated variables using asymptotic and bootstrap distributions: theory and application"] {{Wayback|url=http://ideas.repec.org/a/taf/applec/v38y2006i13p1489-1500.html |date=20110913145837 }}, Applied Economics, Vol. 38(13), pp. 1489–1500.</ref>。這種方法在財金分析上特别實用， 因為許多金融變量不服從常態分佈。<ref>{{cite journal |title=The variation of certain speculative prices |year=1963 |last1=Mandelbrot |first1=Benoit| journal=Journal of Business|volume=36 |issue=1 |pages=394–419 |doi=10.1086/294632}}</ref>

近來，Hacker & Hatemi-J (2012)又進一步改善之，提出一種非對稱的因果關係檢驗模型，據說可以區分正向與負向影響的因果影響。<ref>{{cite journal |doi=10.1007/s00181-011-0484-x |title=Asymmetric causality tests with an application |year=2012 |last1=Hatemi-J |first1=A. |journal=Empirical Economics |volume=42 |issue=6 |pages=forthcoming}}</ref>

时变格兰杰因果关系：将格兰杰因果关系扩展到动态的时变性质，可以更细致地了解时间序列数据中的因果关系是如何随时间演变的。<ref name=":0">{{Cite journal |last=Baum |first=Christopher F. |last2=Hurn |first2=Stan |last3=Otero |first3=Jesús |date=2022-06 |title=Testing for time-varying Granger causality |url=http://journals.sagepub.com/doi/10.1177/1536867X221106403 |journal=The Stata Journal: Promoting communications on statistics and Stata |language=en |volume=22 |issue=2 |doi=10.1177/1536867X221106403 |issn=1536-867X |access-date=2024-01-05 |archive-date=2024-01-05 |archive-url=https://web.archive.org/web/20240105140715/https://journals.sagepub.com/doi/10.1177/1536867X221106403 |dead-url=no }}</ref> 该方法使用递归技术，如前向展（FE）、滚动（RO）和递归演化（RE）窗口，以克服传统格兰杰因果检验的局限性，并理解不同时期因果关系的变化。<ref>{{Cite journal |last=Shojaie |first=Ali |last2=Fox |first2=Emily B. |date=2022-03-07 |title=Granger Causality: A Review and Recent Advances |url=https://www.annualreviews.org/doi/10.1146/annurev-statistics-040120-010930 |journal=Annual Review of Statistics and Its Application |language=en |volume=9 |issue=1 |doi=10.1146/annurev-statistics-040120-010930 |issn=2326-8298 |pmc=10571505 |pmid=37840549 |access-date=2024-01-05 |archive-date=2022-11-09 |archive-url=https://web.archive.org/web/20221109191649/https://www.annualreviews.org/doi/10.1146/annurev-statistics-040120-010930 |dead-url=no }}</ref> 这一方法的核心是 Stata 中的 "tvgc "命令。<ref name=":0" />经验应用，如涉及以太坊交易费用和经济子系统的数据，突出了经济关系随时间变化的动态性质。<ref>{{Cite journal |last=Ante |first=Lennart |last2=Saggu |first2=Aman |date=2024-01 |title=Time-Varying Bidirectional Causal Relationships between Transaction Fees and Economic Activity of Subsystems Utilizing the Ethereum Blockchain Network |url=https://www.mdpi.com/1911-8074/17/1/19 |journal=Journal of Risk and Financial Management |language=en |volume=17 |issue=1 |doi=10.3390/jrfm17010019 |issn=1911-8074}}</ref>

==註腳==
<references/>

==相關條目==
*[[自迴歸模型]]（AR model）
*[[向量自迴歸模型]]（VAR model）
*[[ARIMA模型|差分自回歸滑動平均模型]]（ARIMA模型）

{{DEFAULTSORT:Granger Causality}}

[[Category:时间序列]]
[[Category:因果律]]
[[Category:計量經濟學]]
[[Category:統計學]]
[[Category:迴歸分析]]

[[de:Kausalität#Ökonometrie (Granger)]]