🔖 キーワード索引 — 信頼区間を多角的に理解する

信頼区間（confidence interval）は推測統計の最重要概念です。 関連キーワードを難易度別に整理しました。

🟢 基礎キーワード（まず押さえる）

• 信頼水準（confidence level）：通常 95%、 厳しい場合 99%。 (1 − α) で表現。
• 区間推定（interval estimation）：点推定と対をなす推定法。 値ではなく区間を出す。
• 点推定（point estimation）：標本から 1 つの値を出す推定。 平均、 比率、 回帰係数など。
• 標準誤差（standard error, SE）：推定量のばらつき。 SE = σ/√n。 信頼区間の幅を決める。
• マージン・オブ・エラー（margin of error）：信頼区間の半幅。 z·SE または t·SE。
• 母集団と標本：信頼区間は母数を「区間で挟む」推定。 母集団と標本の区別が前提。

🟡 中級キーワード

• t 分布の信頼区間：母分散未知・小標本の場合、 正規ではなく t 分布の臨界値を使う。
• z 区間と t 区間：母分散既知なら z、 未知なら t（df = n−1）。 N 大なら両者ほぼ同じ。
• 比率の信頼区間：Wald、 Wilson、 Clopper-Pearson、 Agresti-Coull など複数の方法。
• 差の信頼区間：2 群平均差、 比率差、 オッズ比、 リスク比などの区間。
• ブートストラップ信頼区間：再標本化による経験的区間。 分布の仮定なし。
• カバレッジ確率（coverage probability）：信頼区間が真値を含む実際の確率。 理論値（95% など）に近いか。

🔴 上級キーワード

• BCa ブートストラップ：バイアス補正 + 加速。 通常のパーセンタイル法より精度高。
• プロファイル尤度信頼区間：尤度比検定に基づく区間。 最尤推定で頻用。
• 事後信用区間（Bayesian credible interval）：ベイズ統計での区間。 信頼区間と意味が異なる。
• 同時信頼区間（simultaneous CI）：複数パラメータを同時にカバー。 多重比較で重要。
• 非劣性・優越性試験の信頼区間：医薬品開発で重要。 区間が閾値を超えるか否か。
• カバレッジ条件付き信頼区間：選択バイアス補正。 post-selection inference の枠組み。

🎨 直感で掴む

95%CI の例。 n が大きいほど CI は狭く、 推定の確度が高い。 CIが 0 を含むか否かが有意性判定の鍵。

🎨 概念図で押さえる

信頼区間 (CI) は 「同じ手続きを繰り返したとき、 真値を含む区間が 95% できる」幅。 以下 3 図で「点推定 vs 区間推定」「N と CI 幅」「95% の意味」を可視化する。

図 1: 点推定 vs 区間推定の違い

→ 点推定は1 点だけを答えるが、 区間推定は不確かさ込みで答える。 同じデータでも報告の重みが違う。

図 2: 標本サイズ N と CI 幅の関係

→ N が 4 倍 → 幅は √4=2 倍狭くなる (σ/√N の効果)。 N=1000 なら ±0.76 に収まり、 結論に強い根拠を与えられる。

図 3: 「95% 信頼区間」の正しい意味

→ 100 回標本抽出して 100 個の CI を作ると、 約 95 個が真値 μ を含む。 「あなたの 1 つの区間が 95% の確率で真値を含む」ではない (頻度主義の解釈)。

🧩 信頼区間を深める (理解度チェック + 数式読み解き)

🔬 数式を言葉で読み解く (95% 信頼区間の構造)

母平均 $\mu$ の 95% 信頼区間は次式で求めます (正規分布近似、 既知の母分散 $\sigma$ の場合)。

$$ \text{CI}_{0.95} = \bar{x} \pm 1.96 \cdot \frac{\sigma}{\sqrt{n}} $$

記号の意味: $\bar{x}$ は標本平均 (中心)、 $\sigma$ は母標準偏差、 $n$ はサンプルサイズ、 $1.96$ は標準正規分布の 97.5% 点 (両側 95% に対応)。 分子 $\sigma$ は「データのばらつき」、 分母 $\sqrt{n}$ は「標本数の効果」を表します。 $n$ が 4 倍になると区間幅は $\sqrt{4} = 2$ 倍狭くなり、 精度が上がります。 $\sigma$ が未知の場合は標準偏差 $s$ と $t$ 分布を用い、 $1.96$ を $t_{n-1, 0.975}$ に置き換えます。

✅ 理解度チェック (1 分で答える)

1. 「95% 信頼区間 = この区間に真値が 95% の確率で入る」は正しい解釈ですか? 違う場合は正しい言い方を述べてください。
2. サンプルサイズを 100 → 400 に増やすと、 信頼区間の幅はおよそ何倍になりますか?
3. 99% 信頼区間は 95% 信頼区間より広いですか、 狭いですか? その理由は?

このコードでやること: SSDSE-B-2026（2023年度）の47都道府県の消費支出（二人以上の世帯、 月平均、 千円）について、 95% 信頼区間を scipy.stats.t.interval で計算する。

📥 入力データ (SSDSE-B-2026 2023年度抜粋、 消費支出 L3221 (千円)):

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

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

df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1])
x = df.loc[df['SSDSE-B-2026'] == 2023, 'L3221'].values / 1000  # 消費支出(千円)
n = len(x)
mean = np.mean(x)
sem = stats.sem(x)
ci = stats.t.interval(0.95, n-1, loc=mean, scale=sem)
print(f'n = {n}, 標本平均 = {mean:.2f}')
print(f'95% 信頼区間: [{ci[0]:.2f}, {ci[1]:.2f}]')

📤 実行例:

💬 47 都道府県の消費支出の標本平均は 295.86 千円。 95% 信頼区間は [288.77, 302.94]。 「同じ手続きを反復して区間を作ると、 その 95% が母平均を含む」 と読む。

📌 関連用語クイック参照

• 標準誤差 — 信頼区間の幅を決める核
• 仮説検定 — 信頼区間と表裏一体
• t分布 — 母分散未知のときの基礎分布
• 点推定 / 区間推定 — 信頼区間は区間推定の代表例
• ブートストラップ — 分布仮定なしの信頼区間

🎮 触って理解する — 「95%」の意味を体感するシミュレータ

信頼区間の最大の誤解は 「この 1 本の区間に真の値が 95% の確率で入る」 という読み方です。 正しくは 「同じ手続きで標本抽出と区間作成を何度も繰り返すと、 作られた区間の約 95% が真の母平均を含む」。 下のシミュレータでこの頻度論的な被覆を目で確かめてください。

母集団は SSDSE-B-2026・2023 年度の 47 都道府県の消費支出（L3221、 二人以上の世帯、 千円/月）。 この 47 個を「既知の母集団」とみなし、 母平均 μ = 295.86 千円（紫の縦線）を真の値とします。 ここから標本サイズ n の標本を復元抽出し、 各標本の信頼区間を横棒で下へ積み上げます。 緑 = μ を含む、 赤 = μ を外す。 抽出を繰り返すほど、 μ を含む割合が信頼水準に近づくことを観察できます。

※ 乱数を使うため実行のたびに結果は変わります。 t 分布の臨界値は正則化不完全ベータ関数から算出（scipy.stats.t.ppf と小数第 4 位まで一致）。 標本数を数十〜数百と増やすほど、 被覆率は選んだ信頼水準（90 / 95 / 99%）に収束します。

