このページを読み解くキーワードと、 該当セクションへのジャンプリンクです。
🍰 まずはやさしく
知りたい答えがある全体の集まりです。
全体の様子を推測するために使います。
全国民の身長を調べるようなケースです。
この章では母集団の基本について読みます。
推測の対象となる全体集合
🍰 まずはやさしく
統計の土台となる考え方です。
一部から全体を推測するために使います。
クラスの代表で全体の意見を知るような例です。
標本との関係についてこれから読みます。
このページは 推測統計カテゴリの土台となる「母集団」を解説しています。 同じカテゴリには 標本・標本抽出・サンプルサイズ・標準誤差・信頼区間 が並びます。
読み終えると、 「標本から母集団を推測する」という統計の中心テーマが見える視点を獲得できます。 続けて 📚 標本抽出と中心極限定理 を読むと、 母集団とのつながりが一本の流れで理解できます。
🍰 まずはやさしく
全部と一部の関係のようなものです。
現実的に調べられない全体を知るために使います。
スマホのアンケートで全ユーザーを考える例です。
全体をどう定義するかについてこれから読みます。
この概念は「標本から母集団を推測する」考え方の一部です。 標本サイズと不確実性のセットで理解しましょう。
本ページは 母集団 を、 SSDSE-B-2026 47 都道府県を題材に「全数か標本か」の解釈問題として整理する。 母集団を「全国都道府県」と定義した場合と「将来の社会状態」と定義した場合で、 同じ平均値の解釈 (記述統計 vs 推測統計) がどう変わるかを実例で示す。
想像してみてください。 全国 1.24 億人の身長を知りたいとします。 全国民を計測(全数調査)すれば文句なしの真値が出ますが、 時間・予算・物理的に不可能。 そこで、 1,000 人をうまく選んで身長を測り、 そこから全体像を「推し測る」 — これが推測統計の出発点であり、 母集団/標本という言葉が生まれた現場です。
本ページのテーマである 母集団 は、 知りたい「全体」のことです。 SSDSE-B-2026 を例にすると:
同じデータでも「母集団をどう定義するか」で意味が変わります。 47 県の人口を「47 件の母集団」と見れば μ = 2,645,808 人と確定値ですが、 「日本の都道府県という観念的な集団」と見れば、 47 県は 1 つの実現にすぎず、 推測の対象が広がります。
🌳 比喩: 母集団は「森全体」、 標本は「木を 30 本切り出した小山」。 切り出した木の平均高さ x̄ は森の平均高さ μ そのものではないが、 大きな手がかりになる。 標準誤差は「森の平均と切り出し平均がどれだけズレうるか」の見積もり。
🍰 まずはやさしく
推測の対象となる全体の集合です。
分析やレポートで正しく伝えるために使います。
部活の全員をひとつのグループとする例です。
母平均などの詳しい定義についてこれから読みます。
推測の対象となる全体集合
英語名 Population。
この用語を理解・使用するときは、 次のような前提を意識してください:
サイズ $N$ の有限母集団を $\{X_1, X_2, \dots, X_N\}$ とすると、 母平均 $\mu$、 母分散 $\sigma^2$、 母標準偏差 $\sigma$ は次のように定義されます。
$$\mu \;=\; \frac{1}{N}\sum_{i=1}^{N} X_i$$
$$\sigma^2 \;=\; \frac{1}{N}\sum_{i=1}^{N} (X_i - \mu)^2 \qquad \sigma \;=\; \sqrt{\sigma^2}$$
2 値(例: 賛成/反対)の場合、 母比率 $p$ は
$$p \;=\; \frac{\text{該当する個体数}}{N}$$
標本サイズ $n$ の無作為標本 $\{x_1, \dots, x_n\}$ から計算する標本平均 $\bar{x}$ と不偏標本分散 $s^2$ は
$$\bar{x} \;=\; \frac{1}{n}\sum_{i=1}^{n} x_i \qquad s^2 \;=\; \frac{1}{n-1}\sum_{i=1}^{n} (x_i - \bar{x})^2$$
標本平均の標準誤差 SE は
$$\text{SE}(\bar{x}) \;=\; \frac{\sigma}{\sqrt{n}} \;\approx\; \frac{s}{\sqrt{n}}$$
有限母集団 (有限修正) の場合は $\text{SE} = \frac{\sigma}{\sqrt{n}} \sqrt{\frac{N-n}{N-1}}$ と補正します。
| 記号 | 読み方 | 直感的意味 |
|---|---|---|
| $N$ | 「ラージエヌ」 | 母集団全体のサイズ。 47 都道府県なら $N=47$。 全国民なら $N \approx 1.24\times 10^8$。 |
| $X_i$ | 「エックスアイ」 | 母集団の $i$ 番目の個体の値(例: 北海道の人口 5,092,000)。 |
| $\mu$ | 「ミュー」 | 母平均。 母集団全体の平均値。 真値。 |
| $\sigma$ | 「シグマ」 | 母標準偏差。 母集団のバラつきの真値。 |
| $n$ | 「スモールエヌ」 | 標本サイズ(観測した個数)。 |
| $\bar{x}$ | 「エックスバー」 | 標本平均(観測から計算した推定値)。 |
| $s$ | 「エス」 | 標本標準偏差(不偏分散の平方根、 分母 $n-1$)。 |
| SE | 「エスイー」 | 標準誤差。 $\bar{x}$ が $\mu$ からどれだけズレうるかの目安。 |
記号上は $\mu, \sigma$ がギリシア文字(母集団のパラメータ)、 $\bar{x}, s$ がラテン文字(標本統計量)と慣習で書き分けます。 「ギリシア = 真値(普通は知りえない)」「ラテン = 推定値(標本から計算可能)」と覚えると混乱しません。
SSDSE-B-2026 の 2023 年 47 都道府県人口を「母集団」 と見立てて、 母平均 $\mu$ と母分散 $\sigma^2$、 そして「ある標本」 と比較します。
| 統計量 | 記号 | 値 |
|---|---|---|
| サイズ | $N$ | 47 |
| 合計 | $\sum X_i$ | 124,353,000 人(≒ 1.244 億人) |
| 母平均 | $\mu$ | 2,645,808.51 人 |
| 母標準偏差 | $\sigma$ | 2,767,630.21 人 |
| 最小 | $X_\min$ | 537,000 人(鳥取県) |
| 最大 | $X_\max$ | 14,086,000 人(東京都) |
| μ ± σ 範囲内 | — | 41 / 47 県 (87.2%) |
→ σ がほぼ μ と同じ大きさ。 つまり都道府県人口は非常にバラつきが大きい(東京都が外れ値級)。 正規分布ではないので「μ±σ に 68%」 という経験則は機能せず、 実際は 87% が μ±σ 内。
都道府県コード順に並べて 5 番目ごとに 1 県取ると 10 件: 北海道, 山形県, 埼玉県, 富山県, 岐阜県, 京都府, 鳥取県, 徳島県, 佐賀県, 鹿児島県。 これを 1 つの標本と見なします。
| 統計量 | 記号 | 値 |
|---|---|---|
| 標本サイズ | $n$ | 10 |
| 標本平均 | $\bar{x}$ | 2,249,800 人 |
| 標本標準偏差 (不偏) | $s$ | 2,238,068.50 人 |
| 標準誤差 | $s/\sqrt{n}$ | 707,739.40 人 |
| 真値とのズレ | $\mu - \bar{x}$ | 396,008.51 人 |
→ 標本平均 $\bar{x} = 2{,}249{,}800$ は母平均 $\mu = 2{,}645{,}808$ から 約 39.6 万人小さい。 SE = 707,739 のオーダで「これくらいのズレは起こりうる」ことが説明される(実際 $\mu$ は $\bar{x} \pm$ SE には入らないが、 $\bar{x} \pm 2\cdot\text{SE}$ の区間 [834,321, 3,665,279] には収まる)。
合成データで母平均推定の n と信頼区間幅の関係を計算する。
| n | SE | 区間幅 (2·1.96·SE) |
|---|---|---|
| 25 | 2.00 | 7.84 |
| 100 | 1.00 | 3.92 |
| 400 | 0.50 | 1.96 |
| 1600 | 0.25 | 0.98 |
1 2 3 4 5 | import numpy as np sigma = 10 n = np.array([25, 100, 400, 1600]) width = 2 * 1.96 * sigma / np.sqrt(n) print(f"区間幅: {width.round(2)}") |
💬 手計算 (Step 1) と Python 出力が完全一致。
下の図は SSDSE-B-2026 の 2023 年・47 都道府県人口(実データ)を母集団とみなしたものです。 灰色の点が母集団の全 47 県、 緑の点がいま抽出している標本です。 スライダーで標本サイズ $n$ を変えたり、 図の上を左右にドラッグ/スワイプしたり、 「🎲 新しい標本」で抽出し直すと、 標本平均 x̄(青線)が 母平均 μ(赤線)の周りでどう揺れ、 $n$ を大きくするほど μ に近づいて安定していく(大数の法則)ことを体感できます。
↑ 横軸=人口。 灰=母集団 47 県、 緑=標本。 図の上を左右にドラッグ/スワイプしても n を変えられます。
| 量 | 母集団(真値・固定) | 標本(n から毎回計算) |
|---|---|---|
| 平均 | – | – |
| 標準偏差 | – | – |
| |μ − x̄|(推定のズレ) | — | – |
| 標準誤差 SE(理論・有限修正) | — | – |
「📊 1000回抽出」を押すと、 同じ n で 1000 回だけ標本を取り直したときの 標本平均 x̄ の分布(標本分布) を描きます。 n が大きいほど分布は μ(赤線)の周りに細く集中します — これが標準誤差の縮小、 そして中心極限定理の入り口です。
🔑 ポイント: μ と σ は母集団だけで決まる動かない真値。 一方 x̄ と s は標本ごとに変わる推定値です。 「🎲 新しい標本」を何度も押すと x̄ が毎回違う値になる — これが標本誤差。 n を 47 に近づけると標本は母集団そのものに近づき、 x̄ → μ、 SE → 0 に収束します。
① 直感(全体 vs 一部): 母集団は「本当に知りたい全体」(ここでは全 47 県)、 標本は「実際に手元で観測した一部」(緑の点)です。 全数を測れる全数調査(センサス)なら μ が直接わかりますが、 コストや時間の制約で普通は一部しか見られない。 だから標本から母集団を「推し測る」=📚 標本抽出と中心極限定理で扱う推測統計が要ります。
② よくある落とし穴: スライダーを小さな n(例 3〜5)にして「新しい標本」を連打すると、 x̄ が大きく暴れるのが見えます。 このとき手元の標本の x̄ を「母集団の平均そのもの」だと思い込む標本と母集団の混同が典型的な誤りです。 さらに、 抽出が無作為でなく非代表標本(例:大都市だけ、 過疎県だけ)になると x̄ は μ から系統的にズレます — これが選択バイアスで、 n を増やしても解消しません(詳しくは上の「⚠️ よくある落とし穴」節も参照)。
③ 発展: 「📊 1000回抽出」でできる x̄ のヒストグラムが標本分布(sampling distribution)です。 n を上げると分布幅(=標準誤差 $\sigma/\sqrt{n}$)が縮み、 形は母集団が非正規でも釣鐘型に近づく — これが中心極限定理で、 x̄ が μ に集中していく現象が大数の法則です。 この標本分布の幅から信頼区間を作り、 必要なサンプルサイズを逆算できます。
SSDSE-B-2026 のような公的統計データを Python で扱う際の基本パターン:
1 2 3 4 5 6 7 8 9 10 11 12 | import pandas as pd import numpy as np # データ読み込み df = pd.read_csv('data/raw/SSDSE-B-2026.csv', encoding='cp932', skiprows=1) print(df.shape) print(df.dtypes) print(df.describe()) # 「母集団」の文脈で扱う場合の例: # 分野: 推測統計 # 関連手法は同カテゴリの他用語を参照してください。 |
具体的なコードは 標本抽出と中心極限定理 を参照してください。
分析結果を報告するときに含めるべき情報:
🎯 このコードでやること: SSDSE-B-2026 から 2023 年の 47 都道府県人口を抜き出し、 母平均 μ と母分散 σ² を pandas で計算する。 母分散は分母 $N$ (ddof=0)、 標本分散は分母 $n-1$ (ddof=1) — この使い分けが核心。
📥 入力データ (SSDSE-B-2026 抜粋、 2023 年分):
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | import pandas as pd df = pd.read_csv('data/raw/SSDSE-B-2026.csv', skiprows=1, encoding='cp932') latest = df[df['年度'] == 2023].copy() pop = latest['総人口'] N = len(pop) mu = pop.mean() # 算術平均 = 母平均 sigma2 = pop.var(ddof=0) # 母分散 (分母 N) sigma = pop.std(ddof=0) # 母標準偏差 print(f'N = {N}') print(f'μ = {mu:,.2f}') print(f'σ² = {sigma2:,.2f}') print(f'σ = {sigma:,.2f}') print(f'合計 = {pop.sum():,}') |
📤 実行すると次の出力が得られる:
💬 結果の読み方: 47 県を 1 つの母集団と見たときの真値。 μ = 264.6 万人だが、 σ も 276.8 万人とほぼ同じ大きさ — つまり 47 県の人口分布は右に大きく裾を引く非対称で、 東京都 (1,408 万人) が μ + 4σ 付近の外れ値級。 「全国平均」 を語るとき、 県を単位にするか個人を単位にするかで意味が大きく違うことが、 数値からも見えてくる。
🎯 このコードでやること: 47 県から決定論的に 10 県(5 件おき)を抽出して標本とし、 標本平均 $\bar{x}$、 標本標準偏差 $s$、 標準誤差 SE を計算。 母平均 μ との差を見て、 推定誤差を体感する。
📥 入力データ: 上記 latest から 5 件おきに .iloc[::5] で抽出(北海道, 山形, 埼玉, 富山, 岐阜, 京都, 鳥取, 徳島, 佐賀, 鹿児島)。
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 | import numpy as np sample = latest['総人口'].iloc[::5] n = len(sample) x_bar = sample.mean() s = sample.std(ddof=1) # 不偏標本標準偏差 (分母 n-1) se = s / np.sqrt(n) print(f'標本サイズ n = {n}') print(f'標本平均 x̄ = {x_bar:,.2f}') print(f'標本標準偏差 s = {s:,.2f}') print(f'標準誤差 SE = {se:,.2f}') print(f'母平均との差 μ - x̄ = {mu - x_bar:,.2f}') print(f'x̄ ± 2·SE 区間: [{x_bar-2*se:,.0f}, {x_bar+2*se:,.0f}]') print(f'真値 μ はこの区間に含まれるか: {x_bar-2*se <= mu <= x_bar+2*se}') |
📤 実行すると次の出力が得られる:
💬 結果の読み方: 標本平均は母平均から約 39.6 万人小さくズレた。 だが「$\bar{x} \pm 2\cdot$SE 区間」 (≒ 95% 信頼区間の素朴版) を作ると、 真値 μ = 264.6 万人はちゃんと中に入る — これが「標本から母集団を推測する」の感触。 標本は完璧な真値ではないが、 不確実性を込みで報告すれば真値は射程内。
🎯 このコードでやること: 「人口 100 万人以上の県の割合」 を母比率 $p$ として 47 県全体(全数)から計算。 さらに 10 県標本からの推定 $\hat{p}$ と比較。
📥 入力データ: 上記 latest['総人口'] シリーズ(47 行)。
1 2 3 4 5 6 7 8 9 10 11 12 13 | threshold = 1_000_000 # 100万人 # 母集団 (全数) pop_flag = (latest['総人口'] >= threshold) p_true = pop_flag.mean() print(f'母比率 p (>=100万人の県の割合) = {p_true:.4f} ({pop_flag.sum()}/{len(pop_flag)})') # 標本 (n=10) samp_flag = (sample >= threshold) p_hat = samp_flag.mean() se_p = np.sqrt(p_hat*(1-p_hat)/len(samp_flag)) print(f'標本比率 p̂ = {p_hat:.4f}, SE(p̂) = {se_p:.4f}') print(f'母比率との差 p - p̂ = {p_true - p_hat:+.4f}') |
📤 実行すると次の出力が得られる:
💬 結果の読み方: 全数では 47 県のうち 37 県が 100 万人以上 (78.7%)。 一方、 たまたま 10 県標本では 7 県が 100 万人以上で $\hat{p}=0.700$ と母比率をやや下回る — これが標本の偶然による偏り。 標本が小さい (n=10) と SE(p̂)=0.14 と大きくなり、 点推定だけで母比率を語るのは危険。 比率 $p$ が 1 や 0 に近い局面では正規近似式の精度も悪くなるため、 Wilson 区間など区間推定を併用すべき、 という注意点も読み取れる。
🎯 このコードでやること: 47 県人口のヒストグラムを Sturges のビン数で作り、 母集団の形(右に偏った長い尾)を可視化する。 数式 μ, σ だけでは見えない「分布の形」を補う。
📥 入力データ: latest['総人口'](47 行、 単位: 人)。
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 | import math import matplotlib.pyplot as plt values = latest['総人口'].values k = math.ceil(math.log2(len(values)) + 1) # Sturges counts, edges = np.histogram(values, bins=k) print(f'Sturges k = {k}') for i, c in enumerate(counts): print(f' [{edges[i]:>10,.0f}, {edges[i+1]:>10,.0f}) : {c} 県') plt.figure(figsize=(8,4)) plt.hist(values, bins=k, color='#26A69A', edgecolor='#004D40') plt.axvline(mu, color='red', ls='--', label=f'μ={mu:,.0f}') plt.xlabel('総人口 (人)'); plt.ylabel('県の数') plt.title('47都道府県人口の母集団分布 (2023)') plt.legend(); plt.tight_layout() # plt.savefig('population_hist.png', dpi=120) |
📤 実行すると次の出力が得られる:
💬 結果の読み方: 最初のビン(250 万人未満)に 47 県中 34 県 (72%) が集中、 残りの 5 ビンに 12 県、 最後のビン (1,215〜1,409 万人) に東京都 1 県のみ。 母集団は明らかに右に長く裾を引く非正規分布。 平均 μ だけ報告すると東京の影響で「典型的な県」 を誤解しやすい。 メジアン (1,549,000 人) も併記すべき。
🎯 このコードでやること: SSDSE-B-2026 の 47 県を「人口大 (上位 1/3)・中 (中央 1/3)・小 (下位 1/3)」 の 3 層に分け、 各層から 3 県ずつ抜いて層化標本 (n=9) で母平均を推定する。 同じ n の単純無作為と比べて精度がどう変わるかを見る。
📥 入力データ: latest (47 行、 2023 年人口)。 層化基準は pd.qcut(pop, 3, labels=['小','中','大'])。
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 pandas as pd, numpy as np # 3 層に分割 (qcut で人口三分位) latest['層'] = pd.qcut(latest['総人口'], q=3, labels=['小','中','大']) # 各層の母平均・母分散 strat_summary = latest.groupby('層', observed=True)['総人口'].agg(['count','mean','std']) print(strat_summary) # 各層から先頭の 3 県 (決定論的抽出) sample_strat = ( latest.sort_values(['層','地域コード']) .groupby('層', observed=True) .head(3) ) print('--- 層化標本 (n=9) ---') print(sample_strat[['都道府県','層','総人口']].to_string(index=False)) # 層化平均: 各層平均を母集団の層構成比で重み付け weights = latest['層'].value_counts(normalize=True) layer_means = sample_strat.groupby('層', observed=True)['総人口'].mean() mu_strat = (layer_means * weights).sum() print(f'\n層化推定 μ̂ = {mu_strat:,.2f}') print(f'真値 μ = {latest["総人口"].mean():,.2f}') |
📤 実行すると次の出力が得られる:
💬 結果の読み方: 「層化」 と謳っても、 各層内で頭から取ると層内の代表性が壊れる(「大」 層の頭 3 件は北海道・宮城・茨城で、 東京や大阪が入らない)。 真の層化抽出は層内で無作為に取る必要がある。 ここでは決定論的に頭 3 件を取ったため真値より低く推定された — 「機械的な抽出は層化の名を借りた偏りを生む」 という教訓。
🎯 このコードでやること: n が N に近づくと SE がどう縮むかを、 有限母集団修正項 $\sqrt{(N-n)/(N-1)}$ をかけて確認する。 $n/N$ が 5% を超えるあたりから無視できない。
📥 入力データ: $N=47$, $\sigma = 2,767,630$ (上で計算した値)。
1 2 3 4 5 6 7 8 9 10 11 | import numpy as np N = 47 sigma = 2_767_630.21 print(f'{"n":>4} {"n/N":>6} {"SE_naive":>14} {"FPC":>7} {"SE_finite":>14} {"差%":>7}') for n in [5, 10, 15, 20, 30, 40, 45, 47]: se_naive = sigma / np.sqrt(n) fpc = np.sqrt((N-n)/(N-1)) se_finite = se_naive * fpc diff_pct = (1 - se_finite/se_naive) * 100 print(f'{n:>4} {n/N:>6.2%} {se_naive:>14,.0f} {fpc:>7.3f} {se_finite:>14,.0f} {diff_pct:>6.1f}%') |
📤 実行すると次の出力が得られる:
💬 結果の読み方: $n/N$ が 5% 程度なら修正の差は約 5%、 20% 程度で約 10% 縮む。 $n=N$ で SE=0(全数調査)。 一般的な指針として「$n/N > 5\%$ なら有限修正を入れる」 が安全側。 SSDSE-B の 47 県を分析するとき、 10 県以上抜いたらほぼ常に修正が必要なレベル。
分析の最初に必ず通る関門。 ここをスキップすると、 結論が誰にも適用できなくなります。 SSDSE-B-2026 の都道府県人口分析を例に、 7 つの問いに答えてみましょう。
| 問い | SSDSE-B 例での答え |
|---|---|
| ① 単位は何か (人・世帯・県・事業所…) | 県(47 件)。 人を単位にする場合は重み付け平均が必要。 |
| ② 地理範囲はどこまでか | 日本国内 47 都道府県(北方領土・離島は含む)。 |
| ③ 時点はいつか | 2023 年 10 月 1 日時点。 時点が違えば母集団も違う。 |
| ④ 包含基準は何か (誰を入れる) | 住民基本台帳に登録された人(日本人 + 一定要件の外国人)。 短期滞在者・国外在住の日本国籍者は除外。 |
| ⑤ 除外基準は何か | 3 ヶ月未満の在留外国人、 外交官、 国外定住者。 |
| ⑥ 観測 (測定) はどう行うか | 市区町村の住民票登録の集計。 全数で測定誤差ほぼゼロ。 |
| ⑦ 結論はどこまで一般化するか | 2023 年の日本の県別構造についてのみ。 他年・他国・市町村レベルに拡張するときは別途検討。 |
このチェックリストを論文・レポートの「データと方法」 セクションに必ず書くのが、 専門的な統計レポートのお作法。 「母集団は何か」 を明示するだけで、 読者がどの程度結論を信じてよいかの判断ができる。
統計学では「ギリシア文字 = 母集団パラメータ」「ラテン文字 (or ハット付き) = 標本統計量」 という記法の慣習があります。 SSDSE-B-2026 の値を入れて並べると理解が早い。
| 概念 | 母集団 (真値) | 標本 (推定値) | SSDSE-B 例 |
|---|---|---|---|
| サイズ | $N$ | $n$ | $N=47$ vs $n=10$ |
| 平均 | $\mu$ | $\bar{x}$ または $\hat{\mu}$ | $\mu=2{,}645{,}808$ vs $\bar{x}=2{,}249{,}800$ |
| 分散 | $\sigma^2$ | $s^2$ または $\hat{\sigma}^2$ | $\sigma^2 \approx 7.66\times 10^{12}$ |
| 標準偏差 | $\sigma$ | $s$ | $\sigma=2{,}767{,}630$ vs $s=2{,}238{,}068$ |
| 比率 | $p$ または $\pi$ | $\hat{p}$ | 100万人以上の県 $p=0.936$ vs $\hat{p}=1.000$ |
| 中央値 | $\tilde{\mu}$ (= med($X$)) | $\tilde{x}$ | $\tilde{\mu}=1{,}649{,}000$ (奈良県) |
| 相関係数 | $\rho$ | $r$ | 総人口 vs 65歳以上人口 の $r$(強い正の相関) |
| 回帰係数 | $\beta$ | $\hat{\beta}$ または $b$ | 傾き真値 vs OLS 推定値 |
| 分散 (回帰残差) | $\sigma_\varepsilon^2$ | MSE | 残差平均二乗誤差 |
論文や教科書を読むとき、 まず「ギリシア記号 = 真値、 ラテン記号 = 推定値」 の置換ができれば 90% は理解できる。 さらに「ハット (^) は推定値」 という慣習も覚えておく ($\hat{\theta}$ は θ の推定量)。
統計の不確実性は 2 種類あります。 SE が表せるのは前者だけ。 後者は SE では捉えられないので、 設計と運用で潰すしかありません。
| 種類 | 原因 | 対策 | SE で捉えられるか |
|---|---|---|---|
| 標本誤差 | 母集団から一部しか観測しないことによる偶然のズレ | $n$ を増やす、 層化抽出、 適切な推定量 | ✅ できる ($\sigma/\sqrt{n}$) |
| カバレッジ誤差 | 標本枠 (frame) に母集団の一部が含まれない | 最新の住民台帳を使う、 マルチフレーム抽出 | ❌ できない |
| 無回答誤差 | 調査対象が回答しない、 回答内容が偏る | 回収率を上げる、 ウェイト調整、 多重代入 | ❌ できない (本質的にバイアス) |
| 測定誤差 | 質問のあいまいさ、 計測器の偏り、 入力ミス | 質問設計、 再テスト、 計測器校正、 検算 | ❌ できない |
| 処理誤差 | 集計プログラムのバグ、 コード変換ミス | テスト、 ダブルチェック、 再現可能パイプライン | ❌ できない |
よくある勘違い: 「$n$ を 100 万に増やせば SE は限りなく小さくなり、 真値が確定する」 — これは 標本誤差に対してのみ正しい。 オンライン投票の $n$ が 10 万でも、 そもそも特定の層しかアクセスしていなければ、 結論は母集団とずれたまま。 SE 値の小ささを「精度」 と混同しないこと。
1 変数なら母平均 $\mu$ と母分散 $\sigma^2$ だけ気にすれば良いですが、 多変量だと母平均ベクトル $\boldsymbol{\mu}$ と 母共分散行列 $\boldsymbol{\Sigma}$ を考えます。
$$\boldsymbol{\mu} = \begin{pmatrix}\mu_1\\\mu_2\\\vdots\\\mu_p\end{pmatrix}, \quad \boldsymbol{\Sigma} = \begin{pmatrix}\sigma_{11} & \sigma_{12} & \cdots & \sigma_{1p}\\\sigma_{21} & \sigma_{22} & \cdots & \sigma_{2p}\\\vdots & & \ddots & \vdots\\\sigma_{p1} & \sigma_{p2} & \cdots & \sigma_{pp}\end{pmatrix}$$
対角要素 $\sigma_{ii}$ は各変数の母分散、 非対角要素 $\sigma_{ij}$ は変数 $i, j$ の母共分散。 SSDSE-B から「総人口」「15歳未満人口」「65歳以上人口」 の 3 変数を取り出すと、 各々の母平均と 3×3 の母共分散行列が定義できる。 主成分分析・因子分析・線形判別分析はすべて $\boldsymbol{\Sigma}$ を推定して使う手法群。
標本からの推定は標本平均ベクトル $\bar{\mathbf{x}}$ と標本共分散行列 $\mathbf{S}$ で、 やはり分母 $n-1$ の不偏推定を使う ($\bar{\mathbf{x}} = \frac{1}{n}\sum \mathbf{x}_i$, $\mathbf{S} = \frac{1}{n-1}\sum (\mathbf{x}_i-\bar{\mathbf{x}})(\mathbf{x}_i-\bar{\mathbf{x}})^\top$)。 $n < p$ だと $\mathbf{S}$ が特異 (逆行列なし) になるため、 高次元では正則化が必要。
これまでの説明は頻度論 (frequentist) の立場。 母平均 μ は「固定された真値」 とみなし、 標本だけがランダム。 これに対し ベイズ (Bayesian) 流では母平均自体を確率変数として扱います。
| 観点 | 頻度論 | ベイズ |
|---|---|---|
| μ の扱い | 固定された真値(未知) | 確率分布を持つ(事前分布 → 事後分布) |
| 推定の表現 | 点推定 + 信頼区間 | 事後分布全体(信用区間 = credible interval) |
| 「95% 区間」 の意味 | 標本を何度も取れば 95% の区間が真値を含む | μ がこの区間に入る確率 95% |
| 事前知識の扱い | 原則として使わない | 事前分布として明示的に取り込む |
| 代表ツール | scipy.stats, statsmodels | PyMC, Stan, NumPyro |
SSDSE-B の県人口推定でいえば、 頻度論は「47 県の母平均は固定値、 われわれは未知」 と扱う。 ベイズなら「過去の年次データから事前分布を作り、 2023 年の観測でアップデート → 事後分布」 と扱える。 県別の小地域推定 (small area estimation) では、 標本が薄い県のために隣接県の情報を借りる階層ベイズモデルが現代統計の標準ツール。
雑誌 The Literary Digest は 240 万票という当時史上最大級の標本でルーズベルト落選を予測しましたが、 実際はルーズベルトの圧勝。 標本は電話帳・自社購読者・自動車登録者 から作られていて、 1936 年の米国でこれらを持つのは富裕層 — つまり標本枠が共和党寄りに大きく偏っていた。 一方、 ギャラップは僅か 5 万人の無作為抽出標本で正確に予測。 教訓: $n$ より、 母集団を代表する標本枠が大切。
米軍は帰還した爆撃機の弾痕分布を分析し「弾痕の多い部位を装甲しよう」 と考えました。 統計学者 Abraham Wald は逆を提案 — 弾痕のない部位こそ装甲すべき。 母集団は「全爆撃機」 だが、 観測標本は「帰還できた」 爆撃機のみ。 弾痕のない部位(コックピット・エンジン)に被弾した機は帰還できず観測されない。 教訓: 誰が観測から漏れているかを考えないと、 母集団についての結論は逆になる。
米カリフォルニア州サンタクララ郡の住民 3,330 人の抗体検査で陽性率 1.5%、 「市の感染者数は公式の 50-85 倍」 という結論。 だが標本は Facebook 広告で募集。 体調不良の人ほど応募しやすく、 偽陽性率の補正不足もあり、 推定値は過大評価と批判された。 教訓: 標本の集め方が母集団から漏れた瞬間、 推定値は意味を失う。
SSDSE-B-2026 で「47 県の人口平均 = 264.6 万人」 と求めても、 これは県を単位とした平均であって、 日本国民の典型値とは違う。 国民を単位とした平均は「総人口 / 47 県」 ではなく、 そもそも 1 人 = 1 名のスケールで考えるべき。 中央値で県平均は鹿児島県の 155 万人だが、 国民の中央値居住県は人口加重で東京方面に寄る。 教訓: 「単位は何か」 を間違えると、 同じ数値でも結論が逆転する。
🎯 このコードでやること: 母平均 μ は東京都の外れ値に大きく引きずられる。 ロバストな代替指標として中央値、 四分位範囲 (IQR)、 中央絶対偏差 (MAD) を計算し、 「典型的な県」 像を別角度から描く。
📥 入力データ: latest['総人口'] (47 件、 2023 年)。
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 pandas as pd, numpy as np pop = latest['総人口'] # 基本ロケーション指標 mean = pop.mean() median = pop.median() q25, q75 = pop.quantile([0.25, 0.75]) iqr = q75 - q25 mad = (pop - median).abs().median() print(f'平均 μ = {mean:>12,.0f} 人') print(f'中央値 median = {median:>12,.0f} 人') print(f'第1四分位 Q1 = {q25:>12,.0f} 人') print(f'第3四分位 Q3 = {q75:>12,.0f} 人') print(f'IQR (Q3-Q1) = {iqr:>12,.0f} 人') print(f'MAD = {mad:>12,.0f} 人') print(f'平均-中央値の差 = {mean-median:>12,.0f} 人 (右に歪んでいる)') # 外れ値判定 (Tukey 法: Q3 + 1.5*IQR より上) upper = q75 + 1.5*iqr outliers = latest[pop > upper][['都道府県','総人口']].sort_values('総人口', ascending=False) print(f'\n外れ値判定基準 (Q3 + 1.5·IQR) = {upper:,.0f}') print(outliers.to_string(index=False)) |
📤 実行すると次の出力が得られる:
💬 結果の読み方: 平均 264.6 万人と中央値 154.9 万人で約 110 万人もズレている — これが「右に歪んだ分布」 のサイン。 Tukey 基準では 9 県 (上位の大都市圏) が外れ値判定。 「典型的な県」 を 1 つ選ぶなら平均より中央値、 さらに散らばりは σ より IQR や MAD の方が外れ値に頑健。 報告書では μ, σ だけでなく中央値・四分位も併記するのが標準。
🎯 このコードでやること: SSDSE-B-2026 の 2012-2023 年における 47 県人口の母平均 $\mu_t$ の推移を計算し、 母集団のパラメータが時間と共にどう変化しているかを確認する。
📥 入力データ: df(年度 × 47 県 × 多変量)。 年度をキーにグループ化。
1 2 3 4 5 6 7 8 | import pandas as pd annual = (df.groupby('年度')['総人口'] .agg(['mean','std','min','max']) .reset_index()) annual['mean'] = annual['mean'].astype(int) annual['std'] = annual['std'].astype(int) print(annual.to_string(index=False)) |
📤 実行すると次の出力が得られる:
💬 結果の読み方: 母平均 (47 県平均) は 2012 年の 271.6 万人から 2023 年の 264.6 万人へ約 7 万人減少 (-2.6%)。 同時に最小値 (鳥取県人口) も 583 千 → 537 千に減少、 最大値 (東京) は 13.23M → 14.09M に増加 → 分散の変動と東京一極集中が見える。 母集団パラメータは「真値」 とはいえ時間で動く — これが「時点を明示する」 重要性の正体。
レポートや論文の「データ・方法」 セクションで、 母集団を曖昧にせず操作的に書く型を 4 つ示します。
「本分析の母集団は、 2023 年 10 月 1 日時点の住民基本台帳に登録された全 47 都道府県($N=47$)である。 観測単位は県とし、 県別総人口を分析変数とした。 全数調査のため標準誤差は付与しない。」
「本分析の母集団は、 「2020 年代の日本における都道府県人口を生成する潜在過程」 と定義する。 観測された 2023 年 47 県データを、 この母集団からの 1 つの実現 (標本) として扱い、 母平均 μ の点推定と 95% 信頼区間を報告する。」
「本分析の母集団は、 2023 年 10 月 1 日時点で日本に住所を持つ全居住者 ($N \approx 1.243 \times 10^8$ 人) である。 SSDSE-B-2026 から各県の総人口を取得し、 県を 1 ユニットとして集計値を分析した。 個人レベルの分布特性(年齢・性別など)はカバーしない。」
「本研究で構築した県別需要予測モデルの汎化母集団は、 将来 5 年以内の日本国内 47 都道府県とする。 SSDSE-B-2026 の 2012-2022 年データを訓練、 2023 年データをホールドアウト検証に充てた。 母集団分布が大きく変化した場合 (人口構造の急変等)、 モデルの妥当性は別途検証する。」
どの型でも共通するキーワード: 時点・地理範囲・観測単位・包含基準・除外基準。 この 5 項目が明示されていれば、 読者は結論の射程を判断できます。
「東日本と西日本で人口に差があるか」「介入群と対照群で効果があるか」 という比較分析では、 母集団は2 つになります。 それぞれに $\mu_1, \mu_2$ が存在し、 差 $\Delta = \mu_1 - \mu_2$ を推定するのが目標。
$$\hat{\Delta} = \bar{x}_1 - \bar{x}_2, \quad \text{SE}(\hat{\Delta}) = \sqrt{\frac{\sigma_1^2}{n_1} + \frac{\sigma_2^2}{n_2}}$$
SSDSE-B-2026 の例で「東日本 24 県 (北海道〜長野県) vs 西日本 23 県 (岐阜〜沖縄)」 の母平均を比較してみると、 計算上は
| グループ | $n$ | 平均 | SD |
|---|---|---|---|
| 東日本 (北海道〜長野県) | 24 | 2,776,167 人 | 3,257,825 人 |
| 西日本 (岐阜県〜沖縄県) | 23 | 2,510,217 人 | 2,180,889 人 |
差 $\hat{\Delta} = 2,776,167 - 2,510,217 = +265,950$ 人 (東日本の方が若干大きい)。 SE$(\hat{\Delta}) = \sqrt{3{,}257{,}825^2/24 + 2{,}180{,}889^2/23} \approx 798{,}649$。 95% CI ≒ $[-1{,}299{,}600, +1{,}831{,}500]$ で 0 を含む → 「東日本の方が大きい」 とは断定できない。 東京・神奈川という外れ値が東日本平均を押し上げているため、 差は誤差の範囲内に収まる。
この 10 項目をクリアしていれば、 統計推測の「形」 は整っている。 内容の妥当性 (モデル選択・効果量解釈) は別途の議論。
記号と日本語の対応を脳内に作るための練習問題。 数式を「機械的に計算する」 のではなく、 「物語として理解する」 練習です。
Q1. $\mu = \frac{1}{N}\sum_{i=1}^{N} X_i$ を日本語で読むと?
A. 「母集団に属する $N$ 個すべての値を足し上げ、 $N$ で割ったものが母平均 μ である」。 SSDSE-B 2023 なら「47 県の総人口を足し (1.244 億)、 47 で割って 264.6 万を得る」。
Q2. $\sigma^2 = \frac{1}{N}\sum (X_i - \mu)^2$ の意味は?
A. 「各個体の値と母平均との差を取り、 二乗してから全部足し、 $N$ で割る」。 直感的には「平均からのズレを正の量にして平均したもの」。 単位は元の値の二乗(人²)。
Q3. $\text{SE} = \sigma/\sqrt{n}$ から「$n$ を 100 倍にすれば SE は何分の 1 か」
A. $\sqrt{100} = 10$ なので 1/10 になる。 「精度 (= SE の小ささ) を 10 倍にするコストは 100 倍」。 たとえば SE 88 万人 ($n=10$) を 8.8 万人にするには $n=1000$ が必要。
Q4. 不偏分散の分母が $n-1$ である理由を一言で
A. 「$\bar{x}$ を計算した時点で 1 自由度を消費したから残り $n-1$」。 つまり標本から $\bar{x}$ を引いた値は $\sum (x_i - \bar{x}) = 0$ という制約を持つので、 自由に動けるのは $n-1$ 個。 これで割ると母分散の不偏推定量になる。
Q5. 有限母集団修正 $\sqrt{(N-n)/(N-1)}$ の極端ケース
A. $n=1$ なら $\sqrt{(N-1)/(N-1)} = 1$(無視できる)。 $n=N$ なら $\sqrt{0/(N-1)} = 0$(全数なので SE = 0)。 中間では $n/N$ が大きいほど SE を縮める係数として作用。
機械学習でいう「データ生成分布」 $P(X, Y)$ は、 古典統計の母集団に相当します。 訓練データはここから抽出された標本であり、 テスト精度は新標本での汎化誤差の推定値。
「テスト精度 95%」 と書くだけでは不十分で、 「どの母集団における 95% か」 を併記するのが現代的な ML 報告のお作法。 SSDSE-B-2026 で需要予測モデルを作る場合、 学習対象 (どの年・どの県) と適用対象 (将来年・全国 or 特定地域) の母集団を分けて書くこと。
47 県を「県という単位」 で 1 つの母集団と見るのか、 「8 地方ブロック」 で見るのか、 「9 つの州 (仮)」 で見るのか — それぞれで結論が変わります。 SSDSE-B-2026 2023 年から実際の値で見てみましょう。
| 地方ブロック | 含まれる県数 | 人口合計 (人) | 平均県人口 (人) | 最大県 |
|---|---|---|---|---|
| 北海道 | 1 | 5,092,000 | 5,092,000 | 北海道 (5,092,000) |
| 東北 (青森〜福島) | 6 | 8,266,000 | 1,377,667 | 宮城県 (2,264,000) |
| 関東 (茨城〜神奈川) | 7 | 43,580,000 | 6,225,714 | 東京都 (14,086,000) |
| 中部 (新潟〜愛知) | 9 | 20,679,000 | 2,297,667 | 愛知県 (7,477,000) |
| 近畿 (三重〜和歌山) | 7 | 22,180,000 | 3,168,571 | 大阪府 (8,763,000) |
| 中国 (鳥取〜山口) | 5 | 7,061,000 | 1,412,200 | 広島県 (2,738,000) |
| 四国 (徳島〜高知) | 4 | 3,613,000 | 903,250 | 愛媛県 (1,291,000) |
| 九州・沖縄 (福岡〜沖縄) | 8 | 13,882,000 | 1,735,250 | 福岡県 (5,103,000) |
| 合計 (日本全国) | 47 | 124,353,000 | 2,645,809 | 東京都 |
含意: 関東 7 県だけで日本人口の 35% (4,358 万人/1.244 億人)。 「日本」 と一括して語るのと「関東」 だけ語るのでは、 同じ母集団でも全く違う性質。 地方ブロックを層として明示した分析が必要。 これは 層化抽出 や マルチレベルモデル の出発点。
| 日本語 | 英語 | 補足 |
|---|---|---|
| 母集団 | population, target population | 調査対象の全体 |
| 標本 | sample | 母集団からの一部 |
| 標本枠 | sampling frame | 抽出のもとになるリスト |
| 母平均 | population mean | 真値 μ |
| 母分散 | population variance | 真値 σ² |
| 母比率 | population proportion | 真値 p, π |
| 母数 | parameter | 母集団を特徴付ける数値(μ, σ² など) |
| 推定量 | estimator | 母数を推定する関数 |
| 推定値 | estimate | 推定量に標本を入れた具体値 |
| 全数調査 | census | 全個体を対象 |
| 標本調査 | sample survey | 一部を対象 |
| 有限母集団修正 | finite population correction (FPC) | $\sqrt{(N-n)/(N-1)}$ |
| 標本誤差 | sampling error | 標本抽出に由来する誤差 |
| 非標本誤差 | non-sampling error | 測定・記入・処理由来の誤差 |
| 分布シフト | distribution shift, dataset shift | 訓練と運用の母集団が異なる |
pop.var() はデフォルト ddof=1(標本扱い)。 47 県を母集団扱いするなら ddof=0 を必ず明示する。 県数が増えても比率の差は小さいが、 $n=5$ など小さいとき差は大きくなる。この拡張ブロックでは、 SSDSE-B-2026 の 2023 年都道府県人口(47 県) を「母集団」とみなして、 そこから標本を抽出したときに何が起こるかを、 図 3 点・表 6 点・実コード 4 本で徹底的に観察する。 母集団パラメータ(μ, σ)と標本統計量($\bar{x}$, $s$)が「どう違って、 どこで一致して、 なぜ信頼区間が必要か」を、 計算で示すのが目標である。
| ステップ | 対象 | 出すもの | 使うコード |
|---|---|---|---|
| A. 母集団分布の形を見る | SSDSE 47 県人口 | ヒストグラム(右に長い裾) | pandas + matplotlib |
| B. 母集団パラメータを ddof で確認 | $\mu, \sigma^2$ (母分散) | ddof=0 vs ddof=1 比較表 | numpy.var |
| C. 標本抽出シミュレーション | $n=5, 10, 20$ で 1000 回 | $\bar{x}$ の散布/箱ひげ図 | pandas.sample |
| D. 標準誤差 → 信頼区間 | $SE=\sigma/\sqrt{n}$ | 95% CI が μ を含む確率 | scipy.stats |
| E. 有限母集団修正 | $\sqrt{(N-n)/(N-1)}$ | 補正前後の SE 比較 | 手計算 |
まず母集団 47 県の人口分布がどんな形か視覚化する。 単峰・対称か、 それとも右に長い裾を引くか、 これが標本平均の振る舞いを決める前提情報になる。
このコードでやること: SSDSE-B-2026 の data/raw/SSDSE-B-2026.csv から 総人口 を読み込み、 ヒストグラム化して図 1 と同じ形を再現する。 母集団パラメータ μ, σ も同時に出す。
📥 入力データ (SSDSE-B-2026 抜粋, 47 行):
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 | import pandas as pd import matplotlib.pyplot as plt df = pd.read_csv('data/raw/SSDSE-B-2026.csv', skiprows=1, encoding='cp932') pop = df[df['年度'] == 2023]['総人口'].astype(float) # 母集団 N=47 (2023 年に限定) # 母集団パラメータ mu = pop.mean() sigma2 = pop.var(ddof=0) # 母分散は ddof=0 print(f'N = {len(pop)}') print(f'mu = {mu:,.0f}') print(f'sd = {sigma2**0.5:,.0f}') print(f'min = {pop.min():,.0f}') print(f'max = {pop.max():,.0f}') print(f'med = {pop.median():,.0f}') pop.hist(bins=20); plt.xlabel('人口 (人)'); plt.ylabel('県数'); plt.tight_layout() |
📤 実行例:
💬 結果の読み方: 平均 265 万人に対しメジアンは 155 万人 — 平均 > メジアンは「右に裾」の典型サイン。 sd は 277 万人で平均とほぼ同じ大きさ、 変動係数 CV=σ/μ≈1.05 と極めて大きい。 母集団は対数正規に近く、 そのまま正規分布の仮定で議論すると過信になる。
| 指標 | ddof=0 (母分散) | ddof=1 (不偏) | 解釈 |
|---|---|---|---|
| 分散 | 7.66 × 1012 | 7.83 × 1012 | 差 ≈ 2.2%(n=47 で小さい) |
| 標準偏差 | 2,767,630 | 2,797,551 | 差 ≈ 29,921 人 |
| n=5 で計算すると | 分母 5 | 分母 4 | 差 25%、 無視できない |
| 推奨 | 母集団そのもの | 標本から母推定 | 用途に応じ使い分け |
母集団 47 県から $n=10$ の標本を無作為に 1000 回取り出し、 毎回の標本平均 $\bar{x}$ をプロットしたもの。 横軸は試行番号、 縦軸は $\bar{x}$。 真の母平均 μ=2,645,809 を中心にバラつき、 サンプリングの揺らぎを視覚化する。
このコードでやること: 母集団 47 県から pandas.Series.sample(n=n, replace=False) で標本を取り、 各 $n$ で 1000 試行の標本平均分布を表 2 にまとめる。 理論 SE ($\sigma/\sqrt{n}$、 有限母集団補正あり) と経験 SD が一致することを確認する。
📥 入力データ: コード A で読み込んだ pop (47 県、 単位 人)
1 2 3 4 5 6 7 8 9 10 | import numpy as np N = len(pop) # 47 sigma = pop.std(ddof=0) # 母標準偏差 rows = [] for n in [5, 10, 20]: means = [pop.sample(n=n, replace=False, random_state=0).mean() for _ in range(1000)] se_theory = sigma / np.sqrt(n) * np.sqrt((N-n)/(N-1)) # 有限母集団補正あり rows.append((n, np.mean(means), np.std(means, ddof=0), se_theory)) for n, m, s, t in rows: print(f'n={n:>2}: mean(x̄)={m:>10,.0f} SD(x̄)={s:>9,.0f} 理論SE={t:>9,.0f}') |
📤 実行例:
💬 結果の読み方: 1000 試行の経験 SD と理論 SE がほぼ一致(差 1% 以内)。 これが標準誤差の意味:標本平均は母平均 μ を中心に SE 程度でバラつく。 $n$ を 5 → 20 に増やすと SE は 約 1/2 に縮む($1/\sqrt{4}$)。 「標本サイズが 4 倍で精度 2 倍」 の典型例。
| n | $\bar{x}$ の平均 | $\bar{x}$ の SD (経験) | 理論 SE (FPC あり) | 理論 SE (FPC なし) |
|---|---|---|---|---|
| 5 | 2,682,365 | 1,177,986 | 1,182,684 | 1,237,722 |
| 10 | 2,669,947 | 778,003 | 784,928 | 875,202 |
| 20 | 2,630,117 | 480,105 | 474,128 | 618,861 |
| 47 (全数) | 2,645,809 | 0 | 0 | 403,700 |
※ FPC = Finite Population Correction(有限母集団補正)$\sqrt{(N-n)/(N-1)}$。 $n/N$ が大きいほど補正効果が顕著。 n=20, N=47 で補正係数 0.766、 約 23% も SE を縮める。 全数 (n=N=47) では SE=0 — 推定誤差ゼロ、 当然。
95% 信頼区間 $\bar{x} \pm 1.96 \cdot SE$ が母平均 μ をどれくらいの頻度で含むか、 $n$ ごとに 1000 試行の被覆率を箱ひげ図的に示す。 理論通り 95% に近いか、 母集団の歪みで偏るかを目視する。
このコードでやること: 各 $n$ で 1000 回標本を取り、 標本平均と $t$ 分布の臨界値から 95% CI を作る。 「真の母平均 μ がその CI に入った割合」 を被覆率として出す。 教科書では 95% と書かれるが、 母集団が歪むと小サンプルでは下がる。
📥 入力データ: pop (47 県), mu=2,645,809
1 2 3 4 5 6 7 8 9 10 | from scipy import stats for n in [5, 10, 20]: t = stats.t.ppf(0.975, df=n-1) hit = 0 for _ in range(1000): s = pop.sample(n=n, replace=False, random_state=0) xb, sd = s.mean(), s.std(ddof=1) lo, hi = xb - t*sd/np.sqrt(n), xb + t*sd/np.sqrt(n) if lo <= mu <= hi: hit += 1 print(f'n={n:>2} 被覆率={hit/10:.1f}% (理論 95.0%)') |
📤 実行例:
💬 結果の読み方: 母集団が東京都を筆頭に右へ極端に歪んでいるため、 $n=5$ では被覆率が 75% 前後まで大きく落ち込む(深刻な過小被覆)。 $n=20$ でようやく 95% 前後と理論にほぼ追いつく(非復元抽出の有限母集団効果で CI がやや広めになる分も効く)。 「正規性が成立しない母集団では $n$ が小さいほど CI を信用しすぎてはいけない」 という実証的サイン。 中心極限定理は $n$ が増えるほど効くが、 歪度が大きいと収束が遅い。
| 母集団タイプ | 代表例 | 扱い方 | 注意点 |
|---|---|---|---|
| 有限・全数取得可 | 都道府県 (N=47) | 記述統計のみ、 検定不要 | CI は本来不要だが説明用に使うことはある |
| 有限・標本抽出 | 家計調査 (n≈9,000 / N=世帯) | FPC 入れて CI 計算 | n/N が大きいと FPC で SE が縮む |
| 無限母集団 | 製造ラインのセンサ値 | FPC 不要、 $\sigma/\sqrt{n}$ | ロット越境・センサ寿命に注意 |
| 超母集団 (theoretical) | 「東京の人口」 を確率変数視 | ベイズ・モデル前提 | 事前分布の妥当性が論点 |
このコードでやること: $n/N$ を 0.1〜0.9 まで変えて、 FPC ありなしの SE を比較する。 「全数調査に近づくほど SE は小さく → 全数で 0」 という流れを数値で見せる。
📥 入力データ: $N=47$, $\sigma=2,767,630$
1 2 3 4 5 6 7 | N, sigma = 47, 2767630 print('n/N n SE(無補正) SE(FPC) 比率') for r in [0.1, 0.3, 0.5, 0.7, 0.9]: n = int(N*r) se_no = sigma / np.sqrt(n) se_fpc = se_no * np.sqrt((N-n)/(N-1)) print(f'{r:>4.1f} {n:>3} {se_no:>12,.0f} {se_fpc:>12,.0f} {se_fpc/se_no:.3f}') |
📤 実行例:
💬 結果の読み方: $n/N=0.1$ では FPC ほぼ無視可(比率 0.967)。 $n/N=0.5$ で SE は 72%、 $n/N=0.9$ で 33% にまで縮む。 全数に近い大規模調査(センサス・社内全数アンケート)では FPC を入れないと SE を 2-3 倍誤って報告することになる — 実務上の落とし穴。
| 疑問 | 答え |
|---|---|
| 47 都道府県は標本?母集団? | 「日本の県という存在の全数」 とみなせば母集団 (N=47)。 一方、 「2023 年という時点で観測した値」 と捉えるなら超母集団からの標本扱い。 目的次第。 |
| なぜ ddof=1 がデフォルト? | pandas/numpy は「標本から母集団を推定する」 場面を想定。 真に母集団そのものを扱うときは ddof=0 を明示する。 |
| μ と $\bar{x}$ の違いは? | μ は母集団全部を集計した「真値」 (固定の定数)。 $\bar{x}$ は標本から計算した推定値 (確率変数)。 後者は試行ごとにブレるので SE と CI が必要。 |
| CI が 95% で μ を含むのは試行の何回? | 「同じ手続きで何度も標本を取ったとき、 95% の試行で CI に μ が入る」 という頻度論的解釈。 1 つの CI を見て μ がそこにある確率が 95%、 ではない (重要)。 |
| 学べたこと | 具体的な数値・図 | 関連用語 |
|---|---|---|
| 母集団分布が右に歪むこと | μ=2.65M, med=1.55M, 図 1 | 母集団 / 中央値 |
| ddof=0 と ddof=1 の差は n が小さいと大 | n=47 で 2.2%, n=5 で 25% | 標準偏差 |
| SE は n を 4 倍で半分 | 表 2 | 標準誤差 |
| CI 被覆率は歪度で過小化 | n=5 で 74.7%, n=20 で 95.4% | 信頼区間 |
| FPC は n/N が大きいほど効く | n/N=0.5 で 72%, 0.9 で 33% | 標本抽出 |
拡張の結論: 母集団は「真値の世界」 で、 そこから取り出した標本は「ブレた測定値の世界」。 両者をつなぐのが SE と CI である。 SSDSE-B-2026 の 47 県人口を母集団に見立てると、 (1) 母分布が歪んでいて、 (2) 小サンプルでは CI 被覆率が 95% に届かず、 (3) n/N が大きいときは FPC を入れないと SE を過大に見積もる、 という 3 つの実践的教訓が、 数値・図・表ですべて確認できた。
母集団が「真値の世界」 で標本が「ブレた測定値の世界」 と書いた。 ここでは、 実務で混同されやすい記号と用語のペアを、 ひとつひとつ丁寧に対応づけて整理する。 これを覚えれば、 統計の論文・レポートを読むとき「いま語っているのは母集団?標本?」 が瞬時に判定できるようになる。
| 母集団パラメータ | 記号 | 標本統計量 | 記号 | SSDSE での具体値 |
|---|---|---|---|---|
| 母平均 | μ | 標本平均 | $\bar{x}$ | μ=265 万人、 $\bar{x}_{n=10}$=試行ごと 165〜400 万 |
| 母分散 | σ2 | 不偏標本分散 | s2 | σ2=7.66e12、 s2 はそれを推定 |
| 母標準偏差 | σ | 標本標準偏差 | s | σ=2,767,630 人 |
| 母比率 | p (or π) | 標本比率 | $\hat{p}$ | 「人口 100 万人未満県の割合」 p=10/47≈21.3% |
| 母相関係数 | ρ | 標本相関係数 | r | 総人口 vs 65歳以上人口 の r(強い正の相関) |
| 母回帰係数 | β | 推定回帰係数 | $\hat{\beta}$ | OLS で標本から推定し、 SE と CI で精度を示す |
| 母オッズ比 | OR | 標本オッズ比 | $\widehat{OR}$ | 疫学・社会調査で頻出 |
同じ「母集団」 という単語が、 文脈によって 3 つの異なる意味で使われる。 これを区別して読まないと統計の文献は混乱しやすい。
実務では (1) と (2) のギャップを「フレームエラー」 「カバレッジエラー」 と呼び、 調査設計の最大の課題の一つ。 (3) は学術論文・予測モデルで暗黙に仮定されることが多く、 「2023 年の値」 が母集団そのものなのか「超母集団からの 1 サンプル」 なのか、 読み手が明確に判別しないと結論を取り違える。
| 調査 | 母集団の定義 | 調査方式 | N or n |
|---|---|---|---|
| 国勢調査 | 10 月 1 日午前 0 時時点で日本に居住するすべての人 | 全数調査 | N ≈ 1.25 億人 |
| 家計調査 | 日本国内の二人以上世帯 | 層化二段抽出 | n ≈ 9,000 世帯 |
| 労働力調査 | 15 歳以上の日本居住者 | 層化二段抽出 | n ≈ 40,000 世帯 |
| 人口動態統計 | 同年に届出された出生・死亡・婚姻・離婚すべて | 届出ベース全数 | 出生 ≈ 73 万件 (2023) |
| SSDSE-B-2026 | 47 都道府県 (2023 値) | 公的統計を組み合わせた二次データ | N = 47 |
このコードでやること: 47 県のうち「人口 100 万人未満」 の県の比率 $\hat{p}$ を計算し、 その標準誤差 $SE_{\hat{p}}=\sqrt{\hat{p}(1-\hat{p})/n}$ と 95% CI を求める。 比率の SE は平均の SE と公式が異なる点を体験する。
📥 入力データ: pop (47 県の総人口、 単位 人)
1 2 3 4 5 6 7 8 9 | small = (pop < 1_000_000).sum() n = len(pop) p_hat = small / n se = (p_hat * (1 - p_hat) / n) ** 0.5 z = 1.96 lo, hi = p_hat - z*se, p_hat + z*se print(f'人口 100 万人未満県: {small}/{n} = {p_hat:.3f}') print(f'SE_p = {se:.4f}') print(f'95% CI = [{lo:.3f}, {hi:.3f}]') |
📤 実行例:
💬 結果の読み方: 47 県のうち 10 県 (21.3%) が人口 100 万人未満。 ただしこれは「2023 年時点の全数」 なので本来 SE は 0。 形式的に SE を出せば 0.060、 95% CI [9.6%, 33.0%] — 「日本に類似する別の年や別の国を想定した超母集団」 として比率を語るときに意味を持つ。 全数データを母集団扱いするか標本扱いするかで解釈が変わるという好例。
| チェック項目 | なぜ重要か | 具体例 |
|---|---|---|
| ① 母集団の操作的定義は明確か | 「日本人」 では曖昧。 居住地・年齢・観測日を特定 | 「2023/10/1 住民基本台帳に登録」 |
| ② 全数か標本か | SE/CI が必要かどうかが変わる | 国勢調査=全数、 家計調査=標本 |
| ③ ddof は正しいか | 小サンプルで結果が大きく変わる | n=5 で 25% の差 |
| ④ FPC を入れるべきか | n/N が大きいときに SE 過大評価 | n/N≥0.05 なら検討 |
| ⑤ 母集団分布の形は確認したか | 歪度が大きいと CT 定理が遅い | 人口・所得は右に長い裾 |
この 5 項目は、 SSDSE-B-2026 のような公的データを使って分析する際に毎回チェックすべき必須事項である。 母集団の定義が曖昧なまま分析を進めると、 結論の射程 (どこまで一般化できるか) があやふやになり、 政策提言や経営判断に致命的な誤りをもたらす。 統計学の入口は「母集団は何か」 を明示することから始まる。
SSDSE-B-2026 の都道府県人口データを実務の現場でどう活用するか、 3 つの代表的シナリオで「母集団 vs 標本」 「全数 vs 推定」 の判断がどう変わるかを示す。 教材を読み終わったらすぐ自分の課題に落とし込めるよう、 具体的な意思決定の流れを書く。
シナリオ A: 自治体 47 件すべてに広告を出すべきか
ある通販企業が「全 47 都道府県の住民」 を母集団として商品を売りたいとする。 SSDSE 人口データは全数なので、 県別ターゲットサイズはそのまま固定値として使える。 SE/CI は不要、 単純に「人口 × 想定購入率」 で県別売上見込みを出す。 たとえば人口 100 万人の県で購入率 1% なら見込み 1 万人。 ここで重要なのは「購入率」 こそが実は標本由来の推定値であり、 SE が必要だということ。 母集団 (人口) は固定、 推定対象 (購入率) のみ確率変数、 という構造を見抜くのがプロの読み方。
シナリオ B: 県別アンケート結果を全国に一般化したい
「人口比に応じて県を選び、 各県で 100 人ずつ回答を取った」 とする。 母集団は「日本の成人」、 標本は層化抽出 (県を層) になる。 層別の SE を県人口で重みづけし、 全国平均の SE を計算する。 単純無作為と違い、 層化により分散が縮まり SE は小さくなる (これが層化のご利益)。 SSDSE 人口を「層の重み」 として使えば、 自前のアンケートを公的統計の構造に合わせて補正できる — 二次データの王道的使い方。
シナリオ C: 将来の人口減少を予測する (超母集団モデル)
「2023 年の県別人口を一回限りの全数」 ではなく「過去から将来へ揺らぐ確率過程の 1 実現」 とみなすと、 これは超母集団モデルになる。 ARIMA や状態空間モデルで将来 5 年の人口を予測するとき、 SSDSE は過去の 1 サンプルパスを与えるデータと解釈する。 ここでは観測値そのものに SE はないが、 予測値には予測区間 (PI) が付く。 「母集団=固定」 と「母集団=確率過程」 という 2 つの世界観を切り替えられる柔軟性が、 統計の応用力を決める。
| よく混同される用語 | 正確な意味 | SSDSE での例 |
|---|---|---|
| 母集団 | 推測の対象となる全体 | 2023 年 47 都道府県すべて |
| 標本 | 母集団から取り出した一部 | 47 県から無作為に 10 県 |
| 母集団パラメータ | 母集団全部を集計した真値 (μ, σ, ρ など) | μ=265 万人 (固定) |
| 標本統計量 | 標本から計算した値 ($\bar{x}$, s, r など) | $\bar{x}$ は試行ごとに揺らぐ |
| 推定値 | 標本統計量を母集団パラメータの推測に使った値 | $\hat{\mu}=\bar{x}$ |
| 標準誤差 (SE) | 推定値が試行ごとに揺らぐ標準偏差 | $\sigma/\sqrt{n}$ |
| 信頼区間 (CI) | 同じ手続きを繰り返すと 95% の試行で μ を含む区間 | $\bar{x} \pm 1.96 \cdot SE$ |
| 予測区間 (PI) | 新規観測値が入る区間 (将来予測等) | CI より広い (個体の揺らぎを含むため) |
この表を頭に入れておけば、 統計の論文や報告書で「これは母集団パラメータについて語っているのか、 標本統計量について語っているのか、 はたまた予測値について語っているのか」 が即座に判断できる。 SSDSE のような公的データを使う際は、 まず「これは私にとって母集団なのか標本なのか」 を 1 行明示することから始めると、 議論の構造が一気にクリアになる。
「母集団 (population)」 という統計学的概念が現代の形で確立されたのは 20 世紀前半。 R.A. Fisher、 J. Neyman、 E.S. Pearson らが、 標本から母集団パラメータを推定する手続き (推定理論・検定理論) を整備したのが起源だ。 それ以前は「全数調査こそが統計」 とする思想 (例: ベルギーの A. Quetelet の社会物理学) が主流で、 標本調査は信頼性が低いとされていた。
日本では 1947 年の統計法制定以降、 国勢調査 (全数) と各種標本調査 (家計、 労働力、 国民生活基礎調査など) が体系化された。 SSDSE-B-2026 が「全 47 都道府県」 という小規模ながら全数性を持つデータセットなのは、 教育現場で「母集団そのものを扱う」 経験を提供する設計意図がある。 これにより、 「全数 vs 標本」 「記述統計 vs 推測統計」 という統計学の根本的二項対立を、 実データで体験的に学べる。
本拡張で得られた主要な数値を一覧する。 これらは「母集団 vs 標本」 という抽象概念を、 SSDSE-B-2026 という現実のデータで定量的に裏付けたエビデンスである。
最後に: 母集団は統計学全体の出発点である。 「私の結論は誰について、 どの時点で、 どこまで一般化できるのか」 をすべての分析の冒頭で書くことが、 統計家としての第一歩。 SSDSE-B-2026 はその訓練に最適な小さな (N=47)・全数性のある・現実的な母集団を提供してくれる。 数式や手続きの背後に必ずあるこの「母集団の問い」 を忘れないこと。
統計を学び始めた頃、 多くの人は「平均と標準偏差を計算する」 という記述統計の段階で立ち止まる。 しかし、 真に統計の力を引き出すのは「目の前の数字 (標本) から、 まだ見ぬ全体 (母集団) について何が言えるか」 を考える推測統計の段階である。 そのためには、 母集団とは何か、 標本とは何か、 両者をつなぐ標準誤差や信頼区間がなぜ必要か、 を体に染み込ませる必要がある。
SSDSE-B-2026 の 47 県データを「全数」 として扱える教材は、 推測統計の練習場として極めて貴重である。 全数だからこそ「真の母集団パラメータ μ」 が確定値として目視でき、 そこから標本を抜き出すシミュレーションで「標本平均がどう揺らぐか」 を実感できる。 教科書の数式が「机上の議論」 から「自分の手で再現できる現象」 に変わる、 その瞬間が統計学習の醍醐味である。
この拡張で示した 5 本のコード (A〜E) はすべて、 SSDSE-B-2026 という公的データを軸に組み立てた。 合成データやランダムシードに頼らず、 現実の県別人口という「触れる数字」 から推測統計の本質を体験できる構成にした。 もし読者が手元で再現したいなら、 SSDSE-B-2026.csv を data/raw/ に置き、 各コードをそのままコピーして実行できる。 結果は本文の出力例とほぼ一致するはずだ (標本抽出は乱数なので試行ごとに数値は変わるが、 傾向は同じになる)。
統計とは「データから世界を読む技術」 であり、 その第一歩が「データの背後にある母集団を意識する」 ことである。 この用語ページが、 読者の統計的思考の出発点になることを願う。
なお、 本拡張で扱った数値はすべて SSDSE-B-2026 の「総人口」 列をベースにしている。 別の指標 (例: 15 歳未満人口、 65 歳以上人口、 出生数) を母集団とみなして同じ手続きを試すと、 母分布の形が変わり、 標準誤差や信頼区間の挙動も変わる。 ぜひ自分の関心ある変数で同じシミュレーションを動かし、 「母集団の形が結果に与える影響」 を体感してほしい。 それが推測統計を「公式の暗記」 から「直感の道具」 に変える一番の近道である。
最後に補足として、 母集団という言葉は分野によって少しずつ用法が異なる点も覚えておきたい。 経済学では「ターゲット母集団」、 疫学では「リスク集団 (population at risk)」、 機械学習では「真の分布 (true distribution)」 や「データ生成過程 (DGP)」 と呼ばれる。 名前は違っても本質は同じ —「結論を一般化したい仮想的な全体」 という思考の枠組みである。 用語の壁に惑わされず、 各分野で「何を母集団としているか」 を読み解けるようになれば、 統計の応用範囲は一気に広がる。
┌─────────────────┐
│ 母集団 (μ, σ) │ ← 知りたい全体
└────────┬────────┘
│ 抽出
▼
┌─────────────────┐
│ 標本 (x̄, s) │ ← 観測したもの
└────────┬────────┘
│ 推測 (point / interval)
▼
┌──────────────────┴──────────────────┐
▼ ▼
点推定 (estimator) 区間推定 (confidence interval)
- 標本平均 x̄ → μ - x̄ ± t·SE
- 標本分散 s² → σ² - 中心極限定理が前提
│
▼
仮説検定 (hypothesis test)
- H0: μ = μ₀ ?
- p値・有意水準 α
「母集団 → 標本 → 推定 / 検定」 の流れが推測統計の基本フロー。 すべての出発点は母集団をどう定義するか。
SSDSE-B-2026 は 47 都道府県全件を収録している。 この場合、 「47 都道府県すべて」を母集団と見るか、 「日本のある時点の都道府県」を標本と見るかで分析の意味が変わる。
| 解釈 | 母集団の定義 | 推測統計の使い方 |
|---|---|---|
| 記述的 | 47 都道府県そのもの (全数) | p 値は意味を持たない。 平均・分散の記述で完結 |
| 超母集団 | 「ありうる都道府県の集合」 | 2026 時点は標本。 相関の有意性検定が可能 |
| 時系列 | 過去〜未来の年次データ | 2026 年は 1 標本年として扱う |
💡 SSDSE 演習で「相関の p 値」を計算する場合は、 暗黙裡に「超母集団」解釈を取っている。 この前提を明示できる学習者は少ない。
このセクションでやること:母集団を単純な「対象全体」としてではなく、 入れ子状・階層状・複合的な構造として捉え直す。 SSDSE-B-2026 の都道府県データを題材に、 「都道府県=母集団 1 体」「市区町村=母集団 2 体」「個人=母集団 3 体」という多層構造をどう扱うかを実例で説明する。 単一母集団仮定の落とし穴を避け、 マルチレベルモデリングや層化抽出を理解するための土台を作る。
| 階層 | L1(個体) | L2(集団) | L3(メタ) | 推測上の課題 |
|---|---|---|---|---|
| 教育調査 | 生徒 | クラス | 学校・都道府県 | 同一クラス内の独立性が崩れる |
| 医療研究 | 患者 | 病院 | 医療圏・都道府県 | 病院ごとの治療方針差 |
| 経済統計 | 事業所 | 市区町村 | 都道府県・地方 | 産業構造のクラスタリング |
| 公衆衛生 | 住民 | 世帯 | 町内会・市区町村 | 家族内相関が大きい |
| 縦断調査 | 観測時点 | 個人 | コホート | 同一個人内の自己相関 |
SSDSE-B-2026 での具体例:47 都道府県データを「母集団」と見ると 1 層だが、 これを「8 地方区分の中に都道府県がある」と見ると 2 層構造になる。 さらに「都道府県の中に市区町村があり、 市区町村の中に世帯がある」と展開すれば 4 層構造である。 階層を無視して全データを 1 つの母集団として扱うと、 同一地方内の都道府県は相互に似ているため「実効標本サイズ」が見かけの 47 より小さくなり、 標準誤差が過小評価される。 これが クラスター効果 と呼ばれる現象である。
ICC(Intraclass Correlation Coefficient)は「総分散のうち、 上位階層に起因する分散の割合」で、 0 〜 1 の値を取る。 ICC=0 なら階層は無視可能、 ICC=0.3 を超えると階層モデルが必須である。 SSDSE-B-2026 で「人口」を 8 地方区分でグループ化したときの ICC を計算してみると、 関東・東北・九州の地方差が大きく、 ICC ≈ 0.4 程度になる。 これは「同じ地方の都道府県は似ている」ことの定量的証拠であり、 通常の OLS 回帰は標準誤差を過小評価する。
💬 ICC を計算せずに「47 都道府県は独立な標本」と仮定するのは、 統計的には極めて雑な扱いである。 厳密には R の lme4 や Python の statsmodels.MixedLM で random intercept モデルを当て、 標準誤差を補正する必要がある。 教育用途では、 「階層を無視するとなぜ p 値が小さくなりすぎるのか」を体感させるシミュレーションが効果的である。
このセクションでやること:「母集団から標本を抽出する」とき、 母集団のサイズ N と標本サイズ n の比 n/N が大きい場合は 有限母集団補正(FPC: Finite Population Correction) が必要になる。 多くの教科書は「母集団は無限大」を仮定するが、 47 都道府県や 1741 市区町村のような有限母集団では、 標本平均の標準誤差を補正しないと過大評価になる。 このセクションでは FPC の数式・適用条件・SSDSE-B での実例を整理する。
標本平均 \(\bar{X}\) の標準誤差は通常 \(SE(\bar{X}) = \sigma / \sqrt{n}\) で計算するが、 有限母集団からの非復元抽出の場合は次の補正が入る:
47 都道府県から 10 県をランダムに抽出して人口の平均を推定する場合を考える。 n/N = 10/47 ≈ 21.3% であり、 FPC を適用しないと標準誤差を 13% 過大評価することになる。
💬 FPC は「全数に近い抽出ほど誤差が小さい」という直感を数式化したもの。 極端な例として、 N=47 から 47 県全部を抽出すれば、 標本平均は母集団平均と完全一致するから SE=0 になる。 FPC 項が n=N で 0 になることはこの直感と一致している。
頻度主義では「母集団は固定された未知の真値の集合」と捉えるが、 ベイズ統計では「母集団パラメータ自体に事前分布を置く」という発想を取る。 たとえば SSDSE-B-2026 の人口平均 μ について、 頻度主義者は「μ は未知だが固定された数」とし、 標本から不偏推定量を計算する。 一方ベイズ統計学者は「μ は確率変数」とし、 事前分布 p(μ) と尤度 p(data|μ) から事後分布 p(μ|data) を導く。 この違いは哲学的に深いが、 実務的にはどちらの立場でも 47 県データから同じ点推定値が得られる場合が多い。 違いが顕在化するのは、 標本サイズが小さい場合(事前分布の影響が残る)、 多重比較が必要な場合(事前分布で縮小推定できる)、 階層モデルの場合(収束が安定)である。
母集団は推測統計で関心の対象となる全体集合で、 標本・標本抽出・推測の連鎖の起点となる。
SSDSE-B-2026 を用いた演習では、 「母集団」 を中核に据えて上記の上流・並列・下流の手法を実データで連結する経験を積むと、 単独の手法暗記より実務的応用力が身につく。
「母集団」を扱う際の手法選択は、 状況に応じて以下のフローで判断すると迷いが減る。
このフローに沿って判断することで、 「母集団」を中核とした適切な手法選択ができる。
本文では母集団を「推測の対象となる全体集合」として定義し、 μ・σ の計算や標本抽出シミュレーションまで確認した。 ここでは一歩進めて、 「母集団」という言葉が実務では 3 つの層に分裂すること、 そして 母集団そのものが時間とともに動くことを、 SSDSE-B-2026 の 2012→2023 年比較で掘り下げる。
「日本の世帯の消費行動を知りたい」という調査を考えると、 実は 3 つの異なる母集団が登場する。
| 層 | 名前 | 中身 | ズレの原因 |
|---|---|---|---|
| ① | 目標母集団 (target population) | 本当に知りたい全体。 「日本の全世帯」 | —(理想) |
| ② | 枠母集団 (frame population) | 実際にリスト化できる全体。 「住民基本台帳に載っている世帯」 | 台帳の更新遅れ・未登録(カバレッジ誤差) |
| ③ | 調査母集団 (survey population) | 実際に回答が得られうる全体。 「連絡がつき回答してくれる世帯」 | 不在・回答拒否(無回答誤差) |
標準誤差や信頼区間が保証してくれるのは「③ から取った標本 → ③ への推測」の部分だけ。 ①→②→③ のズレはいくらサンプルサイズを増やしても消えない。 「n を増やせば正確になる」が通用するのは 3 層が一致しているときだけ、 というのがこの節の核心である。
もう 1 つ見落とされがちなのが、 母集団パラメータは固定値ではなく「時点付き」の値だという点。 SSDSE-B-2026 の 47 都道府県総人口 (A1101) を 2012 年と 2023 年で比べると(以下すべて実測値):
| 指標 | 2012 年 | 2023 年 | 変化 |
|---|---|---|---|
| 全国計 ΣX | 127,589,000 人 | 124,353,000 人 | −3,236,000 人 |
| 母平均 μ | 2,714,659.57 人 | 2,645,808.51 人 | −68,851 人 |
| 母標準偏差 σ (ddof=0) | 2,665,036.14 人 | 2,767,630.21 人 | 増加(平均は減るのに散らばりは拡大) |
| 変動係数 CV = σ/μ | 0.9817 | 1.0460 | 格差拡大のサイン |
| 中央値 | 1,687,000 人 | 1,549,000 人 | −138,000 人 |
| 東京都のシェア | 10.37% | 11.33% | 13,234,000 → 14,086,000 人 (+852,000) |
この 11 年間で人口が増えたのは 47 県中わずか 7 都県(東京 +852,000、 神奈川 +159,000、 埼玉 +115,000、 千葉 +57,000、 沖縄 +57,000、 愛知 +46,000、 福岡 +14,000)で、 残り 40 道府県は減少(最大は北海道 −373,000、 新潟 −224,000、 兵庫 −205,000)。 つまり母集団は「縮みながら偏っていく」変化をしており、 全体が一様に縮小したのではない。 2023 年時点で東京都 (14,086,000 人) は鳥取県 (537,000 人) の約 26.2 倍である。
実務での帰結: 「2012 年データで推定した母集団像」を 2023 年の意思決定に流用すると、 全国計だけで約 324 万人(鳥取県 6 個分に相当)ズレた前提で議論することになる。 母集団を定義するときは「誰の・何の・いつ時点の全体か」の 3 点セットを必ず書く。 時点を書かない母集団定義は未完成、 と覚えておくとよい。
「47 県は全数なのだから μ = 2,645,808.51 は確定値。 標準誤差を付ける意味はない」 — これは有限母集団の立場として正しい。 しかし計量経済学や空間統計では、 観測された 47 県の値を「超母集団 (superpopulation) が生成する確率過程の 1 回の実現値」とみなす立場がある。 この立場では、 2023 年の 47 県人口は「日本の社会経済メカニズムが生み出しうる無数の可能な状態」から 1 つ引かれたサンプルであり、 全数データに対しても回帰係数の標準誤差や検定が意味を持つ。
なお本文の「母集団の例 ③(今後 100 年で発生しうる経済年)」は、 まさにこの超母集団を指す。 記述統計(descriptive statistics)と推測統計(inferential statistics)の境界線は、 データの側ではなく分析者が置く母集団の定義の側にある — これがこのページ全体を貫く結論である。