🔖 キーワード索引（補強）

Granger 因果性に関連する時系列概念のチップ集。

🔖 キーワード索引 — 完全強化版

「グレンジャー因果」を理解するうえで必要なキーワードを 10 件以上提示します。 各チップから対応セクションへ移動できます。

30 秒結論 文脈 直感 数式 記号読み解き 実値計算 Python 実装 落とし穴 関連手法 関連用語 グループ教材 概念マップ

🎨 直感で掴む — 完全強化版

グレンジャー因果 を一言でいえば「人口 A1101 と出生数 A4101 の時系列でグレンジャー因果検定」。 47 都道府県という小さな母集団でも、 SSDSE-B-2026 の A1101 列に注目すると、 大都市圏と地方の差・人口規模に伴う相対比較など、 様々なパターンが見えてきます。

比喩でいうと、 グレンジャー因果 はデータ分析の「眼鏡」のようなもの。 同じデータでも眼鏡を変えれば、 平均（中心）・分散（ばらつき）・相関（連動）・因果（影響）と、 異なる情報が浮かび上がります。 SSDSE-B-2026 を題材に、 この眼鏡をかけてみるのが本ページの狙いです。

• 具体例 1：47 都道府県の A1101 の中央値 ± 四分位を見ると、 規模感が掴める。
• 具体例 2：A4101 と組み合わせて散布図にすると、 群構造が見える。
• 具体例 3：政令指定都市の有無で群分けすると、 介入研究の素材になる。

🎮 触って理解する

2 つの時系列 X（青） と Y（オレンジ） を表示します。 Y は「自分の過去」と「X の過去」から生成されています。 スライダーで X の過去が Y の将来に効く強さ と ラグ数 p を変えると、 下の 2 本の棒が更新されます。左は Y を自分の過去だけで予測した誤差（RSS_制限）、 右は X の過去も加えて予測した誤差（RSS_非制限）。 後者が 有意に 小さければ「X は Y をグレンジャー予測する」と判定されます。 グラフ上を左右にドラッグ（スマホはタッチ）しても強さを変えられます。

⚠️ 予測改善であって真の因果ではありません。 ここで検定しているのは「X の過去情報が Y の予測誤差を減らすか」だけです。 共通の先行変数（交絡）や 共通トレンド があると、 実際には因果が無くてもグレンジャー予測が有意になります。 「グレンジャー因果あり」は「時間的な先行・予測寄与あり」と読み替えるのが安全です。

💡 直感 — 過去の情報は将来予測を改善するか

グレンジャー因果の核心は「情報」です。今の Y を予測するのに Y 自身の過去 だけで十分なら、 X は Y に対して情報を持ちません。ところが X の過去を足したとたんに予測誤差がガクッと下がるなら、 「X には Y の未来を先取りする情報が含まれている」＝ X は Y をグレンジャー予測する、と考えます。 上の実験で β を上げると RSS_非制限 が RSS_制限 より小さくなり、F 統計量が大きく、p 値が小さくなる様子を確認してください。

⚠️ よくある落とし穴 — 予測 ≠ 因果

• 予測改善 = 因果ではない：グレンジャー因果は「時間的先行と予測寄与」であり、 介入して結果が変わる意味の真の因果とは別物です。
• 交絡（共通原因）：第 3 の変数 Z が X と Y の両方を先に動かすと、 X→Y のグレンジャー因果が偽陽性で出ます。交絡の検討が必須です。
• 共通トレンド・非定常：両系列が上昇トレンドを共有するだけで 見せかけの相関／見せかけの因果が生じます。差分化や定常性の確認が前提です。

🚀 発展 — VAR・インパルス応答・非定常の注意

• VAR（ベクトル自己回帰）：複数系列を同時にモデル化する枠組み。 グレンジャー因果検定は VAR の係数に対するゼロ制約検定として実装されます。
• インパルス応答関数（IRF）：ある系列へのショックが他系列へ時間を追ってどう波及するかを可視化。 「向き」だけでなく「大きさと持続」を読めます。
• 非定常の注意：単位根がある系列では検定統計量の分布が変わり偽陽性が増えます。 ADF/KPSS 検定 →（必要なら）差分化、共和分があれば VECM、という時系列の定石に従います。 ラグ変数の取り方（p の選択）は AIC/BIC で決めます。

📐 数式または定義 — 完全強化版

グレンジャー因果 の代表的な定義式は次のとおりです。

ここで使われる記号や演算の意味は次節で言葉に翻訳します。

🔬 数式を言葉で読み解く — 完全強化版

数式の各記号を、日本語の意味に変換します。

• X, Y — 観測されたデータ系列（SSDSE-B-2026 の A1101, A4101 列など）
• n — 標本サイズ（47 都道府県データなら 47、 多年度パネルなら 47×年数）
• β, θ, α — 推定すべき未知パラメータ
• ε, u — 観測されない誤差項・撹乱項
• E[·] — 期待値（母集団平均）
• Cov, Var — 共分散・分散
• arg min — 「最小化するパラメータ」の意味
• L (ラグ演算子) — 時系列で 1 期前の値を返す操作

🧮 SSDSE-B-2026 で実値計算 — Granger 因果性