💡 この実験で腑に落ちること

• 「95%」は区間ではなく手続きの性質：赤い棒（μ を外した区間）も「間違った計算」ではありません。 正しい手続きでも 約 5% は外れるのが 95% 信頼区間の定義です。 「いま手元にあるこの 1 本が 95% の確率で μ を含む」ではなく、 「この作り方で区間を量産すると 95% が当たる」── 確率は区間の作り方に属し、 固定された真値 μ や個々の区間に属しません（頻度論的解釈）。
• 信頼水準 ↑ → 区間幅 ↑：90 → 95 → 99% とセグメントを上げると、 臨界値 $t_{n-1,\,\alpha/2}$ が大きくなり棒が横に伸びます。 外す本数（赤）は減るが、 その代償に区間は広く（＝情報は曖昧に）なる。 「確実さ」と「精密さ」はトレードオフです。
• 標本サイズ n ↑ → 区間幅 ↓：n を大きくすると全部の棒が一斉に短くなります。 幅は標準誤差 $\text{SE}=s/\sqrt{n}$ に比例するので、 n を 4 倍にすると幅は約 1/2。 精度は $\sqrt{n}$ でしか上がらない（費用対効果の逓減）。 → 標準誤差／標本分布と中心極限定理。
• 被覆率は水準に収束する：数本だけでは被覆率は 100% や 80% に振れますが、 数百本まとめると 90/95/99% の近傍で安定します。 これは 大数の法則の現れ。 「1 回の実験では当たり外れがある／長期頻度としては保証される」という頻度論の二面性そのものです。

🔗 仮説検定・p 値との双対性、 ブートストラップとの関係

• 信頼区間 ⇔ 仮説検定は表裏一体：95% 信頼区間が値 $\theta_0$ を含まない ことと、 帰無仮説 $H_0:\theta=\theta_0$ が有意水準 5% で棄却される ことは同値です。 とくに 区間が 0 を含むか否かは「効果あり／なし」（両側 $p<0.05$）と一致。 → 仮説検定／p 値。 CI は「棄却か否か」の 2 値に加えて効果の大きさと精度まで示すため、 p 値より情報量が多いとされます。
• 分布の仮定が怪しいときはブートストラップ CI：t 区間は「標本平均がほぼ正規」という近似に依存します。 母集団が強く歪む・統計量が中央値や相関係数などで解析式が面倒、 といった場合は、 手元の標本を復元抽出して経験的に区間を作る ブートストラップ信頼区間（パーセンタイル法・BCa 法）が有効。 このページの Python 例（scipy.stats.bootstrap）では消費支出のブートストラップ 95% CI ≈ [288.69, 302.46] と、 t 区間 [288.77, 302.94] がほぼ一致します。

※ ベイズ統計の「信用区間（credible interval）」は「パラメータがこの区間に入る事後確率が 95%」と確率を直接パラメータに与える点で、 ここで扱う頻度論的信頼区間とは意味が異なります（数値は近くても解釈が別物）。

📐 数式

📐 信頼区間の数式

① 平均の信頼区間（σ既知）

$$ \bar{x} \pm z_{\alpha/2} \cdot \frac{\sigma}{\sqrt{n}} $$

② 平均の信頼区間（σ未知、 t分布）

$$ \bar{x} \pm t_{\alpha/2, n-1} \cdot \frac{s}{\sqrt{n}} $$

③ 比率の信頼区間（Wald）

$$ \hat{p} \pm z_{\alpha/2} \sqrt{\frac{\hat{p}(1-\hat{p})}{n}} $$

④ 差の信頼区間（2標本）

$$ (\bar{x}_1 - \bar{x}_2) \pm t_{\alpha/2} \cdot SE $$

典型的な信頼水準と臨界値

🔬 数式を「言葉」で読み解く

🧮 実データでの計算 — SSDSE 食料費

47都道府県の家計食料費（2023）：n = 47、 平均 = 80.598千円、 標準偏差 s = 5.842千円。

1. 標準誤差: SE = s/√n = 5.842/√47 = 0.8521
2. 自由度: df = n-1 = 46
3. t臨界値（α=0.05）: t = 2.013
4. 95%CI: 80.598 ± 2.013 × 0.8521 = [78.883, 82.313]

🧮 SSDSE-B-2026 実値計算例 — 47 都道府県データで信頼区間

合成データではなく公的統計を念頭に、 信頼区間の具体的計算手順を数値で示します。

① 平均寿命の 95% 信頼区間（t 区間）

🎯 このコードでやること：信頼区間 — 母平均が含まれると推定される区間に関連するステップ #4。主要な指標（係数・統計量・スコア）を算出します。

📥 入力例（df.head()）
df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1])
df2023 = df[df['SSDSE-B-2026'] == 2023]
# 2023年度47都道府県の消費支出 L3221（二人以上の世帯、 円/月）を母集団とみなし、 平均の 95% CI を計算
# df2023[['Prefecture','L3221']].head():
#    Prefecture   L3221
# 0  北海道      296888
# 12 青森県      263371
# 24 岩手県      298536
# 36 宮城県      305541
# 48 秋田県      272086

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

# 47 都道府県の平均寿命（仮想）
標本平均 x̄ = 84.55 年
標本標準偏差 s = 0.62 年
n = 47

# 標準誤差
SE = s / √n = 0.62 / √47 ≈ 0.0904

# t 分布の臨界値（df = 46、 α = 0.05）
t(0.025, 46) ≈ 2.013

# 95% 信頼区間
CI_95 = x̄ ± t · SE
      = 84.55 ± 2.013 · 0.0904
      = 84.55 ± 0.182
      = [84.37, 84.73]

→ 「全 47 都道府県の真の平均寿命は 95% 信頼で 84.37〜84.73 年」
→ 幅 0.36 年

📤 このブロックは標準出力には何も出さない

② 99% 信頼区間との比較

🎯 このコードでやること：信頼区間 — 母平均が含まれると推定される区間に関連するステップ #5。仮説検定・モデル評価を行います。

📥 入力例（df.head()）
df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1])
df2023 = df[df['SSDSE-B-2026'] == 2023]
# 2023年度47都道府県の消費支出 L3221（二人以上の世帯、 円/月）を母集団とみなし、 平均の 95% CI を計算
# df2023[['Prefecture','L3221']].head():
#    Prefecture   L3221
# 0  北海道      296888
# 12 青森県      263371
# 24 岩手県      298536
# 36 宮城県      305541
# 48 秋田県      272086

1
2
3
4
5
6
7
8
9

# 99% 信頼水準
t(0.005, 46) ≈ 2.687

CI_99 = 84.55 ± 2.687 · 0.0904
      = 84.55 ± 0.243
      = [84.31, 84.79]

→ 幅 0.49 年（95% より広い）
→ 信頼水準を上げると区間は広がる

📤 このブロックは標準出力には何も出さない

③ 平均差の信頼区間（東日本 vs 西日本）

🎯 このコードでやること：信頼区間 — 母平均が含まれると推定される区間に関連するステップ #6。結果を整形して表示します。

📥 入力例（df.head()）
df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1])
df2023 = df[df['SSDSE-B-2026'] == 2023]
# 2023年度47都道府県の消費支出 L3221（二人以上の世帯、 円/月）を母集団とみなし、 平均の 95% CI を計算
# df2023[['Prefecture','L3221']].head():
#    Prefecture   L3221
# 0  北海道      296888
# 12 青森県      263371
# 24 岩手県      298536
# 36 宮城県      305541
# 48 秋田県      272086

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

