【概率论】大数定律与中心极限定理 & 参数估计

186 阅读1分钟

本文已参与「新人创作礼」活动,一起开启掘金创作之路。

大数定律与中心极限定理

一、切比雪夫不等式

设随机变量XX具有期望E(X)=μE(X)=\mu,方差D(X)=σ2D(X)=\sigma^{2},则对于任意正数ξ\xi,不等式P{XEXξ}DXξ2P\{|X-EX|\geq \xi\}\leq \frac{DX}{\xi^{2}}P{XEX<ξ}1DXξ2P\{|X-EX|< \xi\}\geq1- \frac{DX}{\xi^{2}}成立

 

二、大数定律

1. 依概率收敛定义

设随机变量X1,X2,,Xn,X_{1},X_{2},\cdots,X_{n},\cdots,对于ξ>0\forall \xi>0,都有limnP{Xna<ξ}=1\lim_{n\to \infty}P\{|X_{n}-a|<\xi\}=1,则称XnX_{n}依概率收敛于aa,记作Xn=paX_{n}\overset{p}{\rightarrow=}a

 

2. 辛钦大数定律

设随机变量X1,X2,,Xn,X_{1},X_{2},\cdots,X_{n},\cdots相互独立,且服从于统一分布,数学期望为E(Xi)=μ,i=1,2,E(X_{i})=\mu,i=1,2,\cdots,则1ni=1nXi=pμ\frac{1}{n}\sum\limits^{n}_{i=1}X_{i}\overset{p}{=}\mu。即1ni=1nf(Xi)pE(f(Xi))\frac{1}{n}\sum\limits^{n}_{i=1}f(X_{i})\overset{p}{\rightarrow}E(f(X_{i}))

 

3. 伯努利大数定律

设随机变量X1,X2,,Xn,X_{1},X_{2},\cdots,X_{n},\cdots相互独立,且均服从于B(1,p)B(1,p),则1ni=1nXipp\frac{1}{n}\sum\limits^{n}_{i=1}X_{i}\overset{p}{\rightarrow}p

 

三、中心极限定理

1. 林德贝格-列维中心极限定理

设随机变量X1,X2,,Xn,X_{1},X_{2},\cdots,X_{n},\cdots相互独立,且均服从于同一分布,E(Xi)=μ,D(Xi)=σ2E(X_{i})=\mu,D(X_{i})=\sigma^{2},则i=1nXiN(nμ,nσ2)\sum\limits^{n}_{i=1}X_{i}\overset{\cdot}{\sim}N(n \mu,n \sigma^{2})\overset{\cdot}{\sim}意为近似服从),且对于x\forall x,有limnP{i=1nXinμnσx}=Φ(x)\lim_{n \to \infty}P\Big\{\frac{\sum\limits^{n}_{i=1}X_{i}-n \mu}{\sqrt{n}\sigma}\leq x\Big\}=\Phi (x)

 

2. 棣莫弗-拉普拉斯中心极限定理

设随机变量X1,X2,,Xn,X_{1},X_{2},\cdots,X_{n},\cdots相互独立,均服从于B(1,p)B(1,p),则i=1nXiN[np,np(1p)]\sum\limits^{n}_{i=1}X_{i}\overset{\cdot}{\sim}N[np,np(1-p)]

参数估计

 

 

总体XX的概率密度(概率分布)含有未知数θ\theta

 

一、矩估计量

EX=XˉEX=\bar{X},解得θ\theta的矩估计量

 

二、最大似然估计量

1. 定义

使得似然函数取得最大值的θ\theta

 

