母平均の点推定・区間推定

# 母集団と標本
統計学において、調査対象となる数値・属性等の源泉となる集合全体を**母集団**と言い、母集団から抽出された部分集合のことを**標本**と言います。

統計学において、基本統計量を示す・グラフを作成するなど**「データの特徴を記述する」**学問を**記述統計学**と言い、**「標本から母集団の特徴を推測する」**学問を**推測統計学**と言います。
記述統計学では、大量の要素の標本の存在を前提とするため、母集団と標本はほぼ同一視されます。

![c8856789ec11ab8b1013037cef6929f9-13.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/663460/1cb98db0-7826-5b05-30b6-8e9fc49ac08b.png)

> 画像引用:[16-1. 母集団と標本 | 統計学の時間 | 統計WEB](https://bellcurve.jp/statistics/course/8003.html)

母集団のうち、データの個数が有限の要素から成る母集団を**有限母集団**と言い、データの個数が無限の要素からなる母集団を**無限母集団**と言います。日本に住む女性のデータは非常に多いものの有限であるため有限母集団であり、サイコロを投げて出る目のデータは試行を無限に繰り返すことが可能であるため無限母集団であると言えます。

# 標本抽出
調査を実施する際は、調査対象となる母集団を設定した上でデータを所得します。このとき調査対象となる母集団全てを調べることを**全数調査**と言い、母集団の一部を抽出して調べることを**標本調査**と言います。

標本調査には、母集団を代表する標本を無作為に抽出した**無作為抽出（ランダム・サンプリング）**と、母集団を代表する標本を意図的に抽出した有意抽出があります。標本の抽出には、一度抽出された標本を母集団の中に戻す抽出方法である**「復元抽出法」**と、一度抽出された標本を母集団の中に戻さない抽出方法である**「非復元抽出法」**の2つがあります。

### 単純無作為抽出法
母集団の全ての要素を等確率で抽出し、調査対象とする方法を**「単純無作為抽出法」**と言います。
母集団から調査対象を完全に無作為に抽出するのは、非常に手間と時間がかかるのが問題です。

### 層化抽出法
母集団をその特性に応じていくつかの層に分類することが可能な時、**母集団を層化し、各層からランダムに標本を抽出する**方法を**「層化抽出法」**と言います。特に、層の大きさに比例させて調査対象を抽出する方法を「比例配分法」と言います。

世論調査では、都道府県別・自治体別などに分けてサンプリングする層化抽出法が採られています。母集団内情報の比較が行える、母集団を推測する精度が増すなどのメリットがある一方、層化を行うための母集団の構成情報を事前に知っておく必要があるというデメリットが存在します。

![795316b92fc766b0181f6fef074f03fa-16.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/663460/5c02c743-4f4b-34cf-3410-ce3768701985.png)

> 画像引用： [16-3. 標本の抽出方法 | 統計学の時間 | 統計WEB](https://bellcurve.jp/statistics/course/8007.html)

### 系統抽出法
抽出枠の**「先頭からm番目」の要素を開始点として、そこから「nつ飛ばし」に要素をサンプリング**していく方法を**「系統抽出法」**と言います（m,nは任意の数）。

単純無作為抽出法に比べて手間や時間がかからないなどのメリットがある一方、名簿の並び順に何らかの周期があると標本に偏りが生じる可能性があるというデメリットが存在します。

![3a4f695a458cb0ac0aceaa2eb13ac2dd-10.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/663460/8bfb0519-4865-89b1-9ba2-6c296eabf14b.png)

> 画像引用： [16-3. 標本の抽出方法 | 統計学の時間 | 統計WEB](https://bellcurve.jp/statistics/course/8007.html)

### 集落抽出法（クラスター抽出法）
**母集団をいくつかのクラスター（似た性質を持つ集団）に分け、その中からいくつかのクラスターを無作為抽出し、それぞれのクラスターで全数調査を実施する**という3つの段階をとる方法を**「集落抽出法」（クラスター抽出法）**と言います。

クラスターの情報さえあれば実施することができるため手間や時間がかからないなどのメリットがある一方、同じクラスターに属する調査対象は似た性質を持つため標本に偏りが生じる可能性があるというデメリットが存在します。

![2b530e80c7d0de90885e285c5d798063-14.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/663460/c997f49d-46f5-a3e7-2c64-c0df90846c06.png)

> 画像引用： [16-3. 標本の抽出方法 | 統計学の時間 | 統計WEB](https://bellcurve.jp/statistics/course/8007.html)

### 多段抽出法
母集団をいくつかのグループに分け、そこから下部のグループを無作為抽出し、その中からさらにグループを無作為抽出する…という作業を何度か繰り返し、最後に抽出されたグループから調査対象を無作為抽出するという方法を**「多段抽出法」**と言います。

抽出効率が高いというメリットがある一方、サンプルサイズが小さい場合に標本に偏りが生じる可能性があるというデメリットが存在します。

![c8856789ec11ab8b1013037cef6929f9-15.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/663460/9f0d94ed-3843-37bd-8f64-626aee6b006b.png)

> 画像引用： [16-3. 標本の抽出方法 | 統計学の時間 | 統計WEB](https://bellcurve.jp/statistics/course/8007.html)

# 母平均の点推定
「標本から母集団の特徴を推測する」学問である推測統計学には、推測統計学で求めたいとする**母集団を特徴づける値である「母数」**を統計学的に推測する「推定」と母集団から抽出された標本の統計量に関する仮説が正しいかを統計学的に判定する「検定」の2つがあります。

推定には、標本から求められるただ1つの値によって**母数をピンポイントで推定する**という手法である**「点推定」**と、標本から**母数が含まれるであろう区間を推定する**手法である**「区間推定」**があります。

![16e3bcb46ed7e460f5193e39ee3c7206.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/663460/8959b672-0786-b017-6fec-9f9f3967a8ac.png)

> 画像引用：[統計的推定とは？～点推定と区間推定の違い～ | AVILEN AI Trend](https://ai-trend.jp/basic-study/basic/statistical-estimation/)

※**「母数」**はパラメーターの一種であり、母集団を特徴づける母平均・母分散などが該当します

母数（パラメーター）を推定するために利用する数値の計算方法や計算式を**「推定量」**と言い、実際に試行を行った結果から計算した値を**「推定値」**と言います。

点推定では、「nが大きいとき標本平均は真の平均に近づく」、つまり「標本平均の期待値（平均値）は母平均に一致する」という**大数の法則**から**標本平均を母平均と見なす**ことができます。
また「一致性」と「不偏性」という点推定に関する2つの性質から、標本平均を母平均と見なすこともできます（「大数の法則」≒「一致性」であると言える）。

### 一致性
**「サンプルサイズnが無限大になる時、推定量が母数の真の値に収束する」**という性質を**「一致性」**と言います。
母数θの推定量をθ^とすると次のように表すことができます。

```math
∀_ε>0　n→∞　P(|\hat{θ}-θ|>ε)=0
```
このような推定量θを**一致推定量**と言います。

### 不偏性
**「推定量の期待値が母数（パラメーター）に一致する」**という性質を**「不偏性」**と言います。
母数θの推定量をθ^とすると次のように表すことができます。

```math
E(\hat{θ})=θ
```

このような推定量θを**不偏推定量**と言います。

### 標本分散と不偏分散
偏差平方和をnで割る**標本分散**は、**一致推定量であるものの不偏推定量ではありません**。
nが十分に大きくない場合、標本分散と母分散は一致せず「標準分散＜母分散」となります。

標本分散の期待値が母分散と一致する（不偏性を持つ）ように、標本分散の算出式にn/(n-1)を掛けたものが不偏分散であり、**不偏分散は一致性と不偏性を持つ**ため不偏分散を母分散と見なすことができます。

### 標準誤差
データそのものの散らばりを示す標準偏差とは異なり、平均の散らばりを示し**「標本平均の標準偏差」**である理論的な推定値を**標準誤差**と言います。

「平均μ、分散σ^2を持つあらゆる確率分布から無作為復元抽出した標本平均Xnの分布は、nが十分に大きい時には正規分布 N(μ, σ^2/n) へ近づいていく」という中心極限定理より、標本平均の標準偏差である標準誤差は以下のように表されます。

```math
標準誤差SE=\sqrt{\frac{σ^2}{n}}=\frac{σ}{\sqrt{n}}
```

【点推定のまとめ】

* 母平均の推定量▶︎標本平均を用いる（普遍性・一致性）
* 母分散の推定量▶︎不偏分散を用いる
* 母標準偏差の推定量▶︎標準誤差を用いる（標準偏差をサンプルサイズnの平方根で割る）

# 母平均の区間推定
標本から**母数が含まれるであろう区間を推定する**手法を**「区間推定」**と言います。
この区間のことを**信頼区間（CI）**といい、その両端を**信頼限界**と呼びます。また、信頼限界の小さい端を**下側信頼限界**と言い、大きい端を**上側信頼限界**と言います。

**「母数がある信頼区間に含まれる確率」**を**信頼係数（信頼度）**と言い、一般には90%・95%・99%がよく用いられます。信頼係数が高いほど誤りを犯す確率は小さくなりますが、信頼区間は大きくなります。またサンプルサイズnが大きくなるほど「推定の精度」は高くなるため、信頼区間の幅が狭くなります。

信頼係数が95%であるというのは、標本平均を母平均と見なす場合、**「無作為抽出を繰り返して区間推定を100回行った時、95回は母平均が信頼区間に入るが、5回は信頼区間に入らない可能性がある」**ということです。
「母平均が95%の確率で推定した信頼区間に入る」というのは間違いになります。

![795316b92fc766b0181f6fef074f03fa-1.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/663460/a15b426a-4b87-7132-2580-ed01879f4f55.png)

> 画像引用：[19-3. 95％信頼区間のもつ意味 | 統計学の時間 | 統計WEB](https://bellcurve.jp/statistics/course/8891.html)

母平均の区間推定には、**「母分散既知」**と**「母分散未知」**の2つの場合があります。

* 「母分散既知」▶︎**母分散σ^2**の値を使い、**標準正規分布**を用いて信頼区間を算出する
* 「母分散未知」▶︎**不偏分散s^2**の値を使い、**t分布**を用いて信頼区間を算出する

実際は母平均が分からないのに母分散は分かっているという状況はほとんどないため、後者の「母分散未知」であるt分布を使用した算出法がよく用いられます。

### 母分散既知
母分散既知の母平均の区間推定には、標本平均の分布を正規分布 N(*μ*, σ^2/n) に従わせるという概要の「中心極限定理」を用い、**標本平均を標準化**することで**標準正規分布**を利用します。
標本平均を**標準化した統計量Z**は以下のように表せ、標準正規分布に従います。

```math
Z=\frac{\bar{x}-\mu}{標本平均の標準偏差（標準誤差）}=\frac{\bar{x}-\mu}{\frac{σ}{\sqrt{n}}}
```

信頼区間が95%である場合、統計量Zが標準正規分布の95%の面積（確率）の範囲内にあればいいので、下側確率・上側確率が共に2.5%となるZの値を標準正規分布表を使用して求めます。

![795316b92fc766b0181f6fef074f03fa.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/663460/1ee45a53-9cd0-1a72-4feb-ff781b8d9f57.png)

> 画像引用：[19-2. 母平均の信頼区間の求め方（母分散既知） | 統計学の時間 | 統計WEB](https://bellcurve.jp/statistics/course/8888.html)

標準正規分布表より、Z=1.96であることが分かるため以下の式が成り立ちます。

```math
-1.96 \leq \frac{\bar{x}-\mu}{\frac{σ}{\sqrt{n}}} \leq 1.96
```

求める母平均*μ*について変形すると、95%信頼区間の母分散既知の母平均の区間推定は**標準誤差**を用いて、以下のように表されます。

```math
\bar{x}-1.96・\frac{σ}{\sqrt{n}} \leq \mu \leq \bar{x}+1.96・\frac{σ}{\sqrt{n}}
```

一般化して**信頼係数100(1-α)%**の場合、母平均の区間推定式は次のように表せます。
**Z(α/2)は下側確率（上側確率）**を表し、αの値から**標準正規分布表**を用いてZの値を求めます。

```math
\bar{x}-Z_{(a/2)}・\frac{σ}{\sqrt{n}} \leq \mu \leq \bar{x}+Z_{(a/2)}・\frac{σ}{\sqrt{n}}
```

### 母分散未知
母分散未知の母平均の区間推定には、母分散の代わりに**不偏分散s^2**を使い、標準正規分布の代わりに**t分布**を使います。t分布は**自由度（n-1）の値によって分布の形状が変化する**ことが特徴です。

```math
不偏分散s^2=\frac{1}{n-1} \sum_{i=1}^{n}(x_i-x)^2
```
母分散未知の母平均の区間推定のために標本平均を**標準化した統計量t**を、母分散既知の統計量Zと同じように以下のように設定します。

```math
t=\frac{\bar{x}-\mu}{\frac{s}{\sqrt{n}}}
```

信頼区間が95%である場合、統計量Zが標準正規分布の95%の面積（確率）の範囲内にあればいいので、下側確率・上側確率が共に2.5%となるtの値を、**自由度n-1**に注意しながらt分布表を使用して求めます。

t分布表より、t=2.262であることが分かるため以下の式が成り立ちます。

```math
2.262 \leq \frac{\bar{x}-\mu}{\frac{s}{\sqrt{n}}} \leq 2.262
```
求める母平均*μ*について変形すると、95%信頼区間の母分散未知の母平均の区間推定は**標準誤差**を用いて、以下のように表されます。

```math
\bar{x}-2.262・\frac{s}{\sqrt{n}} \leq \mu \leq \bar{x}+2.262・\frac{s}{\sqrt{n}}
```

一般化して**信頼係数100(1-α)%**の場合、母平均の区間推定式は次のように表せます。
**t(α/2)(n-1)は下側確率（上側確率）**を表し、αの値から**t分布表**を用いてtの値を求めます。

```math
\bar{x}-t_{(a/2)(n-1)}・\frac{s}{\sqrt{n}} \leq \mu \leq \bar{x}+t_{(a/2)(n-1)}・\frac{s}{\sqrt{n}}
```

母分散未知の母平均の区間推定方法は、サンプルサイズ（自由度）が大きくなればt分布が標準正規分布N(0,1)に近づくため、母分散既知の場合と同じになります。

母分散既知・母分散未知のどちらの場合も同じように区間推定を行いますが、母分散未知の時は不偏分散s^2・統計量t・自由度n-1を用いることが母分散既知の場合と異なるポイントです。

# 参考文献
* [例題で学ぶ初歩からの統計学 第2版](https://www.amazon.co.jp/%E4%BE%8B%E9%A1%8C%E3%81%A7%E5%AD%A6%E3%81%B6%E5%88%9D%E6%AD%A9%E3%81%8B%E3%82%89%E3%81%AE%E7%B5%B1%E8%A8%88%E5%AD%A6-%E7%AC%AC2%E7%89%88-%E7%99%BD%E7%A0%82-%E5%A0%A4%E6%B4%A5%E8%80%B6/dp/453555790X)
* [母集団 - Wikipedia](https://ja.wikipedia.org/wiki/%E6%AF%8D%E9%9B%86%E5%9B%A3)
* [標本（統計学） - Wikipedia](https://ja.wikipedia.org/wiki/%E6%A8%99%E6%9C%AC_(%E7%B5%B1%E8%A8%88%E5%AD%A6))
* [16-1. 母集団と標本 | 統計学の時間 | 統計WEB](https://bellcurve.jp/statistics/course/8003.html)
* [16-3. 標本の抽出方法 | 統計学の時間 | 統計WEB](https://bellcurve.jp/statistics/course/8007.html)
* [無作為抽出 - Wikipedia](https://ja.wikipedia.org/wiki/%E7%84%A1%E4%BD%9C%E7%82%BA%E6%8A%BD%E5%87%BA)
* [18-1. 点推定とは | 統計学の時間 | 統計WEB - BellCurve](https://bellcurve.jp/statistics/course/8608.html)
* [18-2. 母平均の点推定と推定量・推定値 | 統計学の時間 | 統計WEB - BellCurve](https://bellcurve.jp/statistics/course/8610.html)
* [18-3. 推定量の性質 | 統計学の時間 | 統計WEB - BellCurve](https://bellcurve.jp/statistics/course/8612.html)
* [18-5. 標準偏差と標準誤差 | 統計学の時間 | 統計WEB - BellCurve](https://bellcurve.jp/statistics/course/8616.html)
* [19-1. 区間推定とは | 統計学の時間 | 統計WEB](https://bellcurve.jp/statistics/course/8885.html)
* [19-2. 母平均の信頼区間の求め方（母分散既知） | 統計学の時間 | 統計WEB](https://bellcurve.jp/statistics/course/8888.html)
* [母平均の区間推定](http://hs-www.hyogo-dai.ac.jp/~kawano/HStat/?plugin=cssj&page=2010%2F11th%2FPopulation_Mean)
* [20-2. t分布表 | 統計学の時間 | 統計WEB](https://bellcurve.jp/statistics/course/8970.html)
* [20-2. 母平均の信頼区間の求め方（母分散未知） | 統計学の時間 | 統計WEB](https://bellcurve.jp/statistics/course/8972.html)