🎯 この記事を読むと何ができるようになるか
- 研究の核心:観光資源の「量」を増やすことが必ずしも観光消費を増やさない、という反直観的な発見の読み解き方
- 分析手法:パネルデータで固定効果/変量効果モデルを使い分け、観察されない「都道府県の固有効果」を統制する考え方
- 分析手法:Hausman検定・F検定で「どのモデルを採用すべきか」を統計的に判断する手順
- 結果の読み方:Pooled OLSと固有効果モデルで係数の符号が反転する現象(欠落変数バイアス)を見抜く力
- 応用:「報告値の可視化」と「実再現」の違いを理解し、再現可能性を意識した分析設計ができる
📥 データ・図の再現について(重要)
このページの図は「報告値の可視化」です
本研究の被説明変数「観光消費額」は、国土交通省観光庁『共通基準による観光入込客統計』から作成されており、SSDSE(教育用標準データセット)には収録されていません。そのため本ページの図1〜4は回帰の再計算ではなく、原論文(藤原2020)の表5・表1の報告値をそのまま可視化したものです。新しい数値は一切生成していません。
1
外部データのダウンロードは不要です
再現スクリプト code/2020_U1_daijin.py には、原論文の報告値(係数・標準誤差・都道府県数)がそのまま埋め込まれています。CSVの配置は不要です。
2
スクリプトをそのまま実行する
ターミナルでプロジェクトルートに移動し、以下を実行します。
python3 code/2020_U1_daijin.py
図は
html/figures/ に自動保存され、コンソールには表5・表1の報告値が出力されます。
3
SSDSEで「実再現」できる部分
観光インフラの代理変数(旅館・ホテルの営業施設数)や統制変数(年平均気温・年間降水日数)は SSDSE-2020B に収録されています。被説明変数さえ入手できれば回帰全体を実再現できますが、観光消費額は公的オープンデータでの都道府県レベル整備が限られる点に注意してください。
観光産業は大きな成長が見込まれる分野である。原論文によれば、2018年の日本国内旅行消費額は20.5兆円、観光がもたらす生産波及効果は46.7兆円(対国民経済計算産出額5.3%)、雇用誘発効果は441万人(全国就業者数の6.4%)とされる。観光は関連消費を誘発するため、日本経済・地域経済の活性化に必要な産業のひとつと位置づけられる。
本研究の出発点は、先行研究尾崎(2016)への批判である。尾崎(2016)は2012年の都道府県データを単年のOLSで分析し、観光資源の充実度が観光消費額に及ぼす影響は限定的だと結論づけた。しかしこの分析には2つの課題が残されていた。
先行研究に残された2つの課題
① 観察されない異質性(固有効果): 各都道府県には「観光資源を発掘・開発する能力(=観光振興力)」のような、データに現れない固有の特性がある。これを統制しない単年OLSでは、説明変数と誤差項が相関し、内生性(欠落変数バイアス)が生じうる。
② トレンドの未考慮: 観光消費は景気・自治体の政策・自然災害の影響を受けやすい。単年(横断)データでは時間を通じた変化を制御できない。
そこで本研究は、2010年から2017年までの都道府県パネルデータを構築し、パネルデータ分析によって都道府県の固有効果と年次効果を統制したうえで、観光資源および観光インフラの充実度が観光消費額に及ぼす影響を検討する。
リサーチクエスチョン
観察されない都道府県の固有効果を統制してもなお、観光資源(自然・歴史・文化…)や観光インフラ(旅館・ホテル)の「量」は観光消費額を押し上げるのか? また、宿泊客と日帰り客で効く資源は異なるのか?
パネルデータ
固定効果/変量効果
Hausman検定
観光経済
SSDSE-2020B
データと変数
使用データ
本研究は、教育用標準データセット(SSDSE-2020B、都道府県別・時系列)と、国土交通省観光庁『共通基準による観光入込客統計』を組み合わせてデータセットを構築した。分析期間は2010〜2017年、対象は大阪府を除く46都道府県である(大阪府は当該統計を導入していないため除外)。導入時期の違い等により一部に欠損があり、データはアンバランスパネル(最終サンプル N=331)となった。
再現可能性の注記(データソースの区別)
記述統計(表3)の脚注どおり、旅館・ホテルの営業施設数、年平均気温、年間降水日数は SSDSE 由来(実再現が可能)である一方、観光消費額と観光地点数(歴史・文化ほか)は『共通基準による観光入込客統計』由来で SSDSE には収録されていない。本ページの図はこのうち報告値の可視化にあたる。
変数の定義
被説明変数は2種類の観光消費額を用いる。観光消費額 = 観光入込客数(千人回)× 1人当たり観光消費額(円/人回)で、当該都道府県を訪れた観光入込客の消費総額と定義される(分析では対数変換)。
| 変数区分 | 変数 | 定義・出典 |
| 被説明変数 | 県外宿泊者の観光消費額 | log変換。共通基準による観光入込客統計 |
| 県外日帰り者の観光消費額 | log変換。共通基準による観光入込客統計 |
| 観光資源の充実度 | 観光地点数(7分類) | 自然/歴史・文化/温泉・健康/スポーツ・レクリエーション/都市型観光/行事・イベント/その他(下表・表2) |
| 観光インフラの充実度 | 旅館営業施設数・ホテル営業施設数 | SSDSE(実再現が可能) |
| 統制変数 | 年平均気温・年間降水日数 | SSDSE(実再現が可能) |
表2 観光地点の定義
| 観光地点 | 定義(例) |
| 自然 | 山岳、高原、湖沼、河川、海岸、海中、島 |
| 歴史・文化 | 史跡、城、神社・仏閣、庭園、歴史的まち並み、博物館、美術館、動植物園、水族館、産業観光 等 |
| 温泉・健康 | 温泉法に基づく温泉施設、スーパー銭湯などの温泉類似施設 |
| スポーツ・レクリエーション | ゴルフ場、テニス場、スキー場、キャンプ場、海水浴場、公園、テーマパーク 等 |
| 都市型観光 | 商業施設、商店街、食の観光拠点、農水産品の直売所、物産館 等 |
| 行事・イベント | 行・祭事、花見、初詣、花火大会、郷土芸能、博覧会、コンサート、スポーツ観戦、国際会議 等 |
| その他 | 上記に分類されない観光地点。道の駅、パーキングエリア 等 |
アンバランスパネルとは
すべての個体(都道府県)がすべての時点で観測されているパネルを「バランスパネル」、一部に欠測があるものを「アンバランスパネル」と呼ぶ。本研究は年ごとに観測県数が37〜43と異なるアンバランスパネルであり、この欠測の非ランダム性はサンプルセレクションバイアスの懸念として後述される。
記述統計と相関
説明変数のGVIF(一般化分散拡大係数)の最大値は4.797(スポーツ・レクリエーション)であり、目安の10を大きく下回るため、多重共線性の問題は深刻ではないと判断された。
表3 記述統計量(原論文の報告値)
| 変数 | 平均 | 標準偏差 | 最小 | 最大 | GVIF |
| log 県外宿泊者の観光消費額 | 11.017 | 0.811 | 9.352 | 12.857 | — |
| log 県外日帰り者の観光消費額 | 10.936 | 1.289 | 3.219 | 14.470 | — |
| 自然 | 34.218 | 20.959 | 3 | 112 | 4.075 |
| 歴史・文化 | 94.885 | 59.916 | 19 | 399 | 3.993 |
| 温泉・健康 | 39.692 | 27.411 | 2 | 149 | 4.362 |
| スポーツ・レクリエーション | 100.517 | 75.629 | 13 | 370 | 4.797 |
| 都市型観光 | 27.378 | 19.962 | 1 | 93 | 2.186 |
| 行事・イベント | 136.641 | 92.563 | 16 | 500 | 3.457 |
| その他 | 27.674 | 18.740 | 1 | 114 | 2.274 |
| 旅館営業施設数 | 1135.109 | 715.357 | 243 | 3602 | 4.127 |
| ホテル営業施設数 | 208.372 | 153.313 | 34 | 718 | 4.455 |
| 年平均気温 | 15.382 | 2.186 | 9.1 | 23.1 | 2.092 |
| 年間降水日数 | 118.740 | 27.422 | 67 | 197 | 1.313 |
📌 相関行列(表4)の要点
- 宿泊者消費額
- 旅館営業施設数(r=.726)・ホテル営業施設数(r=.607)と強く正相関。宿泊需要が施設数と結びつく。
- 日帰り者消費額
- 歴史・文化(r=.541)・スポーツ・レクリエーション(r=.332)と相関。宿泊者ほど施設数と結びつかない。
- 注意
- これらは固有効果を統制する前の単純相関。後述のパネル推定では符号が変わる変数がある(=相関だけで因果を語れない好例)。
パネルデータ分析の利点は、観察不可能な都道府県の固有効果 Fi(観察されない異質性)を識別できる点にある。推定モデルとして固定効果モデルと変量効果モデルを用いる。
固定効果(FE)と 変量効果(RE)の違い
- 固定効果 FE
- 「固有効果 Fi が説明変数と相関してもよい」と仮定するモデル。各都道府県の期間平均を引く within変換 で Fi を消去する。
- 変量効果 RE
- 「固有効果 Fi は説明変数と独立」と仮定し、Fi を誤差項に含めて(誤差の自己相関を考慮したGLMで)推定するより効率的なモデル。(山本2015)
モデルの定式化
基本式:Yit = b·X1,it + b·X2,it + Fi + vit (i = 都道府県, t = 時点)
FE(within変換):(Yit − Ȳi) = b(X1,it − X̄1,i) + b(X2,it − X̄2,i) + (vit − v̄i)
→ 時間で不変な Fi は期間平均との差でゼロになり消去される
RE:Yit = b·X1,it + b·X2,it + (Fi + vit) → 誤差項の自己相関を考慮したGLM推定
どのモデルを採用するか?(2段階の検定)
① F検定(poolability):Pooled OLS と 固定効果 のどちらか。有意なら「固有効果が存在する」→ 固定効果が OLS より妥当。
② Hausman検定:固定効果 と 変量効果 のどちらか。有意なら「Fi と説明変数に相関あり」→ 固定効果を採用。非有意なら相関なし → より効率的な変量効果を採用。
標準誤差は都道府県内の系列相関に頑健なクラスターロバスト標準誤差を用い、係数の検定は片側確率による(原論文の方針)。
① 県外宿泊者の観光消費額 → 変量効果モデルを採択
モデルの全体評価は、Pooled OLS(F=56.011, p<.01)、固定効果(F=7.600, p<.01)、変量効果(χ²=215.772, p<.01)でいずれも1%有意。モデル選択は次のとおり。
モデル選択(宿泊者)
・F検定:F=56.362(1%有意)→ 固有効果が存在 → 固定効果が OLS より妥当。
・Hausman検定:χ²=19.004(非有意)→ 固有効果と説明変数に相関なし → 変量効果モデルを採択。
📌 この回帰係数プロットの読み方
- このグラフは
- 各説明変数の回帰係数(影響の向きと強さ)を水平バーで示したもの。
- 読み方
- 0より右は正の影響、左は負の影響。誤差棒が0をまたがない変数が有意。
- 注意
- これは再計算ではなく、原論文が報告した係数と標準誤差を転記して描いた図。
宿泊者:主な発見(変量効果モデル・原論文の報告値)
- 「スポーツ・レクリエーション」(b=.0033, p<.01)・「その他」(b=.0072, p<.01)が正の影響
- 「歴史・文化」(b=−.0030, p<.01)・「行事・イベント」(b=−.0009, p<.05)が負の影響
- 観光インフラは「ホテル営業施設数」(b=.0033, p<.01)が正、「旅館営業施設数」(b=−.0004, p<.01)が負
- Pooled OLS では非有意だった「歴史・文化」が、固有効果統制後は負に有意化。「旅館営業施設数」は OLS で正有意 → RE で負に反転。
② 県外日帰り者の観光消費額 → 固定効果モデルを採択
Pooled OLS(F=35.735, p<.01)・固定効果(F=3.140, p<.01)・変量効果(χ²=52.502, p<.01)いずれも1%有意。モデル選択は次のとおり。
モデル選択(日帰り者)
・F検定:F=36.809(1%有意)→ 固有効果が存在。
・Hausman検定:χ²=62.320(1%有意)→ 固有効果と説明変数に相関あり → 固定効果モデルを採択。
日帰り者:主な発見(固定効果モデル・原論文の報告値)
- 「自然」(b=.0090, p<.05)・「スポーツ・レクリエーション」(b=.0038, p<.05)・「その他」(b=.0071, p<.10)が正
- 「歴史・文化」(b=−.0036, p<.10)が負
- 観光インフラは「ホテル営業施設数」(b=.0048, p<.05)が正、「旅館営業施設数」(b=−.0007, p<.01)が負
- 宿泊者と異なり「自然」が正に効く一方、「行事・イベント」は非有意。効く資源が客層で異なる。
③ Pooled OLS と 固有効果モデルの「符号反転」
本研究の核心は、固有効果を統制すると単純OLSと係数の符号が変わる変数がある点である。これは尾崎(2016)が固有効果を無視していたために生じた欠落変数バイアスを示唆する。
なぜ符号が反転するのか(教育的ポイント)
「観光振興力の高い県ほど歴史・文化の観光地点も多く登録している」といった見えない共通要因(固有効果)があると、単純OLSではその効果が歴史・文化の係数に紛れ込み、見かけ上プラスに見える。固有効果を差し引くと、この見かけの相関が消え、真の(負の)関係が現れる——これが欠落変数バイアスの典型である。
表5 分析結果(原論文の報告値の転記)
係数の下段( )内はクラスターロバスト標準誤差。***=1%・**=5%・*=10%水準で有意(片側確率)。採択モデルは宿泊者=変量効果、日帰り者=固定効果。
| 説明変数 | 県外宿泊者の観光消費額 | 県外日帰り者の観光消費額 |
| OLS | 固定効果 | 変量効果★ | OLS | 固定効果★ | 変量効果 |
| 切片 | 9.2693*** (.9667) | — | 9.4320*** (.4382) | 14.2439*** (2.0572) | — | 9.2370*** (1.0143) |
| 自然 | −.0024 (.0061) | −.0007 (.0019) | −.0006 (.0019) | −.0171** (.0074) | .0090** (.0088) | .0032 (.0063) |
| 歴史・文化 | .0006 (.0018) | −.0034*** (.0010) | −.0030*** (.0010) | .0188*** (.0039) | −.0036* (.0022) | .0003 (.0024) |
| 温泉・健康 | −.0024 (.0052) | .0028 (.0035) | .0030 (.0033) | −.0188*** (.0068) | −.0048 (.0079) | −.0083* (.0066) |
| スポーツ・レクリエーション | .0040*** (.0028) | .0027*** (.0008) | .0033*** (.0007) | .0010 (.0022) | .0038** (.0018) | .0041** (.0015) |
| 都市型観光 | −.0022 (.0036) | .0006 (.0018) | .0010 (.0018) | −.0039 (.0085) | .0032 (.0031) | .0060* (.0035) |
| 行事・イベント | −.0039*** (.0013) | −.0007** (.0003) | −.0009** (.0003) | .0025* (.0019) | −.0003 (.0006) | −.0007 (.0007) |
| その他 | .0064*** (.0036) | .0060*** (.0022) | .0072*** (.0019) | −.0081** (.0092) | .0071* (.0041) | .0074** (.0046) |
| 旅館営業施設数 | .0007*** (.0002) | −.0006*** (.0001) | −.0004*** (.0001) | .0008*** (.0003) | −.0007*** (.0003) | −.0004* (.0003) |
| ホテル営業施設数 | .0014 (.0009) | .0025*** (.0010) | .0033*** (.0008) | −.0030* (.0023) | .0048** (.0021) | .0023* (.0020) |
| 年平均気温 | .0466*** (.0466) | .0573** (.0286) | .0654*** (.0242) | −.1842*** (.1067) | .0881** (.0469) | .0459 (.0724) |
| 年間降水日数 | .0009 (.0022) | .0013* (.0009) | .0011 (.0009) | −.0083** (.0043) | .0043** (.0021) | .0030* (.0018) |
| サンプルサイズ | 331 | 331 | 331 | 331 | 331 | 331 |
| 自由度調整済み決定係数 | .647 | .077 | .384 | .537 | .118 | .116 |
★=採択モデル。原論文では、宿泊者・日帰り者の両方で有意に正の変数を黄色、それ以外で正を赤、負を青でマーカー表示している。
この表5の報告値を配列に転記して可視化したのが、上の図2〜図4です。再現スクリプトを追ってみましょう。
📝 コード
# 目的: 被説明変数はSSDSE未収録のため、回帰は再計算せず
# 原論文 表5 の「県外宿泊者・変量効果(採択)」の係数・SEを転記する。
labels = ['自然','歴史・文化','温泉・健康','スポーツ・レクリエーション',
'都市型観光','行事・イベント','その他','旅館営業施設数','ホテル営業施設数']
stay_re = [-.0006,-.0030,.0030,.0033,.0010,-.0009,.0072,-.0004,.0033]
stay_re_se = [ .0019, .0010,.0033,.0007,.0018, .0003,.0019, .0001,.0008]
stay_re_s = ['','***','','***','','**','***','***','***'] # ''=n.s.
print("=== 県外宿泊者の観光消費額(変量効果=採択/原論文の報告値) ===")
for lab, c, se, st in zip(labels, stay_re, stay_re_se, stay_re_s):
print(f"{lab:<16}{c:>9.4f} ({se:.4f}) {st or 'n.s.'}") |
▼ 実行結果
=== 県外宿泊者の観光消費額(変量効果=採択/原論文の報告値) ===
自然 -0.0006 (0.0019) n.s.
歴史・文化 -0.0030 (0.0010) ***
温泉・健康 0.0030 (0.0033) n.s.
スポーツ・レクリエーション 0.0033 (0.0007) ***
都市型観光 0.0010 (0.0018) n.s.
行事・イベント -0.0009 (0.0003) **
その他 0.0072 (0.0019) ***
旅館営業施設数 -0.0004 (0.0001) ***
ホテル営業施設数 0.0033 (0.0008) ***
💡 解説
- 目的:回帰の再計算ではなく、原論文が報告した係数と標準誤差を配列に転記する(=報告値の可視化の下ごしらえ)。
- 実行結果の読み取り:「その他」(+.0072)と「スポーツ・レクリエーション」(+.0033)が最も強い正、「歴史・文化」(−.0030)が負。有意記号
*** は1%水準(片側)を示す。
📝 コード
# 目的: どのモデルを採択したかを、原論文が報告した検定統計量で確認する。
print("[県外宿泊者] F検定 F=56.362 → 固有効果あり(FE>OLS)")
print(" Hausman χ²=19.004 → 非有意 → 変量効果(RE)を採択")
print("[県外日帰り者] F検定 F=36.809 → 固有効果あり(FE>OLS)")
print(" Hausman χ²=62.320(p<.01) → 有意 → 固定効果(FE)を採択") |
▼ 実行結果
[県外宿泊者] F検定 F=56.362 → 固有効果あり(FE>OLS)
Hausman χ²=19.004 → 非有意 → 変量効果(RE)を採択
[県外日帰り者] F検定 F=36.809 → 固有効果あり(FE>OLS)
Hausman χ²=62.320(p<.01) → 有意 → 固定効果(FE)を採択
💡 解説
- 目的:2段階の検定(F検定→Hausman検定)の結論を明示し、なぜ客層で採択モデルが違うのかを示す。
- 実行結果の読み取り:両客層とも固有効果は存在する(F検定有意)。宿泊者は Hausman が非有意なので効率的な変量効果、日帰り者は有意なので固定効果を採る。同じデータでも被説明変数が変われば採択モデルが変わる好例。
考察
原論文は考察を3点に整理している。
発見1:観光地点を増やせば消費が増える、とは限らない
固有効果を統制すると、宿泊者消費額に対して「スポーツ・レクリエーション」「その他」は正、「歴史・文化」「行事・イベント」は負に効いた。尾崎(2016)は「歴史・文化」を正と結論づけたが、本研究では逆になる。神社・仏閣などは古くからある地域資源だが、観光地点として登録されても、周辺に飲食店や土産屋がなければ消費喚起にはつながらない。むしろ消費に結びつかない地点への訪問が増えることで、別の場所で使われたはずの支出が失われた可能性がある。資源の「量」を増やすこと自体が負に働きうる点が本研究の大きな貢献である。
発見2:宿泊客と日帰り客で効く資源が異なる
「自然」は日帰り者消費額に正(b=.0090)だが宿泊者には非有意だった。山岳・海岸などの自然は短時間で楽しめるため、日帰り観光と相性がよい。一方、宿泊を促すには長時間の滞在を前提とした施設・サービスが必要で、自然が多いだけでは宿泊客は増えにくい。また「行事・イベント」は宿泊者に負(初詣・花火大会など地域住民向けが多く県外客には魅力が乏しい/交通渋滞で敬遠される可能性)だが、日帰り者には非有意だった。
発見3:ホテルは正、旅館は負
宿泊者消費額に対し「ホテル営業施設数」は正、「旅館営業施設数」は負に効いた。旅館は高級旅館から民宿まで施設・サービスの品質のばらつきが大きく、施設数では観光インフラの充実度をうまく測れなかった可能性がある。品質のばらつきの大きい旅館への宿泊を懸念した結果、旅館数が負に働いたと解釈される。
核心メッセージ
「観光資源を増やせば観光消費が増える」という素朴な想定は、都道府県の固有効果を統制すると崩れる。政策的には、資源の量ではなく、消費に結びつく周辺環境(飲食・宿泊の質)や、狙う客層(宿泊/日帰り)に応じた資源選択が重要になる。
まとめ
本研究の貢献(2点)
| 貢献 | 内容 |
| ① パネルで固有効果・年次効果を統制 | 既存研究(尾崎2016の単年OLS)が抱えた内生性・トレンド未考慮を克服し、「歴史・文化」「行事・イベント」が宿泊者消費額に負という、既存研究では見えなかった知見を得た。 |
| ② 被説明変数を客層で分割 | 「県外宿泊者」と「県外日帰り者」の観光消費額を区別し、消費主体によって喚起する観光資源が異なることを明らかにした。 |
限界(3点)
| 限界 | 内容 |
| ① 資源の「質」未考慮 | 観光資源の充実度を量(地点数)でしか測っておらず、人気の高低など質的側面を捉えていない。世界遺産・世界ジオパーク・日本遺産などを参考に質を測るべき。 |
| ② 都道府県レベルの集計 | 公的集計データが都道府県レベルに限られる。市町村レベルの集計や観光客の個票データの収集が望まれる。 |
| ③ サンプルセレクションバイアス | 一部欠測を含むアンバランスパネルであり、欠損がランダムでない場合はサンプルセレクションバイアスが生じうる。 |
教育的価値(この分析から学べること)
- 固有効果の統制:単純相関やOLSでは見えない関係が、within変換(固定効果)で符号ごと変わりうる。相関≠因果を体感できる好例。
- モデル選択の作法:F検定→Hausman検定という2段階で、Pooled OLS / FE / RE を統計的に選ぶ手順が学べる。
- 再現可能性の意識:被説明変数がSSDSE未収録のため「実再現」ではなく「報告値の可視化」に留まる——データ出所を明示する誠実さが学べる。
SSDSE-2020B
固定効果/変量効果
Hausman検定
観光消費額
パネルデータ
総務大臣賞 2020
⚠️ よくある誤解と注意点
統計分析の解釈で初心者がやりがちな勘違いをまとめます。特に「相関と因果の混同」「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² はモデルの「当てはまりの良さ」を示しますが、R² が高くてもモデルが正しいとは限りません。
R² が高くなる罠:
・説明変数を増やせば R² は自動的に上がる(無関係な変数を追加してもR²は下がらない)
・時系列データでは、共通のトレンド(時間とともに増加)があるだけで R² が 0.9 を超える
・サンプルサイズが小さいとR²が過大評価される
代替指標: 調整済み R²(変数の数でペナルティ)、AIC・BIC(モデル選択基準)を併用してください。予測力の真の評価には交差検証(cross-validation)でテストデータの R² を見ること。
❌ 「ステップワイズで選んだ変数は重要」ではない
ステップワイズ法(バックワード・フォワード選択)は便利ですが、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・空白)。除外するか補完(平均代入・回帰代入など)するかが分析上の重要な判断点。
- VIF
- Variance Inflation Factor(分散拡大係数)。多重共線性の強さを示す指標。VIF > 10 で「強い多重共線性あり」と判断。
- 交絡変数
- 「真の原因」と「結果」の両方に影響する第三の変数。これを統制しないと、見かけ上の関係を真の因果と誤認する。
- 係数(回帰係数)
- 「説明変数 x が1単位増えたとき、目的変数 y が平均でどれだけ変化するか」を示す数値。正の値は正の影響、負の値は負の影響。
- 内生性
- 説明変数と誤差項が相関している状態。逆因果や交絡変数の存在で発生する。これを放置すると係数推定にバイアスが生じる。
- 多重共線性
- 説明変数同士の相関が強すぎる状態。係数推定が不安定になり、解釈を誤る原因になる。VIF > 10 が警告サイン。
- 標準化係数
- 変数の単位の影響を取り除いた係数。複数の変数の影響の大きさを単位に依存せず比較するために使う。
- 決定係数 R²
- 回帰モデルが目的変数のばらつきの何%を説明できるかを示す指標。0〜1の値で、1に近いほどモデルの説明力が高い。
- パネルデータ
- 同じ個体(都道府県など)を複数の時点で観測したデータ。横断面と時系列の両方の情報を持ち、観察されない個体固有の差を統制できる。
- 固定効果モデル(FE)
- 各個体の期間平均を引く within 変換で、時間を通じて不変な固有効果を消去して推定するパネル回帰。固有効果が説明変数と相関していてもよい。
- 変量効果モデル(RE)
- 固有効果を誤差項に含めて推定するパネル回帰。「固有効果は説明変数と独立」を仮定し、FEより効率的。仮定の成否はHausman検定で確認する。
- Hausman検定
- 固定効果と変量効果のどちらを採用すべきかを判断する検定。有意なら固定効果、非有意なら変量効果を選ぶ。
- アンバランスパネル
- 一部の個体・時点に欠測があるパネルデータ。欠損がランダムでないとサンプルセレクションバイアスの懸念が生じる。
📐 使っている手法をわかりやすく解説
この論文で使われている統計手法について「何のためか」「結果をどう読むか」を初心者向けに解説します。注意点は手法ごとに異なる点に絞って示します。
◆ 統計の基本概念(どの論文にも共通)
🔍 p値(有意確率)とは
- 何?
- 「もし本当に効果がなかったとしたら、今回の結果(またはもっと極端な結果)が偶然起きる確率」のこと。
- なぜ必要?
- 帰無仮説(「効果なし」の仮定)のもとで検定統計量の分布から計算する。
- 読み方
- p < 0.05 を「統計的に有意」と判断するのが慣例。本論文は片側確率を用い、*** = 1%、** = 5%、* = 10% 水準を表す。
🏛️ なぜ「パネルデータ」を使うのか
- 何?
- 同じ個体(都道府県)を複数時点で観測したデータ。横断データ(1時点)と時系列(1個体)の両方の情報を持つ。
- なぜ必要?
- 観察できない個体固有の特性(都道府県の「観光振興力」など)を統計的に取り除ける。これが単年のOLSにはできない最大の利点。
- 読み方
- 固有効果を統制すると、単純相関やOLSと係数の符号が変わることがある。これは欠落変数バイアスが是正された証拠と解釈できる。
◆ この論文で使われている手法
🏛️ 固定効果モデル(FE)
- 何?
- 個体ごとの期間平均を引くwithin変換で、時間を通じて不変な固有効果 Fi を消去して推定するパネル回帰。
- どう使う?
- 「固有効果が説明変数と相関してもよい」と仮定するとき使う。本論文では日帰り者消費額の分析で採択された。
- 何がわかる?
- 「東京だから消費が多い」ではなく「この資源が消費を動かした」という、固有特性を除いた効果を推定できる。
- ⚠️ 注意点
- 時間で不変な変数(例:面積・地理)は within 変換で消えて推定できない。決定係数(within R²)は小さくなりがちで、OLSのR²と直接は比べられない(本論文でも宿泊FEのadjR²=.077と低い)。
🎲 変量効果モデル(RE)
- 何?
- 固有効果 Fi を誤差項に含め、誤差の自己相関を考慮したGLS/GLMで推定するパネル回帰。
- どう使う?
- 「固有効果は説明変数と独立」と仮定できるとき使う。FEより効率的(推定の分散が小さい)。本論文では宿泊者消費額で採択。
- 何がわかる?
- 時間で不変な変数の効果も推定でき、係数の標準誤差が小さくなりやすい。
- ⚠️ 注意点
- 「固有効果と説明変数が無相関」という仮定が崩れると係数にバイアスが出る。この仮定の成否をHausman検定で必ず確認する。
⚖️ Hausman検定
- 何?
- 固定効果と変量効果のどちらを採用すべきかを判断する検定。
- どう使う?
- 両モデルの係数の差を検定統計量(χ²)にする。有意(p<.05)なら「固有効果と説明変数に相関あり」→ 固定効果、非有意なら → より効率的な変量効果を採る。
- 何がわかる?
- 本論文では宿泊者 χ²=19.004(非有意)→ RE、日帰り者 χ²=62.320(1%有意)→ FE と、客層で採択が分かれた。
- ⚠️ 注意点
- 非有意は「REで問題ない」ことを示すだけで「REが正しい」ことの証明ではない。検定前に F検定で固有効果の存在自体を確認しておく。
📊 F検定(poolability)
- 何?
- Pooled OLS(全個体をひとまとめ)と固定効果のどちらが妥当かを判断する検定。
- どう使う?
- 「すべての固有効果が0」を帰無仮説にF検定する。有意なら固有効果が存在し、OLSでまとめるのは不適切。
- 何がわかる?
- 本論文では宿泊 F=56.362、日帰り F=36.809 と両方1%有意 → 固有効果は確かに存在する。
- ⚠️ 注意点
- F検定はFE対OLSの判断まで。FE対REの判断はHausman検定が担う(役割分担を混同しない)。
🛡️ クラスターロバスト標準誤差
- 何?
- 同一クラスター(同じ都道府県)内で誤差が相関・不均一でも、正しく推定できる標準誤差。
- どう使う?
- パネルでは同じ県の異なる年の誤差が相関しやすい。都道府県をクラスターに指定して標準誤差を計算する。
- 何がわかる?
- 係数の値は変わらないが、標準誤差・p値がより現実的になる(多くは標準誤差が大きくなる=有意判定が保守的に)。
- ⚠️ 注意点
- クラスター数(本論文は最大46都道府県)が少ないと近似が悪化する。本論文の有意記号は片側確率で付いている点に注意。
🔗 GVIF / 多重共線性
- 何?
- 多重共線性(説明変数同士の強い相関)の程度を測る指標。カテゴリ変数にも使える一般化版が GVIF。
- どう使う?
- 各説明変数についてGVIFを計算し、目安(10、厳しめなら5)を超えないか確認する。
- 何がわかる?
- 本論文の最大GVIFは4.797で、多重共線性は深刻でないと判断された。
- ⚠️ 注意点
- GVIFが低くても、理論的に重複する変数(旅館数とホテル数など)の解釈には注意が必要。
🚀 発展の可能性(結果 X → 新仮説 Y → 課題 Z)
この研究をさらに発展させるための3つの方向性を示します。「今回わかったこと(X)」から「次に検証すべき仮説(Y)」を立て、「具体的に何をするか(Z)」まで考えてみましょう。
① 観光資源の「質」を測る
- 結果 X
- 本論文は観光資源の充実度を量(観光地点数)でしか測っておらず、質的側面を捉えていない。
- 新仮説 Y
- 同じ「歴史・文化」でも、人気の高い資源(世界遺産級)と無名の資源では消費喚起力が違うはず。質を測れば負の効果の一部は説明できるかもしれない。
- 課題 Z
- (1)UNESCO世界遺産・世界ジオパーク、文化庁の日本遺産などのフラグを変数化する。(2)「質×量」の交互作用項を入れて再推定する。(3)質を統制しても「歴史・文化」の負の効果が残るか検証する。
② 分析単位を市町村へ細粒度化
- 結果 X
- 公的集計データが都道府県レベルに限られ、県内の地域差を捉えられていない。
- 新仮説 Y
- 市町村レベルなら、飲食・宿泊が集積した観光地とそうでない地点の違いをより直接に検証できる。
- 課題 Z
- (1)市町村レベルの観光入込・施設データを収集する。(2)観光客の個票データが得られれば、宿泊/日帰りの意思決定を個人単位でモデル化する。(3)都道府県レベルの結論が市町村でも成り立つか比較する。
③ サンプルセレクションへの対処
- 結果 X
- 本論文は欠測を含むアンバランスパネルで、欠損がランダムでない場合のサンプルセレクションバイアスを扱っていない。
- 新仮説 Y
- 統計を早く導入した県ほど観光に積極的、といった偏りがあれば係数にバイアスが生じている可能性がある。
- 課題 Z
- (1)欠測の発生要因を可視化する。(2)Heckmanの2段階推定などセレクション補正を試す。(3)補正あり/なしで結論が変わるかを比較する。
🎯 自分でやってみよう(5つのチャレンジ)
学んだだけでは身につきません。実際に手を動かすのが最強の学習方法です。本論文のスクリプトをベースに、以下のチャレンジに挑戦してみてください。難易度別に5つ用意しました。
★☆☆☆☆ 入門
CH1. 同じデータで分析を再現する
まずは付属の Python スクリプトをそのまま実行し、論文と同じ図を再現してみてください。
ポイント: 各図がどのコード行から生成されているか辿る。エラーが出たら原因を考える。
★★☆☆☆ 初級
CH2. 説明変数を1つ追加・除外して結果を比較
本論文の分析モデルから説明変数を1つ抜いて再実行、あるいは1つ追加して再実行してください。
ポイント: 係数・p値・R² がどう変わったか観察する。多重共線性が原因で結果が変わる例を見つけられたら理想的。
★★★☆☆ 中級
CH3. 別の年度・別の都道府県で同じ分析を試す
SSDSE の別の年度(例:2015年度・2020年度)または特定都道府県のみのデータで同じ分析を実行してください。
ポイント: 時代や地域によって結論が変わるか? 変わるならその理由を考察する。
★★★★☆ 上級
CH4. 別の手法を組み合わせる
本論文の手法 + 1つの追加手法(例:重回帰 + LASSO、相関分析 + 主成分分析)で結果を比較してください。
ポイント: 手法の違いで結論が変わるか? どちらが妥当かを「なぜ」とともに説明できるように。
★★★★★ 発展
CH5. オリジナルの問いを立てて分析する
本論文の手法を借りて、あなた自身の問いを立てて分析してください。
例:「カフェの数と幸福度に関連はあるか」「教育費の高い県は出生率も高いか」など。
ポイント: 問い・データ・手法・結論を1ページのレポートにまとめる。これがデータサイエンスの「実践」。
💡 ヒント: 詰まったら本サイトの他の論文(同じ手法を使っている)のスクリプトをコピーして組み合わせるのが効率的です。手法ガイド・用語集も参考に。
💼 この手法は実社会でこう使われている
本論文で学んだ手法は、研究の世界だけでなく、行政・企業・NPO の現場でも様々に活用されています。具体的なシーンを紹介します。
🏛️
行政の政策立案
都道府県・市区町村の政策担当者は、本論文と同様のデータ分析を用いて「どこに予算を投じれば効果が出るか」を検討します。
例えば医療費削減策、移住促進策、子育て支援策などの効果予測・効果検証に直結します。
🏢
企業のマーケティング・出店戦略
小売チェーン・サービス業の出店戦略では、地域特性(人口構成、所得、ライフスタイル)と売上の関係を本論文と同じ手法で分析します。
ECサイトでも顧客セグメント分析・購買要因分析に類似手法が使われます。
🏥
医療・公衆衛生
感染症の流行予測、医療資源配分の最適化、健康格差の地域要因分析などで、本論文の統計手法は標準的に使われています。
WHO・厚労省レベルの政策評価でも同じ手法が活躍しています。
📊
メディア・ジャーナリズム
新聞・テレビの社会調査記事、選挙予測、世論調査の分析でも、本論文と同じ手法(回帰分析・クラスタリングなど)が使われています。
データジャーナリズムの記事はこの種の分析が中核です。
🎓
学術研究(隣接分野)
経済学・社会学・公衆衛生学・教育学・地理学などの実証研究では、本論文と同じ手法が日常的に使われます。
専門誌に掲載される論文の8割以上が、こうした統計手法に基づいて結論を出しています。
💰
金融・保険業界
与信判断(融資審査)、保険料の地域別設定、不動産価格予測などで、本論文と同様のモデリング手法が広く活用されています。
統計分析の能力は金融業界の必須スキルになっています。
🤔 よくある質問(読者からの想定Q&A)
この論文を読んで初心者が抱きやすい疑問に、教育的観点から答えます。
Q1. この分析、自分でもできますか?
図の再現は誰でもできます。ただし注意が必要です。本研究の被説明変数「観光消費額」は『共通基準による観光入込客統計』由来で SSDSE には収録されていないため、本ページの図は回帰の再計算ではなく原論文の報告値の可視化です。付属スクリプトを実行すれば図と報告値の出力が再現できます。回帰そのものを実再現するには、観光消費額データを別途入手する必要があります。
Q2. 使われている手法は他の分野にも応用できますか?
十分応用可能です。本論文のパネル固定効果/変量効果モデルとHausman検定は、医療・教育・経済・環境など「同じ対象を複数年追跡する」あらゆる分野で標準的に使われます。データの中身(変数)を入れ替えるだけで、別の問いにも適用できます。
Q3. 結論は本当に「因果関係」を示していますか?
本論文は「観察データ」を使った分析であり、厳密な意味での「因果関係」を完全に証明したわけではありません。あくまで「強い関連が見られた」という事実を提示しているにとどまります。真の因果を示すには、無作為化比較試験(RCT)か、自然実験を活用したIV・DiD 等の手法が必要です。
Q4. データの最新版を使うとどうなりますか?
SSDSE は毎年更新されているため、最新版を使えば近年のトレンド(特にコロナ禍以降の変化)も含めて分析できます。ただし、結論が変わる可能性もあります。それ自体が新しい発見につながります。
Q5. もっと深く学ぶには何を読めばいいですか?
「計量経済学」「データサイエンス入門」「統計的因果推論」などのテキストが入門に向いています。Python の場合は『Python ではじめる機械学習』(オライリー)、R の場合は『R で学ぶ統計学』が定番です。本サイトの他の論文も読み比べてみてください。
✅ 理解度チェック(4問)
この論文を理解できたか、クリックで確認しましょう。間違えても解説が出ます。