📊

【統計検定準1級】母分散の区間推定

2024/08/18に公開

はじめに

この記事では、統計検定準1級取得に向けて学習したことをまとめていきます。
工学系の数学ではなく数理あるあるの、論述ゴリゴリな解答になっていると思いますのであらかじめご了承ください。
注意:さらに計算過程は数学文化の『省略の美』を無視してエレファントに書いています。

【リンク紹介】
統計検定準1級のまとめ記事一覧
これまで書いたシリーズ記事一覧

学習書籍について

この記事では「統計学実践ワークブック」を中心に、学んだことをまとめていきます。記事を読んで本格的に勉強してみたいなと思った方は、是非ご購入を検討なさってください。


参考書籍について

統計実践ワークブックは、大量の知識項目と問題が収められている反面、計算過程や知識背景が大きく省略されているため、知識体系をきちんと学ぶ参考書として東京大学から出版されている名著「統計学入門」を使っています。
※ワークブックとしては素晴らしい質だと思いますが、どうしてもその内容量とページ数の都合上、問題のない範囲で削除されているということです。人によっては1冊で問題ない方もおられると思いますが、私には無理でした。


区間推定

データを用いて未知母数\thetaの値の存在範囲を区間として推定する方法を区間推定(interval estimation)という。

母分散の区間推定

X_1, X_2, \cdots, X_nを、母集団からの無作為標本とする。
今、X_1, X_2, \cdots, X_nは互いに独立で、それぞれ同一の(母平均\mu母分散\sigma^2の)正規分布N(\mu, \sigma^2)に従うとする。

このとき、母分散\sigma^2の区間推定を行う。

X_1, X_2, \cdots, X_nにおいて標本分布S^2不偏標本分布U^2とする。
つまり、

S^2 = \cfrac{1}{ \ n \ } \sum_{i = 1}^{n} (X_i - \overline{X})^2
U^2 = \cfrac{1}{n - 1} \sum_{i = 1}^{n} (X_i - \overline{X})^2

である。ここで、統計量

\cfrac{(n - 1) U^2}{\sigma^2}

は、母分散が未知のときの母平均の区間推定の時と同様に自由度n - 1カイ2乗分布\chi^2 (n - 1)に従う。※標本分散S^2も不偏標本分散U^2の代わりに使えます。
ここで、カイ2乗分布表から次の式が成り立つ。

P \left( \chi^2_{0.975} (n - 1) \leqq \cfrac{(n - 1) U^2}{\sigma^2} \leqq \chi^2_{0.025} (n - 1) \right) = 0.95

ただし、自由度n - 1のカイ2乗分布の上側確率\alphaに対する\chi^2_{\alpha} (n - 1)とする。すると、

\begin{alignat*}{2} P \left( \chi^2_{0.975} (n - 1) \leqq \cfrac{(n - 1) U^2}{\sigma^2} \leqq \chi^2_{0.025} (n - 1) \right) &= 0.95 \\ P \left( \cfrac{1}{\chi^2_{0.025} (n - 1)} \leqq \cfrac{\sigma^2}{(n - 1) U^2} \leqq \cfrac{1}{\chi^2_{0.975} (n - 1)} \right) &= 0.95 \\ P \left( \cfrac{(n - 1) U^2}{\chi^2_{0.025} (n - 1)} \leqq \sigma^2 \leqq \cfrac{(n - 1) U^2}{\chi^2_{0.975} (n - 1)} \right) &= 0.95 \\ \end{alignat*}

となる。以上より、母分散\sigma^2信頼係数0.95信頼区間は、

\left[ \cfrac{(n - 1) U^2}{\chi^2_{0.025} (n - 1)}, \ \cfrac{(n - 1) U^2}{\chi^2_{0.975} (n - 1)} \right]

となる。

参考資料

\bf{\textcolor{red}{記事が役に立った方は「いいね」を押していただけると、すごく喜びます \ 笑}}
ご協力のほどよろしくお願いします。

Discussion