# 仮想値
東日本（n=24）：平均 84.3、 標準偏差 0.7
西日本（n=23）：平均 84.8、 標準偏差 0.6

# 平均差の SE（Welch）
SE_diff = √(0.7²/24 + 0.6²/23) = √(0.0204 + 0.0157) ≈ 0.190

# 自由度（Welch–Satterthwaite）
df ≈ 44.5

# 95% 信頼区間（平均差）
CI = (84.3 - 84.8) ± 2.014 · 0.190
   = -0.50 ± 0.383
   = [-0.883, -0.117]

→ 0 を含まないので有意（p < 0.05）
→ 「西日本は東日本より 0.12〜0.88 年長い（95% 信頼）」

📤 このブロックは標準出力には何も出さない

④ 比率の信頼区間（人口減少県の割合）

🎯 このコードでやること：信頼区間 — 母平均が含まれると推定される区間に関連するステップ #7。47都道府県データに当てはめて確認します。

📥 入力例（df.head()）
df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1])
df2023 = df[df['SSDSE-B-2026'] == 2023]
# 2023年度47都道府県の消費支出 L3221（二人以上の世帯、 円/月）を母集団とみなし、 平均の 95% CI を計算
# df2023[['Prefecture','L3221']].head():
#    Prefecture   L3221
# 0  北海道      296888
# 12 青森県      263371
# 24 岩手県      298536
# 36 宮城県      305541
# 48 秋田県      272086

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

# 47 都道府県中 38 県が人口減少（仮想）
p̂ = 38/47 ≈ 0.809
n = 47

# Wald 法
SE_p = √(p̂(1-p̂)/n) = √(0.809·0.191/47) = √0.00329 ≈ 0.0573
CI_Wald = 0.809 ± 1.96 · 0.0573 = [0.696, 0.921]

# Wilson 法（より精度高、 小標本で推奨）
z² = 1.96² = 3.8416
center = (p̂ + z²/(2n)) / (1 + z²/n) = (0.809 + 0.0409) / (1 + 0.0817) ≈ 0.785
margin = (z/(1 + z²/n))·√(p̂(1-p̂)/n + z²/(4n²))
       ≈ 0.1106

CI_Wilson = 0.785 ± 0.1106 ≈ [0.675, 0.896]

# 比較
Wald   :   [0.696, 0.921]（端で不安定）
Wilson :   [0.675, 0.896]（より現実的）

📤 このブロックは標準出力には何も出さない

⑤ 回帰係数の信頼区間

🎯 このコードでやること：信頼区間 — 母平均が含まれると推定される区間に関連するステップ #8。比較・別パターンを検討します。

📥 入力例（df.head()）
df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1])
df2023 = df[df['SSDSE-B-2026'] == 2023]
# 2023年度47都道府県の消費支出 L3221（二人以上の世帯、 円/月）を母集団とみなし、 平均の 95% CI を計算
# df2023[['Prefecture','L3221']].head():
#    Prefecture   L3221
# 0  北海道      296888
# 12 青森県      263371
# 24 岩手県      298536
# 36 宮城県      305541
# 48 秋田県      272086

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11

# 粗死亡率(‰) ~ 高齢化率(%) の単回帰（2023年度、 47 都道府県）
β̂ = 0.610（推定値）
SE(β̂) = 0.022

# 95% 信頼区間
CI = 0.610 ± 2.013 · 0.022
   = 0.610 ± 0.044
   = [0.566, 0.654]

→ 0 を含まないので係数は有意
→ 「高齢化率が 1% 上がると粗死亡率(人口千対)は 0.57〜0.65 ポイント上昇」

📤 このブロックは標準出力には何も出さない

🧮 数式に値を入れて手で計算する — t 分布による 95% 信頼区間

合成データ x = [3, 5, 8, 2, 7]（n=5）を代入し、 平均 x̄ ± t_{n-1, α/2} × s/√n を Step 1〜5 で展開する。 標本平均・標本標準偏差・標準誤差・自由度 4 の t 臨界値・上下限の順に手計算し、 同じ計算を Python（scipy.stats）で再現する。

📐 用語の主要数式（再掲）

$$ \bar{x} \pm t_{n-1,\,\alpha/2} \cdot \frac{s}{\sqrt{n}}, \qquad s = \sqrt{\frac{1}{n-1}\sum_{i=1}^n (x_i - \bar{x})^2} $$

母分散未知・標本サイズ小（n<30）の場合、 標準正規分布ではなく自由度 n-1 の t 分布の臨界値を使う。 α=0.05 で両側 95% 信頼区間を構成する。

Step 1: データ準備（n=5、 合成整数）

Step 2: 標本平均 x̄

Step 3: 標本標準偏差 s（不偏、 分母 n-1）

Step 4: 標準誤差 SE と t 臨界値

Step 5: 95% 信頼区間（最終結果）

n=5 の小サンプルなので、 母平均の不確かさが大きく、 区間幅は約 6.33 と広い。 n が増えれば SE = s/√n が縮み、 区間も狭くなる。

🐍 同じ計算を Python で再現

 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
27
28
29

import numpy as np
from scipy import stats

# Step 1: データ
x = np.array([3, 5, 8, 2, 7])
n = len(x)

# Step 2: 標本平均
mean = x.mean()

# Step 3: 標本標準偏差（ddof=1 で不偏）
s = x.std(ddof=1)

# Step 4: SE と t 臨界値
se = s / np.sqrt(n)
df = n - 1
t_crit = stats.t.ppf(0.975, df)  # 両側 95% → 上側 2.5%
margin = t_crit * se

# Step 5: 95% 信頼区間
lo = mean - margin
hi = mean + margin

print(f"x̄ = {mean:.4f}")
print(f"s  = {s:.4f}")
print(f"SE = {se:.4f}")
print(f"t_(df={df}, 0.025) = {t_crit:.4f}")
print(f"誤差幅 = {margin:.4f}")
print(f"95% CI = [{lo:.4f}, {hi:.4f}]")

📤 実行結果

💬 手計算（Step 2: x̄=5.0、 Step 3: s=2.5495、 Step 4: SE=1.1402, t=2.7764, 誤差幅=3.1656、 Step 5: CI=[1.8344, 8.1656]）と scipy.stats.t.ppf を使った Python 出力が完全一致。 自由度 4 の t 臨界値 2.7764 が、 標準正規の 1.96 より大きいことが、 小サンプルでの保守的な区間構成の根拠となる。

🐍 Python での信頼区間

🎯 このコードでやること：信頼区間 — 母平均が含まれると推定される区間に関連するステップ #1。最初のスニペットです。SSDSE-B-2026 を読み込みます。

📥 入力例（df.head()）
df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1])
df2023 = df[df['SSDSE-B-2026'] == 2023]
# 2023年度47都道府県の消費支出 L3221（二人以上の世帯、 円/月）を母集団とみなし、 平均の 95% CI を計算
# df2023[['Prefecture','L3221']].head():
#    Prefecture   L3221
# 0  北海道      296888
# 12 青森県      263371
# 24 岩手県      298536
# 36 宮城県      305541
# 48 秋田県      272086

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

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

