統計的推測(推定)

最尤推定、ベイズ推定、区間推定など統計検定準1級レベルの統計的推定理論を学習します。

フィッシャー情報量の計算 レベル1

正規分布$N(\mu, \sigma^2)$($\mu$既知)からの標本$X_1, \ldots, X_n$について、$\sigma^2$に関するフィッシャー情報量$I(\sigma^2)$を求めよ。答えを$n$と$\sigma$を用いて表せ。

解説
解答と解説を表示
<h4>フィッシャー情報量の定義と意義</h4><p>フィッシャー情報量は、データから得られるパラメータに関する情報の量を定量化する重要な概念です。</p><h4>フィッシャー情報量の理論的基礎</h4><p class='step'><strong>Step 1: フィッシャー情報量の定義</strong></p><p>パラメータ$\theta$に関するフィッシャー情報量は以下で定義されます:</p><div class='formula'>$I(\theta) = E\left[\left(\frac{\partial \log f(X; \theta)}{\partial \theta}\right)^2\right] = -E\left[\frac{\partial^2 \log f(X; \theta)}{\partial \theta^2}\right]
lt;/div><p>第2の等号は正則条件下で成り立ちます。</p><p class='step'><strong>Step 2: 正規分布の確率密度関数</strong></p><p>正規分布$N(\mu, \sigma^2)$($\mu$既知)の確率密度関数:</p><div class='formula'>$f(x; \sigma^2) = \frac{1}{\sqrt{2\pi\sigma^2}} \exp\left(-\frac{(x-\mu)^2}{2\sigma^2}\right)
lt;/div><h4>対数尤度関数の導出</h4><p class='step'><strong>Step 3: 対数尤度関数</strong></p><p>単一観測値$X$に対する対数尤度:</p><div class='formula'>$\ell(\sigma^2) = \log f(X; \sigma^2) = -\frac{1}{2}\log(2\pi) - \frac{1}{2}\log(\sigma^2) - \frac{(X-\mu)^2}{2\sigma^2}
lt;/div><p class='step'><strong>Step 4: 一階微分(スコア関数)</strong></p><p>$\sigma^2$に関する一階微分:</p><div class='formula'>$\frac{\partial \ell}{\partial \sigma^2} = -\frac{1}{2\sigma^2} + \frac{(X-\mu)^2}{2(\sigma^2)^2}
lt;/div><div class='formula'>$= \frac{1}{2\sigma^2}\left[-1 + \frac{(X-\mu)^2}{\sigma^2}\right]
lt;/div><p class='step'><strong>Step 5: 二階微分</strong></p><p>$\sigma^2$に関する二階微分:</p><div class='formula'>$\frac{\partial^2 \ell}{\partial (\sigma^2)^2} = \frac{1}{2\sigma^4} - \frac{(X-\mu)^2}{(\sigma^2)^3}
lt;/div><div class='formula'>$= \frac{1}{2\sigma^4}\left[1 - \frac{2(X-\mu)^2}{\sigma^2}\right]
lt;/div><h4>期待値の計算</h4><p class='step'><strong>Step 6: 一階微分の期待値</strong></p><p>スコア関数の期待値(正則条件により0):</p><div class='formula'>$E\left[\frac{\partial \ell}{\partial \sigma^2}\right] = \frac{1}{2\sigma^2}\left[-1 + \frac{E[(X-\mu)^2]}{\sigma^2}\right] = \frac{1}{2\sigma^2}\left[-1 + \frac{\sigma^2}{\sigma^2}\right] = 0
lt;/div><p class='step'><strong>Step 7: 一階微分の二乗の期待値</strong></p><div class='formula'>$E\left[\left(\frac{\partial \ell}{\partial \sigma^2}\right)^2\right] = E\left[\left(\frac{1}{2\sigma^2}\left[-1 + \frac{(X-\mu)^2}{\sigma^2}\right]\right)^2\right]
lt;/div><div class='formula'>$= \frac{1}{4\sigma^4} E\left[\left(-1 + \frac{(X-\mu)^2}{\sigma^2}\right)^2\right]
lt;/div><p>$(X-\mu)^2/\sigma^2 \sim \chi^2_1$なので、$E[(X-\mu)^2/\sigma^2] = 1$、$\text{Var}[(X-\mu)^2/\sigma^2] = 2
lt;/p><p>したがって:</p><div class='formula'>$E\left[\left(-1 + \frac{(X-\mu)^2}{\sigma^2}\right)^2\right] = \text{Var}\left[\frac{(X-\mu)^2}{\sigma^2}\right] = 2
lt;/div><div class='formula'>$E\left[\left(\frac{\partial \ell}{\partial \sigma^2}\right)^2\right] = \frac{1}{4\sigma^4} \cdot 2 = \frac{1}{2\sigma^4}
lt;/div><h4>代替計算:二階微分による方法</h4><p class='step'><strong>Step 8: 二階微分の期待値</strong></p><div class='formula'>$E\left[\frac{\partial^2 \ell}{\partial (\sigma^2)^2}\right] = \frac{1}{2\sigma^4}\left[1 - \frac{2E[(X-\mu)^2]}{\sigma^2}\right]
lt;/div><div class='formula'>$= \frac{1}{2\sigma^4}\left[1 - \frac{2\sigma^2}{\sigma^2}\right] = \frac{1}{2\sigma^4}(1 - 2) = -\frac{1}{2\sigma^4}
lt;/div><p>したがって:</p><div class='formula'>$I(\sigma^2) = -E\left[\frac{\partial^2 \ell}{\partial (\sigma^2)^2}\right] = -\left(-\frac{1}{2\sigma^4}\right) = \frac{1}{2\sigma^4}
lt;/div><h4>標本に対するフィッシャー情報量</h4><p class='step'><strong>Step 9: $n$個の独立標本</strong></p><p>独立な標本$X_1, \ldots, X_n$に対して:</p><div class='formula'>$I_n(\sigma^2) = n \cdot I(\sigma^2) = n \cdot \frac{1}{2\sigma^4} = \frac{n}{2\sigma^4}
lt;/div><div class='key-point'><div class='key-point-title'>フィッシャー情報量の解釈</div><ul><li><strong>情報の加法性</strong>:独立観測値の情報量は加算される</li><li><strong>パラメータの精度</strong>:$I(\theta)$が大きいほど$\theta$の推定精度が高い</li><li><strong>クラメール・ラオ下界</strong>:不偏推定量の分散の下界は$1/I(\theta)
lt;/li><li><strong>漸近有効性</strong>:最尤推定量は漸近的に$1/I(\theta)$の分散を持つ</li></ul></div><h4>実用的応用と検証</h4><p class='step'><strong>Step 10: 結果の検証</strong></p><p>$\sigma^2$の最尤推定量$\hat{\sigma}^2 = \frac{1}{n}\sum_{i=1}^n (X_i - \mu)^2$について:</p><div class='formula'>$\text{Var}(\hat{\sigma}^2) = \text{Var}\left(\frac{\sigma^2}{n} \chi^2_n\right) = \frac{\sigma^4}{n^2} \cdot 2n = \frac{2\sigma^4}{n}
lt;/div><p>クラメール・ラオ下界:$\frac{1}{I_n(\sigma^2)} = \frac{2\sigma^4}{n}
lt;/p><p>一致しており、最尤推定量が漸近有効であることが確認できます。</p><p>したがって、$\sigma^2$に関するフィッシャー情報量は$I(\sigma^2) = \frac{n}{2\sigma^4}$です。</p>
問題 1/10
カテゴリ一覧に戻る