はじめに
この記事では、統計検定準1級取得に向けて学習したことをまとめていきます。
工学系の数学ではなく数理あるあるの、論述ゴリゴリな解答になっていると思いますのであらかじめご了承ください。
注意:さらに計算過程は数学文化の『省略の美』を無視してエレファントに書いています。
【リンク紹介】
・統計検定準1級のまとめ記事一覧
・これまで書いたシリーズ記事一覧
学習書籍について
この記事では「統計学実践ワークブック」を中心に、学んだことをまとめていきます。記事を読んで本格的に勉強してみたいなと思った方は、是非ご購入を検討なさってください。
![](https://storage.googleapis.com/zenn-user-upload/e8c0414ab477-20240702.jpg)
参考書籍について
統計実践ワークブックは、大量の知識項目と問題が収められている反面、計算過程や知識背景が大きく省略されているため、知識体系をきちんと学ぶ参考書として東京大学から出版されている名著「統計学入門」を使っています。
※ワークブックとしては素晴らしい質だと思いますが、どうしてもその内容量とページ数の都合上、問題のない範囲で削除されているということです。人によっては1冊で問題ない方もおられると思いますが、私には無理でした。
![](https://storage.googleapis.com/zenn-user-upload/37d6f140a49a-20240702.jpg)
カイ2乗分布
確率変数X_i \ (1 \leqq i \leqq n)が独立でそれぞれ標準正規分布N(0, 1)に従うとする。このとき、統計量Yを
Y = {X_{1}}^{2} + {X_{2}}^{2} + \cdots + {X_{n}}^{2}
とおくとき、Yの確率密度関数fが
\begin{alignat*}{2}
f(y) &= \cfrac{1}{2^{\frac{n}{2}} \ \Gamma \left( \frac{n}{2} \right)} \
y^{\frac{n}{2} - 1} e^{- \frac{y}{2}}
\hspace{5mm}
(y > 0)
\end{alignat*}
であるとき、この確率密度関数fを自由度nカイ2乗分布(chi-square distribution with n degrees of freedom)といい、\chi^2 (n)と表す。ただし、\Gammaはガンマ分布のガンマ関数を表し、
\Gamma \left( \frac{n}{2} \right) =
\int_{0}^{\infty} x^{\frac{n}{2} - 1} e^{-x} dx
である。
カイ2乗分布の期待値
E[Y]:Yの期待値
カイ2乗分布の分散
V[Y]:Yの分散
カイ2乗分布のモーメント母関数
m(\theta) = E[e^{\theta Y}]:Yのモーメント母関数
\begin{alignat*}{2}
m(\theta) &= E[e^{\theta Y}] \\
&= (1 - 2 \theta)^{- \frac{n}{2}}
\hspace{5mm}
\left( \theta < \cfrac{1}{ \ 2 \ } \right) \\
\end{alignat*}
非心カイ2乗分布
確率変数X_i \ (1 \leqq i \leqq n)が独立でそれぞれ正規分布N(\mu_i, 1)に従うとする。
統計量Yを
Y = {X_{1}}^{2} + {X_{2}}^{2} + \cdots + {X_{n}}^{2}
とし、
\lambda = {\mu_{1}}^{2} + {\mu_{2}}^{2} + \cdots + {\mu_{n}}^{2}
とする。このときYの確率密度関数fが
\begin{alignat*}{2}
f(y) &= \sum_{k = 0}^{\infty}
\cfrac{e^{- \frac{\lambda}{2}} \ (\frac{\lambda}{2})^k}{k!} \cdot
\cfrac{1}{2^{\frac{n}{2} + k} \ \Gamma(\frac{n}{2} + k)}
\ y^{\frac{n}{2} + k - 1}
\ e^{- \frac{y}{2}}
\hspace{5mm}
(y > 0)
\end{alignat*}
であるとき、この確率密度関数fを自由度n非心度\lambdaの非心カイ2乗分布といい、\chi^2 (n, \lambda)と表す。ただし、\Gammaはガンマ分布のガンマ関数を表し、
\Gamma \left( \frac{n}{2} + k \right) =
\int_{0}^{\infty} x^{\frac{n}{2} + k - 1} e^{-x} dx
である。
※確率密度関数は以下のサイトを参考にしました。
https://pystyle.info/statistics-noncentral-chi-square-distribution/
非心カイ2乗分布の期待値
E[Y]:Yの期待値
非心カイ2乗分布の分散
V[Y]:Yの分散
非心カイ2乗分布のモーメント母関数
m(\theta) = E[e^{\theta Y}]:Yのモーメント母関数
\begin{alignat*}{2}
m(\theta) &= E[e^{\theta Y}] \\
&= (1 - 2 \theta)^{- \frac{n}{2}}
\exp \left( \cfrac{\lambda \theta}{1 - 2 \theta} \right)
\hspace{5mm}
\left( \theta < \cfrac{1}{ \ 2 \ } \right) \\
\end{alignat*}
参考資料
\bf{\textcolor{red}{記事が役に立った方は「いいね」を押していただけると、すごく喜びます \ 笑}}
ご協力のほどよろしくお願いします。
Discussion