df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1])
data = df.loc[df['SSDSE-B-2026'] == 2023, 'L3221'].values / 1000  # 消費支出(千円)

# 平均の95%信頼区間
def mean_ci(data, alpha=0.05):
    mean = np.mean(data)
    se = stats.sem(data)  # 標準誤差
    t_crit = stats.t.ppf(1 - alpha/2, df=len(data)-1)
    return mean - t_crit * se, mean + t_crit * se

ci_low, ci_high = mean_ci(data)

# scipy で一発
ci = stats.t.interval(0.95, df=len(data)-1, loc=np.mean(data), scale=stats.sem(data))

# ブートストラップ信頼区間
from scipy.stats import bootstrap
res = bootstrap((data,), np.mean, confidence_level=0.95, random_state=0)
print(f'ブートストラップ 95% CI: [{res.confidence_interval.low:.2f}, {res.confidence_interval.high:.2f}]')

📤 実行例（実行時の標準出力）
ブートストラップ 95% CI: [288.69, 302.46]

比率の信頼区間（Wilson、 Clopper-Pearson）

🎯 このコードでやること：信頼区間 — 母平均が含まれると推定される区間に関連するステップ #2。基本統計量を計算します。

📥 入力例（df.head()）
df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1])
df2023 = df[df['SSDSE-B-2026'] == 2023]
# 2023年度47都道府県の消費支出 L3221（二人以上の世帯、 円/月）を母集団とみなし、 平均の 95% CI を計算
# df2023[['Prefecture','L3221']].head():
#    Prefecture   L3221
# 0  北海道      296888
# 12 青森県      263371
# 24 岩手県      298536
# 36 宮城県      305541
# 48 秋田県      272086

1
2
3
4
5

from statsmodels.stats.proportion import proportion_confint

# 成功 30 / 試行 100
ci_low, ci_high = proportion_confint(count=30, nobs=100, alpha=0.05, method='wilson')
print(f'Wilson 95% CI: [{ci_low:.4f}, {ci_high:.4f}]')

📤 実行例（実行時の標準出力）
Wilson 95% CI: [0.2189, 0.3958]

🎲 ブートストラップ信頼区間

分布の仮定なしに信頼区間を計算する手法。 元データから復元抽出を繰り返し、 統計量の分布を直接シミュレーション：

1. 元データから n 個を復元抽出 → ブートストラップ標本
2. その平均（または他の統計量）を計算
3. 1〜2を 1000〜10000 回繰り返す → 統計量の分布
4. 2.5%〜97.5%パーセンタイルが 95% CI

正規分布の仮定が崩れる場合、 メディアンや分位数の CI を計算する場合に強力。 現代の標準ツール。

🧠 信頼区間の正しい解釈

📏 CIの幅は何で決まる？

CIの幅 = 2 × z × σ / √n。 これを左右する3要素：

• 信頼水準：95% → 1.96、 99% → 2.58。 高くするほど広く
• 標本サイズ n：n が4倍で幅が半分。 大きいほど狭く
• データのばらつき σ：σ が大きいほど広く

「CIが広い」→ 推定の精度が低い、 サンプル不足の可能性。 「CIが狭い」→ 精度が高い推定。

🎲 信用区間（credible interval）— ベイズ版

ベイズ統計では「信用区間（credible interval）」が信頼区間に対応。 解釈がより直感的：

💡 ベイズの95%信用区間：「真値がこの範囲に入る確率が95%」。 頻度主義の信頼区間と違って、 直接的な確率表現が可能。

事前分布 + 尤度 → 事後分布 → 事後分布の2.5%〜97.5%が95%信用区間。

🔗 「信頼区間」と他の概念の関係を深掘り

① CI幅と仮説検定の関係

95%CIが特定値（例えば0）を含まない ⇔ 両側5%水準でその値が棄却される。 つまりCIと検定は表裏一体。

• 95%CIが [2.5, 5.0] → 0 を含まない → p < 0.05 で「ゼロと異なる」
• 95%CIが [-1.0, 3.0] → 0 を含む → p ≥ 0.05 で「ゼロと異なるとは言えない」

② 推定の精度評価

CIの幅は推定精度の指標。 半分の幅にしたければ、 サンプルサイズを4倍にする（√n則）。

③ 効果量との併記

近代的な統計報告では点推定 + 95%CI + 効果量を併記するのが標準（APA 7版マニュアル、 ASA声明）。

④ 種類別CIの公式まとめ

⑤ 落とし穴 — CIの誤用

• CIを「中の値はどれも同じ確率」と解釈しがちだが、 中心の方がより確からしい
• 2つのCIが少しでも重なる ≠ 差がない。 厳密には差のCIで判断
• n が大きいと CI は狭くなるが、 小さな差まで「有意」になりやすい
• CIは正規分布近似に依存することが多い → 小標本ではブートストラップを

✅ 実務チェックリスト — 推測統計を使う前に

1. データの確認

• サンプルサイズ n は十分か？（n < 30 なら t分布、 非正規性に注意）
• 独立同分布の仮定は妥当か？（時系列、 クラスター構造に注意）
• 外れ値の影響を確認したか？（box plot で）
• 正規性は確認したか？（QQプロット、 Shapiro-Wilk検定）
• 欠損値の扱いは適切か？

2. 検定・推定の設計

• 仮説（H₀、 H₁）は事前に定義したか？
• 片側 vs 両側を選択しているか？
• 有意水準 α は事前に設定したか？
• 多重比較の補正は必要か？
• サンプルサイズの事前計算（power analysis）したか？

3. 結果の報告

• 点推定 + 信頼区間を併記しているか？
• p値だけでなく効果量も報告したか？
• サンプルサイズを明記したか？
• 仮定の確認結果を述べたか？
• 「統計的有意 = 実質的に重要」と混同していないか？

4. 解釈の注意

• 「相関 ≠ 因果」を意識
• シンプソンのパラドックスを警戒
• 探索的分析と確認的分析を区別
• 結果を再現できるか

🐍 Python 実装のバリエーション — scipy / statsmodels / pingouin

① scipy.stats による平均の信頼区間

🎯 このコードでやること：信頼区間 — 母平均が含まれると推定される区間に関連するステップ #9。ハイパーパラメータを変えて再計算します。

📥 入力例（df.head()）
df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1])
df2023 = df[df['SSDSE-B-2026'] == 2023]
# 2023年度47都道府県の消費支出 L3221（二人以上の世帯、 円/月）を母集団とみなし、 平均の 95% CI を計算
# df2023[['Prefecture','L3221']].head():
#    Prefecture   L3221
# 0  北海道      296888
# 12 青森県      263371
# 24 岩手県      298536
# 36 宮城県      305541
# 48 秋田県      272086

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

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

df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1])
x = df.loc[df['SSDSE-B-2026'] == 2023, 'L3221'].dropna().values / 1000  # 消費支出(千円)

# t 区間
mean = x.mean()
sem = stats.sem(x)  # 標準誤差
ci_low, ci_high = stats.t.interval(0.95, df=len(x)-1, loc=mean, scale=sem)
print(f'平均 = {mean:.2f}, 95% CI = [{ci_low:.2f}, {ci_high:.2f}]')

# z 区間（母分散既知の場合）
z_low, z_high = stats.norm.interval(0.95, loc=mean, scale=sem)
print(f'z CI = [{z_low:.2f}, {z_high:.2f}]')

