青の統計学-DS Playground-

極限・漸近理論

中心極限定理、漸近分布、デルタ法、U統計量、経験分布関数、極値理論、M推定量、ノンパラメトリック推定など、統計的漸近理論の理論と応用

デルタ法による分散近似 レベル2

確率変数Xがn個の独立同分布の標本平均で、$E[X] = μ = 4、Var(X) = σ²/n = 1/100$とする。\n\n$Y = 1/X$の漸近分散をデルタ法により求めよ。小数第5位まで答えよ(小数第6位を四捨五入すること)。

解説

デルタ法(Delta Method)による非線形変換の漸近理論

デルタ法は確率変数の非線形変換の漸近分布を求める強力な手法です。統計学や計量経済学で推定量の標準誤差を計算する際に頻繁に使用されます。

デルタ法の重要性

非線形推定:複雑なパラメータの関数の分散を推定できます。実用性:最尤推定量や回帰係数の変換における標準誤差計算に不可欠です。

Step 1: デルタ法の基本理論

確率変数Xₙが漸近的に正規分布に従う場合:

$$\sqrt{n}(X_n - \mu) \xrightarrow{d} N(0, \sigma^2)$$

関数g(x)が点μで微分可能で g'(μ) ≠ 0 ならば:

$$\sqrt{n}(g(X_n) - g(\mu)) \xrightarrow{d} N(0, [g'(\mu)]^2 \sigma^2)$$

デルタ法の条件

条件内容本例での確認
漸近正規性Xₙ が漸近正規標本平均なので満足
微分可能性g(x)がμで微分可能g(x) = 1/x、μ = 4で微分可能
微分非零g'(μ) ≠ 0g'(4) = -1/16 ≠ 0

Step 2: 問題設定の整理

与えられた情報:

  • 標本平均:X ~ 漸近的に N(μ, σ²/n)
  • 平均:μ = 4
  • 分散:Var(X) = σ²/n = 1/100
  • 変換関数:g(x) = 1/x

求めるもの:Y = g(X) = 1/X の漸近分散

本問の漸近分散は、Y そのものの分散の近似値です。√n(Y − 1/μ) の漸近分散 [g'(μ)]²σ² を n で割った値にあたります。

Step 3: 微分の計算

変換関数の微分:

$$g(x) = \frac{1}{x} = x^{-1}$$
$$g'(x) = -x^{-2} = -\frac{1}{x^2}$$

μ = 4での微分値:

$$g'(4) = -\frac{1}{4^2} = -\frac{1}{16}$$

Step 4: デルタ法の適用

デルタ法により、Y = 1/X の漸近分散は:

$$Var(Y) \approx [g'(\mu)]^2 \times Var(X)$$
$$= \left(-\frac{1}{16}\right)^2 \times \frac{1}{100}$$
$$= \frac{1}{256} \times \frac{1}{100} = \frac{1}{25600}$$
$$= 0.0000390625$$

小数第5位まで:0.00004

計算の検証

各ステップの確認:

  • 微分:(1/x)' = -1/x²
  • 点での値:g'(4) = -1/16
  • 二乗:(-1/16)² = 1/256
  • 分散計算:(1/256) × (1/100) = 1/25600
  • 小数変換:1/25600 ≈ 0.0000390625

多変量デルタ法への拡張

Step 5: 一般的なデルタ法の公式

k次元確率変数ベクトル X = (X₁, ..., Xₖ)' に対して:

$$\sqrt{n}(\mathbf{X}_n - \boldsymbol{\mu}) \xrightarrow{d} N(\mathbf{0}, \boldsymbol{\Sigma})$$

関数 g: ℝᵏ → ℝ に対して:

$$\sqrt{n}(g(\mathbf{X}_n) - g(\boldsymbol{\mu})) \xrightarrow{d} N(0, \nabla g(\boldsymbol{\mu})' \boldsymbol{\Sigma} \nabla g(\boldsymbol{\mu}))$$

ここで ∇g(μ) は勾配ベクトルです。

実際の統計解析での応用例

分野変換例g(x)用途
比率推定率の逆数1/p平均到着時間
回帰分析対数オッズlog(p/(1-p))ロジスティック回帰
信頼性工学累積ハザード H(t)(瞬間ハザードは h(t)=f(t)/(1-F(t)))-log(1-F(t))故障率分析
経済学弾性値β₁ × (x̄/ȳ)需要弾性

標準誤差の計算例

Step 6: 標準誤差の導出

Y = 1/X の標準誤差:

$$SE(Y) = \sqrt{Var(Y)} = \sqrt{0.0000390625} \approx 0.00625$$

1/μ の95%信頼区間は、推定値 Y = 1/X を中心に作ります。実際のデータでは μ は未知なので、SE も g'(μ) の μ に観測値 X を代入して SE = (1/X²) × √(1/100) と評価します。

$$\frac{1}{X} \pm 1.96 \times SE(Y)$$

たとえば X = 4 が観測されたとすると、SE = (1/16) × 0.1 = 0.00625 なので:

$$0.25 \pm 1.96 \times 0.00625 = 0.25 \pm 0.01225 = [0.23775, 0.26225]$$

区間の幅は 2 × 0.01225 = 0.0245 です。

デルタ法の適用条件と限界

条件満足時違反時の対処
微分可能性通常の変換関数数値微分、ブートストラップ
微分非零単調変換高次デルタ法
漸近正規性CLTが成立変換後の分布確認
十分なサンプルサイズn ≥ 30程度(目安)有限標本調整
問題 3/20
カテゴリ一覧に戻る
問題検索