現代数理統計学の基礎 7章問2 (1) - 脳内ライブラリアン

検定が来年になったため非常にゆっくりなペースとなっていますが、今回は7章の問2。

複合仮説の尤度比検定の問題ですね。

なかなか解答も理解しづらく、説明いれてくとかなり長くなってしまいました、、。似たような問題が今後も続くため今回は丁寧めに説明を考えています。

目次：

１、尤度比検定の棄却域を求める
２、確率上界の証明をする

１、尤度比検定の棄却域を求める

まずは帰無仮説 $H_0:\mu\leq\mu_0$ における尤度比検定の棄却域を求めます。

尤度比検定は

尤度比＝(帰無仮説下でのMLE＜最尤推定量＞)/(特に指定されない場合のMLE)を用いる検定でした。

尤度比検定についての説明記事はこちら↓

medibook.hatenablog.com

そこで今回の複合仮説において帰無仮説下のMLEを考えると、 $\mu_0$ もしくは $\hat\mu$ （制限をつけない通常の最尤推定量）のどちらか小さいほうとなります。

よって式で表現するとしたら帰無仮説下のMLEを $\hat\mu_0$ とすると、 $\hat\mu_0=min(\mu_0, \hat\mu)$ ですね。どちらが帰無仮説下のMLEになるかは場合分けすればよいと言えます。

肝心の尤度比の式は、尤度関数を $L(\mu)$ とすると

$\frac{L(\hat\mu_0)}{L(\hat\mu)}=exp\{-\frac{n}{2\sigma^2}(\bar X-\hat\mu_0)^2\}$

となります。

先ほどの場合分けをして代入してみると

① $\hat\mu\leq\mu_0$ の場合

$\hat\mu_0=\hat\mu$ となります。これを尤度比の式に代入すると

$exp\{-\frac{n}{2\sigma^2}(\bar X-\hat\mu)^2\}=1$

これは $\bar X=\hat\mu$ であることによります。

これはつまり、帰無仮説下のMLEとそうでない場合のMLEが一致してしまう（尤度比が１）ので、そもそも帰無仮説の棄却ができません。

公式の解答では詳しく書いてはいないですが、この場合はそもそも棄却域が存在しないと言うことで良いのではないでしょうか。

② $\hat\mu\geq\mu_0$ の場合

この場合、 $\hat\mu_0=\mu_0$ となります。これを尤度比の式に代入すると

$exp\{-\frac{n}{2\sigma^2}(\bar X-\mu_0)^2\}$ となるので、-2をかけて対数をとるとカイ二乗分布に従うことを利用して

$\frac{n}{\sigma^2}(\bar X-\mu_0)^2\sim\chi_1^2$

となります。

解答にあるように $\hat\mu=\bar X\geq\mu_0$ なので

尤度比検定の式は、有意水準がα％となる分位点 $z_\alpha$ を用いて

$\sqrt{\frac{n}{\sigma^2}}(\bar X-\mu_0)\gt z_\alpha$ となります。

ここまででようやく尤度比検定の式が求まりました。

２、確率上界の証明をする

あとは複合仮説の場合、Ｐ値の定義として、確率上界が

$sup_{\mu\leq\mu_0}P(\sqrt{\frac{n}{2\sigma^2}}(\bar X-\mu_0)\gt z_\alpha)$

となることを示さないといけません。

この辺は個人的にとてもややこしかったので以前に一度記事を書きました↓

medibook.hatenablog.com

過去記事と同様に進めてみると

$\sqrt{\frac{n}{\sigma^2}}(\bar X-\mu_0)\gt z_\alpha$ をまず $\mu$ を用いた式に変形します。

$\sqrt{\frac{n}{\sigma^2}}(\bar X-\mu)+\sqrt{\frac{n}{\sigma^2}}(\mu-\mu_0)\gt z_\alpha\\\sqrt{\frac{n}{\sigma^2}}(\bar X-\mu)\gt z_\alpha-\sqrt{\frac{n}{\sigma^2}}(\mu-\mu_0)$

ここで

$\sqrt{\frac{n}{\sigma^2}}(\bar X-\mu)\sim N(0,1)$

となることと

$-\mu+\mu_0\geq 0$

と利用して、上記の式の確率をとると

$P\{\sqrt{\frac{n}{\sigma^2}}(\bar X-\mu)\gt z_\alpha-\sqrt{\frac{n}{\sigma^2}}(\mu-\mu_0)\}\\\leq P\{\sqrt{\frac{n}{\sigma^2}}(\bar X-\mu)\gt z_\alpha=\alpha$

となり、確率上界がαとなるような検定であることが証明できました。

5件のコメント

確率上界の証明での式でμ0をμの式で書き直すところですが、μ-μ0に係数(n/σ^2)^1/2が書かれていません。
最終的に定数倍を掛けるか掛けないかの差になるので、結果は変わらないと思いますが。

返信

>nuka137さん
ご指摘ありがとうございます。最終的に消えていくものなので、気付いてなかったみたいです。助かりました。

脳内ライブラリアン様
いつも大変助かっております。
独学で周りに質問できない身としてはありがたい限りです。
帰無仮説下のMLEがμ0もしくはμ^なのは尤度関数が-∞からμ^まで単調増加だからという理由でよろしいでしょうか。
ご教示願います。
おしょ

＞おしょさん
コメント有難うございます。自分もほぼ独学なのでお役立ていただけますと幸いです。
MLEは書いていただいていますように、尤度関数がμ^まで単調増加であるため、パラメータの範囲にμ^が含まれる場合は、μ^で、範囲にμ^が含まれず端がμ0になる場合はμ0ということになっていると考えています。