カイ2乗分布の再生性

公開日: 更新日:

【2023年4月1週】 【B000】数理統計学 【B040】連続型の確率分布

この記事をシェアする
  • B!
サムネイル画像

本稿では、①確率変数のたたみこみによる方法、②モーメント母関数を用いる方法、③ガンマ分布との関係を用いる方法の3通りの方法で、カイ2乗分布の再生性を証明しています。①の方法は計算に骨がある問題ですが、②・③の方法はとても簡単です。

なお、閲覧にあたっては、以下の点にご注意ください。

  • スマートフォンやタブレット端末でご覧の際、数式が見切れている場合は、横にスクロールすることができます。

【定理】カイ2乗分布の再生性

【定理】
カイ2乗分布の再生性
Reproductive Property of Chi-squared Distribution

確率変数 $X,Y$ がそれぞれ独立に $\chi^2$分布 \begin{align} X \sim \chi^2 \left(n_1\right) \quad Y \sim \chi^2 \left(n_2\right) \end{align} に従うとき、 確率変数 $X$ と $Y$ の和を \begin{align} Z=X+Y \end{align} とすると、 新たな確率変数 $Z$ は、$\chi^2$分布 \begin{align} Z \sim \chi^2 \left(n_1+n_2\right) \end{align} に従う。

証明法①:確率変数のたたみこみによる方法

証明

確率変数のたたみこみの公式 $k \left(z\right)=\int_{0}^{z}{g \left(x\right) \cdot h \left(s-x\right)}$ より、 \begin{align} k \left(z\right)&=\int_{0}^{z}{ \left\{\frac{1}{2^\frac{n_1}{2} \cdot \Gamma \left(\frac{n_1}{2}\right)}x^{\frac{n_1}{2}-1}e^{-\frac{x}{2}}\right\} \cdot \left\{\frac{1}{2^\frac{n_2}{2} \cdot \Gamma \left(\frac{n_2}{2}\right)} \left(z-x\right)^{\frac{n_2}{2}-1}e^{-\frac{z-x}{2}}\right\}dx}\\ &=\frac{z^{\frac{n_2}{2}-1} \cdot e^{-\frac{z}{2}}}{2^\frac{n_1+n_2}{2} \cdot \Gamma \left(\frac{n_1}{2}\right)\Gamma \left(\frac{n_2}{2}\right)}\int_{0}^{z}{x^{\frac{n_1}{2}-1} \cdot \left\{z \left(1-\frac{x}{z}\right)\right\}^{\frac{n_2}{2}-1}dx}\\ &=\frac{z^{\frac{n_2}{2}-1} \cdot e^{-\frac{z}{2}}}{2^\frac{n_1+n_2}{2} \cdot \Gamma \left(\frac{n_1}{2}\right)\Gamma \left(\frac{n_2}{2}\right)}\int_{0}^{z}{x^{\frac{n_1}{2}-1} \left(1-\frac{x}{z}\right)^{\frac{n_2}{2}-1}dx} \end{align} ここで、以下のように変数変換すると、 \begin{gather} u=\frac{x}{z}\Leftrightarrow x=uz\\ \frac{dx}{du}=z\Rightarrow dx=zdu\\ x:0\rightarrow z \quad \Rightarrow \quad u:0\rightarrow1 \end{gather} となるので、 置換積分法により、 \begin{align} k \left(z\right)&=\frac{z^{\frac{n_2}{2}-1} \cdot e^{-\frac{z}{2}}}{2^\frac{n_1+n_2}{2} \cdot \Gamma \left(\frac{n_1}{2}\right)\Gamma \left(\frac{n_2}{2}\right)}\int_{0}^{1}{ \left(uz\right)^{\frac{n_1}{2}-1} \left(1-u\right)^{\frac{n_2}{2}-1} \cdot z d u}\\ &=\frac{z^{\frac{n_1+n_2}{2}-1} \cdot e^{-\frac{z}{2}}}{2^\frac{n_1+n_2}{2} \cdot \Gamma \left(\frac{n_1}{2}\right)\Gamma \left(\frac{n_2}{2}\right)}\int_{0}^{1}{u^{\frac{n_1}{2}-1} \left(1-u\right)^{\frac{n_2}{2}-1}du} \end{align} ベータ関数の定義式 $B \left(s,t\right)=\int_{0}^{1}{u^{s-1} \left(1-u\right)^{t-1}du}$ より、 \begin{align} k \left(z\right)=\frac{z^{\frac{n_1+n_2}{2}-1} \cdot e^{-\frac{z}{2}}}{2^\frac{n_1+n_2}{2} \cdot \Gamma \left(\frac{n_1}{2}\right)\Gamma \left(\frac{n_2}{2}\right)} \cdot B \left(\frac{n_1}{2},\frac{n_2}{2}\right) \end{align} ベータ関数の性質 $B \left(s,t\right)=\frac{\Gamma \left(s\right)\Gamma \left(t\right)}{\Gamma \left(s+t\right)}$ より、 \begin{align} k \left(z\right)&=\frac{z^{\frac{n_1+n_2}{2}-1} \cdot e^{-\frac{z}{2}}}{2^\frac{n_1+n_2}{2} \cdot \Gamma \left(\frac{n_1}{2}\right)\Gamma \left(\frac{n_2}{2}\right)} \cdot \frac{\Gamma \left(\frac{n_1}{2}\right)\Gamma \left(\frac{n_2}{2}\right)}{\Gamma \left(\frac{n_1+n_2}{2}\right)}\\ &=\frac{1}{2^\frac{n_1+n_2}{2} \cdot \Gamma \left(\frac{n_1+n_2}{2}\right)} \cdot z^{\frac{n_1+n_2}{2}-1} \cdot e^{-\frac{z}{2}} \end{align} これは、$\chi^2$分布の確率密度関数 \begin{align} f \left(x\right)=\frac{1}{2^\frac{n}{2} \cdot \Gamma \left(\frac{n}{2}\right)}x^{\frac{n}{2}-1}e^{-\frac{x}{2}} \end{align} において、 \begin{align} n\rightarrow n_1+n_2 \quad x\rightarrow z \end{align} と置き換えたものとみなすことができる。 したがって、確率関数の一意性により、確率変数 $Z$ は、$\chi^2$分布 \begin{align} Z \sim \chi^2 \left(n_1+n_2\right) \end{align} に従う。 $\blacksquare$

