跳到主要内容

法政大学 理工学研究科 システム理工学専攻 経営システム系 2020年8月実施 概率统计

Author

思齐塾, 祭音Myyura

Description

(X1,,Xn)(X_1, \dots, X_n) は正規母集団からの無作為標本で、平均は同じだが分散が異なる。 i=1,,ni=1, \dots, n に対し XiN(μ,σi2)X_i \sim N(\mu, \sigma_i^2) で分散を既知としたとき、以下の問いに答えよ。

(1) μ\mu の最尤推定量 μ^\hat{\mu}

μ^=1i=1nσi2i=1nXiσi2\hat{\mu} = \frac{1}{\sum_{i=1}^n \sigma_i^{-2}} \sum_{i=1}^n \frac{X_i}{\sigma_i^2}

で表されることを示せ。

(2) 前問で求めた最尤推定量が不偏性を満たすか否か示せ。

题目描述

(X1,,Xn)(X_1,\dots,X_n) 是取自正态总体的随机样本,各变量均值相同而方差不同。对 i=1,,ni=1,\dots,n,设

XiN(μ,σi2),X_i\sim N(\mu,\sigma_i^2),

其中各方差 σi2\sigma_i^2 已知。回答下列问题。

(1)证明 μ\mu 的最大似然估计量 μ^\hat\mu 可表示为

μ^=1i=1nσi2i=1nXiσi2.\hat\mu =\frac{1}{\sum_{i=1}^n\sigma_i^{-2}} \sum_{i=1}^n\frac{X_i}{\sigma_i^2}.

(2)说明上一问所得最大似然估计量是否具有无偏性。

Kai

(1) 尤度関数は

L(μ)=i=1n12πσi2exp((Xiμ)22σi2)L(\mu) = \prod_{i=1}^n \frac{1}{\sqrt{2\pi\sigma_i^2}} \exp \left( -\frac{(X_i - \mu)^2}{2\sigma_i^2} \right)

対数尤度関数は

logL(μ)=i=1n[12log(2πσi2)(Xiμ)22σi2]\log L(\mu) = \sum_{i=1}^n \left[ -\frac{1}{2} \log(2\pi\sigma_i^2) - \frac{(X_i - \mu)^2}{2\sigma_i^2} \right]

μ\mu で微分して0とおくと、

logL(μ)μ=i=1nXiμσi2=0\frac{\partial \log L(\mu)}{\partial \mu} = \sum_{i=1}^n \frac{X_i - \mu}{\sigma_i^2} = 0
i=1nXiσi2μi=1n1σi2=0\sum_{i=1}^n \frac{X_i}{\sigma_i^2} - \mu \sum_{i=1}^n \frac{1}{\sigma_i^2} = 0
μ^=i=1nXiσi2i=1n1σi2=1i=1nσi2i=1nXiσi2\hat{\mu} = \frac{\sum_{i=1}^n \frac{X_i}{\sigma_i^2}}{\sum_{i=1}^n \frac{1}{\sigma_i^2}} = \frac{1}{\sum_{i=1}^n \sigma_i^{-2}} \sum_{i=1}^n \frac{X_i}{\sigma_i^2}

(2) 不偏性について

E[μ^]=E[1i=1nσi2i=1nXiσi2]E[\hat{\mu}] = E\left[ \frac{1}{\sum_{i=1}^n \sigma_i^{-2}} \sum_{i=1}^n \frac{X_i}{\sigma_i^2} \right]
=1i=1nσi2i=1nE[Xi]σi2= \frac{1}{\sum_{i=1}^n \sigma_i^{-2}} \sum_{i=1}^n \frac{E[X_i]}{\sigma_i^2}
=1i=1nσi2i=1nμσi2= \frac{1}{\sum_{i=1}^n \sigma_i^{-2}} \sum_{i=1}^n \frac{\mu}{\sigma_i^2}
=μi=1nσi2i=1nσi2=μ= \mu \frac{\sum_{i=1}^n \sigma_i^{-2}}{\sum_{i=1}^n \sigma_i^{-2}} = \mu

よって、 μ^\hat{\mu}μ\mu の不偏推定量である。