正态分布

常态分布
机率密度函数
Probability density function for the Normal distribtion
红线代表标准常态分布
累积分布函数
Cumulative distribution function for the Normal distribution
颜色与机率密度函数同
母数 数学期望(实数)
变异数(实数)
值域
机率密度函数
累积分布函数
期望值
中位数
众数
变异数
偏度 0
峰度 0
动差母函数
特征函数

常态分布(中国大陆作正态分布,香港作常态分布,澳门作常态分布,英语:Normal distribution),又名高斯分布(英语:Gaussian distribution)、正规分布,是一个非常常见的连续机率分布。常态分布在统计学上十分重要,经常用在自然社会科学来代表一个不明的随机变数。

随机变数服从一个位置母数为、尺度母数为的常态分布,记为:

[3]

则其机率密度函数[3]

常态分布的数学期望值或期望值等于位置母数,决定了分布的位置;其变异数的开平方或标准差等于尺度母数,决定了分布的幅度。

常态分布的机率密度函数曲线呈钟形,因此人们又经常称之为钟形曲线(类似于寺庙里的大钟,因此得名)。我们通常所说的标准常态分布是位置母数,尺度母数的常态分布[3](见右图中红色曲线)。

概要

常态分布是自然科学行为科学中的定量现象的一个方便模型。各种各样的心理学测试分数和物理现象比如光子计数都被发现近似地服从常态分布。尽管这些现象的根本原因经常是未知的,理论上可以证明如果把许多小作用加起来看做一个变数,那么这个变数服从常态分布(在R.N.Bracewell的Fourier transform and its application中可以找到一种简单的证明)。常态分布出现在许多区域统计:例如,采样分布均值是近似地常态的,即使被采样的样本的原始群体分布并不服从常态分布。另外,常态分布资讯熵在所有的已知均值及变异数的分布中最大,这使得它作为一种均值以及变异数已知的分布的自然选择。常态分布是在统计以及许多统计测试中最广泛应用的一类分布。在机率论,常态分布是几种连续以及离散分布的极限分布。

历史

常态分布最早是棣美弗在1718年著作的书籍的(Doctrine of Change),及1734年发表的一篇关于二项分布文章中提出的,当二项随机变数的位置母数n很大及形状母数p为1/2时,则所推导出二项分布的近似分布函数就是常态分布。拉普拉斯在1812年发表的《分析机率论》(Theorie Analytique des Probabilites)中对棣莫佛的结论作了扩展到二项分布的位置母数为n及形状母数为1>p>0时。现在这一结论通常被称为棣莫佛-拉普拉斯定理

拉普拉斯在误差分析试验中使用了常态分布。勒让德于1805年引入最小平方法这一重要方法;而高斯则宣称他早在1794年就使用了该方法,并通过假设误差服从常态分布给出了严格的证明。

「钟形曲线」这个名字可以追溯到Jouffret他在1872年首次提出这个术语「钟形曲面」,用来指代二元常态分布。常态分布这个名字还被查尔斯·皮尔士法兰西斯·高尔顿威尔赫姆·莱克希斯在1875分别独立地使用。这个术语是不幸的,因为它反映和鼓励了一种谬误,即很多机率分布都是常态的。(请参考下面的「实例」)

这个分布被称为「常态」或者「高斯」正好是Stigler名字由来法则的一个例子,这个法则说「没有科学发现是以它最初的发现者命名的」。

常态分布的定义

有几种不同的方法用来说明一个随机变数。最直观的方法是机率密度函数,这种方法能够表示随机变数每个取值有多大的可能性。累积分布函数是一种机率上更加清楚的方法,请看下边的例子。还有一些其他的等价方法,例如cumulant、特征函数动差生成函数以及cumulant-生成函数。这些方法中有一些对于理论工作非常有用,但是不够直观。请参考关于机率分布的讨论。

机率密度函数

四个不同母数集的机率密度函数(红色线代表标准常态分布)

常态分布机率密度函数均值为 变异数 (或标准差)是高斯函数的一个实例:

(请看指数函数以及.)

如果一个随机变数服从这个分布,我们写作 ~ . 如果并且,这个分布被称为标准常态分布,这个分布能够简化为

右边是给出了不同母数的常态分布的函数图。

常态分布中一些值得注意的量:

  • 密度函数关于平均值对称
  • 平均值与它的众数(statistical mode)以及中位数(median)同一数值。
  • 函数曲线下68.268949%的面积在平均数左右的一个标准差范围内。
  • 95.449974%的面积在平均数左右两个标准差的范围内。
  • 99.730020%的面积在平均数左右三个标准差的范围内。
  • 99.993666%的面积在平均数左右四个标准差的范围内。
  • 函数曲线的反曲点(inflection point)为离平均数一个标准差距离的位置。