SSDSE-B-2026 を年度パネルとして並べ替えれば、 都道府県の時系列を作ることができます。 ここでは年次データを擬似時系列として Granger 因果性を計算します。

例：2変量で Granger 因果性検定

📤 出力の読み方
# grangercausalitytests は maxlag までの各ラグについて
# ssr based F test / likelihood ratio test の p 値を返す。
# あるラグで p < 0.05 なら「X の過去は Y の予測を改善する（Granger 因果あり）」。
# 逆向き（Y→X）も検定し、因果の非対称性を確認する。

💬 読み方：算出された統計量を判定基準と比較し、有意性／効果量を評価する。

🎯 このコードでやること：Granger 因果性 — 過去値の予測寄与に基づく因果推論に関連するステップ #4。主要な指標（係数・統計量・スコア）を算出します。

📥 入力の形
# Granger 検定に渡すのは「定常な2系列」を並べた DataFrame。
# 列は説明変数候補 X（原因側）と被説明変数 Y（結果側）の2列で足りる。
# 実データが非定常なら df.diff().dropna() で差分化してから渡す。
#      X        Y
# 0  ...      ...

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13

import pandas as pd
from statsmodels.tsa.stattools import grangercausalitytests

df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', header=1)
# 東京都の時系列を抽出（年度列が存在する想定）
tokyo = df[df['都道府県'] == '東京都'].sort_values(df.columns[0])
# 年度列は使わず、人口 A1101（総人口）と出生数 A4101 の2系列を使う
ts = tokyo[['総人口', '出生数']].dropna()

# Granger 因果性：列1 が 列2 を予測できるか（最大 3 ラグ）
res = grangercausalitytests(ts, maxlag=3, verbose=False)
for lag, r in res.items():
    print(f'lag={lag}: F={r[0]["ssr_ftest"][0]:.3f}, p={r[0]["ssr_ftest"][1]:.4f}')

例：VAR モデルでの因果性

📤 出力の読み方
# grangercausalitytests は maxlag までの各ラグについて
# ssr based F test / likelihood ratio test の p 値を返す。
# あるラグで p < 0.05 なら「X の過去は Y の予測を改善する（Granger 因果あり）」。
# 逆向き（Y→X）も検定し、因果の非対称性を確認する。

💬 読み方：p 値や信頼区間と合わせて読み、効果の有無＋大きさを両輪で判断する。

🎯 このコードでやること：Granger 因果性 — 過去値の予測寄与に基づく因果推論に関連するステップ #5。仮説検定・モデル評価を行います。

📥 入力の形
# Granger 検定に渡すのは「定常な2系列」を並べた DataFrame。
# 列は説明変数候補 X（原因側）と被説明変数 Y（結果側）の2列で足りる。
# 実データが非定常なら df.diff().dropna() で差分化してから渡す。
#      X        Y
# 0  ...      ...

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13

import pandas as pd
from statsmodels.tsa.api import VAR

df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', header=1)
tokyo = df[df['都道府県'] == '東京都'].sort_values(df.columns[0])
cols = ['総人口', '出生数']  # 人口 A1101 と出生数 A4101
ts = tokyo[cols].dropna().diff().dropna()  # 1階差で定常化

model = VAR(ts)
selection = model.select_order(maxlags=2)  # 年12点の短系列なので小さめに
print('推奨ラグ (AIC/BIC):', selection.aic, selection.bic)
fit = model.fit(maxlags=2)
print(fit.test_causality(cols[0], [cols[1]], kind='f'))

🧮 実値で計算してみる — SSDSE-B-2026 で グレンジャー因果（完全強化版）

SSDSE-B-2026（公的統計の社会・教育系データセット、 47 都道府県 × 10 年分超 × 100 以上の列）を用いて、 「グレンジャー因果」を体感します。 ファイル名は SSDSE-B-2026.csv、 読み込みは下記の Python コードで行います。

1
2
3
4
5
6
7
8

import pandas as pd

# SSDSE-B-2026 を読み込む（cp932 / Shift_JIS）
df = pd.read_csv('data/raw/SSDSE-B-2026.csv', skiprows=[1], encoding='cp932')
print(df.shape)          # (564, 112)
print(df['SSDSE-B-2026'].unique())  # 含まれる年度
latest = df[df['SSDSE-B-2026'] == df['SSDSE-B-2026'].max()].copy()
print(latest[['Prefecture', 'A1101', 'A4101']].head())

ここで使った中心列 A1101 は SSDSE-B-2026 における 人口 A1101 と出生数 A4101 の時系列でグレンジャー因果検定 に関連する指標です。 算出例：

• 47 都道府県・最新年度の A1101 平均と標準偏差を求める
• A1101 と A4101 の相関（線形・順位）を比較する
• 群分け（例：人口上位 10 県 vs それ以外）で平均差を見る

🧮 数式に値を入れて手で計算する: AR と AR+X モデルの RSS 比較

合成データで X→Y の Granger 検定統計量を概算する。

Step 1: モデル別残差平方和

n=100 サンプル

Step 2: F 統計量

🐍 Python で再現

 1
 2
 3
 4
 5
 6
 7
 8
 9
10