2. 做法

  1. 对样本值x1,x2,,xnx_{1},x_{2},\cdots,x_{n},则似然函数为L(θ)={i=1nP(xi;θ)i=1nf(xi;θ)\begin{aligned}L(\theta)=\left\{\begin{aligned}&\prod^{n}_{i=1}P(x_{i};\theta)\\&\prod^{n}_{i=1}f(x_{i};\theta)\end{aligned}\right.\end{aligned}

  2. 似然函数两端取对数,求导数

  3. dln(L(θ))dθ=0\begin{aligned}\frac{d\ln(L(\theta))}{d \theta}=0\end{aligned},得θ\theta的最大似然估计量

 

例1:设XB(1,p)X\sim B(1,p)X1,X2,,XnX_{1},X_{2},\cdots,X_{n}是来自总体XX的样本,试求参数pp的最大似然估计量

x1,x2,,xnx_{1},x_{2},\cdots,x_{n}X1,X2,,XnX_{1},X_{2},\cdots,X_{n}中的样本值

XX的分布律为P{X=x}=px(1p)1x,x=0,1P\{X=x\}=p^{x}(1-p)^{1-x},x=0,1

则似然函数L(p)=i=1np(xi;p)=i=1npxi(1p)1xi=pi=1nxi(1p)i=1n(1xi)L(p)=\prod^{n}_{i=1}p(x_{i};p)=\prod^{n}_{i=1}p^{x_{i}}\cdot (1-p)^{1-x_{i}}=p^{\sum\limits^{n}_{i=1}x_{i}}\cdot (1-p)^{\sum\limits^{n}_{i=1}(1-x_{i})}

取对数ln(L(p))=i=1nxilnp+(ni=1n)ln(1p)\ln(L(p))=\sum\limits^{n}_{i=1}x_{i}\ln p+(n-\sum\limits^{n}_{i=1})\ln(1-p)

dln(L(p))dp=i=1nxipni=1nxi1p=0\frac{d\ln(L(p))}{dp}=\frac{\sum\limits^{n}_{i=1}x_{i}}{p}-\frac{n-\sum\limits^{n}_{i=1}x_{i}}{1-p}=0

解得p^=1ni=1nxi=xˉ\hat{p}=\frac{1}{n}\sum\limits^{n}_{i=1}x_{i}=\bar{x}

pp的最大似然估计量为p^=xˉ\hat{p}=\bar{x}

 

例2:设XN(μ,σ2)X\sim N(\mu,\sigma^{2})μ,σ2\mu,\sigma^{2}为未知参数,x1,x2,,xnx_{1},x_{2},\cdots,x_{n}是来自XX的一个样本值,求μ,σ2\mu,\sigma^{2}的最大似然估计量

f(x;μ,σ2)=12πσe(xμ)22σ2f(x;\mu,\sigma^{2})=\frac{1}{\sqrt{2\pi}\sigma}e^{\frac{-(x-\mu)^{2}}{2\sigma^{2}}}

L(μ,σ2)=i=1nf(xi;μ,σ2)=2πn2(σ2)n2ei=1n(xiμ)22σ2L(\mu,\sigma^{2})=\prod^{n}_{i=1}f(x_{i};\mu,\sigma^{2})=2\pi^{-\frac{n}{2}}\cdots (\sigma^{2})^{-\frac{n}{2}}\cdot e^{-\frac{\sum\limits^{n}_{i=1}(x_{i}-\mu)^{2}}{2\sigma^{2}}}

lnL(μ,σ2)=n2ln(2π)n2ln(σ2)i=1n(xiμ)22σ2\ln L(\mu,\sigma^{2})=-\frac{n}{2}\ln (2\pi)-\frac{n}{2}\ln(\sigma^{2})--\frac{\sum\limits^{n}_{i=1}(x_{i}-\mu)^{2}}{2\sigma^{2}}

lnLμ=0,lnLσ2=0\frac{\partial\ln L}{\partial\mu}=0,\frac{\partial\ln L}{\partial\sigma^{2}}=0

解得μ^=1ni=1n=xˉ,σ2^=1ni=1n(xixˉ)2\hat{\mu}=\frac{1}{n}\sum\limits^{n}_{i=1}=\bar{x},\hat{\sigma^{2}}=\frac{1}{n}\sum\limits^{n}_{i=1}(x_{i}-\bar{x})^{2}

故解得μ,σ2\mu,\sigma^{2}的最大似然估计量分别为μ^=1ni=1n=xˉ,σ2^=1ni=1n(xixˉ)2\hat{\mu}=\frac{1}{n}\sum\limits^{n}_{i=1}=\bar{x},\hat{\sigma^{2}}=\frac{1}{n}\sum\limits^{n}_{i=1}(x_{i}-\bar{x})^{2}