論文一覧に戻る 🕸 論文ネットワーク 🗺 用語ネットワーク 統計データ分析コンペ 教育用再現集
2024年 統計データ分析コンペティション | 統計活用奨励賞 [大学生・一般の部]

金融資産購入経験の要因分析
―金融教育・損失回避傾向・Digital Capability Indexに注目して―

⏱️ 推定読了時間: 約27分
NGUYEN THI NGOC ANH・NGUYEN THI MINH QUY 青森中央学院大学
🔬 VIF/多重共線性🔬 パネルデータ分析🔬 ロジスティック回帰🔬 重回帰🏷 金融🏷 IT・デジタル
🔬 手法タグ / 🏷 応用タグ(クリックで同タグの論文をネットワーク検索)
🔬 再現カルテ — 原論文と同じ種類のデータで再現(細部は簡略化)

この教材は原論文と同じ種類のデータで計算し直しています。ただし工程の一部は、学習しやすさのために簡略化しています(下の「できないこと」を参照)。

原論文が使ったデータ金融リテラシー調査・デジタル・ケイパビリティ・インデックス
分析単位:都道府県
中核手法:重回帰分析
この教材が使うデータ
原論文(PDF)金融資産購入経験の要因分析 ―金融教育、損失回避傾向、Digital Capability Index に注目して―
統計活用奨励賞/NGUYEN THI NGOC ANH、NGUYEN THI MINH QUY (青森中央学院大学経営法学部経営法学科)
✅ この教材でできること
  • 原論文の中核手法(重回帰分析)を実データで実行できる
  • 分析に使うデータは同梱済みで、ブラウザ上でそのまま読み込める
  • 数値・図を最後まで再計算できる(重回帰分析(OLS)・VIF(多重共線性の確認)・相関分析)
⚠️ この教材ではできないこと(原論文との違い)
  • 手順の一部が原論文と違う:論文=重回帰(OLS、2022年クロスセクション)。教材=ロジスティック回帰+オッズ比(別分析)
📄 原論文との違いを、もっと詳しく
📄 原論文:金融リテラシー調査 2022(都道府県別)+野村総合研究所 DCI で、株式・投信・外貨預金の購入経験者割合を OLS 重回帰(8 説明変数・頑健標準誤差)により分析。
📘 本教材:原論文と同じ金融リテラシー調査 2022 の公開統計表から 8 指標を抽出し、同じ OLS 重回帰 3 モデルを再現。DCI(非公開)は除外し、年収は SSDSE-E の 1 人当たり県民所得で代理。
⚠️ 注意:DCI と年収の扱いが違うため係数の値は原論文と完全には一致しないが、中心的発見(金融教育の経験は非有意・正答率が 1%有意で正)は同じ結果が得られている。

▶ ここに書いてある分析は、ページ下部の「🐍 ブラウザで動かす」でそのまま実行できます(インストール不要)。動かしているのは code/2024_U4_katsuyo.py(124 行)そのものです。

🏅 論文審査会コメント(審査員はここを評価した)
「分析テーマや目的変数、説明変数の設定や分析手順が的確であり、重回帰分析の結果の解釈も手堅くよくまとまっている。政策的知見もあり好ましい。変数の本当の意味や実態を理解し分析するとなおよい。」
出典:統計センター「統計データ分析コンペティション」受賞論文PDF掲載の審査講評。プロの審査員が何を評価し、何を注文したかは論文の読み方の最良の手本になる。
📝 3行で分かる要約

目次

  1. 研究概要と背景
  2. 原論文の分析枠組み(OLS 重回帰)
  3. データの入手と結合
  4. VIF(多重共線性チェック)
  5. 重回帰 3 モデルの推定(表3の再現)
  6. 結果の解釈
  7. まとめ
  8. 📥 データの準備
  9. 💼 実社会での応用
  10. ⚠️ よくある誤解
  11. 📖 用語集
  12. 📐 手法ガイド
  13. 🚀 発展の可能性
  14. 🎯 自分でやってみよう
  15. 🤔 Q&A
  16. 🐍 ブラウザで動かす

🎯 この記事を読むと何ができるようになるか

📥 データの準備(再現コードを動かす前に)

このページの分析を自分で再現するには、以下の手順でデータを準備してください。コードの編集は不要です。

1
データをダウンロードする 統計センターの SSDSE 配布ページから、以下のファイルをダウンロードします。
kinyu_literacy_2022_pref.csv ← 金融リテラシー調査 2022(都道府県別・data/2024_U4/ に配置)📥 直接DL
SSDSE-B-2026.csv ← SSDSE-B(都道府県データ)📥 直接DL
SSDSE-E-2026.csv ← SSDSE-E(都道府県の指標2)📥 直接DL
⬇ SSDSEダウンロードページを開く
2
ファイルを所定のフォルダに配置する ダウンロードしたCSVを、プロジェクトの data/raw/ フォルダに入れます。
2026 統計・データ解析コンペ/ ├── code/ │ └── 2024_U4_katsuyo.py ← 実行するスクリプト └── data/ └── raw/ SSDSE-B-2026.csv ← ここに置く SSDSE-E-2026.csv ← ここに置く
3
スクリプトをそのまま実行する ターミナルでプロジェクトルートに移動し、以下を実行します。
python3 code/2024_U4_katsuyo.py
図は html/figures/ に自動保存されます。
研究概要と背景