from scipy import stats
rss_r, rss_u = 50, 40
n = 100
k_u = 4
q = 2
F = ((rss_r - rss_u)/q) / (rss_u/(n - k_u))
crit = stats.f.ppf(0.95, q, n - k_u)
print(f"F = {F:.2f}")
print(f"臨界値: {crit:.3f}")
print(f"Granger 因果あり: {F > crit}")

📤 実行結果

💬 手計算 (Step 2) F=12.0 と Python 出力が完全一致。

🐍 Python 実装バリエーション（statsmodels / scipy / arch / pmdarima）

1. statsmodels.tsa — Granger / VAR 標準実装

📤 出力の読み方
# grangercausalitytests は maxlag までの各ラグについて
# ssr based F test / likelihood ratio test の p 値を返す。
# あるラグで p < 0.05 なら「X の過去は Y の予測を改善する（Granger 因果あり）」。
# 逆向き（Y→X）も検定し、因果の非対称性を確認する。

💬 読み方：表示された数値テーブルから個別の都道府県の位置づけを読み取る。

🎯 このコードでやること：Granger 因果性 — 過去値の予測寄与に基づく因果推論に関連するステップ #6。結果を整形して表示します。

📥 入力の形
# Granger 検定に渡すのは「定常な2系列」を並べた DataFrame。
# 列は説明変数候補 X（原因側）と被説明変数 Y（結果側）の2列で足りる。
# 実データが非定常なら df.diff().dropna() で差分化してから渡す。
#      X        Y
# 0  ...      ...

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15

import pandas as pd
from statsmodels.tsa.stattools import grangercausalitytests, adfuller
from statsmodels.tsa.api import VAR

df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', header=1)
tokyo = df[df['都道府県'] == '東京都'].sort_values(df.columns[0])
ts = tokyo[['総人口', '出生数']].dropna()

# まず定常性チェック
for c in ts.columns:
    stat, p = adfuller(ts[c])[:2]
    print(f'{c}: ADF stat={stat:.3f}, p={p:.4f}')

# Granger
res = grangercausalitytests(ts.diff().dropna(), maxlag=3, verbose=False)

2. scipy.stats — F検定で簡易 Granger

📤 出力の読み方
# grangercausalitytests は maxlag までの各ラグについて
# ssr based F test / likelihood ratio test の p 値を返す。
# あるラグで p < 0.05 なら「X の過去は Y の予測を改善する（Granger 因果あり）」。
# 逆向き（Y→X）も検定し、因果の非対称性を確認する。

💬 読み方：SSDSE-B-2026 の実値に当てはめると教科書例より分散が大きいことに注意。

🎯 このコードでやること：Granger 因果性 — 過去値の予測寄与に基づく因果推論に関連するステップ #7。47都道府県データに当てはめて確認します。

📥 入力の形
# Granger 検定に渡すのは「定常な2系列」を並べた DataFrame。
# 列は説明変数候補 X（原因側）と被説明変数 Y（結果側）の2列で足りる。
# 実データが非定常なら df.diff().dropna() で差分化してから渡す。
#      X        Y
# 0  ...      ...

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24

import numpy as np
import pandas as pd
from scipy import stats
import statsmodels.api as sm

df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', header=1)
tokyo = df[df['都道府県'] == '東京都'].sort_values(df.columns[0])
ts = tokyo[['総人口', '出生数']].dropna().values

# 制約モデルと非制約モデルの SSR を比較
def granger_f(y, x, lag=2):
    Y = y[lag:]
    X_rest = np.column_stack([y[lag-i-1:-i-1] for i in range(lag)])
    X_full = np.column_stack([Y]*0 + [y[lag-i-1:-i-1] for i in range(lag)] +
                             [x[lag-i-1:-i-1] for i in range(lag)])
    X_rest = sm.add_constant(X_rest); X_full = sm.add_constant(X_full)
    ssr_r = sm.OLS(Y, X_rest).fit().ssr
    ssr_f = sm.OLS(Y, X_full).fit().ssr
    n, k = len(Y), lag
    F = ((ssr_r - ssr_f) / k) / (ssr_f / (n - 2*k - 1))
    p = 1 - stats.f.cdf(F, k, n-2*k-1)
    return F, p

print(granger_f(ts[:, 0], ts[:, 1], lag=2))

3. arch — ブートストラップ Granger

📤 出力の読み方
# grangercausalitytests は maxlag までの各ラグについて
# ssr based F test / likelihood ratio test の p 値を返す。
# あるラグで p < 0.05 なら「X の過去は Y の予測を改善する（Granger 因果あり）」。
# 逆向き（Y→X）も検定し、因果の非対称性を確認する。

💬 読み方：別パターンと比べることで、手法選択の感度を体感できる。

🎯 このコードでやること：Granger 因果性 — 過去値の予測寄与に基づく因果推論に関連するステップ #8。比較・別パターンを検討します。

📥 入力の形
# Granger 検定に渡すのは「定常な2系列」を並べた DataFrame。
# 列は説明変数候補 X（原因側）と被説明変数 Y（結果側）の2列で足りる。
# 実データが非定常なら df.diff().dropna() で差分化してから渡す。
#      X        Y
# 0  ...      ...

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15

from arch.bootstrap import StationaryBootstrap
import numpy as np
import pandas as pd

df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', header=1)
tokyo = df[df['都道府県'] == '東京都'].sort_values(df.columns[0])
ts = tokyo[['総人口', '出生数']].dropna().values

