正态总体参数的点估计和区间估计

本文详细介绍了区间估计的基本步骤及应用,包括在正态分布下参数μ和σ²的置信区间计算,以及最大似然估计法在两点分布和泊松分布参数估计中的应用。

区间估计求解步骤:
1.构造统计量u(X1,X2⋯ ,Xn;θ)u(X_1,X_2\cdots ,X_n;\theta)u(X1,X2,Xn;θ)
2.P{c<u(X1,X2⋯ ,Xn;θ)<d}=1−αP\left \{ c< u(X_1,X_2\cdots ,X_n;\theta)<d \right \}=1-\alphaP{c<u(X1,X2,Xn;θ)<d}=1α
3.求得θ\thetaθ1−α1-\alpha1α置信区间

eg1:
设总体X∼N(μ,σ2)X\sim N(\mu ,\sigma ^2)XN(μ,σ2), 则μ\muμ的区间估计:
σ2\sigma ^2σ2已知:
由于 X‾∼N(μ,σ2n)\overline{X}\sim N(\mu ,\frac{\sigma ^2}{n})XN(μ,nσ2),故 X‾−μσ/n∼N(0,1)\frac{\overline{X}-\mu }{\sigma/\sqrt{n}}\sim N(0,1)σ/nXμN(0,1)
对于给定的α\alphaα,uα2u_\frac{\alpha}{2}u2α是标准正态分布的α2\frac{\alpha}{2}2α上侧分位数,于是
P{−uα2<X‾−μσ/n<uα2}=1−αP\left \{ -u_\frac{\alpha}{2}< \frac{\overline{X}-\mu }{\sigma/\sqrt{n}}<u_\frac{\alpha}{2} \right \}=1-\alphaP{u2α<σ/nXμ<u2α}=1α

P{X‾−uα2σn<μ<X‾+uα2σn}=1−αP\left \{ \overline{X}-u_\frac{\alpha}{2}\frac{\sigma}{\sqrt{n}}< \mu< \overline{X}+u_\frac{\alpha}{2}\frac{\sigma}{\sqrt{n}} \right \}=1-\alphaP{Xu2αnσ<μ<X+u2αnσ}=1α
μ\muμ的置信度为1−α1-\alpha1α的置信区间为(X‾−uα2σn,X‾+uα2σn)\left (\overline{X}-u_\frac{\alpha}{2}\frac{\sigma}{\sqrt{n}},\overline{X}+u_\frac{\alpha}{2}\frac{\sigma}{\sqrt{n}} \right )(Xu2αnσ,X+u2αnσ)
α=0.05\alpha=0.05α=0.05时,uα2u_\frac{\alpha}{2}u2α=1.96

σ2\sigma ^2σ2未知:
μ\muμ的置信度为1−α1-\alpha1α的置信区间为(X‾−tα2(n−1)Sn∗n,X‾+tα2(n−1)Sn∗n)\left (\overline{X}-t_\frac{\alpha}{2}(n-1)\frac{S_n^*}{\sqrt{n}},\overline{X}+t_\frac{\alpha}{2}(n-1)\frac{S_n^*}{\sqrt{n}} \right )(Xt2α(n1)nSn,X+t2α(n1)nSn)
其中:uα2σnu_\frac{\alpha}{2}\frac{\sigma}{\sqrt{n}}u2αnσ是抽样误差.

eg2:
设总体X∼N(μ,σ2)X\sim N(\mu ,\sigma ^2)XN(μ,σ2), μ\muμ, σ2\sigma^2σ2 均未知,则总体方差σ2\sigma^2σ2或标准差σ\sigmaσ的区间估计(讨论方差的区间估计可用于分析生产的稳定性与估计精度等问题):

σ2\sigma^2σ2置信度为1−α1-\alpha1α的置信区间为
((n−1)Sn∗2χα22(n−1),(n−1)Sn∗2χ1−α22(n−1))\left(\frac{(n-1){S_n^*}^2}{\chi_\frac{\alpha}{2}^2(n-1)} , \frac{(n-1){S_n^*}^2}{\chi_{1-\frac{\alpha}{2}}^2(n-1)}\right)(χ2α2(n1)(n1)Sn2,χ12α2(n1)(n1)Sn2)

最大似然估计
eg:已知X1,X2...,XnX_1,X_2...,X_nX1,X2...,Xn是从两点分布Bernoulli(1,p)中抽取出来的独立同分布样本,参数p的似然估计:
L=∏i=1npxi(1−p)1−xi=p∑xi(1−p)n−∑xiL =\displaystyle\prod_{i=1}^{n}p^{x_i}(1-p)^{1-x_i}=p^{\sum x_i}(1-p)^{n-\sum x_i}L=i=1npxi(1p)1xi=pxi(1p)nxi
ln⁡(L)=∑xiln⁡p+(n−∑xi)ln⁡(1−p)\ln(L)=\sum x_i \ln p+(n-\sum x_i)\ln(1-p)ln(L)=xilnp+(nxi)ln(1p)
∂ln⁡(L)∂p=0\frac{\partial\ln(L)}{\partial p} =0pln(L)=0
p^=Xˉ\hat p=\bar{X}p^=Xˉ

泊松分布
Possion(λ)Possion(\lambda)Possion(λ)分布参数λ\lambdaλ的极大似然估计的渐近分布求置信区间
P(X=x)=λxe−λx!P(X=x)=\frac{\lambda^x e^{-\lambda}}{x!}P(X=x)=x!λxeλ
L(λ)=∏i=1nλxie−λxi!=e−nλλ∑xi/∏i=1nxi!L(\lambda)=\displaystyle\prod_{i=1}^{n}\frac{\lambda^{x_i}e^{-\lambda}}{x_i!}=e^{-n\lambda}\lambda^{\sum x_i}/\displaystyle\prod_{i=1}^{n}x_i!L(λ)=i=1nxi!λxieλ=enλλxi/i=1nxi!
ln⁡(λ)=\ln(\lambda)=ln(λ)=
∂ln⁡(λ)∂λ=0\frac{\partial\ln(\lambda)}{\partial \lambda} =0λln(λ)=0
λ^=X‾\hat\lambda =\overline{X}λ^=X
根据CLT: n→∞,λ^∼N(λ,λ/n)n\to\infty,\hat\lambda\sim N(\lambda ,\lambda/n)n,λ^N(λ,λ/n)
λ^−λλ/n∼N(0,1)\frac{\hat\lambda-\lambda}{\sqrt{\lambda/n}}\sim N(0 ,1)λ/nλ^λN(0,1)
λ^−uα/2λ^/n<λ<λ^+uα/2λ^/n\hat\lambda-u_{\alpha/2}\sqrt{\hat\lambda/n}<\lambda<\hat\lambda+u_{\alpha/2}\sqrt{\hat\lambda/n}λ^uα/2λ^/n<λ<λ^+uα/2λ^/n

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值