【過去問解答】2022年統計検定1級＜数理統計問5＞

統計検定1級の過去問解答解説を行います。目次は以下をご覧ください。

不適切な内容があれば，記事下のコメント欄またはお問い合わせフォームよりご連絡下さい。

問題

統計検定1級の過去問からの出題になります。統計検定の問題の著作権は日本統計学会に帰属していますので，本稿にて記載することはできません。「演習問題を俯瞰する」で詳しく紹介している公式の過去問題集をご購入いただきますようお願い致します。

解答

分散分析モデルの検定と欠損値に関する出題でした。

(1)

正規分布の再生性より$D_i$の従う分布は正規分布であり，平均と分散は下記のように計算されます。

$$\begin{array}{}\text{(2)}& \{\begin{array}{l}E[D_i]=E[Y_i]-E[X_i]=\theta \\ V[D_i]=V[Y_i]+V[X_i]=2{\sigma }^2\end{array}\end{array}$$

(2)

与式(1)の両辺の期待値をとることにより，

$$\begin{array}{}\text{(4)}& \{\begin{array}{l}E[Z_{i1}]=E[X_i]=\nu +a_i+b_1={\mu }_i\\ E[Z_{i2}]=E[Y_i]=\nu +a_i+b_2={\mu }_i+\theta \end{array}\end{array}$$

が得られます。制約条件を使うために両辺の$i=1,\dots ,n$に関する総和を考えると，

$$\begin{array}{}\text{(5)}& \{\begin{array}{l}{b}_1=\overline{\mu }-\nu \\ b_2=\overline{\mu }+\theta -\nu \end{array}\end{array}$$

が得られます。ただし，$\overline{\mu }=\mu /n$とおきました。二つの式の和を考えることにより，

が得られます。これを元の式に代入することにより，

$$\begin{array}{}\text{(7)}& \{\begin{array}{l}{\displaystyle b_1=-\frac{\theta }{2}}\\ {\displaystyle b_2=\frac{\theta }{2}}\\ {\displaystyle a_i={\mu }_i-(\nu +b_1)={\mu }_i-\overline{\mu }}\end{array}\end{array}$$

が得られます。

(3)

$m=\sum _{i,j}{Z}_{ij}/(2n)$，パラメータ$a_1,\dots ,a_n$の因子を$A$，パラメータ$b_1,b_2$の因子を$B$とおき，

とおくとき，総平方和$S$は$A$間平方和$S_A$，$B$間平方和$S_B$，残差平方和$S_E$を用いて

と分解されます。

自分自身の二乗和以外の項は$0$の項が出現して消えてしまうからです。

$S_A$，$S_B$，$S_E$はそれぞれ

$$\begin{array}{}\text{(10)}& \{\begin{array}{l}{\displaystyle S_A=2\sum _{i=1}^n(m_{Ai}-m{)}^2}\\ {\displaystyle S_B=n\sum _{j=1}^2(m_{Bj}-m{)}^2}\\ {\displaystyle S_E=\sum _{i=1}^n\sum _{j=1}^2(Z_{ij}-m_{Ai}-m_{Bj}+m{)}^2}\end{array}\end{array}$$

と表され，自由度はそれぞれ${\varphi }_A=n-1$，${\varphi }_B=2-1=1$，${\varphi }_E=2n-n-2+1=n-1$となります。

${\varphi }_A=n-1$で「$-1$」されているのは，$n-1$群の平均が分かれば残りの$1$群の平均も分かるからです。同様に，${\varphi }_B=2-1$で「$-1$」されているのは，$2-1$群の平均が分かれば残りの$1$群の平均も分かるからです。同様に，${\varphi }_E=2n-n-2+1=n-1$で「$-(n+2-1)$」されているのは，$n+2-1$群の平均が分かれば残りの$1$群の平均も分かるからです。

帰無仮説(a)について，$\theta =0$は$b_1=b_2=0$と同値になりますので，統計量$(S_B/{\varphi }_B)/(S_E/{\varphi }_E)$で要因$B$に関する検定を行えばよく，帰無分布は$\mathrm{F}({\varphi }_B,{\varphi }_E)$となります。同様にして，帰無仮説(b)について，${\mu }_1=\cdots ={\mu }_n=0$は$a_1=\cdots =a_n=0$と同値になりますので，統計量$(S_A/{\varphi }_A)/(S_E/{\varphi }_E)$で要因$A$に関する検定を行えばよく，帰無分布は$\mathrm{F}({\varphi }_A,{\varphi }_E)$となります。

(4)

具体的に最尤推定量を求めて議論できますが，too muchな印象を受けます。

(5)

前問同様，具体的に最尤推定量を求めて議論できますが，too muchな印象を受けます。