📤 実行例（実行時の標準出力）
平均 = 295.86, 95% CI = [288.77, 302.94]
z CI = [288.95, 302.76]

② 比率の信頼区間（statsmodels）

🎯 このコードでやること：信頼区間 — 母平均が含まれると推定される区間に関連するステップ #10。最終結果のまとめ・保存を行います。

📥 入力例（df.head()）
df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1])
df2023 = df[df['SSDSE-B-2026'] == 2023]
# 2023年度47都道府県の消費支出 L3221（二人以上の世帯、 円/月）を母集団とみなし、 平均の 95% CI を計算
# df2023[['Prefecture','L3221']].head():
#    Prefecture   L3221
# 0  北海道      296888
# 12 青森県      263371
# 24 岩手県      298536
# 36 宮城県      305541
# 48 秋田県      272086

1
2
3
4
5
6
7
8
9

from statsmodels.stats.proportion import proportion_confint

# 2023年度、 47 都道府県中 45 県が前年（2022年度）より総人口減少
count = 45
n = 47

for method in ['normal', 'wilson', 'beta', 'agresti_coull']:
    ci = proportion_confint(count, n, alpha=0.05, method=method)
    print(f'{method:15s}: [{ci[0]:.4f}, {ci[1]:.4f}]')

📤 実行例（実行時の標準出力）
normal         : [0.8997, 1.0000]
wilson         : [0.8575, 0.9883]
beta           : [0.8546, 0.9948]
agresti_coull  : [0.8496, 0.9962]

③ 平均差の信頼区間（Welch）

🎯 このコードでやること：信頼区間 — 母平均が含まれると推定される区間に関連するステップ #11。可視化を仕上げ、レポートに統合します。

📥 入力例（df.head()）
df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1])
df2023 = df[df['SSDSE-B-2026'] == 2023]
# 2023年度47都道府県の消費支出 L3221（二人以上の世帯、 円/月）を母集団とみなし、 平均の 95% CI を計算
# df2023[['Prefecture','L3221']].head():
#    Prefecture   L3221
# 0  北海道      296888
# 12 青森県      263371
# 24 岩手県      298536
# 36 宮城県      305541
# 48 秋田県      272086

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

import pandas as pd

# df2023 はこのあとのブロックで作っているので、ここでも用意しておく
df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1])
df2023 = df[df['SSDSE-B-2026'] == 2023]

from statsmodels.stats.weightstats import CompareMeans, DescrStatsW

# 便宜上、 地域コード R24000（三重県）以東を東日本とする（東 24 県・西 23 県）
east = df2023.loc[df2023['Code'] <= 'R24000', 'L3221'] / 1000
west = df2023.loc[df2023['Code'] >  'R24000', 'L3221'] / 1000

cm = CompareMeans(DescrStatsW(east), DescrStatsW(west))
ci = cm.tconfint_diff(usevar='unequal', alpha=0.05)
print(f'平均差の 95% CI = [{ci[0]:.3f}, {ci[1]:.3f}]')

📤 実行例（実行時の標準出力）
平均差の 95% CI = [9.337, 34.867]

④ ブートストラップ信頼区間

🎯 このコードでやること：信頼区間 — 母平均が含まれると推定される区間に関連するステップ #12。追加検証・感度分析を実行します。

📥 入力例（df.head()）
df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1])
df2023 = df[df['SSDSE-B-2026'] == 2023]
# 2023年度47都道府県の消費支出 L3221（二人以上の世帯、 円/月）を母集団とみなし、 平均の 95% CI を計算
# df2023[['Prefecture','L3221']].head():
#    Prefecture   L3221
# 0  北海道      296888
# 12 青森県      263371
# 24 岩手県      298536
# 36 宮城県      305541
# 48 秋田県      272086

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11

from scipy.stats import bootstrap

# 平均値のブートストラップ CI（パーセンタイル法）
res = bootstrap((x,), np.mean, confidence_level=0.95, n_resamples=10000,
                method='percentile', random_state=0)
print(f'パーセンタイル: [{res.confidence_interval.low:.3f}, {res.confidence_interval.high:.3f}]')

# BCa 法（バイアス補正 + 加速）
res_bca = bootstrap((x,), np.mean, confidence_level=0.95, n_resamples=10000,
                    method='BCa', random_state=0)
print(f'BCa: [{res_bca.confidence_interval.low:.3f}, {res_bca.confidence_interval.high:.3f}]')

📤 実行例（実行時の標準出力）
パーセンタイル: [288.873, 302.577]
BCa: [288.695, 302.457]

⑤ 回帰係数の信頼区間（statsmodels）

🎯 このコードでやること：信頼区間 — 母平均が含まれると推定される区間に関連するステップ #13。応用パターン（別データ・別手法）に拡張します。

📥 入力例（df.head()）
df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1])
df2023 = df[df['SSDSE-B-2026'] == 2023]
# 2023年度47都道府県の消費支出 L3221（二人以上の世帯、 円/月）を母集団とみなし、 平均の 95% CI を計算
# df2023[['Prefecture','L3221']].head():
#    Prefecture   L3221
# 0  北海道      296888
# 12 青森県      263371
# 24 岩手県      298536
# 36 宮城県      305541
# 48 秋田県      272086

 1
 2
 3
 4
 5
 6
 7
 8
 9
10

import statsmodels.api as sm

d = df2023.assign(高齢化率 = df2023['A1303'] / df2023['A1101'] * 100,
                  粗死亡率 = df2023['A4200'] / df2023['A1101'] * 1000)
X = sm.add_constant(d[['高齢化率']])
y = d['粗死亡率']

model = sm.OLS(y, X).fit()
print('係数 CI:')
print(model.conf_int(alpha=0.05))

📤 実行例（実行時の標準出力）
係数 CI:
              0         1
const -6.573784 -3.762576
高齢化率   0.565729  0.654249

⑥ 中央値の信頼区間（順序統計量）

🎯 このコードでやること：信頼区間 — 母平均が含まれると推定される区間に関連するステップ #14。実務での落とし穴対策を実装します。

📥 入力例（df.head()）
df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1])
df2023 = df[df['SSDSE-B-2026'] == 2023]
# 2023年度47都道府県の消費支出 L3221（二人以上の世帯、 円/月）を母集団とみなし、 平均の 95% CI を計算
# df2023[['Prefecture','L3221']].head():
#    Prefecture   L3221
# 0  北海道      296888
# 12 青森県      263371
# 24 岩手県      298536
# 36 宮城県      305541
# 48 秋田県      272086

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

import numpy as np
from scipy import stats

x_sorted = np.sort(x)
n = len(x_sorted)

# 中央値の漸近 95% CI（n 大の場合）
median = np.median(x_sorted)
ci_idx_low = int(np.floor(n/2 - 1.96*np.sqrt(n)/2))
ci_idx_high = int(np.ceil(n/2 + 1.96*np.sqrt(n)/2))
print(f'中央値 = {median:.2f}')
print(f'95% CI ≈ [{x_sorted[ci_idx_low]:.2f}, {x_sorted[ci_idx_high]:.2f}]')

# ブートストラップでの中央値 CI
res = stats.bootstrap((x,), np.median, confidence_level=0.95,
                      n_resamples=10000, method='percentile', random_state=0)
print(f'ブートストラップ中央値 CI: [{res.confidence_interval.low:.2f}, {res.confidence_interval.high:.2f}]')

