【概率论】6-3:中心极限定理(The Central Limit Theorem)
原文地址:https://www.face2ai.com/Math-Probability-6-3-The-Central-Limit-Theorem转载请标明出处
Abstract: 本文介绍中心极限定理
Keywords: The Central Limit Theorem,The Normal distribution,The Delta Method
中心极限定理
读书的一个重要用途就是建立自己对事情的理解方法,在数学领域,尤其是概率和数理统计,学习这两门课程,可以让你对世界上所有的事情的理解改变一个角度,甚至统计最后可以解释哲学,那么这样解释自然的三种方法——神学,哲学,科学,就会被改成“神学”和“科学”了,如果哪天神学也被建模了,哈哈哈,世界大一统,这里并不是对宗教或者哲学家的任何不尊重,只是谈一个小想法😆
概率论快接近尾声了,本文讲完基本也就剩下一篇了,学了微积分,学了线性代数,学了数学分析,学了概率论,但是感觉自己还是没什么长进,这就是数学基础的困难,变现速度慢,但是作为长久投资,我相信是值得的。
本文我们介绍中心极限定理,上一篇的大数定理,围绕的一个核心观点就是,样本均值概率极限是分布均值。而今天的中心极限定理是描述样本均值的分布的。一个数量足够多的随机变量样本的样本期望(Sample Mean,一个随机变量)的期望是 μ \mu μ 以及有限的方差 σ 2 \sigma^2 σ2 那么他的分布近似于均值为 μ \mu μ 方差为 σ 2 / n \sigma^2/n σ2/n 的正态分布。这个结论可以帮助证明正态分布可以用来建模一些随机变量,当然这些随机变量是多个独立的部分组成的。
我们本文会提出Delta方法来帮助我们一组随机变量样本的期望的函数变化后产生的新随机变量。
中心极限定理和大数定理放在后面是有目的的,因为我们马上要过渡到数理统计,而前面我们说了,我们的主要应用时根据数据(已知样本)找模型,也就是我们在数理统计里面要研究的,而中心极限定理和大数定理给出的就是样本和模型之间的理论关系。
中心极限定理 Statement of the Theorem
先举个🌰 :
临床试验,100个患者将接受治疗,患者有 0.5的可能性接受治疗,而其本人不知道自己是否接受治疗,我们这里假设所有患者之间相互独立,本实验的目的是观察新的治疗能否提高患者的治愈率,随机变量
X
X
X 表示 100 个人中治愈的人数,如果接受治疗的治愈率是0.5和没有接受治疗治愈率的一样,那么
X
X
X 就是
n
=
100
,
p
=
0.5
n=100,p=0.5
n=100,p=0.5 的二项分布,那么根据二项分布的随机变量的数值,可以描绘出如下的图片
如果我们绘制一个均值为
μ
=
50
\mu=50
μ=50 方差为
σ
2
=
25
\sigma^2=25
σ2=25 的正态分布,刚好可以把这个二项分布包裹起来
在Poisson分布的学习中,当我们遇到
n
n
n 很大
p
p
p 很小,而
n
p
np
np 有不大不小(这个标准有点含糊) 的时候可以用泊松分布来近似二项分布,上面的例子我们看出,当
n
n
n 很大而
p
p
p 又不太小的时候,正态分布可以用来近似二项分布。
中心极限定理就是一个形式化的描述,来描述正态分布是如何近似i.i.d的随机变量的一般和(随机变量)或者均值(样本均值,随机变量)的分布的。
在正态分布的那一课证明了
n
n
n 个独立同正态分布的的随机变量,如果其分布均值是
μ
\mu
μ 分布方差是
σ
2
\sigma^2
σ2 那么样本均值
X
ˉ
n
\bar{X}_n
Xˉn 也是正态分布,均值是
μ
\mu
μ 方差是
σ
2
/
n
\sigma^2/n
σ2/n 而我们本篇就会介绍根据中心极限定理,无论原始分布是否是正态分布,其样本均值分布都可以用正态分布来近似,样本均值。
而对于样本是独立的不同分布的随机变量的话,我们也有一套中心极限定理能够建模这些样本的和。
注意,中心极限定理,有两个版本,就像微积分基本定理有两个一样,中心极限定理,一个针对独立同分布的样本,一个针对独立不同分布的样本,大家要注意区分。
首先来看同分布的,这个定理由 Lindeberg 和Levy 在1920s初期证明:
Theorem Central Limit Theorem(Lindeberg and Levy) If the random variables X 1 , … , X n X_1,\dots,X_n X1,…,Xn form a random sample of size n n n from a given distribution with mean μ \mu μ and variance σ 2 \sigma^2 σ2 ( 0 < σ 2 < ∞ 0<\sigma^2<\infty 0<σ2<∞ ) ,then for each fixed number x x x ,
l i m n → ∞ P r [ X ˉ n − μ σ n 1 / 2 ≤ x ] = Φ ( x ) lim_{n\to \infty}Pr[\frac{\bar{X}_n-\mu}{\frac{\sigma}{n^{1/2}}}\leq x]=\Phi(x) limn→∞Pr[n1/2σXˉn−μ≤x]=Φ(x)
where Φ \Phi Φ denotes the c.d.f. of the standard normal distribution.
书上并没有给出上述定理的证明,而是给出了定义的解释说明(我觉得还是证明了的定理用的踏实,但是有些定理证明过于复杂,所以理解定理含义,也算是站在巨人的肩膀上,但还是证明一下用的踏实。),定理说明是,从任何一个分布中选取大量的随机变量,这个分布的均值是 μ \mu μ 方差是 σ 2 \sigma^2 σ2 不需要考虑这个分布式离散的还是连续的,则 n 1 / 2 ( X ˉ n − μ ) / σ n^{1/2}(\bar{X}_n-\mu)/\sigma n1/2(Xˉn−μ)/σ 近似于标准正太分布;那么根据正态分布的性质,样本均值 X ˉ \bar{X} Xˉ 近似于 μ X ˉ = μ , σ X ˉ 2 = σ 2 / n \mu_{\bar{X}}=\mu,\sigma^2_{\bar{X}}=\sigma^2/n μXˉ=μ,σXˉ2=σ2/n 的正态分布; ∑ i = 1 n X i \sum^{n}_{i=1}X_i ∑i=1nXi 的分布近似于 μ Σ = n μ , σ Σ 2 = n σ 2 \mu_{\Sigma}=n\mu,\sigma^2_{\Sigma}=n\sigma^2 μΣ=nμ,σΣ2=nσ2 的正态分布
关于中心极限定理的证明,我们在后面的高级课程中给出,这里我们需要做到的是理解定理内容。
又一个🌰 :
从一个均匀分布中采样,假设从
[
0
,
1
]
[0,1]
[0,1] 的均匀分布重采样得到
n
=
12
n=12
n=12 个样本,求
P
r
(
∣
X
ˉ
n
−
1
2
∣
≤
0.1
)
Pr(|\bar{X}_n-\frac{1}{2}|\leq 0.1)
Pr(∣Xˉn−21∣≤0.1) 的大概值。
解:
- [ 0 , 1 ] [0,1] [0,1] 之间的均匀分布的期望是 1 2 \frac{1}{2} 21
- 方差是 σ 2 = 1 12 \sigma^2=\frac{1}{12} σ2=121
- 这里 n = 12 n=12 n=12
- X ˉ n \bar{X}_n Xˉn 根据中心极限定理,其分布近似为 μ X ˉ n = 1 2 , σ X ˉ n 2 = 1 144 \mu_{\bar{X}_n}=\frac{1}{2},\sigma_{\bar{X}_n}^2=\frac{1}{144} μXˉn=21,σXˉn2=1441
- 根据中心极限定理,设随机变量
Z
=
12
(
X
ˉ
n
−
1
2
)
Z=12(\bar{X}_n-\frac{1}{2})
Z=12(Xˉn−21) 近似于标准正太分布,于是有:
P r ( ∣ X ˉ n − 1 2 ∣ ≤ 0.1 ) = P r [ 12 ∣ X ˉ n − 1 2 ∣ ≤ 1.2 ] = P r ( ∣ Z ∣ ≤ 1.2 ) ≈ 2 Φ ( 1.2 ) − 1 = 0.7698 \begin{aligned} Pr(|\bar{X}_n-\frac{1}{2}|\leq 0.1)&=Pr[12|\bar{X}_n-\frac{1}{2}|\leq 1.2]\\ &=Pr(|Z|\leq 1.2)\\ &\approx 2\Phi(1.2)-1\\ &=0.7698 \end{aligned} Pr(∣Xˉn−21∣≤0.1)=Pr[12∣Xˉn−21∣≤1.2]=Pr(∣Z∣≤1.2)≈2Φ(1.2)−1=0.7698
上面是书上的答案,但是我看有点问题就是第5步,根据中心极限定理 Z = 1 2 1 / 2 ( X ˉ n − 1 2 ) 1 12 = 24 3 ( X ˉ n − 1 2 ) Z=\frac{12^{1/2}(\bar{X}_n-\frac{1}{2})}{\frac{1}{12}}=24\sqrt{3}(\bar{X}_n-\frac{1}{2}) Z=121121/2(Xˉn−21)=243(Xˉn−21) 但不知道为啥直接被变成12了。
下面我们来介绍在中心极限定理里面出现的收敛在其他过程中的应用,以及他有自己的名字。
Definition Convergence in Distribution/Asymptotic Distribution.Let X 1 , X 2 , … X_1,X_2,\dots X1,X2,… be a sequence of random variables,and for n = 1 , 2 , … n=1,2,\dots n=1,2,… ,let F n F_n Fn denote the c.d.f. of X n X_n Xn .Also,let F ′ F' F′ be a c.d.f. Then it is said that the sequence X 1 , X 2 , … X_1,X_2,\dots X1,X2,… converges in distribution to F ′ F' F′ if
l i m n → ∞ F n ( x ) = F ′ ( x ) lim_{n\to \infty}F_n(x)=F'(x) limn→∞Fn(x)=F′(x)
for all x at which F ′ ( x ) F'(x) F′(x) is continuous.Sometimes,it is simply said that X n X_n Xn converges in distribution to F ′ F' F′ ,and F ′ F' F′ is called the asymptotic distribution of X_n,If F ′ F' F′ has a name ,then we say that X_n converges in distribution to that name
定义收敛到分布,或者渐进分布:一个随机变量序列 X 1 , X 2 , … X_1,X_2,\dots X1,X2,… 对应的c.d.f. 为 F n F_n Fn 其中 n = 1 , 2 … n=1,2\dots n=1,2… 并且 F ′ F' F′ 是一个c.d.f. 那么当 l i m n → ∞ F n ( x ) = F ′ ( x ) lim_{n\to \infty}F_n(x)=F'(x) limn→∞Fn(x)=F′(x) 时,我们说随机变量序列 X 1 , X 2 , … X_1,X_2,\dots X1,X2,… 的分布收敛到 F ′ F' F′ 。 x x x 对于所有 F ′ ( x ) F'(x) F′(x) 都是连续的,有时候简单的说 X n X_n Xn 收敛到 F ′ F' F′ ,如果这个 F ′ F' F′ 有自己的名字,比如正态分布,我们就说 X n X_n Xn 收敛到正态分布。也可以说 X n X_n Xn 的渐进分布是正态分布。
根据上面的定义,我们可以说 n 1 / 2 ( X ˉ n ) − μ σ \frac{n^{1/2}(\bar{X}_n)-\mu}{\sigma} σn1/2(Xˉn)−μ 收敛到标准正态分布。
中心极限定理也解释了为什么自然界中的大量数据都收敛到正态分布,比如人的身高,体重,收集足够的数据后会发现其收敛到正态分布。
Delta 方法 The Delta Method
🌰 :
服务时间问题,顾客们在排队,假设第
i
i
i 个顾客到达柜台至完成服务用时
X
i
X_i
Xi ,我们会得到
X
1
,
X
2
,
…
X_1,X_2,\dots
X1,X2,… 共
n
n
n 个随机变量,我们假设这些随机变量独立同分布,并且分布均值为
μ
\mu
μ 方差为
σ
2
\sigma^2
σ2 那么我们可以根据中心极限定理给出
X
ˉ
n
\bar{X}_n
Xˉn 的分布,但是如果我们关心服务速度,也就是
1
X
ˉ
n
\frac{1}{\bar{X}_n}
Xˉn1 那么我们改如何建模呢?
例子引出的问题,独立同分布随机变量样本期望可以用中心极限定理建模,如果这个期望经过某个函数变换呢?比如上面例子中的取倒数,如何对这个新的随机变量建模呢?
统计中一个已知的方法叫做 Delta方法,能够帮我们解决这个问题。
Theorem Delta Method. Let Y 1 , Y 2 , … Y_1,Y_2,\dots Y1,Y2,… be a sequence of random variables,and let F ′ F' F′ be a continuous c.d.f. Let θ \theta θ be a real number,and let α 1 , α 2 , … \alpha_1,\alpha_2,\dots α1,α2,… be a sequence of positive numbers that increase to ∞ \infty ∞ .Suppose that a n ( Y n − θ ) a_n(Y_n-\theta) an(Yn−θ) converges in distribution to F ′ F' F′ . Let α \alpha α be a function with contimuous derivative such that α ′ ( θ ) ≠ 0 \alpha'(\theta)\neq 0 α′(θ)=0 .Then a n [ α ( Y n ) − α ( θ ) ] / α ′ ( θ ) a_n[\alpha(Y_n)-\alpha(\theta)]/\alpha'(\theta) an[α(Yn)−α(θ)]/α′(θ) converges in distribution to F ′ F' F′
一个随机变量序列的线性变换收敛到某分布,那么这个随机变量序列的某函数映射后的随机变量的线性变换依旧收敛到这个分布。
证明:
这个过程相对粗糙相当于给出了证明框架。
-
a n → ∞ a_n\to\infty an→∞ 的时候 Y n Y_n Yn 会概率收敛到 θ \theta θ 否则 ∣ a n ( Y n − θ ) ∣ → ∞ |a_n(Y_n-\theta)|\to \infty ∣an(Yn−θ)∣→∞ ,也就是不收敛。
-
因为 α \alpha α 是连续函数, α ( Y n ) \alpha(Y_n) α(Yn) 将会非常接近 α ( θ ) \alpha(\theta) α(θ)
-
使用泰勒级数(忽略更高阶项):
α ( Y n ) ≈ α ( θ ) + α ′ ( θ ) ( Y n − θ ) \alpha(Y_n)\approx \alpha(\theta)+\alpha'(\theta)(Y_n-\theta) α(Yn)≈α(θ)+α′(θ)(Yn−θ) -
两边同时减去 α ( θ ) \alpha(\theta) α(θ) 后乘以 a n α ′ ( θ ) \frac{a_n}{\alpha'(\theta)} α′(θ)an 得到:
a n α ′ ( θ ) ( Y n − θ ) ≈ a n ( Y n − θ ) \frac{a_n}{\alpha'(\theta)}(Y_n-\theta)\approx a_n(Y_n-\theta) α′(θ)an(Yn−θ)≈an(Yn−θ) -
4中结论左右两边都收敛到 F ′ F' F′
这个定理给出的方法最常见的一个就是 Y n Y_n Yn 是随机样本的平均数,随机样本的方差有限,可以得出以下推论。
Corollary Delta Method for Average of a Ramdom Sample.Let X 1 , X 2 … X_1,X_2\dots X1,X2… be a sequence of i.i.d. random variables from a distribution with mean μ \mu μ and finite variance σ 2 \sigma^2 σ2 .Let α \alpha α be a function with continuous derivative such that α ′ ( μ ) ≠ 0 \alpha'(\mu)\neq 0 α′(μ)=0 .Then the asymptotic distribution of
n 1 / 2 σ α ′ ( μ ) [ α ( X ˉ n ) − α ( μ ) ] \frac{n^{1/2}}{\sigma\alpha'(\mu)}[\alpha(\bar{X}_n)-\alpha(\mu)] σα′(μ)n1/2[α(Xˉn)−α(μ)]
is the standard normal distribution.
证明方法就是使用Delta法则,就能得到结果。
总结
本文介绍了中心极限定理的基础知识,因为本定理在概率论和数理统计中的超级低位,希望大家能仔细阅读文章内容。
下一篇我们对本章进行收尾,同时结束概率论基础部分。
更多推荐
所有评论(0)