跳到主要内容

千葉大学 理学研究科 基盤理学専攻 数学・情報数理学コース 2013年8月実施 専門 B9

Author​

祭音Myyura (co-authored with GPT 6 Astra)

Description​

独立な確率変数 Xi∼N(μ,σi2)X_i\sim N(\mu,\sigma_i^2) (1≤i≤n1\le i\le n) を考える。μ∈R\mu\in\mathbb R は未知、σi>0\sigma_i>0 は既知とする。

(1) S=∑i=1nXi/σi2S=\sum_{i=1}^nX_i/\sigma_i^2 が μ\mu の十分統計量であることを示せ。

(2) SS を用いて μ\mu の不偏推定量 μ^\widehat\mu を作れ。(3) その分布を求めよ。

(4) μ\mu の不偏推定量の分散に関する Cramér–Rao 下限と μ^\widehat\mu の効率を求めよ。

(5) μ‾=n−1∑iXi\overline\mu=n^{-1}\sum_iX_i とする。相対効率 r=Var⁡(μ^)/Var⁡(μ‾)r=\operatorname{Var}(\widehat\mu)/\operatorname{Var}(\overline\mu) を求め、r≤1r\le1 を示せ。

题目描述​

独立变量 Xi∼N(μ,σi2)X_i\sim N(\mu,\sigma_i^2) 具有共同未知均值和已知正标准差。(1) 证明 S=∑iXi/σi2S=\sum_iX_i/\sigma_i^2 是充分统计量。(2) 用它构造无偏估计 μ^\widehat\mu。(3) 求估计量分布。(4) 求 Cramér–Rao 方差下界和该估计量效率。(5) 对算术平均 μ‾\overline\mu,求上述相对效率并证明不超过一。

Kai​

T=∑iσi−2T=\sum_i\sigma_i^{-2} とおく。

(1)​

同時密度は

pμ(x)=[∏i12πσiexp⁡(−12∑ixi2σi2)]exp⁡(μS−μ2T2).p_\mu(x)=\left[\prod_i\frac1{\sqrt{2\pi}\sigma_i} \exp\left(-\frac12\sum_i\frac{x_i^2}{\sigma_i^2}\right)\right] \exp\left(\mu S-\frac{\mu^2T}{2}\right).

因子分解定理により SS は十分統計量である。

(2), (3)​

E[S]=μTE[S]=\mu T、Var⁡(S)=T\operatorname{Var}(S)=T であり、独立正規変数の線形結合は正規分布に従う。よって

μ^=ST∼N(μ,1T)\boxed{\widehat\mu=\frac ST\sim N\left(\mu,\frac1T\right)}

は不偏推定量である。

(4)​

スコアは S−μTS-\mu T、Fisher 情報量は I(μ)=TI(\mu)=T である。したがって下限は 1/T1/T であり、μ^\widehat\mu はこれを達成するので効率は 11 である。

(5)​

Var⁡(μ‾)=1n2∑iσi2,r=n2(∑iσi2)(∑iσi−2).\operatorname{Var}(\overline\mu)=\frac1{n^2}\sum_i\sigma_i^2, \qquad \boxed{r=\frac{n^2}{(\sum_i\sigma_i^2)(\sum_i\sigma_i^{-2})}}.

Cauchy–Schwarz の不等式から分母は (∑iσiσi−1)2=n2(\sum_i\sigma_i\sigma_i^{-1})^2=n^2 以上である。ゆえに r≤1r\le1。等号はすべての σi\sigma_i が等しい場合に限る。