def stat_func(data):
    # 簡易: 共分散を統計量に
    return np.cov(data[:, 0], data[:, 1])[0, 1]

bs = StationaryBootstrap(3, ts)
results = bs.apply(stat_func, 200)
print('ブートストラップ平均:', np.mean(results), '95%CI:', np.percentile(results, [2.5, 97.5]))

4. pmdarima — 自動 ARIMA で予測誤差比較

📤 出力の読み方
# grangercausalitytests は maxlag までの各ラグについて
# ssr based F test / likelihood ratio test の p 値を返す。
# あるラグで p < 0.05 なら「X の過去は Y の予測を改善する（Granger 因果あり）」。
# 逆向き（Y→X）も検定し、因果の非対称性を確認する。

💬 読み方：ハイパーパラメータで結果が大きく変わる場合は安定性を疑う。

🎯 このコードでやること：Granger 因果性 — 過去値の予測寄与に基づく因果推論に関連するステップ #9。ハイパーパラメータを変えて再計算します。

📥 入力の形
# Granger 検定に渡すのは「定常な2系列」を並べた DataFrame。
# 列は説明変数候補 X（原因側）と被説明変数 Y（結果側）の2列で足りる。
# 実データが非定常なら df.diff().dropna() で差分化してから渡す。
#      X        Y
# 0  ...      ...

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13

import pmdarima as pm
import pandas as pd

df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', header=1)
tokyo = df[df['都道府県'] == '東京都'].sort_values(df.columns[0])
ts = tokyo[['総人口', '出生数']].dropna()
y, x = ts.iloc[:, 0], ts.iloc[:, 1]

# 単独 ARIMA vs 外生変数つき ARIMA
mod1 = pm.auto_arima(y, seasonal=False)
mod2 = pm.auto_arima(y, exogenous=x.values.reshape(-1,1), seasonal=False)
print(f'BIC 単独: {mod1.bic():.2f}, 外生あり: {mod2.bic():.2f}')
print('外生変数で予測改善 →', mod2.bic() < mod1.bic())

5. 多変量 VAR と インパルス応答

📤 出力の読み方
# grangercausalitytests は maxlag までの各ラグについて
# ssr based F test / likelihood ratio test の p 値を返す。
# あるラグで p < 0.05 なら「X の過去は Y の予測を改善する（Granger 因果あり）」。
# 逆向き（Y→X）も検定し、因果の非対称性を確認する。

💬 読み方：最終結果は CSV／プロットとして保存しておくと後続分析で再利用できる。

🎯 このコードでやること：Granger 因果性 — 過去値の予測寄与に基づく因果推論に関連するステップ #10。最終結果のまとめ・保存を行います。

📥 入力の形
# Granger 検定に渡すのは「定常な2系列」を並べた DataFrame。
# 列は説明変数候補 X（原因側）と被説明変数 Y（結果側）の2列で足りる。
# 実データが非定常なら df.diff().dropna() で差分化してから渡す。
#      X        Y
# 0  ...      ...

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12

import pandas as pd
from statsmodels.tsa.api import VAR
import matplotlib.pyplot as plt

df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', header=1)
tokyo = df[df['都道府県'] == '東京都'].sort_values(df.columns[0])
ts = tokyo[['総人口', '出生数', '死亡数']].dropna().diff().dropna()

model = VAR(ts).fit(maxlags=2)
irf = model.irf(10)
irf.plot()
plt.savefig('var_irf.png', dpi=120, bbox_inches='tight')

🐍 Python 実装 — 完全強化版

scipy / pandas / scikit-learn / statsmodels を中心とした標準的な実装例です。 まず CSV を読み込み、 次に グレンジャー因果 の解析を行います。

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20

import pandas as pd
import numpy as np
from scipy import stats

df = pd.read_csv('data/raw/SSDSE-B-2026.csv', skiprows=[1], encoding='cp932')
df = df[df['SSDSE-B-2026'] == df['SSDSE-B-2026'].max()].copy()

x = df['A1101'].astype(float).values
y = df['A4101'].astype(float).values

# 基本統計量
print('n            =', len(x))
print('mean(x)      =', np.mean(x))
print('std(x)       =', np.std(x, ddof=1))

# グレンジャー因果 の代表的計算（用途に応じて scipy/statsmodels を切替える）
r, p = stats.pearsonr(x, y)
print(f'Pearson r = {r:.4f}, p = {p:.4g}')
rs, ps = stats.spearmanr(x, y)
print(f'Spearman rho = {rs:.4f}, p = {ps:.4g}')

用途別の追加実装：

1
2
3
4
5
6
7
8
9

# 標準化と簡易クラスタリングの例
from sklearn.preprocessing import StandardScaler
from sklearn.cluster import KMeans

X = df[['A1101', 'A4101']].astype(float).values
Xs = StandardScaler().fit_transform(X)
km = KMeans(n_clusters=4, n_init=10, random_state=0).fit(Xs)
df['cluster'] = km.labels_
print(df[['Prefecture', 'A1101', 'A4101', 'cluster']].head(10))

1
2
3
4
5
6
7
8

# 時系列（年度別・47都道府県平均の A1101）— 例として ARIMA 系の前処理
import statsmodels.api as sm

