方針

カイ二乗分布について書く．

Z_1,...,Z_n:iid\sim\mathcal{N}(0,1)

の時，

X=Z_1^2+...+Z_n^2\sim\chi_n^2=Ga(\frac{n}{2},\frac{1}{2})

で，確率密度関数は，

f_{X}(x)=\frac{(\frac{1}{2})^{n/2}}{\Gamma(\frac{n}{2})}x^{\frac{n}{2}-1}\exp(-\frac{1}{2}x)

である．積率母関数は

\begin{align}
M_X(t)&= \mathbb{E}_{X\sim\chi_n^2}[e^{tX}]\\
&= \int_0^\infty\frac{(\frac{1}{2})^{n/2}}{\Gamma(\frac{n}{2})}x^{\frac{n}{2}-1}\exp\{-(\frac{1}{2}-t)x\}dx\\
&= \frac{(\frac{1}{2})^{n/2}}{\Gamma(\frac{n}{2})}\frac{\Gamma(\frac{n}{2})}{(\frac{1}{2}-t)^{n/2}}\int_0^\infty\frac{(\frac{1}{2}-t)^{n/2}}{\Gamma(\frac{n}{2})}x^{\frac{n}{2}-1}\exp\{-(\frac{1}{2}-t)x\}dx\\
&= \frac{1}{(1-2t)^{n/2}}
\end{align}

である．まずは

Z_1^2\sim\chi_1^2

を示す．

Y_1=Z_1^2,\ Y_2=Z_2^2,\ ...,\ Y_n=Z_n^2

と置くと

\begin{align}
f_{Y_1}(y)&= f_{Z_1}(\sqrt{y})\frac{1}{2\sqrt{y}}-f_{Z_1}(-\sqrt{y})(-\frac{1}{2\sqrt{y}})\\
&= \frac{1}{\sqrt{y}}\frac{1}{\sqrt{2\pi}}\exp(-\frac{y}{2})\\
&= \frac{(\frac{1}{2})^{1/2}}{\Gamma(\frac{1}{2})}y^{\frac{1}{2}-1}\exp(-\frac{1}{2}y)
\end{align}

ゆえ，

Y_1,...,Y_n:iid\sim\chi_1^2

が示された．次に，

X=Y_1+...+Y_n

と置くと

\begin{align}
M_X(t)&= \mathbb{E}_X[e^{tX}]\\
&= \mathbb{E}_{Y_1,...,Y_n:iid\sim\chi_1^2}[e^{t(Y_1+...+Y_n)}]\\
&= \mathbb{E}_{Y_1\sim\chi_1^2}[e^{tY_1}]...\mathbb{E}_{Y_n\sim\chi_1^2}[e^{tY_n}]\\
&= \frac{1}{(1-2t)^{1/2}}...\frac{1}{(1-2t)^{1/2}}\\
&= \frac{1}{(1-2t)^{n/2}}
\end{align}

より，

X=Y_1+...+Y_n=Z_1^2+...+Z_n^2\sim\chi_n^2

が示された．問題はカイ二乗分布の期待値演算なので，ガンマ分布の時と同様に計算すれば良い．

1.答案

\begin{align}
\mathbb{E}_{X\sim\chi_n^2}[X^\nu]&= \int_0^\infty\frac{(\frac{1}{2})^{n/2}}{\Gamma(\frac{n}{2})}x^{\nu+\frac{n}{2}-1}\exp(-\frac{1}{2}x)dx\\
&= \frac{(\frac{1}{2})^{n/2}}{\Gamma(\frac{n}{2})}\frac{\Gamma(\nu+\frac{n}{2})}{(\frac{1}{2})^{\nu+\frac{n}{2}}}\int_0^\infty\frac{(\frac{1}{2})^{\nu+\frac{n}{2}}}{\Gamma(\nu+\frac{n}{2})}x^{\nu+\frac{n}{2}-1}\exp(-\frac{1}{2}x)dx\\
&= \frac{2^\nu\Gamma(\nu+\frac{n}{2})}{\Gamma(\frac{n}{2})}
\end{align}

2.答案

自由度nのカイ二乗分布に従う確率変数をX_nと置くことにする．

X_n\sim\chi_n^2

この時，X_{n+2}の確率密度関数は

\begin{align}
f_{X_{n+2}}(x)&= \frac{(\frac{1}{2})^{\frac{n}{2}+1}}{\Gamma(\frac{n}{2}+1)}x^{n/2}\exp(-\frac{1}{2}x)\\
&= \frac{\frac{1}{2}x}{\frac{n}{2}}\frac{(\frac{1}{2})^{\frac{n}{2}}}{\Gamma(\frac{n}{2})}x^{\frac{n}{2}-1}\exp(-\frac{1}{2}x)\\
&= \frac{x}{n}f_{X_n}(x)
\end{align}

と表せ，

\begin{align}
f_{X_n}(x)&= \frac{n}{x}f_{X_{n+2}}(x)\\
f_{X_n}(x)&= \frac{x}{n-2}f_{X_{n-2}}(x)
\end{align}

が成り立つ．ゆえ，

\begin{align}
\mathbb{E}_{X_n\sim\chi_n^2}[h(X_n)]&= \int_0^\infty h(x)f_{X_n}(x)dx\\
&= \int_0^\infty h(x)\frac{n}{x}f_{X_{n+2}}(x)dx\\
&= n\mathbb{E}_{X_{n+2}\sim\chi_{n+2}^2}[\frac{h(X_n)}{X_n}]\\
\mathbb{E}_{X_n\sim\chi_n^2}[h(X_n)]&= \int_0^\infty h(x)f_{X_n}(x)dx\\
&= \int_0^\infty h(x)\frac{x}{n-2}f_{X_{n-2}}(x)dx\\
&= \frac{1}{n-2}\mathbb{E}_{X_{n-2}\sim\chi_{n-2}^2}[X_{n-2}h(X_{n-2})]
\end{align}

参考文献

『現代数理統計学の基礎』(久保川達也　著)