区间估计求解步骤:
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)X∼N(μ,σ2), 则μ\muμ的区间估计:
σ2\sigma ^2σ2已知:
由于 X‾∼N(μ,σ2n)\overline{X}\sim N(\mu ,\frac{\sigma ^2}{n})X∼N(μ,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{X−u2α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 )(X−u2α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 )(X−t2α(n−1)nSn∗,X+t2α(n−1)nSn∗)
其中:uα2σnu_\frac{\alpha}{2}\frac{\sigma}{\sqrt{n}}u2αnσ是抽样误差.
eg2:
设总体X∼N(μ,σ2)X\sim N(\mu ,\sigma ^2)X∼N(μ,σ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(n−1)(n−1)Sn∗2,χ1−2α2(n−1)(n−1)Sn∗2)
最大似然估计
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=1∏npxi(1−p)1−xi=p∑xi(1−p)n−∑xi
ln(L)=∑xilnp+(n−∑xi)ln(1−p)\ln(L)=\sum x_i \ln p+(n-\sum x_i)\ln(1-p)ln(L)=∑xilnp+(n−∑xi)ln(1−p)
∂ln(L)∂p=0\frac{\partial\ln(L)}{\partial p} =0∂p∂ln(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=1∏nxi!λxie−λ=e−nλλ∑xi/i=1∏nxi!
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
本文详细介绍了区间估计的基本步骤及应用,包括在正态分布下参数μ和σ²的置信区间计算,以及最大似然估计法在两点分布和泊松分布参数估计中的应用。

286

被折叠的 条评论
为什么被折叠?



