跳到主要内容

大阪大学 基礎工学研究科 数理科学 (システム創成専攻) 2023年度 数理科学 I [5]

Author​

祭音Myyura (co-authored with GPT 5.6 SOL)

Description​

α,λ\alpha,\lambda は正の実数、nn は nα>2n\alpha>2 を満たす整数とする。独立な X1,…,XnX_1,\ldots,X_n は密度

f(x;λ,α)=λαΓ(α)xα−1e−λx(x>0)f(x;\lambda,\alpha)=\frac{\lambda^\alpha}{\Gamma(\alpha)}x^{\alpha-1}e^{-\lambda x}\quad(x>0)

をもち、x≤0x\le0 では密度は0とする。ここで Γ(α)=∫0∞xα−1e−x dx\Gamma(\alpha)=\int_0^\infty x^{\alpha-1}e^{-x}\,dx である。α\alpha は既知とし、λ^=(nα−1)/∑iXi\widehat\lambda=(n\alpha-1)/\sum_iX_i とおく。S=∑iXiS=\sum_iX_i の密度が f(s;λ,nα)f(s;\lambda,n\alpha) であることを用いてよい。

(1) λ^\widehat\lambda は λ\lambda の不偏推定量か。理由を述べよ。

(2) V[λ^]V[\widehat\lambda] と {nE[(∂λlog⁡f(X1;λ,α))2]}−1\{nE[(\partial_\lambda\log f(X_1;\lambda,\alpha))^2]\}^{-1} を比較せよ。

Kai​

(1)​

k=nα>2k=n\alpha>2 とする。ガンマ積分より

E[S−1]=λk−1,E[S−2]=λ2(k−1)(k−2).E[S^{-1}]=\frac\lambda{k-1},\qquad E[S^{-2}]=\frac{\lambda^2}{(k-1)(k-2)}.

したがって E[λ^]=(k−1)E[S−1]=λE[\widehat\lambda]=(k-1)E[S^{-1}]=\lambda であり、不偏である。

(2)​

V[λ^]=(k−1)2E[S−2]−λ2=λ2k−2.V[\widehat\lambda]=(k-1)^2E[S^{-2}]-\lambda^2=\boxed{\frac{\lambda^2}{k-2}}.

一方、スコアは α/λ−X1\alpha/\lambda-X_1 で平均 00、分散 α/λ2\alpha/\lambda^2。よって

{nE[(∂λlog⁡f)2]}−1=λ2nα=λ2k.\left\{nE[(\partial_\lambda\log f)^2]\right\}^{-1}=\frac{\lambda^2}{n\alpha}=\frac{\lambda^2}k.

k>2k>2 なので V[λ^]>λ2/(nα)\boxed{V[\widehat\lambda]>\lambda^2/(n\alpha)}。