概率论与数理统计教程(六)-参数估计02:矩估计及相合性

本文介绍了概率论与数理统计中的矩估计方法,通过实例展示了如何利用样本矩来估计总体参数,如均值、方差和中位数。讨论了矩估计的相合性,并提供了相关分布的矩估计示例,包括指数分布和均匀分布。此外,还探讨了在不同情况下估计量的唯一性和相合性的性质。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

§ 6.2 矩估计及相合性
6.2.1 替换原理和矩法估计
1900 年 K.皮尔逊提出了一个替换原理, 后来人们称此方法为矩法.
替换原理常指如下两句话:
- 用样本矩去替换总体矩, 这里的矩可以是原点矩也可以是中心矩.
- 用样本矩的函数去替换相应的总体矩的函数.
根据这个替换原理,在总体分布形式未知场合也可对各种参数作出估计, 譬如:
- 用样本均值 xˉ\bar{x}xˉ 估计总体均值 E(X)E(X)E(X).
- 用样本方差 s2s^{2}s2 估计总体方差 Var⁡(X)\operatorname{Var}(X)Var(X).
- 用事件 AAA 出现的频率估计事件 AAA 发生的概率.
- 用样本的 ppp 分位数估计总体的 ppp 分位数, 特别,
用样本中位数估计总体中位数.
例 6.2.1 对某型号的 20 辆汽车记录其每 5 L5 \mathrm{~L}5 L 汽油的行驶里程
(单位: km\mathrm{km}km ), 观测数据如下:
29.827.628.327.930.128.729.928.027.928.7\begin{array}{llllllllll}29.8 & 27.6 & 28.3 & 27.9 & 30.1 & 28.7 & 29.9 & 28.0 & 27.9 & 28.7\end{array}29.827.628.327.930.128.729.928.027.928.7外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传{width=“420px”}
这是一个容量为 20 的样本观测值, 对应总体是该型号汽车每 5 L5 \mathrm{~L}5 L
汽油的行驶里程, 其分布形式尚不清楚, 可用矩法估计其均值、方差和中位数等.
本例中经计算有
xˉ=28.695,s2=0.9668,m0.5=28.6,\bar{x}=28.695, \quad s^{2}=0.9668, \quad m_{0.5}=28.6,xˉ=28.695,s2=0.9668,m0.5=28.6,
由此给出总体均值、方差和中位数的估计分别为 28.695,0.966828.695,0.966828.695,0.9668 和 28.6 .
矩法估计的统计思想 (替换原理) 十分简单明确, 众人都能接受, 使用场合甚广.
它的实质是用经验分布函数去替换总体分布,其理论基础是格利文科定理.
6.2.2 概率函数已知时未知参数的矩估计
设总体具有已知的概率函数
p(x;θ1,θ2,⋯ ,θk),(θ1,θ2,⋯ ,θk)∈Θp\left(x ; \theta_{1}, \theta_{2}, \cdots, \theta_{k}\right),\left(\theta_{1}, \theta_{2}, \cdots, \theta_{k}\right) \in \Thetap(x;θ1,θ2,,θk),(θ1,θ2,,θk)Θ
是未知参数或参数向量, x1,x2,⋯ ,xnx_{1}, x_{2}, \cdots, x_{n}x1,x2,,xn 是样本, 假定总体的
kkk 阶原点矩 μk\mu_{k}μk 存在, 则对所有的 jjj, 0<j<k,μj0<j<k, \mu_{j}0<j<k,μj 都存在,
若假设 θ1,θ2,⋯ ,θk\theta_{1}, \theta_{2}, \cdots, \theta_{k}θ1,θ2,,θk 能够表示成
μ1,μ2,⋯ ,μk\mu_{1}, \mu_{2}, \cdots, \mu_{k}μ1,μ2,,μk 的函数
θj=θj(μ1\theta_{j}=\theta_{j}\left(\mu_{1}\right.θj=θj(μ1,
μ2,⋯ ,μk)\left.\mu_{2}, \cdots, \mu_{k}\right)μ2,,μk), 则可给出诸 θj\theta_{j}θj
的矩估计:
θ^j=θj(a1,a2,⋯ ,ak),j=1,2,⋯ ,k,\hat{\theta}_{j}=\theta_{j}\left(a_{1}, a_{2}, \cdots, a_{k}\right), \quad j=1,2, \cdots, k,θ^j=θj(a1,a2,,ak),j=1,2,,k,
其中 a1,a2,⋯ ,aka_{1}, a_{2}, \cdots, a_{k}a1,a2,,ak 是前 kkk 阶样本原点矩
aj=1n∑i=1nxija_{j}=\frac{1}{n} \sum_{i=1}^{n} x_{i}^{j}aj=n1i=1nxij. 进一步, 如果我们要估计
θ1,θ2,⋯ ,θk\theta_{1}, \theta_{2}, \cdots, \theta_{k}θ1,θ2,,θk 的函数
η=g(θ1,θ2,⋯ ,θk)\eta=g\left(\theta_{1}, \theta_{2}, \cdots, \theta_{k}\right)η=g(θ1,θ2,,θk),
则可直接得到 η\etaη 的矩估计
η^=g(θ^1,θ^2,⋯ ,θ^k),\hat{\eta}=g\left(\hat{\theta}_{1}, \hat{\theta}_{2}, \cdots, \hat{\theta}_{k}\right),η^=g(θ^1,θ^2,,θ^k),
k=1k=1k=1 时, 我们通常可以由样本均值出发对未知参数进行估计; 如果 k=2k=2k=2,
我们可以由一阶、二阶原点矩 (或二阶中心矩) 出发估计未知参数.
例 6.2.2 设总体为指数分布, 其密度函数为
p(x;λ)=λe−λx,x⩾0,p(x ; \lambda)=\lambda \mathrm{e}^{-\lambda x}, \quad x \geqslant 0,p(x;λ)=λeλx,x0,
x1,x2,⋯ ,xnx_{1}, x_{2}, \cdots, x_{n}x1,x2,,xn 是样本, 此处 k=1k=1k=1, 由于
E(X)=1/λE(X)=1 / \lambdaE(X)=1/λ, 亦即 λ=1/E(X)\lambda=1 / E(X)λ=1/E(X), 故 λ\lambdaλ 的矩估计为
λ^=1/xˉ. \hat{\lambda}=1 / \bar{x} \text {. }λ^=1/xˉ
另外, 由于 Var⁡(X)=1/λ2\operatorname{Var}(X)=1 / \lambda^{2}Var(X)=1/λ2, 其反函数为
λ=1/Var⁡(X)\lambda=1 / \sqrt{\operatorname{Var}(X)}λ=1/Var(X) , 因此, 从替换原理来看,
λ\lambdaλ的矩估计也可取为
λ^1=1/s,\hat{\lambda}_{1}=1 / s,λ^1=1/s,
sss 为样本标准差. 这说明矩估计可能是不唯一的,
此时通常应该尽量采用低阶矩给出未知参数的估计.
例 6.2.3 设 x1,x2,⋯ ,xnx_{1}, x_{2}, \cdots, x_{n}x1,x2,

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值