はじめに
代表的な標本分布についてまとめています。正規母集団からの標本分布から始まり、t-分布と、F分布についてもまとめていきます。
標本分布
平均が μ、分散が σ2 の確率分布に従う母集団からランダムサンプルされた標本を用いて計算される標本平均、標本分散はそれぞれ
Xˉ=n1i=1∑nXi
S2=n1i=1∑n(Xi−Xˉ)2
で定義されます。このように母数(この場合だと μ や σ2)を含んでいないものは統計量と呼ばれ、それらの従う確率分布を標本分布と呼びます。
ここで問題となってくるのが、標本分布は果たしてどのような分布になっているのか?ということです。母集団が正規分布である場合と、母集団以外の一般的な確率分布に従う場合とで考え方が異なってきます。
正規母集団からの標本分布
平均が μ、分散が σ2 の正規分布に従う母集団からランダムサンプルされた標本
X1,…,Xn,i.i.d.∼N(μ,σ2)
を用いた標本平均と不偏分散 が従う確率分布(標本分布)について考えてみます。
標本平均について
標本平均は
Xˉ=n1i=1∑nXi=n1(X1+…+Xn)
のように正規分布に従う確率分布 Xi の和であるため、正規分布の再生性より Xˉ も正規分布に従うことが分かります。あとは平均と分散を求めると、
E[Xˉ]=…=μ
V[Xˉ]=…=nσ2
であるので、Xˉ は N(μ,σ2/n) に従うことが分かります。
不偏分散について
不偏分散は
V2=n−11i=1∑n(Xi−Xˉ)2
で、これを式変形していきます。かなり天下り的なため、変形の雰囲気だけでも覚えておいた方が良いと思います。
(n−1)V2(n−1)V2/σ2=i=1∑n(Xi−μ+μ−Xˉ)2=i=1∑n(Xi−μ)2−2i=1∑n(Xi−μ)(Xˉ−μ)+i=1∑n(Xˉ−μ)2=i=1∑n(Xi−μ)2−2n(Xˉ−μ)2+n(Xˉ−μ)2=i=1∑n(Xi−μ)2−n(Xˉ−μ)2=i=1∑n(σXi−μ)2−(σ/nXˉ−μ)2
ここで右辺第一項の中身は
σXi−μ∼N(0,1)
であり、右辺第二項の中身は
σ/nXˉ−μ∼N(0,1)
です。標準正規分布に従う確率変数の二乗和がカイ二乗分布に従うことを思い出すと、
(n−1)V2/σ2=χn2−χ12=χn−12
となり、不偏分散は自由度 n−1 のカイ二乗分布に従うことが分かりました。
補足
標本平均、不偏分散についての標本分布はヘルマート行列を用いて計算すると証明することができるのですが、毎回計算しているのも面倒なので、上記の流れで覚えておいたほうが忘れたときにぱっと計算し直せるかと思います。
t-分布
分布の定義
ここで、まず t-分布と呼ばれる確率分布を定義します。互いに独立な確率変数 Z∼N(0,1)、U∼χm2 とするとき、
T=U/mZ
で定義される確率変数 T は自由度 m の t-分布と呼ばれる確率分布に従います。
標本を用いた定義
標本平均は
Xˉ∼N(μ,σ2/n)
のように正規分布に従い、標準化すると
Z=σ/nXˉ−μ∼N(0,1)
として標準正規分布に従います。また不偏分散は
σ2n−1V2∼χn−12
として自由度 n−1 のカイ二乗分布に従います。そのためこれらを用いた確率変数
σ2n−1V2/n−11σ/nXˉ−μ
は t-分布に従うことが分かります。式変形していくと
V/nXˉ−μ
となり、標本平均を標準化した式において母分散を不偏分散に置き換えた式が自由度 n−1 の t-分布に従うという性質を導くことができました。
σ/nXˉ−μ∼N(0,1)→V/nXˉ−μ∼T
F分布
分布の定義
互いに独立な確率変数 S∼χm2、T∼χn2 とするとき、
Y=T/nS/m
で定義される確率変数 Y は自由度 (m,n) の F分布と呼ばれる確率分布に従います。
標本を用いた定義
不偏分散が
σ2n−1V2∼χn−12
として自由度 n−1 のカイ二乗分布に従うことを用いると
σ2m−1V22/(m−1)σ2n−1V12/(n−1)=V22/σ22V12/σ12
となり、不偏分散を用いて表すことができます。