累积分布函数

上图所示的机率密度函数的累积分布函数

累积分布函数是指随机变数小于或等于的机率,用机率密度函数表示为

常态分布的累积分布函数能够由一个叫做误差函数特殊函数表示:

标准常态分布的累积分布函数习惯上记为,它仅仅是指的值,

将一般常态分布用误差函数表示的公式简化,可得:

它的反函数被称为反误差函数,为:

该分位数函数有时也被称为probit函数。probit函数已被证明没有初等原函数。

常态分布的分布函数没有解析表达式,它的值可以通过数值积分泰勒级数或者渐进序列近似得到。

生成函数

动差母函数

动差生成函数,或称动差母函数被定义为的期望值。

常态分布的动差产生函数如下:


可以通过在指数函数内配平方得到。

特征函数

特征函数被定义为期望值,其中是虚数单位. 对于一个常态分布来讲,特征函数是:

把动差生成函数中的换成就能得到特征函数。

性质

常态分布的一些性质:

  1. 如果实数,那么 (参见期望值变异数).
  2. 如果统计独立的常态随机变数,那么:
    • 它们的和也满足常态分布 (proof).
    • 它们的差也满足常态分布.
    • 两者是相互独立的。(要求X与Y的变异数相等)
  3. 如果是独立常态随机变数,那么:
    • 它们的积服从机率密度函数为的分布
      其中是修正贝塞尔函数(modified Bessel function)
    • 它们的比符合柯西分布,满足.
  4. 如果为独立标准常态随机变数,那么服从自由度为n卡方分布

标准化常态随机变数