「貯蓄から投資へ」が政策目標とされる中、金融資産への投資行動(株式・投資信託等の購入)を阻む・促す要因の解明が重要である。本研究は、金融リテラシー調査 2022 の都道府県別集計データ(N=47)を用い、金融教育・金融リテラシー(正答率)・損失回避傾向・デジタル能力指数(DCI)などが株式・投資信託・外貨預金の購入経験者割合に与える効果を OLS 重回帰分析(頑健標準誤差)で分析した。

原論文の概要(統計センター公式):「都道府県別に人々の金融に関連する行動の違いを、金融に関する教育や理解度により説明するモデルにより推定し、金融教育の経験の有無よりも、実際に金融に関する知識の深さや判断力の高さが、金融資産の購入経験者割合を高めることを示すとともに、デジタル機器の利用可能性の高さが、株式購入等の経験者割合を高める可能性があることを示した。」
本ページではこの分析の流れを実データでたどりながら、使われた統計手法を一つずつ学んでいく。

3つの主要変数
  • 金融教育受講経験:学校・職場・セミナーなどで金融の知識・スキルを学んだことがある(0/1)
  • 損失回避傾向:「損失は利益より2倍以上辛い」という行動経済学の概念。スコアが高いほど損失をより強く嫌う
  • Digital Capability Index(DCI):オンラインサービスの利用・情報収集能力を数値化した指標
分析フロー
金融リテラシー
調査 2022年
N=25,000
VIF確認
(多重共線
性排除)
OLS
重回帰×3
(頑健SE)
係数の符号
・有意性の解釈
Forest Plot

重回帰(OLS) 頑健標準誤差 VIF DCI・損失回避

1
原論文の分析枠組み — OLS 重回帰(47都道府県クロスセクション)

原論文は、都道府県別の集計データ(N=47)を使った最小二乗法(OLS)による重回帰分析です。金融資産の購入経験の程度を 3 つの被説明変数で捉え、それぞれについて同じ 8 つの説明変数で回帰します。

Yij = α0 + α1×金融教育i + α2×正答率i + α3×損失回避i + α4×近視眼i + α5×横並びi + α6×DCIi + α7×年収i + α8×高齢化率i + μi (i = 1, …, 47) # j=1: 株式購入経験者割合, j=2: 投資信託購入経験者割合, j=3: 外貨預金経験者割合

変数の定義と出典(原論文 表1)

変数定義出典符号予想
❶〜❸ 被説明変数株式/投資信託/外貨預金を購入した経験がある人の割合(%)金融リテラシー調査 2022
1 金融教育を受けた人の割合学校等で金融教育を受けた人の割合(%)金融リテラシー調査 2022
2 正誤問題25問の正答率金融リテラシーの客観的指標(%)金融リテラシー調査 2022
3 損失回避傾向が強い人の割合行動経済学的バイアス(%)金融リテラシー調査 2022
4 近視眼的行動バイアスが強い人の割合行動経済学的バイアス(%)金融リテラシー調査 2022
5 横並び行動バイアスが強い人の割合行動経済学的バイアス(%)金融リテラシー調査 2022±
6 DCIデジタル・ケイパビリティ・インデックス野村総合研究所(2022)
7 平均年収都道府県別の平均年収(千円)(本教材では県民所得で代理)
8 高齢化率65歳以上人口割合(%)SSDSE 等±
本教材の再現範囲 説明変数 1〜5 と被説明変数 ❶〜❸ は、原論文と同じ出典(金融広報中央委員会「金融リテラシー調査 2022」の都道府県別統計表)から取得した実データを使う。6 DCI は野村総合研究所の非公開データのため除外し、7 年収は SSDSE-E の「1人当たり県民所得」で代理する。この 2 点を除き、原論文と同じモデルを推定する。

DS LEARNING POINT 1

頑健標準誤差(Robust SE)

原論文の表3の括弧内は「頑健な標準誤差」。誤差項の分散が観測ごとに異なる(不均一分散)場合でも標準誤差を正しく推定するための方法で、statsmodels では fit(cov_type='HC1') と書くだけで使える。

2
データの入手と結合(原論文と同じ出典の実データ)

金融リテラシー調査 2022 の統計表(22lite_toukeir.xlsx)は知るぽると(金融広報中央委員会)で公開されています。本教材ではそのうち「86 都道府県比較表 (1)(3)(4)」から 8 指標を抽出して CSV 化したもの(data/2024_U4/kinyu_literacy_2022_pref.csv)を使います。

やってみようデータ読み込み — 金融リテラシー調査 + SSDSE の結合
📝 コード
 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
import matplotlib
matplotlib.use('Agg')
import matplotlib.pyplot as plt
import numpy as np
import pandas as pd
import statsmodels.api as sm

plt.rcParams['font.family'] = 'Hiragino Sans'
plt.rcParams['axes.unicode_minus'] = False

# 1) 金融リテラシー調査 2022(都道府県別、原論文と同じ出典)
lit = pd.read_csv('data/2024_U4/kinyu_literacy_2022_pref.csv', index_col='都道府県')

# 2) SSDSE から高齢化率(B、2022年度)と1人当たり県民所得(E)
b_raw = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', header=1)
b = b_raw[(b_raw['年度'] == 2022)
          & b_raw['地域コード'].str.match(r'^R\d{5}$', na=False)].set_index('都道府県')