df_all = pd.read_csv('data/raw/SSDSE-B-2026.csv', skiprows=[1], encoding='cp932')  # 全年度パネルを再読込
ts = df_all.sort_values('SSDSE-B-2026').groupby('SSDSE-B-2026')['A1101'].mean()
print(ts.tail())
res = sm.tsa.stattools.adfuller(ts)
print('ADF stat:', res[0], 'p:', res[1])

⚠️ Granger 因果性の落とし穴（補強・各 100 文字以上）

⚠️ 落とし穴 — 完全強化版

グレンジャー因果 を実務で扱う際に踏みやすい落とし穴を 5 件挙げます。

• 因果ではなく予測可能性: グレンジャー因果は「先行情報」の意味で、構造的因果ではない。
• ラグ長の選定: 不適切なラグで結果が変わる。AIC/BIC で選定。
• 非定常性: 単位根があると検定統計量の分布が変わる。差分か共和分。
• 共通要因: 第 3 変数が真の駆動なら見かけのグレンジャー因果が出る。
• 非線形関係: 線形 VAR では非線形依存を検出できない。

🌐 関連手法・派生 — 完全強化版

• 時系列解析 — グレンジャー因果が前提とする枠組み
• 定常性 — 非定常系列にそのまま当てると見せかけの因果が出る
• VAR モデル — 複数系列を同時に扱う枠組み。グレンジャー因果はその上の検定
• 自己相関 — ラグ次数を決める手がかり
• 同時性 — 同じ期に互いに決まる関係は、ラグでは切り分けられない

🔗 関連用語ネットワーク（前提・並列・発展）

📚 前提となる用語

• 時系列の基礎
• 定常性 — Granger の前提
• 自己相関
• OLS — VAR の根幹
• F検定 — Granger の核

🔀 並列の概念

• クロス相関
• 共和分
• VECM（誤差修正モデル）
• インパルス応答
• 分散分解

🚀 発展的トピック

• ベイズ VAR
• DCC-GARCH — 時変共分散
• 転送エントロピー — 非線形 Granger
• 因果推論 — 真の因果へ
• 状態空間モデル

🔗 直近の隣人

• ARIMA
• VAR
• 単位根検定
• パネル因果推論

🔗 関連用語（前提・並列・発展） — 完全強化版

本ページ「グレンジャー因果」に関連する用語を 前提（基礎となる概念）／並列（同カテゴリ）／発展（応用先） の 3 区分で整理します。

📥 前提となる用語

• 平均 — 中心傾向の最も基本的な要約値
• 分散 — ばらつきの定量化
• 標準偏差 — 分散の平方根、原データと同じ単位
• 確率 — 不確実性の数学的基礎
• 正規分布 — 連続確率分布の基盤

↔️ 並列・関連の用語

• 相関 — 2 変数の線形関係
• 共分散 — 相関の規格化前の量
• 線形回帰 — 関係性のモデル化
• 仮説検定 — 統計的判断の枠組み
• 信頼区間 — 推定の不確実性
• p 値 — 帰無仮説下での観測の極端さ
• 効果量 — 実質的な大きさ
• 標準誤差 — 推定量のばらつき

🚀 発展・応用の用語

• t 検定 — 平均差の検定
• ANOVA — 多群比較
• カイ二乗検定 — カテゴリ変数の独立性
• ロジスティック回帰 — 二値結果のモデリング
• クラスタリング — 教師なし分類
• 時系列解析 — 時間構造の解析
• パネルデータ — 時間と個体の二次元構造

📚 統計学習の総合ガイド

📚 統計学習の総合ガイド

🎯 学習目標

このページの概念をマスターすることで、 以下のスキルが身につきます：

• 定義と公式を正確に理解
• 適切な使用場面を判断
• Python で実装し、 結果を可視化
• 仮定の確認と診断
• 結果の解釈と報告
• 限界と注意点の理解
• 関連手法との使い分け

📊 SSDSE-B-2026 データの構造

このコンペの主要データセット（SSDSE-B-2026）の構造：

• 47都道府県 × 過去複数年（パネル形式）
• 112列の社会経済指標
• 人口、 出生、 死亡、 婚姻、 経済、 教育、 環境、 家計など多次元
• 政府統計を統合した信頼性の高いデータ

🔍 主要な変数群

💡 ジャストインタイム型学習

このガイドは「必要なときに必要な知識」を提供する設計：

• 論文中の用語をクリック → 該当の用語解説へジャンプ（ポップアップ）
• 概念マップで関連用語を辿る
• 包含マップで体系を把握
• ツリーマップで全体を俯瞰
• Python コードをコピーして実行
• SSDSE データで実際に試す

🛠️ Python データサイエンス環境

🎯 このコードでやること：Granger 因果性 — 過去値の予測寄与に基づく因果推論に関連するステップ #2。基本統計量を計算します。

📥 入力の形
# Granger 検定に渡すのは「定常な2系列」を並べた DataFrame。
# 列は説明変数候補 X（原因側）と被説明変数 Y（結果側）の2列で足りる。
# 実データが非定常なら df.diff().dropna() で差分化してから渡す。
#      X        Y
# 0  ...      ...

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22

