跳到主内容
数理基础统计推断(Statistical Inference) · 统计第1–10讲8/10 · 配套练习模块导览 →
练习

统计第8讲:最大似然估计量的大样本分布 | Large-Sample Distribution of MLEs · 练习

练习计算单参数与多参数模型的费希尔信息,推导最大似然估计量的渐近方差,并据此构造置信区间。

← 返回本讲讲义
完成状态保存在当前浏览器,可随时取消。

练习题#

1. 伯努利分布的 Fisher 信息量#

XBernoulli(p)X \sim \text{Bernoulli}(p),证明其 Fisher 信息量为 I(p)=1p(1p)I(p) = \frac{1}{p(1-p)}


2. 几何分布的 MLE 渐近方差#

X1,,XniidGeometric(p)X_1,\dots,X_n \overset{\text{iid}}{\sim} \text{Geometric}(p),最大似然估计量为 p^=1/Xˉn\hat{p} = 1/\bar{X}_n

求其渐近方差的理论表达式。


3. 正态分布均值参数的渐近分布#

X1,,XniidN(μ,σ2)X_1,\dots,X_n \overset{\text{iid}}{\sim} N(\mu, \sigma^2)σ2\sigma^2 已知。

证明样本均值 Xˉn\bar{X}_n 的渐近分布为 N(μ,σ2/n)N\left( \mu,\, \sigma^2/n \right)


4. Fisher 信息矩阵计算(二元参数)#

XN(μ,σ2)X \sim N(\mu, \sigma^2),参数向量 θ=(μ,σ)\theta = (\mu, \sigma)

计算其 Fisher 信息矩阵 I(θ)I(\theta)


5. 渐近正态性分析#

对于几何分布 Geometric(p)\text{Geometric}(p),MLE p^=1/Xˉn\hat{p} = 1/\bar{X}_n 是否满足渐近无偏性?

通过计算 p^\hat{p} 的期望说明原因。


6. 实际应用分析#

在正态分布中,当样本量 n=100n=100 时,测得样本均值 xˉ=5.2\bar{x} = 5.2,样本标准差 σ^=2.1\hat{\sigma} = 2.1

给出 μ\muσ\sigma 的近似 95% 置信区间。


查看参考答案

1. 伯努利分布的 Fisher 信息量#

解答

对数似然函数为 logf(Xp)=Xlogp+(1X)log(1p)\log f(X|p) = X\log p + (1-X)\log(1-p)

一阶导数:ddplogf(Xp)=Xp1X1p\frac{d}{dp} \log f(X|p) = \frac{X}{p} - \frac{1-X}{1-p}

二阶导数:d2dp2logf(Xp)=Xp21X(1p)2\frac{d^2}{dp^2} \log f(X|p) = -\frac{X}{p^2} - \frac{1-X}{(1-p)^2}

取期望:

E[Xp21X(1p)2]=pp21p(1p)2=1p11p=1p(1p)E\left[ -\frac{X}{p^2} - \frac{1-X}{(1-p)^2} \right] = -\frac{p}{p^2} - \frac{1-p}{(1-p)^2} = -\frac{1}{p} - \frac{1}{1-p} = -\frac{1}{p(1-p)}

因此 I(p)=1p(1p)I(p) = \frac{1}{p(1-p)}


2. 几何分布的 MLE 渐近方差#

解答

由讲义中结果,I(p)=1p2(1p)I(p) = \frac{1}{p^2(1-p)}

渐近方差为 1nI(p)=p2(1p)n\frac{1}{nI(p)} = \frac{p^2(1-p)}{n}


3. 正态分布均值参数的渐近分布#

解答

样本均值 XˉnN(μ,σ2/n)\bar{X}_n \sim N(\mu, \sigma^2/n)

根据中心极限定理,即使原分布非正态,当 nn 足够大时 Xˉn\bar{X}_n 仍近似正态分布。

此处原分布为正态,故渐近分布为精确分布。


4. Fisher 信息矩阵计算#

解答

对数似然函数:

logf(Xθ)=logσ(Xμ)22σ2+常数\log f(X|\theta) = -\log \sigma - \frac{(X-\mu)^2}{2\sigma^2} + \text{常数}

二阶偏导数:

  • 2μ2logf=1σ2\frac{\partial^2}{\partial \mu^2} \log f = -\frac{1}{\sigma^2}

  • 2μσlogf=2(Xμ)σ3\frac{\partial^2}{\partial \mu \partial \sigma} \log f = -\frac{2(X-\mu)}{\sigma^3}

  • 2σ2logf=1σ23(Xμ)2σ4\frac{\partial^2}{\partial \sigma^2} \log f = \frac{1}{\sigma^2} - \frac{3(X-\mu)^2}{\sigma^4} 取期望后:

I(θ)=[1σ2002σ2]I(\theta) = \begin{bmatrix} \frac{1}{\sigma^2} & 0 \\ 0 & \frac{2}{\sigma^2} \end{bmatrix}

5. 渐近无偏性分析#

解答

E(p^)=E(1Xˉn)>1E(Xˉn)=11/p=pE(\hat{p}) = E\left( \frac{1}{\bar{X}_n} \right) > \frac{1}{E(\bar{X}_n)} = \frac{1}{1/p} = p(由 Jensen 不等式)

因此 p^\hat{p} 是渐近无偏的(当 nn \to \infty 时偏差消失),但小样本下有正偏差。


6. 置信区间计算#

解答

  • μ\mu 的 95% CI:

    5.2±1.96×2.1100=5.2±0.412(4.788,5.612)5.2 \pm 1.96 \times \frac{2.1}{\sqrt{100}} = 5.2 \pm 0.412 \Rightarrow (4.788, 5.612)

  • σ\sigma 的近似 95% CI:

    2.1±1.96×2.1200=2.1±0.291(1.809,2.391)2.1 \pm 1.96 \times \frac{2.1}{\sqrt{200}} = 2.1 \pm 0.291 \Rightarrow (1.809, 2.391)

完成练习后,自评一次

建议先独立作答,再展开参考答案进行核对。