証明法②:モーメント母関数を用いる方法

証明

$\boldsymbol{\chi}^\boldsymbol{2}$分布のモーメント母関数の公式より、 \begin{gather} M_X \left(\theta\right)= \left(\frac{1}{1-2\theta}\right)^\frac{n_1}{2}\\ M_Y \left(\theta\right)= \left(\frac{1}{1-2\theta}\right)^\frac{n_2}{2} \end{gather} モーメント母関数の性質 $M_Z \left(\theta\right)=M_X \left(\theta\right) \cdot M_Y \left(\theta\right)$ より、 \begin{align} M_Z \left(\theta\right)&= \left(\frac{1}{1-2\theta}\right)^\frac{n_1}{2} \cdot \left(\frac{1}{1-2\theta}\right)^\frac{n_2}{2}\\ &= \left(\frac{1}{1-2\theta}\right)^\frac{n_1+n_2}{2} \end{align} これは、$\boldsymbol{\chi}^\boldsymbol{2}$分布のモーメント母関数 \begin{align} M_X \left(\theta\right)= \left(\frac{1}{1-2\theta}\right)^\frac{n}{2} \end{align} において、 \begin{align} n\rightarrow n_1+n_2 \quad X\rightarrow Z \end{align} と置き換えたものとみなすことができる。 したがって、モーメント母関数の一意性により、確率変数 $Z$ は、$\boldsymbol{\chi}^\boldsymbol{2}$分布 \begin{align} Z \sim \chi^2 \left(n_1+n_2\right) \end{align} に従う。 $\blacksquare$

証明法③:ガンマ分布との関係を用いる方法

証明

$\chi^2$分布とガンマ分布の関係より、 \begin{gather} \chi^2 \left(n_1\right)=\mathrm{Ga} \left(\frac{n_1}{2},\frac{1}{2}\right)\\ \chi^2 \left(n_2\right)=\mathrm{Ga} \left(\frac{n_2}{2},\frac{1}{2}\right) \end{gather} ガンマ分布の再生性より、 \begin{align} \chi^2 \left(n_1\right)+\chi^2 \left(n_2\right)=\mathrm{Ga} \left(\frac{n_1}{2}+\frac{n_2}{2},\frac{1}{2}\right)=\chi^2 \left(n_1+n_2\right) \end{align} $\blacksquare$

参考文献

  • 野田 一雄, 宮岡 悦良 著. 入門・演習数理統計. 共立出版, 1990, p.145

関連記事

自己紹介

自分の写真

yama

大学時代に読書の面白さに気づいて以来、読書や勉強を通じて、興味をもったことや新しいことを学ぶことが生きる原動力。そんな人間が、その時々に学んだことを備忘録兼人生の軌跡として記録しているブログです。

このブログを検索

ブログ アーカイブ

QooQ