📤 実行例（実行時の標準出力）
中央値 = 300.65
95% CI ≈ [294.54, 306.54]
ブートストラップ中央値 CI: [294.84, 305.54]

⑦ オッズ比の信頼区間（ロジスティック回帰）

🎯 このコードでやること：信頼区間 — 母平均が含まれると推定される区間に関連するステップ #15。出力結果を解釈・考察します。

📥 入力例（df.head()）
df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1])
df2023 = df[df['SSDSE-B-2026'] == 2023]
# 2023年度47都道府県の消費支出 L3221（二人以上の世帯、 円/月）を母集団とみなし、 平均の 95% CI を計算
# df2023[['Prefecture','L3221']].head():
#    Prefecture   L3221
# 0  北海道      296888
# 12 青森県      263371
# 24 岩手県      298536
# 36 宮城県      305541
# 48 秋田県      272086

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

import statsmodels.api as sm

# 転出超過（転出者数 A5102 > 転入者数 A5101）か否かを 2 値の目的変数に
X = sm.add_constant(pd.DataFrame({
    '高齢化率': df2023['A1303'] / df2023['A1101'] * 100,
    '消費支出': df2023['L3221'] / 1000}))
y = (df2023['A5102'] > df2023['A5101']).astype(int)

model = sm.Logit(y, X).fit()

# 係数と CI
params = model.params
conf = model.conf_int()

# 指数変換してオッズ比とその CI
or_table = pd.DataFrame({
    'OR': np.exp(params),
    'OR_CI_low': np.exp(conf[0]),
    'OR_CI_high': np.exp(conf[1])
})
print(or_table)

📤 実行例（実行時の標準出力）
Optimization terminated successfully.
         Current function value: 0.190567
         Iterations 9
             OR     OR_CI_low    OR_CI_high
const  0.002118  5.087428e-11  88138.979013
高齢化率   2.117222  1.191540e+00      3.762045
消費支出   0.955091  8.980132e-01      1.015797

⚠️ よくある落とし穴

👁️ 直感 — 信頼区間は「真値の入る範囲」

信頼区間（confidence interval, CI）は、 「母数（真の平均など）が一定の確率で含まれる範囲」。 推測統計の中心概念。

20回の試行で計算した95%CI。 約95%（19/20）が真の平均（破線）を含むのが正しい解釈。

⚠️ よくある誤解：「95%CI に真値が95%の確率で入る」→ 不正確。 正しくは「同じ方法で多くの標本から CI を計算すれば、 95%の CI が真値を含む」。 1つの CI に対しては入っているか入っていないかどちらかであり、 確率ではない。

⚠️ 信頼区間の適用条件・限界・誤解回避

信頼区間は推測統計の基礎ツールですが、 解釈を誤ると科学的判断を歪める危険があります。 適用前に下表のチェックリストで前提を確認しましょう。

📋 適用条件チェックリスト

🚨 よくある誤解 TOP 5

🔄 典型ワークフロー

1. 研究設計: 信頼水準 (95%/99%) と必要検出効果サイズを事前決定 → 必要 n を power analysis で算出
2. データ収集: ランダム抽出 + 欠損理由の記録 + 外れ値検査 (箱ひげ図)
3. 分布確認: ヒストグラム/QQプロットで正規性チェック → 大きく外れたら t 分布またはブートストラップへ
4. 区間計算: scipy.stats.t.interval() か statsmodels の関数。 ブートストラップなら scipy.stats.bootstrap
5. 解釈・報告: 「95%CI [a, b]」と区間幅 (b-a) を併記、 効果サイズ (Cohen's d など) も提示
6. 追試・再現性: 別データでも同じ CI 構築手順で結果が安定か確認

📊 信頼区間 vs 関連手法の使い分け

📤 読み方: 信頼区間は「母数 (定数) を含むかどうか」、 予測区間は「次の観測値が入るかどうか」と区別する。 予測区間は CI より常に広い (個別観測の散らばりが加わるため)。

⚠️ 信頼区間の落とし穴 — 実務で必ず引っかかるポイント 7 選

🌐 関連手法・派生

信頼区間は単独で読むものではなく、 点推定 (標本平均など) と標準誤差・p 値・効果量・検出力という「推測統計の 5 点セット」の一翼を担う。 とくに p 値だけ報告して CI を併記しない論文は、 効果の「大きさ」と「不確実性」が読み取れず再現性危機の温床になってきた。

• 上位概念: 推測統計 / 推定理論。 標本から母集団パラメータを推定する大枠の中で、 CI は区間推定の中核を担う。
• 並列手法 (同カテゴリ): 予測区間 (未来の 1 観測の範囲)・許容区間 (母集団の特定割合を含む範囲)・信用区間 (ベイズ的等価物) — 「何の不確実性を表現するか」が異なる。
• 前段処理: サンプルサイズ の設計、 正規性検定 や Q-Q プロット (t 区間の前提)、 外れ値 確認 (平均ベースの CI は外れ値に弱い)。
• 後段処理: 仮説検定 との対応 (CI が 0 を含むか否か = 5% 有意か否か)、 効果量 との併記、 forest plot による複数研究の統合 (メタ解析)。
• 派生・拡張: ブートストラップ CI (パーセンタイル法 / BCa)、 プロファイル尤度 CI、 ベイズの信用区間、 Fieller の比率 CI、 条件付き CI (post-selection inference)。
• 競合・代替: 単一の点推定 + 標準誤差表記、 ベイズ事後分布の可視化、 p 値のみ報告 (推奨されない)。 「区間で語る」文化を組織に根付かせることが目的。

信頼区間は「次に同じ実験を繰り返したらどの幅にパラメータが収まるか」の保証ではなく、 「同じ手順で CI を作る試行を 100 回繰り返したら 95 回はその区間が真値を含む」という頻度論的保証である。 この区別が曖昧なまま使うと、 結論の誤解釈につながる。

🔗 関連用語

🔗 関連用語ネットワーク — 前提・並列・発展で 12 項目

📚 前提となる用語

• 標準誤差（SE） — 信頼区間の幅を決める基本量。 SE = σ/√n。
• t 分布・正規分布 — 信頼区間の臨界値を決める分布。
• 標本サイズ — 区間の幅を決める鍵。 √n に反比例。
• 点推定 — 信頼区間の中心。 まず点推定を理解。

🤝 並列に学ぶ用語

• 仮説検定（p 値） — 信頼区間と表裏一体。 0 を含むか / α を超えるか。
• p 値 — 区間と検定の橋渡し。
• 効果量 — 区間で表現すると不確実性も伝わる。
• 予測区間 — 信頼区間とは別物。 個別観測の不確実性。

🚀 発展先

• ブートストラップ法 — 分布の仮定なしに経験的区間。 BCa などの先進手法。
• ベイズ信用区間（credible interval） — 解釈が直感的な区間。 事後分布から。
• 同時信頼区間 — Tukey HSD、 Scheffé、 Bonferroni 補正（ANOVA 関連）。
• 非劣性試験の信頼区間 — 医薬品開発の中核。 マージン設定。

🔗 さらに関連する基本用語

• t 検定 — 平均差の検定。 信頼区間の双子概念
• 分散分析（ANOVA） — 群間比較
• 平均 — 信頼区間が囲む中心
• 分散 / 標準偏差 — 区間幅の源
• 相関係数 — 信頼区間付き解釈が定番
• 正規分布 — 標本平均の漸近分布

📚 推測統計を学ぶための重要文献

• Fisher (1925) "Statistical Methods for Research Workers" — 古典中の古典
• Neyman & Pearson (1933) — 仮説検定の理論的基礎
• Cohen (1988) "Statistical Power Analysis" — 効果量とサンプルサイズの教科書
• ASA Statement (2016) — p値の正しい使い方
• Gelman & Hill (2007) "Data Analysis Using Regression and Multilevel/Hierarchical Models"
• Wasserman (2004) "All of Statistics" — 現代統計学の総括

🆚 推測統計の主要用語 — 一目で分かる対比表

📖 包括的解説 — この概念を完全マスター

📍 学習の3ステップ

1. 定義を理解する：この概念は何か？ 数式や条件を確認
2. 具体例を見る：実データ（SSDSE 等）で計算してみる
3. 応用する：自分のデータに適用、 結果を解釈

🔧 Python実装パターン

🎯 このコードでやること：信頼区間 — 母平均が含まれると推定される区間に関連するステップ #3。可視化（散布図／樹形図／時系列プロット）を描きます。

📥 入力例（df.head()）
df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1])
df2023 = df[df['SSDSE-B-2026'] == 2023]
# 2023年度47都道府県の消費支出 L3221（二人以上の世帯、 円/月）を母集団とみなし、 平均の 95% CI を計算
# df2023[['Prefecture','L3221']].head():
#    Prefecture   L3221
# 0  北海道      296888
# 12 青森県      263371
# 24 岩手県      298536
# 36 宮城県      305541
# 48 秋田県      272086

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