lit['aging'] = b['65歳以上人口'] / b['総人口'] * 100

e_raw = pd.read_csv('data/raw/SSDSE-E-2026.csv', encoding='cp932', header=2)
e = e_raw[e_raw['地域コード'].str.match(r'^R\d{5}$', na=False)].set_index('都道府県')
lit['income'] = pd.to_numeric(e['1人当たり県民所得(平成27年基準)'], errors='coerce') / 1000

print('=== データ確認(先頭5県) ===')
print(lit.head().round(2))
▼ 実行結果
=== データ確認(先頭5県) === seitou stock toshin gaika ... kinshigan yokonarabi aging income 都道府県 ... 島根県 58.83 29.22 33.12 19.48 ... 47.40 18.18 34.80 2.91 奈良県 58.80 34.58 33.96 15.26 ... 43.30 14.33 32.39 2.55 香川県 58.55 34.53 35.43 18.83 ... 43.95 18.83 32.33 2.85 千葉県 58.28 36.50 32.74 20.55 ... 45.45 16.93 27.98 3.06 愛媛県 57.35 33.33 33.97 16.03 ... 46.47 17.31 33.92 2.67 [5 rows x 10 columns]
💡 解説
  • lit が金融リテラシー調査の 47 都道府県 × 8 指標。stock(株式)・toshin(投資信託)・gaika(外貨預金)が被説明変数です。
  • SSDSE-B から高齢化率(2022年度)、SSDSE-E から 1 人当たり県民所得を結合。原論文の説明変数 7・8 に対応します。
💡 Python TIPS set_index("都道府県") どうしで割り算・代入すると、pandas が県名で自動整列(アラインメント)してくれます。
3
VIF による多重共線性チェック(原論文 表2)

原論文は重回帰の前に説明変数間の VIF(分散拡大要因)を確認し、最大値 約3.6 < 5 であることから「多重共線性の発生可能性は低い」と判断しています。本教材のデータでも同じチェックをします。

やってみようVIF の計算(補助回帰の R² から)
📝 コード
1
2
3
4
5
6
7
8
9
X_cols = ['kyoiku', 'seitou', 'sonshitsu', 'kinshigan', 'yokonarabi', 'income', 'aging']

# VIF(多重共線性チェック、原論文の表2に対応)
print('=== VIF(分散拡大要因) ===')
X = sm.add_constant(lit[X_cols])
for col in X_cols:
    others = [c for c in X_cols if c != col]
    r2 = sm.OLS(lit[col], sm.add_constant(lit[others])).fit().rsquared
    print(f'  {col:12s} VIF = {1 / (1 - r2):.2f}')
▼ 実行結果
=== VIF(分散拡大要因) === kyoiku VIF = 1.29 seitou VIF = 1.55 sonshitsu VIF = 1.90 kinshigan VIF = 1.36 yokonarabi VIF = 1.02 income VIF = 1.61 aging VIF = 1.53
💡 解説
  • VIFj = 1/(1−R²j)。R²j は「説明変数 j を他の説明変数で回帰したときの決定係数」です。
  • 本教材の最大 VIF は 1.90(損失回避)。原論文(最大約3.6)と同様、目安の 5 を大きく下回り、多重共線性の心配はありません。
💡 Python TIPS 専用関数(variance_inflation_factor)を使わなくても、VIF の定義どおり OLS を 1 回まわせば計算できます。
4
OLS 重回帰 3 モデルの推定(原論文 表3の再現)
やってみよう3つの被説明変数それぞれに OLS を当てる(頑健標準誤差)
📝 コード
 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
# 重回帰(OLS)3モデル、原論文の表3に対応
# cov_type='HC1' は頑健標準誤差(原論文の表3注1と同じ)
for y_col, y_name in [('stock', '❶株式購入経験者割合'),
                      ('toshin', '❷投資信託購入経験者割合'),
                      ('gaika', '❸外貨預金経験者割合')]:
    model = sm.OLS(lit[y_col], X).fit(cov_type='HC1')
    print(f'=== 被説明変数: {y_name} ===')
    print(f'  N=47  自由度調整済みR2={model.rsquared_adj:.3f}')
    for name in ['const'] + X_cols:
        c, se, p = model.params[name], model.bse[name], model.pvalues[name]
        sig = '***' if p < 0.01 else ('**' if p < 0.05 else ('*' if p < 0.10 else ''))
        print(f'  {name:12s} {c:+8.3f} ({se:.3f}) {sig}')
