※ブログ記事の商品・サービスリンクにはアフィリエイトリンクが含まれます。

現代数理統計学の基礎 7章 問10(3)

さて、(3)は一様最強力不偏検定を示す問題です。

公式の解答を見ると

\(\sqrt{n}|x| > z_{\frac{\alpha}{2}}\) が一様最強力不偏検定である

から唐突に始まってます。

それを導出していくというより、これが一様最強力不偏検定であることを示す感じになってますね。問題の意図がそれでいいならまあ仕方ないのですが、、、。

なので同様に、一様最強力不偏検定の示し方と先程の検定がそれで良いのかどうかを書きます。

昨日書いた記事はμの値を一般化したものなので、やり方は本当にほぼ一緒です。なので、コピペしまくります。

medibook.hatenablog.com

まず棄却域をRとして検定関数を以下で定義します。

\(\phi(x)=1 (x\in R)\\\phi(x)=0 (x\notin R)\)

一般性を失わないので、\(\bar{X}=x, n=1\)とします。

次に帰無仮説において、有意水準αを満たし、かつ対立仮説の検出力はα以上となるのでμが連続であることから以下の式が成り立ちます。

\(\int_{-\infty}^{\infty}\phi(x)\frac{1}{\sqrt{2\pi}}exp\{-\frac{x^2}{2}\}=\alpha\)

また、検出力関数は帰無仮説において最小値を取るため以下が成り立ちます。

\(\frac{d}{d\mu} \beta(\mu_0) \\ = \int_{-\infty}^{\infty} \phi(x) \frac{1}{\sqrt{2\pi}} x \exp\left( -\frac{x^2}{2} \right) = 0\)

よって検出力関数を最大化するための式は、上の条件と合わせると

\(\int_{-\infty}^{\infty} \phi(x) \frac{1}{\sqrt{2\pi}} \exp\left( -\frac{(x-\mu)^2}{2} \right) + a\alpha + b \cdot 0 \\ = \int_{-\infty}^{\infty} \phi(x) \frac{1}{\sqrt{2\pi}} \exp\left( -\frac{(x-\mu)^2}{2} \right) – a \int_{-\infty}^{\infty} \phi(x) \frac{1}{\sqrt{2\pi}} \exp\left( -\frac{x^2}{2} \right) – b \int_{-\infty}^{\infty} \phi(x) \frac{1}{\sqrt{2\pi}} x \exp\left( -\frac{x^2}{2} \right) \\ = \int_{-\infty}^{\infty} \phi(x) \frac{1}{\sqrt{2\pi}} \left[ \exp\left( -\frac{(x-\mu)^2}{2} \right) – a \exp\left( -\frac{x^2}{2} \right) – bx \exp\left( -\frac{x^2}{2} \right) \right]\)

となるので

\(\exp\left( -\frac{(x-\mu)^2}{2} \right) \geq a \exp\left( -\frac{x^2}{2} \right) + bx \exp\left( -\frac{x^2}{2} \right)\) のとき \(\phi(x)=1\)

\(\exp\left( -\frac{(x-\mu)^2}{2} \right) < a \exp\left( -\frac{x^2}{2} \right) + bx \exp\left( -\frac{x^2}{2} \right)\) のとき \(\phi(x)=0\)

となります。

これを整理すると

\(\exp(\mu x) \geq a’ + b’x\) のとき \(\phi(x)=1\)

\(\exp(\mu x) < a' + b'x\) のとき \(\phi(x)=0\)

となり

\(x \leq c, x \geq d\)と言えます。

そうすると

\(\sqrt{n} |x| > z_{\frac{\alpha}{2}}\)

は有意水準αを満たし、かつμで微分すると0になります。

また、両側検定となっているので、これが一様最強力不偏検定であることが分かります。

次の問11は適合度検定の話でちょっとまた内容が違いそうなので、5章に戻ろうかと思います。

コメントを残す

メールアドレスが公開されることはありません。 が付いている欄は必須項目です

CAPTCHA


日本語が含まれない投稿は無視されますのでご注意ください。(スパム対策)