格蘭傑因果關係檢驗是一種假設檢定的統計方法，檢驗一組時間序列 x 是否為另一組時間序列 y 的原因。它的基礎是迴歸分析當中的自迴歸模型。迴歸分析通常只能得出不同 變數間的同期 相關性；自迴歸模型只能得出同一 變數前後期 的相關性；但諾貝爾經濟學獎得主克萊夫格蘭傑（Clive Granger）於 1969 年論證 ，在自迴歸模型中透過一系列的檢定進而揭示不同變數之間的時間落差相關性是可行的。 格蘭傑本人在其 2003 年獲獎演說中強調了其引用的侷限性，以及 很多荒謬論文的出現（Of course many ridiculous papers appeared）。格蘭傑因果關係檢驗的結論只是一種統計估計，不是真正意義上的因果關係，不能作為肯定或否定因果關係的根據。同時，格蘭傑因果關係檢驗也有一些不足之處，如並未考慮幹擾因素的影響，也未考慮時間序列間非線性的相互關係。一些基於格蘭傑因果關係檢驗的方法一定程度上解決了這些問題。 核心概念 過去值（lag value，或稱滯後值）：同一變項比當期時間上更早的值。例如：當期為 y_10，它的滯後值為 y_i10。 格蘭傑因果關係檢驗的基本觀念在於：未來的事件不會對目前與過去產生因果影響，而過去的事件才可能對現在及未來產生影響。 也就是說，如果我們試圖探討變數 x 是否對變數 y 有因果影響，那麼只需要估計 x 的滯後值是否會影響 y 的現在值，因為 x 的未來值不可能影響 y 的現在值。假如在控制了 y 變數的過去值以後，x 變數的過去值仍能對 Y 變數有顯著的解釋能力，我們就可以稱 x 能「Granger 影響」（Grangercause）y。 侷限性和改進 最初版的格蘭傑因果測試，有時候無法發現真正的因果關係。因為雖然對於認定因果關係而言，理論上還必須控制其他可能的幹擾因素，但在 Granger 最初提出這套因果測試的版本中，並未納入幹擾變數的分析，而是假設其他可能解釋變數的資訊包含在 y 的落後值中。如果事實上帶來因果關係的是第三變數（幹擾變數），亦即若事實上操控 x 並無法改變 y ，格蘭傑因果關係的零假設仍然可能被拒絕。因此標準版的格蘭傑因果測試結果可能會產生誤導性。 1980 年代由其他的計量經濟學家對 Granger 測試加以修改、擴充，將可能的第三（以上）變數納入測試，成為使用的向量自迴歸模型。相較於最初版的 Granger 測試，擴充版可以產生更有效的估計結果。 步驟描述 準備工作：一開始要用幾個落後期來建立模型，需要研究者的評估，通常使用 赤池資訊量準則 或 來判斷。 格蘭傑因果關係檢驗的第一步是建立用 y 的滯後值來預測 y 的自迴歸模型。此際，如果時間序列 y 是廣義平穩的，則可以直接使用滯後值。如果不平穩，就必須對不平穩的時間序列先做（一階或更多階）差分，直到得出平穩時間數列。 如果發現 y 的某期滯後值 1 在迴歸分析中具有顯著性（根據 t 檢定的 p 值來判斷），且 2 這期滯後值加入模型後可提高迴歸模型的解釋力（根據迴歸分析的 F 檢定），這個滯後值便被留在模型中。 然後進一步加入 x（或 Δx）的滯後期來擴充迴歸模型。關於平穩時間序列的要求、某期滯後值留在模型中的條件，同上述 y 的處理。 若且唯若（充分必要）沒有任何解釋變項 x（或 Δx）的滯後值被留在模型中，便無法拒絕無格蘭傑因果關係的零假設。 研究人員希望發現明顯的證據，比如 x 是 y 的格蘭傑原因但反之不成立，便能做出因果關係的推論。然而在實際操作中也可能會發現沒有變數是對方的格蘭傑原因，或者 x 和 y 兩個變數互為格蘭傑原因。 數學定義 1 令 x 和 y 為廣義平穩序列。如要檢測 x 非 y 的格蘭傑原因之零假設，首先引入 y 的滯後值建立 y 的自迴歸模型（AR model on y）： y_t a_0 a_1y_t1 a_2y_t2 a_my_tm residual_t 所有的 y 滯後值中：1 在迴歸分析中具有顯著性（根據 t 統計值的 p 值來判斷）的，且 2 這期滯後值加入模型後可提高迴歸模型的解釋力（根據迴歸分析的 F 檢定）的－－將被留在模型中。m 表示的是 y 變數滯後期中檢定為顯著的時間上最早一個。 2 接著，引入 x 的滯後值建立增廣迴歸模型： y_t a_0 a_1y_t1 a_2y_t2 a_my_tm b_px_tp b_qx_tq residual_t 所有的 x 滯後值中：1 在迴歸分析中具有顯著性（根據 t 檢驗的 p 值來判斷）的，且 2 這期滯後值加入模型後可提高迴歸模型的解釋力（根據迴歸分析的 F 檢定）的－－將被留在模型中。在以上增廣迴歸模型中，p 代表 x 變數滯後值中檢定為顯著的時間上最早一個，q 則是 x 變數滯後值中檢定為顯著的時間上最近一個。 3 如果沒有任何 x 的滯後值被留在模型中，無格蘭傑因果關係的零假設就成立。 軟體執行 一些統計軟體可以執行 Granger causality test。例如：Stata、SPSS、EViews、R 語言。 這裡舉個 R 語言中 lmtest 程式庫裡 grangertest 指令的例子： Granger causality test Model 1 fii Lags fii 15 Lags rM 15 Model 2 fii Lags fii 15 ResDf Df F Pr F 1 629 2 634 5 25115 002896 Signif codes 0 0001 001 005 01 1 Granger causality test Model 1 rM Lags rM 15 Lags fii 15 Model 2 rM Lags rM 15 ResDf Df F Pr F 1 629 2 634 5 11804 03172 模型 1 檢驗將滯後的 rM 從解釋 FII 的回歸模型中移除是否可行，答案是不可行的（因為 p 值 002896）。但由模型 1 和模型 2 的組合可發現從解釋 rM 的模型中移除 FII 的滯後值是可能的。我們可以由此斷定 rM 是 FII 的格蘭傑原因，反之則不成立。 延伸 承繼著回歸模型的基本性質，格蘭傑因果關係分析也假設實際值與預測值之間的誤差呈常態分佈，若實際現象不呈常態分佈將嚴重影響推論的有效性。 Hacker HatemiJ 2006 發展出一種不必在乎誤差項是否呈正態分佈的格蘭傑因果關係研究方法 Hacker RS and HatemiJ A 2006 Tests for causality between integrated variables using asymptotic and bootstrap distributions theory and application Applied Economics Vol 38 13 pp 14891500。這種方法在財金分析上特別實用， 因為許多金融變數不服從常態分佈。 近來，Hacker HatemiJ 2012 又進一步改善之，提出一種非對稱的因果關係檢驗模型，據說可以區分正向與負向影響的因果影響。 時變格蘭傑因果關係：將格蘭傑因果關係擴充套件到動態的時變性質，可以更細緻地瞭解時間序列資料中的因果關係是如何隨時間演變的。 該方法使用遞迴技術，如前向展（FE）、滾動（RO）和遞迴演化（RE）視窗，以克服傳統格蘭傑因果檢驗的侷限性，並理解不同時期因果關係的變化。 這一方法的核心是 Stata 中的 tvgc 命令。經驗應用，如涉及以太坊交易費用和經濟子系統的資料，突出了經濟關係隨時間變化的動態性質。 註腳 相關條目 自迴歸模型（AR model） 向量自迴歸模型（VAR model） 差分自回歸滑動平均模型（ARIMA 模型） Category 時間序列 Category 因果律 Category 計量經濟學 Category 統計學 Category 迴歸分析 deKausalitätÖkonometrie Granger