▼ 実行結果
=== 被説明変数: ❶株式購入経験者割合 === N=47 自由度調整済みR2=0.512 const +22.209 (27.609) kyoiku -0.141 (0.332) seitou +0.869 (0.263) *** sonshitsu -0.554 (0.221) ** kinshigan +0.043 (0.186) yokonarabi +0.410 (0.243) * income +1.550 (1.247) aging -0.302 (0.181) * === 被説明変数: ❷投資信託購入経験者割合 === N=47 自由度調整済みR2=0.709 const +23.591 (18.528) kyoiku -0.298 (0.210) seitou +1.025 (0.189) *** sonshitsu -0.599 (0.136) *** kinshigan -0.153 (0.173) yokonarabi +0.283 (0.148) * income +0.250 (0.438) aging -0.050 (0.113) === 被説明変数: ❸外貨預金経験者割合 === N=47 自由度調整済みR2=0.667 const -12.595 (24.473) kyoiku -0.185 (0.387) seitou +0.769 (0.197) *** sonshitsu -0.164 (0.185) kinshigan +0.087 (0.118) yokonarabi +0.070 (0.126) income +2.092 (0.669) *** aging -0.372 (0.126) ***
💡 解説
  • 3 モデルとも 2 正答率が 1%水準で有意に正1 金融教育を受けた人の割合はすべて非有意 — 原論文の中心的発見がそのまま再現できました。
  • 損失回避傾向は投資信託に対して負で 1%有意(原論文と同じ)。株式にも 5%有意の負(原論文では非有意)で、DCI・年収の扱いの違いによる差と考えられます。
💡 Python TIPS fit(cov_type="HC1") — OLS の係数はそのままで、標準誤差だけ頑健版に置き換わります。
重回帰係数のForest Plot(3モデル比較)
図1:3 モデルの回帰係数と95%信頼区間(赤 = 5%水準で有意)。正答率だけが 3 モデル共通で右側(正)に有意。

原論文(表3)と本教材の再現結果の対比

説明変数❶株式❷投資信託❸外貨預金
原論文本教材原論文本教材原論文本教材
1 金融教育を受けた人の割合−0.112−0.141−0.308−0.298−0.069−0.185
2 正誤問題25問の正答率+0.658***+0.869***+0.992***+1.025***+0.725***+0.769***
3 損失回避傾向が強い人の割合−0.298−0.554**−0.607***−0.599***−0.060−0.164
4 近視眼的行動バイアス+0.258+0.043−0.147−0.153+0.146+0.087
5 横並び行動バイアス+0.269+0.410*+0.289*+0.283*+0.070
6 DCI+0.180*(除外)+0.048(除外)+0.119*(除外)
7 年収/県民所得+0.005+1.550+0.001+0.250+0.002**+2.092***
8 高齢化率−0.016−0.302*+0.080−0.050−0.239−0.372***
自由度調整済み R²0.6340.5120.7020.7090.6740.667

注:***, **, * は順に 1%, 5%, 10% 水準で有意(両者とも頑健標準誤差)。原論文の年収は千円単位、本教材の県民所得は百万円単位なので係数の大きさは直接比較できない。

5
結果の解釈 — 「教育を受けたか」より「リテラシーの高さ」
原論文の中心的発見(本教材でも再現) 「学校等で金融教育を受けた人の割合」はどのモデルでも有意でないのに対し、「正誤問題25問の正答率」は3 モデルすべてで 1%水準の有意な正。つまり金融教育を受けた経験そのものではなく、実際に身についた金融リテラシーの高さが金融資産の購入経験者割合を高めている。
正答率と株式購入経験者割合の散布図
図2:正誤問題25問の正答率(金融リテラシー)と株式購入経験者割合。r = +0.57 の明確な正の関係。

その他の変数の解釈(原論文の議論)

集計データの限界(原論文 5節) 47 都道府県の集計データによるクロスセクション分析なので、個人レベルの因果関係は特定できない(生態学的誤謬に注意)。また DCI が 2020 年からしか公開されていないため、固定効果を使ったパネル分析はできず「固定効果を有する説明変数の欠落から生じるバイアスを回避できていない可能性は残されている」と原論文も明記している。

まとめ

主要な発見(原論文・本教材の再現の両方で確認)

  1. 金融教育を受けた経験は非有意:「学校等で金融教育を受けた人の割合」は 3 モデルすべてで有意でない。
  2. 金融リテラシー(正答率)が最重要:正誤問題25問の正答率は 3 モデルすべてで 1%水準の有意な正。「受けたかどうか」より「身についているか」。
  3. 損失回避傾向は投資信託購入を抑制:係数 −0.6 前後で 1%有意(原論文・再現とも)。
  4. 政策示唆:金融教育は「実施したか」ではなく「リテラシーが身についたか」で評価すべき。教育内容の質と定着の測定が重要。
OLS 重回帰+頑健標準誤差という選択 被説明変数が「割合(連続値)」なので、原論文は 2 値データ用のロジスティック回帰ではなくOLS 重回帰を使い、不均一分散に備えて頑健標準誤差で推論している。N=47 と小さいので、説明変数の数(8個)とのバランスにも注意が必要。
教育的価値(この分析から学べること)
  • 符号条件の事前予想:推定前に各説明変数の符号を理論・先行研究から予想し、結果と突き合わせる。
  • VIF による多重共線性診断:重回帰の前の定番チェック(目安 VIF < 5)。
  • 複数の被説明変数で頑健性確認:株式・投信・外貨預金の 3 つで同じ結論(正答率が有意)が出ることが説得力を生む。
  • 集計データの限界:都道府県集計の相関は個人レベルの因果と別物(生態学的誤謬)。

データ・コードのダウンロード

分析スクリプト(2024_U4_katsuyo.py)
データ出典
金融リテラシー調査 2022(都道府県別 8 指標)金融広報中央委員会「金融リテラシー調査 2022」統計表(22lite_toukeir.xlsx)より抽出
高齢化率(2022年度)SSDSE-B-2026(統計センター)
1人当たり県民所得SSDSE-E-2026(統計センター)