# 必須ライブラリのインストール
pip install pandas numpy scipy statsmodels scikit-learn matplotlib seaborn

# 標準的なインポート
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
import seaborn as sns
from scipy import stats
from sklearn.preprocessing import StandardScaler
from sklearn.model_selection import train_test_split
from sklearn.metrics import r2_score, mean_squared_error

# 日本語表示の設定（matplotlib）
plt.rcParams['font.family'] = 'Hiragino Sans'
plt.rcParams['axes.unicode_minus'] = False

# データ読み込み（SSDSE は cp932 エンコーディング）
df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932')
print(df.shape)
print(df.head())
print(df.describe())

🌟 効果的なEDAテンプレート

🎯 このコードでやること：Granger 因果性 — 過去値の予測寄与に基づく因果推論に関連するステップ #3。可視化（散布図／樹形図／時系列プロット）を描きます。

📥 入力の形
# Granger 検定に渡すのは「定常な2系列」を並べた DataFrame。
# 列は説明変数候補 X（原因側）と被説明変数 Y（結果側）の2列で足りる。
# 実データが非定常なら df.diff().dropna() で差分化してから渡す。
#      X        Y
# 0  ...      ...

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26

import seaborn as sns

def quick_eda(df, target=None):
    """探索的データ分析の基本テンプレート"""
    print(f"Shape: {df.shape}")
    print(f"\nColumn types:\n{df.dtypes}")
    print(f"\nMissing values:\n{df.isnull().sum()}")
    print(f"\nBasic stats:\n{df.describe()}")

    # 数値列の可視化
    numeric_cols = df.select_dtypes(include=[np.number]).columns
    df[numeric_cols].hist(bins=20, figsize=(15, 10))
    plt.tight_layout()
    plt.show()

    # 相関ヒートマップ
    if len(numeric_cols) > 1:
        plt.figure(figsize=(12, 10))
        sns.heatmap(df[numeric_cols].corr(), annot=True, fmt='.2f',
                    cmap='RdBu_r', center=0)
        plt.show()

    # ターゲットがあれば散布図行列
    if target and target in df.columns:
        sns.pairplot(df[numeric_cols[:5]], hue=target if df[target].dtype == 'O' else None)
        plt.show()

📈 報告書テンプレート

分析結果を報告する際の標準的な構成：

1. 背景・目的：なぜこの分析が必要か
2. データ：出所、 サンプルサイズ、 期間
3. 方法：使用した統計手法、 仮定
4. 結果：図表、 統計量、 検定結果
5. 解釈：結果が何を意味するか
6. 限界：分析の制約
7. 結論：要点まとめ、 今後の課題

• クラスタリング教材
• 時系列解析教材
• 処置群・対照群教材
• 対照群教材
• 自然言語処理教材

🎯 まとめ — 完全強化版

本ページでは「グレンジャー因果」を 12 セクション（🔖 キーワード索引／💡 30 秒結論／📍 文脈／🎨 直感／📐 数式／🔬 数式を言葉で読み解く／🧮 実値計算／🐍 Python 実装／⚠️ 落とし穴／🌐 関連手法／🔗 関連用語／📚 グループ教材）で完結に整理しました。 SSDSE-B-2026 を素材に、 概念の輪郭・式の意味・実装手順・典型的な失敗パターンの 4 点を最低限押さえれば、 統計データ分析コンペの現場で迷わず使えるはずです。

🖼 図で振り返るGranger因果

Granger 因果は「過去の X が Y の将来を予測する追加情報を持つか」を時系列の文脈で判定する。 ここでは時系列・自己相関・パネル構造の 3 図で直感を掴む。

読み方: 時系列 (図 1) で先行・後続の関係を視察 → 自己相関 (図 2) でラグ次数を決定 → パネル (図 3) で個体ごとの異質性を考慮、 という 3 ステップが Granger 検定の実務フローと一致する。

🗺️ 統計手法選択フローチャート

Q1: 何を知りたい？

• 記述したい → 平均、 分散、 ヒストグラム
• 比較したい → t検定、 ANOVA、 χ²検定
• 関係を見たい → 相関、 回帰
• 予測したい → 回帰、 機械学習
• 分類したい → ロジスティック回帰、 SVM、 RF
• グループ分けしたい → クラスタリング
• 次元を減らしたい → PCA、 因子分析
• 因果関係を知りたい → RCT、 IV、 DiD、 PSM

Q2: データの種類は？

• 連続値 → t検定、 ANOVA、 線形回帰
• カテゴリ → χ²検定、 ロジスティック回帰
• 順序 → ノンパラ検定、 順位回帰
• カウント → ポアソン回帰、 負の二項回帰
• 時系列 → ARIMA、 VAR、 状態空間
• パネル → 固定効果、 ランダム効果

Q3: サンプルサイズは？

• n < 30：ノンパラ、 ベイズ、 ブートストラップ
• 30 ≤ n < 200：古典的検定、 単純な回帰
• n ≥ 200：複雑なモデル、 機械学習
• n ≥ 10000：深層学習も可能

Q4: 仮定は？

• 正規性：満たす → パラメトリック / 満たさない → ノンパラ
• 独立性：必須 / 違反 → クラスター調整、 時系列モデル
• 等分散性：満たす → OLS / 違反 → WLS、 ロバスト