# 基本パターン
import pandas as pd
import numpy as np
from scipy import stats
import matplotlib.pyplot as plt
import seaborn as sns

# データ読み込み
df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=[1])
df2023 = df[df['SSDSE-B-2026'] == 2023]

# 基本統計量
df2023.describe()

# 可視化（食料費・住居費・教育費）
sns.pairplot(df2023[['L322101', 'L322102', 'L322108']])
plt.show()

📤 このブロックは標準出力には何も出さない

📚 統計概念マップでの位置

このページの上にある3つの概念マップ（関係マップ、 包含マップ、 ツリーマップ）でこの概念の位置づけが視覚的に分かります。 関連手法を辿って学習を進めましょう。

🎯 SSDSE-B-2026 で挑戦

統計データ活用コンペティションのSSDSE-B-2026データは、 47都道府県の社会経済データ。 この概念を使って以下のような分析ができます：

• 地域別の特徴抽出
• 家計支出パターンの解析
• 人口動態と社会経済指標の関連
• 気候要因の影響評価

💡 よく使うコマンド集

🚧 一般的な落とし穴と対策

• 外れ値の影響：散布図・ 箱ひげ図で確認、 ロバスト手法も検討
• サンプルサイズ不足：power analysis で事前に確認
• 仮定の違反：正規性、 独立性、 等分散性をチェック
• 多重比較問題：補正（Bonferroni、 FDR）を適用
• p-hacking：事前登録（pre-registration）で防ぐ
• 因果と相関の混同：観察データから因果結論を出さない

📊 結果報告の標準フォーマット

• 点推定：得られた値
• 不確実性：信頼区間または標準誤差
• サンプルサイズ：n を明記
• 効果量：実質的な意義
• p値：統計的有意性
• 仮定の確認：診断プロット

🌐 関連分野での応用

• マーケティング：A/Bテスト、 顧客分析
• 医療：臨床試験、 疫学研究
• 金融：リスク管理、 ポートフォリオ
• 製造：品質管理、 工程最適化
• 公共政策：効果評価、 計画立案
• 研究：仮説検証、 探索的解析

🎓 さらに学ぶための文献

• Wasserman "All of Statistics"
• Hastie, Tibshirani & Friedman "The Elements of Statistical Learning"
• Gelman & Hill "Data Analysis Using Regression"
• VanderPlas "Python Data Science Handbook"

🔗 統計用語ネットワーク

この概念は、 他の多くの統計概念と密接に関連しています。 ジャストインタイム型学習では、 必要に応じて関連用語へジャンプしながら全体像を構築します。

主要な関連概念のグループ

学習順序の推奨

1. 記述統計（平均、 分散、 標準偏差）
2. 可視化（ヒストグラム、 散布図）
3. 確率分布（正規分布）
4. 推測統計（標準誤差、 信頼区間、 p値）
5. 仮説検定（t検定、 χ²検定）
6. 相関と回帰（単回帰、 重回帰）
7. 多変量解析（PCA、 クラスタリング）
8. 機械学習（決定木、 RF、 NN）
9. 時系列・因果推論（応用）

📝 実践練習 — SSDSE-B-2026 で挑戦

初級課題

1. 東北6県の家計食料費の基本統計量を計算
2. 食料費のヒストグラムを描く
3. 食料費と教育費の散布図を描く
4. 都道府県を「東日本/西日本」に分け、 平均を比較

中級課題

1. 家計支出 5項目で相関行列を作成、 ヒートマップ可視化
2. 食料費 → 教育費の単回帰を実行、 残差分析
3. 家計5項目で PCA を実施、 バイプロット表示
4. k-means (k=3) で都道府県をクラスタリング、 解釈

上級課題

1. 地域別の家計パターンに有意差があるか ANOVA で検定
2. 重回帰で教育費を予測、 多重共線性を VIF で確認
3. Ridge/LASSO で正則化、 CV で α を最適化
4. 階層クラスタリングと Ward 法で都道府県を分類、 デンドログラム作成

🖼️ 信頼区間の追加図解 — 標本サイズと幅の関係

信頼区間は 標本サイズ n が大きくなるほど狭くなります。 下図は SSDSE-B-2026 都道府県人口に対し、 n=10/20/30/47 でリサンプルした際の 95%CI 幅を可視化したものです。 n が増えると CI が縮み、 推定の精度が上がる様子が直感的に分かります。

図: n を増やすほど 95%CI は √n に反比例して狭くなる。 n=100 なら n=25 の半分の幅。

📊 表: 標本サイズ別 95%CI 幅の理論値 (σ=15 と仮定)

📤 読み方: n を 10 倍にすると CI 幅は約 1/3 に縮む。 「精度を倍にしたい」 → 「n を 4 倍に増やす」が経験則。

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

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

📍 体系階層のパス

🌐 統計・データサイエンス › 推測統計 › 推定 › 信頼区間

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

中心に 信頼区間 を置き、 そこから p値・標準誤差・正規分布・標準偏差・標本サイズ・平均推定 など 計 16 個の用語へ、 前提・兄弟・発展形といった関係を矢印で結んでいます。 線がどちら向きかを見れば、 先に読むべき用語とあとから読む用語が分かります。 ノードはドラッグで動かせ、 ホイールでズーム、 クリックでその用語のページへ遷移します。

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

大きな円が小さな円を包含する Circle Packing 図。 「信頼区間」は緑色でハイライト。

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

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

長方形を入れ子に分割した Treemap 図。 各分野の規模感を面積で比較。 「信頼区間」は緑色でハイライト。

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

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

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

🔗 隣接手法への橋渡し

「信頼区間」は推定値の不確実性を範囲で示す指標で、 上流の標準誤差計算と下流の仮説検定 (p 値) との整合を取らないと、 「95% 信頼区間に 0 を含むのに有意」のような矛盾報告が出る。