本教育用コードは実データのみを使用(合成データ・np.random 不使用)。DCI(野村総合研究所)は非公開のため本再現には含まれない。

教育用再現コード | 2024年 統計データ分析コンペティション 統計活用奨励賞 [大学生・一般の部] | 青森中央学院大学

⚠️ よくある誤解と注意点

統計分析の解釈で初心者がやりがちな勘違いをまとめます。特に「相関因果の混同」「p値の過信」は研究現場でもよく起きる落とし穴です。本文を読む前にも、読んだ後にも、目を通してみてください。

❌ 「相関がある=因果関係がある」ではない
疑似相関spurious correlationとは、見かけ上は関係があるように見えるが、実際は無関係、または第三の変数(交絡変数)が両方に影響しているだけの現象です。

古典例: アイスクリームの売上 と 水難事故件数 は強く相関するが、片方が他方を引き起こしているわけではない。両者とも「夏の暑さ」という第三の変数に引きずられているだけ。

論文を読むときの心構え: 「○○と△△に強い相関が見られた」だけで終わっている主張は、本当に因果関係があるのか、それとも第三の変数(人口・所得・地理など)が共通要因として効いているだけではないかを必ず疑ってください。
❌ 「p値が小さい=重要な発見」ではない
p値が小さい(例えば p < 0.001)ことは「統計的に偶然とは考えにくい」という意味であって、「実用的に大きな効果がある」という意味ではありません。

例: 巨大なサンプルサイズ(n=100,000)では、相関係数 r=0.02 でも p < 0.001 になります。しかし r=0.02 は実用上ほぼ無視できる関係です。

正しい読み方: p値効果量係数の大きさ、相関係数の値)の両方をセットで判断してください。p値だけで「重要な発見」と結論づけるのは誤りです。
❌ 「回帰係数が大きい=重要な変数」ではない
回帰係数の絶対値は、説明変数単位に強く依存します。「年収(万円)」と「失業率(%)」の係数を直接比較しても意味がありません。

正しい比較方法: (1) 標準化係数(各変数を平均0・分散1に変換した上での係数)を使う、(2) 限界効果(変数を1標準偏差動かしたときのyの変化)で比較する。

また、係数の大きさが「因果関係の強さ」を意味するわけでもありません。あくまで「相関的な関連の強さ」です。
❌ 「外れ値を除外すれば正しい結果」ではない
外れ値(極端な値)を「目障りだから」「結果が綺麗にならないから」という理由で除外するのは分析の改ざんに近い行為です。

外れ値が示すもの: 本当に重要な情報(東京の超高密度、北海道の超低密度など)であることが多い。外れ値を取り除くと「日本全体の傾向」を見誤る原因になります。

正しい対処: (1) 外れ値の出現要因を調査する(なぜ東京だけ突出するのか)、(2) ノンパラメトリック手法(Spearman相関Kruskal-Wallis)を使う、(3) 外れ値を含む結果と除外した結果の両方を提示し、解釈を読者に委ねる。
❌ 「サンプルサイズが大きい=信頼できる」ではない
サンプルサイズ(n)が大きいと統計的検定の検出力は上がりますが、それは「偶然による誤差を減らす効果」にすぎません。

nが大きくても解消されない問題:
選択バイアス標本が偏っている)
測定誤差(変数の定義が曖昧)
欠損値のパターン(欠損がランダムでない)
交絡変数の見落とし

例: 1万人にWeb調査して「ネット利用と幸福度は強く相関」と言っても、そもそも回答者がネットユーザー寄りに偏っているため、母集団全体の結論にはなりません。
❌ 「複雑なモデル=より良い分析」ではない
ランダムフォレストニューラルネット・複雑な階層モデルなど、高度な手法を使えば「良い分析」と感じがちですが、必ずしもそうではありません。

過学習(overfitting)の罠: モデルが複雑すぎると、訓練データ偶然のパターンまで学習してしまい、新しいデータでは予測精度が落ちます。

シンプルさの価値: 重回帰分析相関分析は「結果が解釈しやすい」「再現性が高い」という大きな利点があります。複雑な手法はシンプルな手法で答えが出ない時の最後の手段です。
❌ 「多重共線性は気にしなくていい」ではない
多重共線性とは、説明変数同士の相関が極めて強い状態のこと。これを放置すると、回帰係数符号や大きさが入れ替わる異常事態が起こります。

典型例: 「総人口」と「労働力人口」を同時に投入すると、両者の相関が r=0.99 になり、係数推定が極端に不安定になります。「総人口は正だが、労働力人口は負」のような解釈不能な結果になりがちです。

診断と対処:
VIF(分散拡大係数)を計算し、VIF > 10 の変数を確認
相関行列で |r| > 0.8 のペアをチェック
・対処法:一方を除外、合成変数(PCA)に変換、Ridge回帰で安定化
❌ 「R²が高い=良いモデル」ではない
決定係数 R² はモデルの「当てはまりの良さ」を示しますが、 が高くてもモデルが正しいとは限りません

が高くなる罠:
説明変数を増やせば は自動的に上がる(無関係な変数を追加してもは下がらない)
時系列データでは、共通のトレンド(時間とともに増加)があるだけで が 0.9 を超える
サンプルサイズが小さいとが過大評価される