📏 効果量の参照表

p値だけでなく効果量も併記するのが現代統計の標準。 主要な指標と Cohen の解釈基準：

🚀 実務応用の深掘り

典型的なプロジェクトの流れ

1. 問題理解：ステークホルダーとの対話、 KGI/KPI 設定
2. データ収集：内部DB、 公的データ（SSDSE等）、 API
3. EDA：データの全体像把握、 異常検出
4. 仮説立案：ドメイン知識からの仮説
5. モデリング：シンプルから複雑へ段階的に
6. 検証：CV、 ホールドアウト、 A/Bテスト
7. 解釈：可視化、 SHAP、 部分依存プロット
8. 展開：本番デプロイ、 監視

ベストプラクティス

• シンプルなモデルから始める（線形回帰、 単純ルール）
• 必ずベースラインと比較
• 過学習を防ぐ（CV、 正則化、 早期停止）
• 解釈可能性を重視
• 再現可能なコード・ノートブック
• バージョン管理（Git）と環境管理（venv, conda）
• ドキュメント化を怠らない

論文・コンペでよく使う言い回し

統計データ活用コンペでのコツ

• SSDSE データの構造を理解し、 適切なテーブルを選ぶ
• 地域別・年度別の比較で時空間的視点を入れる
• 1つの分析で多角的に切り口を変える
• 仮説と発見の両方を持つ
• ストーリーラインを明確に
• 図表を1枚1枚作り込む
• 政策提言や実務的意義に繋げる

🗺️ 概念マップ — 3つの視点で体系を理解する

Granger因果検定 がデータサイエンスの体系の中でどこに位置するかを、 3つの異なる視点で可視化します。 同じ情報でも見方を変えると気付きが変わります。

📍 体系階層のパス

🌐 統計・データサイエンス › 因果推論 › 因果手法 › Granger因果検定

① 🔗 関係マップ — 「他の手法とどう繋がっているか」

中心に Granger因果検定 を置き、 そこから 時系列分析・ARIMA・VAR 計 3 個の用語へ、 前提・兄弟・発展形といった関係を矢印で結んでいます。 線がどちら向きかを見れば、 先に読むべき用語とあとから読む用語が分かります。 ノードはドラッグで動かせ、 ホイールでズーム、 クリックでその用語のページへ遷移します。

② ⭕ 包含マップ — 「どのカテゴリに含まれているか」

大きな円が小さな円を包含する Circle Packing 図。 「Granger因果検定」は緑色でハイライト。

• カテゴリ円をクリック：その内部にズームイン
• 白背景クリック：1階層戻る
• 用語円をクリック：詳細ページへ遷移
• マウスホバー：階層パス表示

③ 🌳 ツリーマップ — 「面積で見るボリューム比較」

長方形を入れ子に分割した Treemap 図。 各分野の規模感を面積で比較。 「Granger因果検定」は緑色でハイライト。

• カテゴリ矩形をクリック：その内部にドリルダウン
• パンくず（上のリンク）クリック：その階層に戻る
• 用語矩形をクリック：詳細ページへ遷移
• マウスホバー：階層パスと値を表示

🎯 3つのマップの使い分け

💡 3 つの図は同じ位置関係を別の角度から描いています。 関係マップは Granger因果検定 の隣に何が並ぶかを、 包含マップとツリーマップは 統計・データサイエンス → 因果手法 → Granger因果検定 という入れ子の位置を示します。 「因果手法には他に何があったか」を思い出したいときは包含マップを、 「次に何を読むか」を決めたいときは関係マップを開いてください。

🗺 概念マップ — 完全強化版

• 基礎：平均／分散／標準偏差／確率
• 分布：正規分布／標本／中心極限定理
• 関係性：相関／共分散／線形回帰
• 推測：仮説検定／信頼区間／p 値／効果量
• 因果：内生性／IV／DID／平行トレンド／RDD
• 分類：クラスタリング／階層クラスタリング／Ward 法
• 時系列：ARIMA／VAR／グレンジャー因果／自己相関

🧭 解説深化 — 直感・落とし穴・発展（追記）

ここまでの本文・ウィジェットに加えて、グレンジャー因果を「予測因果」として腹落ちさせるための追補です。既存の説明と重複する部分もありますが、初学者が最短で誤解を避けられるよう、直感 → 落とし穴 → 発展の順に整理し直しています。

🎨 直感 — 「過去の情報が将来予測を助けるか」

グレンジャー因果の本質は「時間的先行」と「予測寄与」の2点に尽きます。
今の Y を当てるのに Y 自身の過去 だけで十分か、それとも X の過去 を足すと予測誤差が有意に減るか——これだけを見ています。X を足して初めて予測が改善するなら「X の過去には Y の未来を先取りする情報が含まれる」＝「X は Y をグレンジャー予測する」と読みます。

• キーワードは「情報」：因果の矢印そのものではなく、「X を知っていると Y の予測が得か？」という情報量の話。
• 非対称性が要点：X→Y が有意でも Y→X は非有意、ということが起こり得る。向きで結果が変わるのが相関との決定的な違いです（相関は対称）。
• 「予測因果」であって真の因果ではない：予測が当たることと、介入して Y を動かせることは別問題です（因果関係／逆の因果）。

