中心極限定理、漸近分布、デルタ法、U統計量、経験分布関数、極値理論、M推定量、ノンパラメトリック推定など、統計的漸近理論の理論と応用
確率変数Xがn個の独立同分布の標本平均で、$E[X] = μ = 4、Var(X) = σ²/n = 1/100$とする。\n\n$Y = 1/X$の漸近分散をデルタ法により求めよ。小数第5位まで答えよ(小数第6位を四捨五入すること)。
デルタ法は確率変数の非線形変換の漸近分布を求める強力な手法です。統計学や計量経済学で推定量の標準誤差を計算する際に頻繁に使用されます。
非線形推定:複雑なパラメータの関数の分散を推定できます。実用性:最尤推定量や回帰係数の変換における標準誤差計算に不可欠です。
Step 1: デルタ法の基本理論
確率変数Xₙが漸近的に正規分布に従う場合:
関数g(x)が点μで微分可能で g'(μ) ≠ 0 ならば:
| 条件 | 内容 | 本例での確認 |
|---|---|---|
| 漸近正規性 | Xₙ が漸近正規 | 標本平均なので満足 |
| 微分可能性 | g(x)がμで微分可能 | g(x) = 1/x、μ = 4で微分可能 |
| 微分非零 | g'(μ) ≠ 0 | g'(4) = -1/16 ≠ 0 |
Step 2: 問題設定の整理
与えられた情報:
求めるもの:Y = g(X) = 1/X の漸近分散
本問の漸近分散は、Y そのものの分散の近似値です。√n(Y − 1/μ) の漸近分散 [g'(μ)]²σ² を n で割った値にあたります。
Step 3: 微分の計算
変換関数の微分:
μ = 4での微分値:
Step 4: デルタ法の適用
デルタ法により、Y = 1/X の漸近分散は:
小数第5位まで:0.00004
各ステップの確認:
Step 5: 一般的なデルタ法の公式
k次元確率変数ベクトル X = (X₁, ..., Xₖ)' に対して:
関数 g: ℝᵏ → ℝ に対して:
ここで ∇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 の標準誤差:
1/μ の95%信頼区間は、推定値 Y = 1/X を中心に作ります。実際のデータでは μ は未知なので、SE も g'(μ) の μ に観測値 X を代入して SE = (1/X²) × √(1/100) と評価します。
たとえば X = 4 が観測されたとすると、SE = (1/16) × 0.1 = 0.00625 なので:
区間の幅は 2 × 0.01225 = 0.0245 です。
| 条件 | 満足時 | 違反時の対処 |
|---|---|---|
| 微分可能性 | 通常の変換関数 | 数値微分、ブートストラップ |
| 微分非零 | 単調変換 | 高次デルタ法 |
| 漸近正規性 | CLTが成立 | 変換後の分布確認 |
| 十分なサンプルサイズ | n ≥ 30程度(目安) | 有限標本調整 |