代替指標: 調整済み (変数の数でペナルティ)AICBICモデル選択基準)を併用してください。予測力の真の評価には交差検証(cross-validation)テストデータ を見ること。
❌ 「ステップワイズで選んだ変数は重要」ではない
ステップワイズ法(バックワード・フォワード選択)は便利ですが、p値ベースの変数選択は再現性に問題があると批判されています。

問題点:
同じデータでも実行順序によって最終モデルが変わる
p値を繰り返し見ることで「偶然に有意な変数」を拾ってしまう(p-hacking
係数標準誤差が過小評価され、信頼区間が嘘っぽくなる

より良い方法:
事前に変数を理論で絞る(先行研究から候補を選ぶ)
LASSO回帰(自動かつ統計的に正当化された変数選択)を使う
交差検証AIC/BIC 最小モデルを選ぶ
❌ 「線形回帰なら線形関係を前提にすべき」
重回帰分析線形関係を前提とします。実際の関係が非線形なのに線形モデルで分析すると、本当の関係を見逃します

非線形の例:
U字型関係: 失業率と物価上昇率(フィリップス曲線)
逓減効果: 所得と幸福度(年収 800万円までは強い正の効果、それ以上は飽和)
閾値効果: 高齢化率と医療費(ある水準を超えると急激に上がる)

診断と対処:
残差プロット残差が0周辺に均等に分布しているか確認
変数の対数変換・二乗項追加で非線形性を取り込む
・どうしても線形では捉えられないなら、機械学習RF・GBM)を併用する
❌ 「データに当てはまった=予測に使える」ではない
「過去のデータでフィットしたから将来も予測できる」と思うのは危険です。

過学習(overfitting)の例: 47都道府県のデータに10個の説明変数を投入すれば、ほぼ完璧にフィットします(自由度がほぼゼロ)。でもそのモデルを新しい年度に適用すると、予測精度はほぼランダム並みに落ちることがあります。

正しい予測力の評価:
・データを訓練用 70%テスト用 30%に分割し、テスト用での予測精度を見る
k分割交差検証(k-fold CV)で予測の安定性を確認
・「説明変数の数 ≪ サンプルサイズ」のバランスを意識(目安:n > 10 × 変数数)

📖 用語集(この記事に出てくる統計用語)

統計の基本用語を初心者向けに解説します。本文中で見慣れない言葉が出てきたら、ここに戻って確認してください。

p値
「効果がない」と仮定したときに、観察されたデータ(またはより極端なデータ)が得られる確率。0〜1の値で、慣例的に 0.05(5%)未満を「有意」と判断する。
有意水準
「偶然」と「意味のある違い」を分ける基準。通常 α=0.05(5%)を使う。p値 < α なら「有意」と判定。
信頼区間
「真の値はこの範囲にあるだろう」という幅。95%信頼区間 = 同じ実験を100回繰り返したら95回はこの範囲に真の値が入る。
サンプルサイズ
分析に使ったデータ点の数(n)。一般にnが大きいほど推定が安定し、わずかな差も検出できるようになる。
標準誤差
推定値(係数など)のばらつきの目安。標準誤差が小さいほど推定値が安定している。
正規分布
釣鐘型の左右対称な分布。多くのパラメトリック検定(t検定F検定など)は「データが正規分布に従う」ことを仮定する。
因果相関
相関がある」と「原因と結果の関係(因果)」は別物。アイスクリームの売上と水難事故は相関するが、原因は両者とも「夏の暑さ」。
外れ値
他のデータから極端に離れた値。分析結果を歪める原因になるため、検出して除外するか別途扱う必要がある。
欠損値
データが取得できなかった部分(NaN・空白)。除外するか補完(平均代入・回帰代入など)するかが分析上の重要な判断点。
係数回帰係数
説明変数 x が1単位増えたとき、目的変数 y が平均でどれだけ変化するか」を示す数値。正の値は正の影響、負の値は負の影響。

📐 使っている手法をわかりやすく解説

統計手法について「何のためか」「結果をどう読むか」を初心者向けに解説します。

◆ 統計の基本概念(どの論文にも共通)

🔍 p値有意確率)とは
何?
「もし本当に効果がなかったとしたら、今回の結果(またはもっと極端な結果)が偶然起きる確率」のこと。
なぜ必要?
帰無仮説(「効果なし」の仮定)のもとで検定統計量の分布から計算する。
何がわかる?
「この関係は偶然ではなく、統計的に意味がある」と主張するための客観的な根拠になる。
読み方
p < 0.05(5%未満)を「統計的に有意」と判断するのが慣例。ただし「p値が小さい=効果が大きい」ではない。効果量係数の大きさ)とセットで判断する。
🗂️ ノンパラメトリック検定とは(なぜ使うのか)
何?
「データが正規分布に従う」という仮定を置かない検定手法の総称。Kruskal-Wallis検定・Mann-Whitney U検定などが代表例。
なぜ必要?
データの値ではなく「順位」に変換して検定統計量を計算する。外れ値や偏った分布に対しても安定して機能する。
何がわかる?
サンプルサイズが小さい・データが歪んでいる・外れ値がある場合でも、グループ差の有無を検定できる。
読み方
「なぜノンパラメトリックを選ぶのか」の理由を示すには、正規性検定(Shapiro-Wilk)の結果を添えるのが望ましい。結果の解釈は対応するパラメトリック検定と同様(p < 0.05 で有意差あり)。

