跳到主要内容

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

Author

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

Description

α,λ\alpha,\lambda は正の実数、nnnα>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)

をもち、x0x\le0 では密度は0とする。ここで Γ(α)=0xα1exdx\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[(λlogf(X1;λ,α))2]}1\{nE[(\partial_\lambda\log f(X_1;\lambda,\alpha))^2]\}^{-1} を比較せよ。

Kai

(1)

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

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

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

(2)

V[λ^]=(k1)2E[S2]λ2=λ2k2.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[(λlogf)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)}