动差(moment

一些常态分布的一阶动差如下:

阶数 原动差 主动差 累积量
0 1 0
1 0
2
3 0 0
4 0

标准常态的所有二阶以上的累积量为零。

生成常态随机变数

中央极限定理

常态分布的机率密度函数,母数为μ = 12,σ = 3,趋近于n = 48、p = 1/4的二项分布的机率质量函数。

常态分布有一个非常重要的性质:在特定条件下,大量统计独立的随机变数的平均值的分布趋于常态分布,这就是中央极限定理。中央极限定理的重要意义在于,根据这一定理的结论,其他机率分布可以用常态分布作为近似。

  • 母数为二项分布,在相当大而且接近0.5时近似于常态分布(有的参考书建议仅在至少为5时才能使用这一近似)。

近似常态分布平均数为且变异数为.

  • 卜瓦松分布带有母数当取样样本数很大时将近似常态分布.

近似常态分布平均数为且变异数为.

这些近似值是否完全充分正确取决于使用者的使用需求

无限可分性

常态分布是无限可分的机率分布。

稳定性

常态分布是严格稳定的机率分布。

标准偏误

深蓝色区域是距平均值小于一个标准差之内的数值范围。在常态分布中,此范围所占比率为全部数值之68%,根据常态分布,两个标准差之内的比率合起来为95%;三个标准差之内的比率合起来为99%

在实际应用上,常考虑一组数据具有近似于常态分布的机率分布。若其假设正确,则约68.3%数值分布在距离平均值有1个标准差之内的范围,约95.4%数值分布在距离平均值有2个标准差之内的范围,以及约99.7%数值分布在距离平均值有3个标准差之内的范围。称为「68-95-99.7法则」或「经验法则」。

数字比率
标准差值
机率 包含之外比例
百分比 百分比 比例
0.318639σ 25% 75% 3 / 4
0.674490σ 50% 50% 1 / 2
0.994458σ 68% 32% 1 / 3.125
1σ 68.2689492% 31.7310508% 1 / 3.1514872
1.281552σ 80% 20% 1 / 5
1.644854σ 90% 10% 1 / 10
1.959964σ 95% 5% 1 / 20
2σ 95.4499736% 4.5500264% 1 / 21.977895
2.575829σ 99% 1% 1 / 100
3σ 99.7300204% 0.2699796% 1 / 370.398
3.290527σ 99.9% 0.1% 1 / 1000
3.890592σ 99.99% 0.01% 1 / 10000
4σ 99.993666% 0.006334% 1 / 15787
4.417173σ 99.999% 0.001% 1 / 100000
4.5σ 99.9993204653751% 0.0006795346249% 1 / 147159.5358
3.4 / 1000000 (每一边)
4.891638σ 99.9999% 0.0001% 1 / 1000000
5σ 99.9999426697% 0.0000573303% 1 / 1744278
5.326724σ 99.99999% 0.00001% 1 / 10000000
5.730729σ 99.999999% 0.000001% 1 / 100000000
6σ 99.9999998027% 0.0000001973% 1 / 506797346
6.109410σ 99.9999999% 0.0000001% 1 / 1000000000
6.466951σ 99.99999999% 0.00000001% 1 / 10000000000
6.806502σ 99.999999999% 0.000000001% 1 / 100000000000
7σ 99.9999999997440% 0.000000000256% 1 / 390682215445

相关分布

估计

母数的最大概似估计

概念一般化

多元常态分布共变异数矩阵的估计的推导是比较难于理解的。它需要了解谱原理(spectral theorem)以及为什么把一个标量看做一个1×1矩阵的迹(trace)而不仅仅是一个标量更合理的原因。请参考共变异数矩阵的估计(estimation of covariance matrices)。

母数的动差估计

常见实例

光子计数

计量误差

饮料装填量不足与超量的机率

某饮料公司装瓶流程严谨,每罐饮料装填量符合平均600毫升,标准差3毫升的常态分配法则。随机选取一罐,求(1)容量超过605毫升的机率;(2)容量小于590毫升的机率。

容量超过605毫升的机率 = p ( X > 605)= p ( ((X-μ) /σ) > ( (605 – 600) / 3) )= p ( Z > 5/3) = p( Z > 1.67) = 1 - 0.9525 = 0.0475

容量小于590毫升的机率 = p (X < 590) = p ( ((X-μ) /σ) < ( (590 – 600) / 3) )= p ( Z < -10/3) = p( Z < -3.33) = 0.0004

6-标准差(6-sigma或6-σ)的品质管制标准

6-标准差(6-sigma或6-σ),是制造业流行的品质管制标准。在这个标准之下,一个标准常态分配的变数值出现在正负三个标准差之外,只有2* 0.0013= 0.0026 (p (Z < -3) = 0.0013以及p(Z > 3) = 0.0013)。也就是说,这种品质管制标准的产品不良率只有万分之二十六。假设例中的饮料公司装瓶流程采用这个标准,而每罐饮料装填量符合平均600毫升,标准差3毫升的常态分配。那么预期装填容量的范围应该多少?

6-标准差的范围 = p ( -3 < Z < 3)= p ( - 3 < (X-μ) /σ < 3) = p ( -3 < (X- 600) / 3 < 3)= p ( -9 < X – 600 < 9) = p (591 < X < 609) 因此,预期装填容量应该介于591至609毫升之间。

生物标本的物理特性

金融变数

寿命

测试和智力分布

计算学生智商高低的机率

假设某校入学新生的智力测验平均分数与标准差分别为100与12。那么随机抽取50个学生,他们智力测验平均分数大于105的机率?小于90的机率?

本例没有常态分配的假设,还好中央极限定理提供一个可行解,那就是当随机样本长度超过30,样本平均数近似于一个常态变数,

因此标准常态变数

平均分数大于105的机率

平均分数小于90的机率

计算统计应用

生成常态分布随机变数

在计算机模拟中,经常需要生成常态分布的数值。最基本的一个方法是使用标准的常态累积分布函数的反函数。除此之外还有其他更加高效的方法,Box-Muller转换就是其中之一。另一个更加快捷的方法是ziggurat算法。下面将介绍这两种方法。一个简单可行的并且容易编程的方法是:求12个在(0,1)上均匀分布的和,然后减6(12的一半)。这种方法可以用在很多应用中。这12个数的和是Irwin-Hall分布;选择一个变异数12。这个随即推导的结果限制在(-6,6)之间,并且密度为12,是用11次多项式估计常态分布。

Box-Muller方法是以两组独立的随机数U和V,这两组数在(0,1]上均匀分布,用U和V生成两组独立的标准常态分布随机变数X和Y:

这个方程式的提出是因为二自由度的卡方分布(见性质4)很容易由指数随机变数(方程式中的lnU)生成。因而通过随机变数V可以选择一个均匀环绕圆圈的角度,用指数分布选择半径然后转换成(常态分布的)x,y坐标。

参考文献

  1. ^ Normal Distribution页面存档备份,存于网际网路档案馆), Gale Encyclopedia of Psychology
  2. ^ Casella & Berger (2001, p. 102)
  3. ^ 3.0 3.1 3.2 Shaou-Gang Miaou; Jin-Syan Chou. 《Fundamentals of probability and statistics》. 高立图书. 2012: 第147页. ISBN 9789864128990. 

外部连结

参见