◆ この論文で使われている手法

📐 OLS 重回帰分析+頑健標準誤差
何?
複数の説明変数で連続的な目的変数(本論文では「購入経験者割合」)を同時に説明する最も基本的な回帰手法。
どう使う?
sm.OLS(y, X).fit(cov_type='HC1')。係数は最小二乗法で推定し、標準誤差だけを不均一分散に頑健な推定量(HC1)に置き換える。
何がわかる?
「正答率が1%ポイント高い県は株式購入経験者割合が約0.9%ポイント高い」のように、他の変数を一定とした各要因の影響を読み取れる。
結果の読み方
係数の符号(+/−)と有意性(***: 1%、**: 5%、*: 10%)をセットで見る。事前の符号予想と一致するかの確認が重要。
⚠️ 注意点
(1) N=47 は小標本—説明変数を増やしすぎると自由度が失われる。(2) 集計データ—県レベルの関係を個人に当てはめると生態学的誤謬。(3) 割合が0-100%に有界—予測値が範囲を超えうる点は限界(分数ロジットが発展形)。
↔️ VAR(ベクトル自己回帰)/ Granger因果検定
何?
複数の時系列変数が互いに影響し合う関係を分析する手法(VAR)と、「AがBの予測に役立つか」を検定する手法(Granger因果)。
どう使う?
VARは全変数を互いに説明変数として同時回帰Granger因果F検定でAのラグ変数がBの予測精度を向上させるかを確認する。
何がわかる?
「女性就業率と出生率はどちらが先に動くか」「リード・ラグ関係」を特定できる。
結果の読み方
Granger因果 p < 0.05 → 「Aの過去値はBの予測に役立つ」(ただし真の因果とは限らない)。
⚠️ 注意点
(1) 「予測に役立つ」=真の因果ではない—Granger因果はあくまで予測的先行関係。(2) ラグ次数の選択に敏感—AIC/BICで選ぶ。(3) 非定常な時系列では見せかけの関係が出やすい—単位根検定・差分化が前提。

🚀 発展の可能性(結果 X → 新仮説 Y → 課題 Z)

この研究をさらに発展させるための3つの方向性を示します。「今回わかったこと(X)」から「次に検証すべき仮説(Y)」を立て、「具体的に何をするか(Z)」まで考えてみましょう。

① データ・時間的拡張
結果 X
本論文は特定の年度・地域の断面データ(または限られた時系列)で分析を行った。
新仮説 Y
より新しい年度のデータや市区町村レベルの細粒度データを使えば、知見の時間的頑健性や地域内格差を検証できる。
課題 Z
(1)統計センターから最新の SSDSE をダウンロードし、同じ分析を再実行する。(2)結果が変わった場合、その要因(コロナ・政策変化など)を考察する。(3)市区町村データ(SSDSE-A/C/F)で分析単位を細かくした場合の結果と比較する。
② 手法の発展:OLS 重回帰の次のステップ
結果 X
本論文は 47 都道府県の単年クロスセクションで OLS 重回帰(頑健標準誤差)を推定した。
新仮説 Y
DCI が複数年分公開されれば、パネルデータ化して固定効果モデルを使うことで、都道府県固有の未観測要因(金融機関の立地・県民性など)によるバイアスを除去できる(原論文自身が挙げている限界)。
課題 Z
(1)金融リテラシー調査 2016・2019・2022 の3時点をつないでパネルを構築する。(2)固定効果モデルと変量効果モデルを推定し Hausman 検定で選択する。(3)割合データの有界性に対応する分数ロジット(fractional logit)でも推定し、OLS の結果と比較する。
③ 政策提言・実践への応用
結果 X
本論文は分析結果から特定の変数が目的変数に影響することを示した。
新仮説 Y
分析対象を日本全国から特定地域に絞ること、または逆に国際比較に拡張することで、政策の移転可能性と文脈依存性を検証できる。
課題 Z
(1)有意な変数を「政策で変えられるもの」と「変えにくいもの」に分類する。(2)政策で変えられる変数について、係数の大きさから「どれだけ変えればどれだけ効果があるか」を試算する。(3)自治体・政策立案者への提言として、実現可能なアクションプランを1枚にまとめる。

🎯 自分でやってみよう(5つのチャレンジ)

学んだだけでは身につきません。実際に手を動かすのが最強の学習方法です。本論文のスクリプトをベースに、以下のチャレンジに挑戦してみてください。難易度別に5つ用意しました。

