統計検定1級の過去問解答解説を行います。目次は以下をご覧ください。
不適切な内容があれば，記事下のコメント欄またはお問い合わせフォームよりご連絡下さい。
問題
統計検定1級の過去問からの出題になります。統計検定の問題の著作権は日本統計学会に帰属していますので，本稿にて記載することはできません。「演習問題を俯瞰する」で詳しく紹介している公式の過去問題集をご購入いただきますようお願い致します。
解答
順序統計量の確率密度関数に関する出題でした。
(1)
f_{1}(y)=3(1-y)^{2},\quad f_{3}(y)=3y^{2},\quad E[Y_{1}]=\frac{1}{4},\quad E[Y_{3}]=\frac{3}{4}
\end{align}
一様分布$\U(0,1)$の確率密度関数と累積密度関数は
f(x) &=
\begin{cases}
1&(0<x<1)\\[0.7em]
0&(\text{その他})
\end{cases},\quad
F(x) =
\begin{cases}
0&(x\leq 0)\\[0.7em]
x&(0<x\leq 1)\\[0.7em]
1&(1<x)
\end{cases}
\end{align}
となります。順序統計量の確率密度関数は累積分布関数から求めます。$0{<}y{<}1$に対し，$Y_{1}$の累積分布関数は
F_{1}(y) &= P(Y_{1}\leq y) = 1-P(Y_{1}>y,Y_{2}>y,Y_{3}>y) = 1-P(X>y)^{3} = 1-(1-y)^{3}
\end{align}
となるため，確率密度関数は$F_{1}^{\prime}(y)=3(1-y)^{2}$と求められます。したがって，
E[Y_{1}] &= \int_{0}^{1}3(1-y)^{2}dy = \frac{1}{4}
\end{align}
が得られます。同様に，$0{<}y{<}1$に対し$Y_{3}$の累積分布関数は
F_{3}(y) &= P(Y_{3}\leq y) = P(Y_{1}<y,Y_{2}<y,Y_{3}<y) = P(X<y)^{3} = y^{3}
\end{align}
となるため，確率密度関数は$F_{3}^{\prime}(y)=3y^{2}$と求められます。したがって，
E[Y_{3}] &= \int_{0}^{1}3y^{2}dy = \frac{3}{4}
\end{align}
が得られます。
補足
一般に，順序統計量$X_{(i)}$の累積分布関数と確率密度関数は
F_{X_{(i)}}(x) &= \sum_{k=i}^{n} {}_{n} C{k}~F(x)^{k}\left\{1-F(x)\right\}^{n-k}\label{F_i} \\[0.7em]
f_{X_{(i)}}(x) &= {}_{n-1} C_{i-1}~ nf(x)F(x)^{i-1}\left\{1-F(x)\right\}^{n-i}\label{f_i}
\end{align}
と表されることを利用すると，
\begin{cases}
f_{1}(y) = {}_{2}C_{0}~3\cdot 1\cdot y^{0}(1-y)^{2} = 3(1-y)^{2}\\[0.7em]
f_{3}(y) = {}_{2}C_{2}~3\cdot 1\cdot y^{2}(1-y)^{0} = 3y^{2}
\end{cases}
\end{align}
と求められます。
(2)
f_{2}(y) &= 6y(1-y),\quad P(Y_{2}\leq 1/2) = 1/2
\end{align}
$0{<}y{<}1$に対し，$Y_{1}$の累積分布関数は
F_{2}(y) &= P(Y_{2}\leq y)\\[0.7em]
&= P(Y_{1}\leq y,Y_{2}\leq y,Y_{3}\leq y)+P(Y_{i}\leq y,Y_{j}\leq y,Y_{k}>y)\\[0.7em]
&= y^{3}+3y^{2}(1-y)
\end{align}
となります。ただし，第二項目については$(i,j)$の選び方として${}_{3}C_{2}$通り，$Y_{k}$の選び方として残りの$1$通りであることを利用しています。確率密度関数は$F_{3}^{\prime}(y)=6y(1-y)$と求められます。したがって，
P(Y_{2}\leq 1/2) &= \int_{0}^{1/2}6y(1-y)dy = \frac{1}{2}
\end{align}
が得られます。
補足
式($\ref{f_i}$)を利用すると，
f_{2}(y) = {}_{2}C_{1}~3\cdot 1\cdot y^{1}(1-y)^{1} = 6y(1-y)
\end{align}
と求められます。
(3)
E[Z] &= \frac{1}{2},\quad V[Z] = \frac{1}{20}
\end{align}
期待値の線形性より，
E[Y_{3}-Y_{1}] &= E[Y_{3}]-E[Y_{1}] = \frac{3}{4}-\frac{1}{4} = \frac{1}{2}
\end{align}
が得られます。続いて，$(Y_{1},Y_{3})$は独立ではないため，同時確率密度関数を求めることで分散を定義から計算することを目指します。累積密度関数$P(Y_{1}\leq y_{1},Y_{3}\leq y_{3})$は，
- $X_{1},X_{2},X_{3}$のすべてが$y_{1}$以下
- $X_{1},X_{2},X_{3}$のうち$2$つが$y_{1}$以下で$1$つが$y_{1}$以上$y_{3}$以下
- $X_{1},X_{2},X_{3}$のうち$1$つが$y_{1}$以下で$2$つが$y_{1}$以上$y_{3}$以下
となる確率と等しくなるため，
F_{13}(y_{1},y_{3})
= P(Y_{1}\leq y_{1},Y_{3}\leq y_{3})
&= y_{1}^{3}+{}_{3}C_{2}~y_{1}^{2}(y_{3}-y_{1})+{}_{3}C_{1}~y_{1}(y_{3}-y_{1})^{2}\\[0.7em]
&= y_{1}^{3}-3y_{1}^{2}y_{3}+3y_{1}y_{3}^{2}
\end{align}
と計算できます。
累積密度関数を$P(Y_{1}\leq y,Y_{3}\leq y)$としてしまうと，同時累積密度関数ではなく一変数の累積密度関数となってしまうため，注意してください。
したがって，$0{<}y{<}1$に対し同時確率密度関数は
f_{13}(y_{1},y_{3}) &= \frac{\partial^{2}}{\partial y_{1}\partial y_{3}} = 6(y_{3}-y_{1})
\end{align}
となります。よって，原点周りの二次モーメントは
E[Z^{2}] = E[(Y_{3}-Y_{1})^{2}]
&= 6\int_{0}^{1}\int_{0}^{y_{3}}(y_{3}-y_{1})^{3}dy_{1}dy_{3}\\[0.7em]
&= 6\int_{0}^{1}\left[\frac{1}{4}(y_{3}-y_{1})^{3}\right]_{0}^{y_{3}}dy_{3}\\[0.7em]
&= \frac{3}{2}\int_{y_{1}}^{1}y_{3}^{4}~dy_{3}
= \frac{3}{10}\left[y_{3}^{5}\right]_{0}^{1} = \frac{3}{10}
\end{align}
が得られます。以上より，
V[Z] &= E[Z^{2}]-E[Z]^{2} = \frac{3}{10}-\frac{1}{4} = \frac{1}{20}
\end{align}
となります。
補足
f_{Y_{1},Y_{3}}(y_{1},y_{3}) = 6(y_{3}-y_{1})
\end{align}
に対し，
\begin{cases}
z=y_{3}-y_{1}\\[0.7em]
w=y_{1}
\end{cases}
\end{align}
とおくと，逆変換は
\begin{cases}
y_{1}=w\\[0.7em]
y_{3}=w+z
\end{cases}
\end{align}
と表されるため，ヤコビアンは$1$となります。これより，
f_{Z,W}(z,w) = f_{Y_{1},Y_{3}}(w,w+z) = 6((w+z)-w) = 6z
\end{align}
が得られます。ただし，積分領域は$0<y_{1}<y_{3}<1$より$0<w<w+z<1$となるため，
\begin{cases}
0 < z < 1\\[0.7em]
0 < w < 1-z
\end{cases}
\end{align}
となります。$W$について周辺化すると，
\int_{0}^{1-z}f_{Z,W}(z,w)~dw = \int_{0}^{1-z}6z~dw = 6z(1-z)
\end{align}
となるため，
\begin{cases}
\displaystyle
E[Z] = \int_{0}^{1}z\cdot 6z(1-z)dz = 6\int_{0}^{1}z^{2}-z^{3}dz = 6\cdot \frac{1}{12} = \frac{1}{2}\\[0.7em]
\displaystyle
E[Z^{2}] = \int_{0}^{1}z^{2}\cdot 6z(1-z)dz = 6\int_{0}^{1}z^{3}-z^{4}dz = 6\cdot \frac{1}{20} = \frac{3}{10}\\[0.7em]
\displaystyle
V[Z] = \frac{3}{10} - \frac{1}{4} = \frac{1}{20}
\end{cases}
\end{align}
が得られます。
コメント
コメント一覧 （2件）
問(3)につきまして、 式(20)から Z=Y3-Y1, W=Y1 と変数変換することで、ZとWの同時密度関数 f(z,w)=6z を算出し、それをWについて周辺化することでZの密度関数を求めようとしたのですが、うまくいきませんでした。 この方法でこの問に解答することは可能でしょうか？また、可能であれば周辺化のためのWの積分範囲はどうなるでしょうか。 ご教示いただけますと幸いです。（※変数変換が誤っていたので、再コメントしました。）
ホワイト様
ご返信遅れましてすみません。おそらく試験終わってしまったかと思いますが，本文に追記してみました。