⚠️ 落とし穴（重要）— グレンジャー因果 ≠ 真の因果

グレンジャー因果は誤用が非常に多い概念です。以下は「有意に出たのに結論を誤る」典型パターンです。

⚠️ 合言葉：「グレンジャー因果あり」は「時間的先行・予測寄与あり」と読み替える。介入で結果が変わる真の因果（因果関係）とは常に区別して報告します。

🚀 発展 — VAR・条件付き・非線形・真の因果推論

• VAR での検定：グレンジャー因果検定は VAR の係数に対する「ゼロ制約の F/Wald 検定」として実装される。多変量で複数方向を一度に検定できる。
• 定常化とラグ選択：ADF/KPSS →（非定常なら）差分化、共和分があれば VECM。ラグ次数は AIC/BIC/HQIC で選ぶ（定常性・自己相関）。
• 条件付きグレンジャー因果：第3変数 Z を回帰に含めた上で X→Y を検定し、交絡を統制する。多変量 VAR での部分因果に相当。
• 伝達エントロピー（Transfer Entropy）：情報理論に基づく非線形版グレンジャー因果。分布全体の依存を測るため線形 VAR が見落とす関係も検出できる（情報エントロピー）。ガウス変数では線形グレンジャーと等価になることが知られる。
• インパルス応答関数（IRF）・分散分解：「向き」だけでなく「ショックの大きさと持続」を可視化。有無（検定）と効果量（IRF）を両輪で読む。
• 真の因果推論との違い：介入・反実仮想に基づく枠組み（パネル因果推論・因果関係）は「もし X を動かしたら Y はどうなるか」を問う。グレンジャー因果は観測データの予測改善にとどまり、これに置き換わるものではない。

🧮 実データで確認 — 東京都の総人口 A1101 × 出生数 A4101（SSDSE-B-2026 実測）

SSDSE-B-2026 の東京都を 2012–2023 年の12点の時系列として取り出し、実際にグレンジャー因果検定を回した結果です（encoding='cp932', skiprows=[1] で読み込み）。原系列は明確に非定常なので1階差で定常化してから検定しました。

💬 読み方（重要な注意つき）：どの向き・どのラグでも p>0.05 で、統計的にはグレンジャー因果は検出されませんでした。ただしこれは「因果が無い証明」ではありません。わずか12年（差分後11点）という極端な短系列では検定の検出力が非常に低く、真に効果があっても見逃しやすい（第2種の過誤）。都道府県×年度パネル（47県×12年）をパネル手法でプールする、より長期の系列を使う、といった対処が必要です。短系列での「有意でない」も「有意」も、そのまま因果の結論にしないことが最大の教訓です。

🔗 関連ページ（このサイト内・実在ページのみ）

• VAR（ベクトル自己回帰） — グレンジャー因果検定の土台となる多変量枠組み
• 時系列 — 定常化・ラグ・差分などの前提知識
• 定常性 — 検定の前提。ADF/KPSS と差分化
• 自己相関 — ラグ構造の把握
• ラグ変数 — ラグ次数 p の選択（AIC/BIC）
• 因果関係 — 予測因果と真の因果の違い
• パネル因果推論 — 介入・反実仮想に基づく真の因果へ
• 交絡 — 共通原因による偽陽性
• 見せかけの相関 — 共通トレンド・見せかけの回帰
• 逆の因果 — 向きの取り違え
• 情報エントロピー — 非線形版（伝達エントロピー）の基礎

※ 「伝達エントロピー」「VECM」「Toda-Yamamoto 検定」「インパルス応答」等は本サイトに独立ページが無いため、本文中では用語のみの言及にとどめています。

🔗 隣接手法への橋渡し

「グレンジャー因果」は単独で完結する手法ではなく、 隣接領域と連携することで真価を発揮する。

• 上流: 入力データの前処理・整備 — 時系列分析 ・ VAR モデル (Vector AutoRegression)
• 並列: 類似目的を持つ代替手法 — 転移エントロピー (Schreiber 2000) ・ CCM (Sugihara 2012)
• 下流: 出力結果の評価・解釈 — 因果推論 (Pearl 2009) ・ 構造方程式モデル

上流の VAR モデル・定常性検定 (ADF) で前提を確認し、 並列の Transfer Entropy・Convergent Cross Mapping と非線形版を比較し、 下流の構造 VAR・インパルス応答関数で因果効果を解釈する。 グレンジャー因果は「予測改善」という限定的因果概念で、 必ず真の介入因果 (Pearl 流) と区別して報告する。

🌳 手法選択フロー

「granger causality」を含む手法選択は、 データの性質・分析目的・運用制約の 3 軸で決まる。 以下に典型シナリオと対応する手法の組合せを示す。

典型シナリオ別の手法選択

選んだ後の検証ステップ

1. 前処理確認: 標準化 / 欠損処理 / カテゴリ変数のエンコード が適切か
2. ハイパラ調整: 交差検証で安定する値を選ぶ
3. 性能評価: タスクに応じた指標 (RMSE / Accuracy / F1 / AUC / シルエット等) を複数併用
4. 頑健性チェック: 別手法 / 別シードで結果が大きく違わないか確認
5. 解釈: 結果を分野知識と照らし合わせ、 意味のある発見になっているか検討