• 上流: 信頼区間 の入力となるA1101 人口の平均に対する 95% 信頼区間 を整える前処理 (データクレンジング、 欠損値処理)
• 並列: 類似目的を持つ代替手法 (仮説検定、 信用区間) と比較・併用
• 下流: 出力結果の評価・解釈 (p 値、 効果量)

上流で標本サイズと標準誤差を整え、 並列のブートストラップ CI と理論 CI を比較し、 下流の仮説検定・効果量と同時に報告すれば、 SSDSE 都道府県の平均値推定や差の検定を統一的に解釈できる。

🌳 手法選択フロー

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

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

選んだ後の検証ステップ

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

🧭 解説を深める — 「47都道府県の信頼区間」は何を推定しているのか

このページでは 47 都道府県の消費支出（L3221、 2023 年度）から 95% 信頼区間 [288.77, 302.94] を計算しました。 ここで一歩立ち止まって考えたいのが、 「47 都道府県は標本なのか、 それとも母集団そのものなのか」という問いです。 日本の都道府県は 47 しかなく、 SSDSE にはその全部が入っています。 「48 番目の都道府県」を将来抽出することはできません。 全数データに信頼区間を付けるとき、 その区間は何を推定しているのかを意識すると、 信頼区間の理解が一段深まります。

👁️ 直感 — 信頼区間の意味は「問いの立て方」で変わる

同じ数値 [288.77, 302.94] でも、 分析の枠組みによって意味がまったく違います。

• 枠組み A: 47 県を「母集団の全数調査」とみなす — 知りたいのが「2023 年の 47 都道府県平均」そのものなら、 標本平均 295.86 千円は推定値ではなく母平均そのものです。 抽出誤差は存在せず、 論理的には信頼区間は不要（幅ゼロ）。 記述統計として平均と散らばり（s = 24.14 千円、 最小 223.4〔愛媛県〕〜最大 344.1〔埼玉県〕）を報告すれば足ります。
• 枠組み B: 47 県を「超母集団からの 1 実現値」とみなす — 「日本の社会経済構造が生み出す都道府県消費支出のデータ生成過程」という仮想的な母集団（superpopulation、 超母集団）を考え、 観測された 2023 年の 47 値はその 1 回の実現とみなす立場です。 このとき [288.77, 302.94] は「データ生成過程の期待値」に対する区間として意味を持ちます。 経済学・地域分析で全数データに CI を付けるときの標準的な正当化がこれです。
• 枠組み C: 47 県から一部を抽出する — 予算の都合で 20 県しか調査できない、 という状況なら通常の標本抽出の話になります。 ただし母集団が N = 47 と小さいので、 次の落とし穴で述べる有限母集団修正が必要になります。

上の 🎮 シミュレータは、 47 県を「既知の母集団」と宣言したうえで復元抽出する枠組み C の練習装置でした。 論文やレポートで CI を報告するときは、 自分がどの枠組みで語っているのかを明示するのが誠実な書き方です。

⚠️ 落とし穴（重要）— 全数データ・小さな母集団で CI が狂う 3 パターン

• ① 全数データに機械的に CI を付けて「推定」と言ってしまう — 枠組み A のつもり（2023 年の 47 県平均を知りたいだけ）なのに stats.t.interval を回して「母平均の 95% CI は [288.77, 302.94]」と書くと、 存在しない抽出誤差を捏造したことになります。 CI を付けるなら「超母集団の期待値を推定している」（枠組み B）と目的を書き換える必要があります。
• ② 小さな母集団からの非復元抽出で有限母集団修正（FPC）を忘れる — N = 47 県から n = 20 県を非復元で選ぶと、 母集団のかなりの部分（約 43%）を見てしまうため、 標準誤差は無限母集団の式より小さくなります。 修正係数は $$ \text{FPC} = \sqrt{\frac{N-n}{N-1}} = \sqrt{\frac{47-20}{46}} \approx 0.766 $$ で、 SE をこの分だけ縮めます。 FPC を忘れると信頼区間は約 1.31 倍（= 1/0.766）広すぎになり、 「精度が低い」と過小評価してしまいます。 逆に、 シミュレータのような復元抽出や N ≫ n の場合は FPC ≈ 1 となり不要です（n/N < 5% が実務の目安）。
• ③ 「値そのものの誤差」を無視する（二段階の不確実性） — SSDSE-B の消費支出は総務省「家計調査」に基づく値で、 各県の数値自体が世帯標本から推定された統計量です。 つまり本当は「家計調査の抽出誤差（県ごと）」と「県 → 全体をまとめる際の不確実性」の二段構えなのに、 このページの計算は各県の値を誤差ゼロの真値として扱っています。 公的統計の二次利用ではこの「見えない第一段の誤差」がしばしば無視されることを覚えておいてください。

🚀 発展 — 実データで見る「信頼水準と幅」の値段表、 そして推測の 2 学派

枠組み B（超母集団）を採用した場合、 信頼水準を上げる「値段」がどれくらいかを、 同じ実データ（L3221、 2023 年、 n = 47、 平均 295.86 千円、 SE = 3.52 千円）で並べると次のとおりです。

95% → 99% で外す確率を 1/5（5% → 1%）にする代償は幅 +34%（14.18 → 18.93 千円）。 「確実さの購入コスト」が具体的な金額の幅で見えます（値はすべて SSDSE-B-2026 実測に基づく計算値）。

• デザインベース推測 vs モデルベース推測 — 標本調査論では、 不確実性の源を「どの単位を抽出したかの無作為性」に置くデザインベース（枠組み C、 FPC はここに登場）と、 「値を生む確率モデル」に置くモデルベース（枠組み B、 超母集団）の 2 学派があります。 全数データの CI はモデルベースでしか正当化できない、 という整理を知っておくと公的統計の分析論文が読みやすくなります。
• 等重み平均と人口重み平均は別の推定対象 — 47 県の単純平均は東京都と鳥取県を同じ重み 1/47 で扱います。 「日本の平均的な世帯の消費支出」を知りたいなら人口（世帯数）で重み付けた平均が対象であり、 CI もその推定量に対して作り直す必要があります。 推定対象（estimand）を先に決めてから区間を作る、 という順序が本質です。
• クラスター構造と実効サンプルサイズ — 県単位のデータは県内世帯の集計値なので、 「n = 47」は世帯レベルの情報量よりはるかに小さい実効サンプルサイズです。 集計データの CI は個票データの CI と情報量が違う、 という視点は生態学的誤謬（集団の関係を個人に当てはめる誤り）とも接続します。

🔗 関連ページ

• 標準誤差 — FPC は SE に掛かる補正。 まず SE の構造を押さえる。
• 標本抽出 — 復元/非復元、 デザインベース推測の出発点。
• 標本分布と中心極限定理 — 「繰り返し抽出」という思考実験の基礎。
• 区間推定 — 信頼区間・予測区間・信用区間の使い分け。
• サンプルサイズ — 「n = 47 の実効情報量」を考える入口。
• ブートストラップ — 超母集団の考え方を計算機で体感する方法。
• 仮説検定 / p値 — 全数データの検定にも同じ「何を推測するのか」問題が現れる。
• 有限母集団修正（FPC）・超母集団モデル・デザインベース推測 — 個別ページは未整備のため、 本節の説明を参照。