★☆☆☆☆ 入門
CH1. 同じデータで分析を再現する
まずは付属の Python スクリプトをそのまま実行し、論文と同じ図を再現してみてください。
ポイント: 各図がどのコード行から生成されているか辿る。エラーが出たら原因を考える。
★★☆☆☆ 初級
CH2. 説明変数を1つ追加・除外して結果を比較
本論文の分析モデルから説明変数を1つ抜いて再実行、あるいは1つ追加して再実行してください。
ポイント: 係数p値 がどう変わったか観察する。多重共線性が原因で結果が変わる例を見つけられたら理想的。
★★★☆☆ 中級
CH3. 別の年度・別の都道府県で同じ分析を試す
SSDSE の別の年度(例:2015年度・2020年度)または特定都道府県のみのデータで同じ分析を実行してください。
ポイント: 時代や地域によって結論が変わるか? 変わるならその理由を考察する。
★★★★☆ 上級
CH4. 別の手法を組み合わせる
本論文の手法 + 1つの追加手法(例:重回帰 + LASSO相関分析 + 主成分分析)で結果を比較してください。
ポイント: 手法の違いで結論が変わるか? どちらが妥当かを「なぜ」とともに説明できるように。
★★★★★ 発展
CH5. オリジナルの問いを立てて分析する
本論文の手法を借りて、あなた自身の問いを立てて分析してください。 例:「カフェの数と幸福度に関連はあるか」「教育費の高い県は出生率も高いか」など。
ポイント: 問い・データ・手法・結論を1ページのレポートにまとめる。これがデータサイエンスの「実践」。
💡 ヒント: 詰まったら本サイトの他の論文(同じ手法を使っている)のスクリプトをコピーして組み合わせるのが効率的です。手法ガイド・用語集も参考に。

💼 この手法は実社会でこう使われている

本論文で学んだ手法は、研究の世界だけでなく、行政・企業・NPO の現場でも様々に活用されています。具体的なシーンを紹介します。

🏛️
行政の政策立案
都道府県・市区町村の政策担当者は、本論文と同様のデータ分析を用いて「どこに予算を投じれば効果が出るか」を検討します。 例えば医療費削減策、移住促進策、子育て支援策などの効果予測・効果検証に直結します。
🏢
企業のマーケティング・出店戦略
小売チェーン・サービス業の出店戦略では、地域特性(人口構成、所得、ライフスタイル)と売上の関係を本論文と同じ手法で分析します。 ECサイトでも顧客セグメント分析・購買要因分析に類似手法が使われます。
🏥
医療・公衆衛生
感染症の流行予測、医療資源配分の最適化、健康格差の地域要因分析などで、本論文の統計手法は標準的に使われています。 WHO・厚労省レベルの政策評価でも同じ手法が活躍しています。
📊
メディア・ジャーナリズム
新聞・テレビの社会調査記事、選挙予測、世論調査の分析でも、本論文と同じ手法(回帰分析・クラスタリングなど)が使われています。 データジャーナリズムの記事はこの種の分析が中核です。
🎓
学術研究(隣接分野)
経済学・社会学・公衆衛生学・教育学・地理学などの実証研究では、本論文と同じ手法が日常的に使われます。 専門誌に掲載される論文の8割以上が、こうした統計手法に基づいて結論を出しています。
💰
金融・保険業界
与信判断(融資審査)、保険料の地域別設定、不動産価格予測などで、本論文と同様のモデリング手法が広く活用されています。 統計分析の能力は金融業界の必須スキルになっています。

🤔 よくある質問(読者からの想定Q&A)

この論文を読んで初心者が抱きやすい疑問に、教育的観点から答えます。

Q1. この分析、自分でもできますか?
はい、できます。SSDSE データは無料で公開されており、Python の pandas, scikit-learn, statsmodels を使えば全く同じ手順で再現可能です。本ページ下部のスクリプトを実行するだけで結果が得られます。
Q2. 使われている手法は他の分野にも応用できますか?
十分応用可能です。本論文の手法は、医療・教育・経済・環境など他のドメインでも標準的に使われる手法です。データの中身(変数)を入れ替えるだけで、別の問いにも適用できます。
Q3. 結論は本当に「因果関係」を示していますか?
本論文は「観察データ」を使った分析であり、厳密な意味での「因果関係」を完全に証明したわけではありません。あくまで「強い関連が見られた」という事実を提示しているにとどまります。真の因果を示すには、無作為化比較試験(RCT)か、自然実験を活用したIVDiD 等の手法が必要です。
Q4. データの最新版を使うとどうなりますか?
SSDSE は毎年更新されているため、最新版を使えば近年のトレンド(特にコロナ禍以降の変化)も含めて分析できます。ただし、結論が変わる可能性もあります。それ自体が新しい発見につながります。
Q5. もっと深く学ぶには何を読めばいいですか?
「計量経済学」「データサイエンス入門」「統計的因果推論」などのテキストが入門に向いています。Python の場合は『Python ではじめる機械学習』(オライリー)、R の場合は『R で学ぶ統計学』が定番です。本サイトの他の論文も読み比べてみてください。

✅ 理解度チェック(4問)

この論文を理解できたか、クリックで確認しましょう。間違えても解説が出ます。

🐍 ブラウザで動かす — インストール不要でこの論文の分析を再現する

このページの分析は、この画面の中でそのまま実行できます。 Python をインストールする必要も、CSV をダウンロードする必要もありません。 下のセルの 「▶ ブラウザで実行」 を上から順に押すか、 「▶ 最初から全部実行」 で一気に流してください。 表示されるのは、本文の図表とまったく同じ計算の結果です (動かしているのは再現スクリプト code/2024_U4_katsuyo.py そのもの)。

コードは書き換えられます。「✏️ 書き換える」で数字や変数を変えて ▶ を押すと、 結果がどう変わるかをその場で確かめられます(元に戻すボタンつき)。 ページを開いた時点で裏で実行環境の準備が始まっているので、待ち時間はほとんどありません。 初回だけ実行環境の取得にインターネット接続が必要です。