正規分布

連続的変数に対する確率分布の1つ

正規分布(せいきぶんぷ、: normal distribution)またはガウス分布: Gaussian distribution)は、確率論統計学で用いられる連続的な変数に関する確率分布の一つである[1]。データが平均の付近に集積するような分布を表す。主な特徴としては平均値と最頻値中央値が一致する事や平均値を中心にして左右対称である事などが挙げられる[1][2]

正規分布
確率密度関数
正規分布の確率密度関数。赤は標準正規分布
累積分布関数
正規分布の累積分布関数:色は確率密度関数と同じ
母数(位置)
σ2 > 0 スケールの2乗(実数)
確率密度関数
累積分布関数
期待値μ
中央値μ
最頻値μ
分散σ2
歪度0
尖度0(定義によっては3)
エントロピー
モーメント母関数
特性関数
テンプレートを表示

中心極限定理により、独立な多数の因子の和として表される確率変数は正規分布に従う。このことによって正規分布は統計学や自然科学、社会科学の様々な場面で複雑な現象を簡単に表すモデルとして用いられている[1]

たとえば、実験における測定の誤差は正規分布に従って分布すると仮定され、不確かさの評価が計算されている。

正規分布の確率密度関数フーリエ変換は再び正規分布の密度関数になることから、フーリエ解析および派生した様々な数学物理の理論の体系において、正規分布は基本的な役割を果たしている。

確率変数 X が1次元正規分布に従う場合は と表記し、確率変数 Xn 次元正規分布に従う場合は などと表記する。

概要

平均μ, 分散σ2 > 0 とする(1次元)正規分布とは、確率密度関数が次の形(ガウス関数と呼ばれる)

で与えられる確率分布のことである[1][3][4]。この分布を N(μ, σ2) と表す[1][3]。(N は「正規分布」を表す英語 "normal distribution" の頭文字から取られている)[1]

標準正規分布

特に μ = 0, σ2 = 1 のとき、この分布は(1次元)標準正規分布(または基準正規分布)と呼ばれる[5]。つまり標準正規分布 N(0, 1)

なる確率密度関数を持つ確率分布として与えられる[1]

再生性

正規分布は再生性を持つ[6] —— つまり確率変数 X1, …, Xn が独立にそれぞれ正規分布 N(μ1, σ12), …, N(μn, σn2) に従うとき、線型結合 aiXi は正規分布 N(∑aiμi, ∑ai2σi2) に従う。

確率密度関数

正規分布の確率密度関数をグラフ化した正規分布曲線は左右対称な釣鐘状の曲線であり、の形に似ていることからベル・カーブ(鐘形曲線)とも呼ばれる。直線 x = μ に関して対称であり、x 軸は漸近線である。なお、曲線は σ の値が大きいほど扁平になる。

なお、中心極限定理により、巨大な n に対する二項分布とも考えることができる。

平均値の周辺の nモーメントは、各次数 n に対して

となることが知られている[7]。ただし(2n − 1)!! ≔ (2n − 1) ⋅ (2n − 3) ⋅ … ⋅ 3 ⋅ 1。 (odd: 奇数/even: 偶数)

多変量正規分布

また、多変量の統計として共分散まで込めた多次元の正規分布も定義され、平均 μ = (μ1, μ2, …, μn)n 次元正規分布の同時密度関数は次の式で与えられる。

ここで、 = (σij)分散共分散行列と呼ばれる正定値対称行列である。|Σ| は Σ の行列式。なお、A[x] は(対称)行列 A とベクトル x に対して二次形式 xTAx を意味するもの(ジーゲルの記号)とすると (xμ)T−1(xμ) = −1[xμ] と書くこともできる。

この n 次元正規分布を Nn(μ, ) と表す[8]。特に1次元の場合、平均 (μ) と分散共分散行列 = (σ2) は共に1次元の平均と分散を意味する1つの実数値であり、記号 N1((μ), ) = N1((μ), (σ2)) は単に N(μ, σ2) と書かれる(先に述べた1次元の場合の記号と同じものと理解してよい)。

歪正規分布

歪正規分布の確率密度関数

正規分布の拡張としては、上で示した多次元化を施した多変量正規分布の他に、歪正規分布 (Skew-Normal (SN) distribution) がある。これは三変数で表現され、そのうち1つの変数について α = 0 のときに正規分布となることから、分布を平均と分散の二変数で表現する正規分布の拡張であるといえる。φ(x) を標準正規分布の確率密度関数とする。

その累積確率密度関数は次で与えられる。

ここに "erf" は誤差関数シグモイド関数)である。このとき、標準正規分布に対応する歪正規分布 SN(0, 1, α) の確率密度関数は次で与えられる。

これに平均のようなもの相当する変数と分散のようなものに相当する変数を加えるためにZ変換(標準化)の逆 y = ξ + ωx を施す。すると歪正規分布は一般の形になり、以下の関係が成り立つ。

歴史

正規分布はアブラーム・ド・モアブルによって1733年に導入された[9]。この論文はド・モアブル自身による1738年出版の The Doctrine of Chances 第二版の中で、高い次数に関する二項分布の近似の文脈において再掲されている。ド・モアブルの結果はピエール=シモン・ラプラスによる『確率論の解析理論』(1812年)において拡張され、いまではド・モアブル–ラプラスの定理と呼ばれている。

ラプラスは正規分布を実験の誤差の解析に用いた。その後アドリアン=マリ・ルジャンドルによって1805年に最小二乗法が導入され[10]、1809年のカール・フリードリヒ・ガウスによる誤差論で詳細に論じられた(ガウスは1794年から最小二乗法を知っていたと主張していた)。

「ベル・カーブ」という名前は、1872年に2変数正規分布に対して「鐘形曲面」という言葉を用いた Esprit Jouffret英語版 にさかのぼる。「正規分布」という言葉はチャールズ・サンダース・パースフランシス・ゴルトンヴィルヘルム・レキシスの3人によって1875年頃に独立に導入された。

統計的な意味

標準正規分布がもつ確率密度関数のグラフ

正規分布 N(μ, σ2) からの無作為標本 x を取ると、平均 μ からのずれが ±1σ 以下の範囲に x が含まれる確率は 68.27%、±2σ 以下だと 95.45%、さらに ±3σ だと 99.73% となる[1]

正規分布は、t分布F分布といった種々の分布の考え方の基礎になっているだけでなく、実際の推測統計学においても、仮説検定区間推定など、様々な場面で利用される。

正規分布 N(μ, σ) に従う確率変数 X が与えられたとき Z = Xμ/σ と標準化すれば確率変数 Z は標準正規分布に従う。大学レベルの統計入門のクラスでは必ず行われているが、Z 値を求めることで標準正規分布表と呼ばれる変量に対応した確率を表す一覧表を用いて、コンピュータを使うことなく正規分布に従った事象の確率を求めることができる。

不連続値をとる確率変数についての検定の場合でも、連続変数と同様の考え方で正規分布を近似的に用いることがある。これは標本の大きさ n が大きく、かつデータの階級幅が狭いほど、近似の精度が高い。

標準正規分布における信頼度の推移
標準正規分布におけるσ区間の推移
信頼区間に対する信頼度の推移
信頼区間信頼度危険率
百分率百分率
0.318 639σ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/10,0000
4.5σ99.9993204653751%0.0006795346249%1/14,7159.5358
4.891638σ99.9999%0.0001%1/100,0000
5σ99.9999426697%0.0000573303%1/174,4278
5.326724σ99.99999%0.00001%1/1000,0000
5.730729σ99.999999%0.000001%1/1,0000,0000
6σ99.9999998027%0.0000001973%1/5,0679,7346
6.109410σ99.9999999%0.0000001%1/10,0000,0000
6.466951σ99.99999999%0.00000001%1/100,0000,0000
6.806502σ99.999999999%0.000000001%1/1000,0000,0000
7σ99.9999999997440%0.000000000256%1/3906,8221,5445

正規分布の適用

自然界の事象の中には正規分布に従う数量の分布をとるものがあることが知られている[11]。また、そのままでは変数が正規分布に従わない場合もその対数をとると正規分布に従う場合がある。

正規分布が統計学上特別な地位を持つのは中心極限定理が存在するためである。中心極限定理は、「独立な同一の分布に従う確率変数の算術平均(確率変数の合計を変数の数で割ったもの)の分布は、元の確率変数に標準偏差が存在するならば、元の分布の形状に関係なく、変数の数が多数になったとき、正規分布に収束する」というものである[1]。このため大標本の平均値の統計には、正規分布が仮定されることが非常に多い。

前述のごとく自然界の事象の中には、正規分布に従う数量の分布をとるものがあることが知られている。しかしそれは必ずしも多数派というわけではない。19世紀ではさながら「正規分布万能主義」といったものがまかり通っていたが、20世紀以降そういった考え方に修正が見られた。今日においては社会現象、生物集団の現象等々、種別から言えば、正規分布に従うものはむしろ少数派であることが確認されている。例えば、フラクタルな性質を持つ物は正規分布よりも、パレート分布になることが多い。人間は自然界の事象とは違って自分の意思をもっているため、たとえば、子供の成績などは決して正規分布にはならない[11]。しかし、そもそも理論上、正規分布の x の値は負の無限大から正の無限大まで取れるのに対して、多くの事象は最小値(例えば比例尺度におけるゼロ)と最大値(例えばテストにおける100点満点)が予め定まっている場合があり、そのような事象が完全な正規分布に従うとするには無理がある(その際はcensoringつまり打ち切りを考慮したり、対数正規分布を用いたりするとより正確な確率を求めることが出来る場合がある)。また、0 および自然数しかとらない離散確率分布、例えばポアソン分布二項分布連続確率分布である正規分布で近似することも一般的に行われている。

検定

正規Q-Qプロット

何らかの事象について法則性を捜したり理論を構築しようとしたりする際、その確率分布がまだ分かっていない場合にはそれが正規分布であると仮定して推論することは珍しくないが、誤った結論にたどりついてしまう可能性がある。標本データが正規分布に近似しているかどうを判断するためには、尖度歪度を調べる、ヒストグラムを見る、正規Q-Qプロットをチェックする、あるいはシャピロ–ウィルク検定コルモゴロフ–スミルノフ検定(正規分布)を利用する方法などが一般的に行われている。

点推定

平均や分散が未知の正規分布に従うデータから、母数 θ = (μ, σ2) を推定したいことがある。これには次の推定量 がよく用いられる。正規分布 N(μ, σ2) からの無作為標本 x1, …, xn が与えられたとき、

最小分散不偏推定量英語版である[12]

区間推定

点推定が1つの値を用いて母数の推定を行うのに対し、一定の区間を設けて推定することを区間推定という。

例えば、

「2022年6月の岸田内閣の支持率は59%である[13]

という推定が点推定であるのに対し、

「2022年1月から12月まで支持率は33%から59%である[13]

という推定は区間推定に分類される[14]

また、推定する区間を信頼区間と呼び、水準に応じて「90%信頼区間」「95%信頼区間」「99%信頼区間」などとも呼ばれる[15]

正規分布表

引用元:成美 清松、坂井 忠次『数理統計学要説』培風館、1952年。doi:10.11501/1371195 

標準正規分布 における確率 の値をまとめた。

Z00.010.020.030.040.050.060.070.080.09
0.0.0000.0040.0080.0120.0160.0199.0239.0279.0319.0359
0.1.0398.0438.0478.0517.0557.0596.0636.0675.0714.0753
0.2.0793.0832.0871.0910.0948.0987.1026.1064.1103.1141
0.3.1179.1217.1255.1293.1331.1368.1406.1443.1480.1517
0.4.1554.1591.1628.1664.1700.1736.1772.1808.1844.1879
0.5.1915.1950.1985.2019.2054.2088.2123.2157.2190.2224
0.6.2257.2291.2324.2357.2389.2422.2454.2486.2517.2549
0.7.2580.2611.2642.2673.2704.2734.2764.2794.2823.2852
0.8.2881.2910.2939.2967.2995.3023.3051.3078.3106.3133
0.9.3159.3186.3212.3238.3264.3289.3315.3340.3365.3389
1.0.3413.3438.3461.3485.3508.3531.3554.3577.3599.3621
1.1.3643.3665.3686.3708.3729.3749.3770.3790.3810.3830
1.2.3849.3869.3888.3907.3925.3944.3962.3980.3997.4015
1.3.4032.4049.4066.4082.4099.4115.4131.4147.4162.4177
1.4.4192.4207.4222.4236.4251.4265.4279.4292.4306.4319
1.5.4332.4345.4357.4370.4382.4394.4406.4418.4429.4441
1.6.4452.4463.4474.4484.4495.4505.4515.4525.4535.4545
1.7.4554.4564.4573.4582.4591.4599.4608.4616.4625.4633
1.8.4641.4649.4656.4664.4671.4678.4686.4693.4699.4706
1.9.4713.4719.4726.4732.4738.4744.4750.4756.4761.4767
2.0.4772.4778.4783.4788.4793.4798.4803.4808.4812.4817
2.1.4821.4826.4830.4834.4838.4842.4846.4850.4854.4857
2.2.4861.4864.4868.4871.4875.4878.4881.4884.4887.4890
2.3.4893.4896.4898.4901.4904.4906.4909.4911.4913.4916
2.4.4918.4920.4922.4925.4927.4929.4931.4932.4934.4936
2.5.4938.4940.4941.4943.4945.4946.4948.4949.4951.4952
2.6.4953.4955.4956.4957.4959.4960.4961.4962.4963.4964
2.7.4965.4966.4967.4968.4969.4970.4971.4972.4973.4974
2.8.4974.4975.4976.4977.4977.4978.4979.4979.4980.4981
2.9.4981.4982.4982.4983.4984.4984.4985.4985.4986.4986
3.0.4987.4987.4987.4988.4988.4989.4989.4989.4990.4990
3.1.4990.4991.4991.4991.4992.4992.4992.4992.4993.4993
3.2.4993.4993.4994.4994.4994.4994.4994.4995.4995.4995
3.3.4995.4995.4995.4996.4996.4996.4996.4996.4996.4997
3.4.4997.4997.4997.4997.4997.4997.4997.4997.4997.4998
3.5.4998.4998.4998.4998.4998.4998.4998.4998.4998.4998
3.6.4998.4998.4999.4999.4999.4999.4999.4999.4999.4999
3.7.4999.4999.4999.4999.49991.49992.49992.49992.49992.49992
3.8.49993.49993.49993.49994.49994.49994.49994.49995.49995.49995
3.9.49995.49995.49996.49996.49996.49996.49996.49996.49997.49997
4.0.49997.49997.49997.49997.49997.49997.49997.49997.49997.49997
4.1.49998.49998.49998.49998.49998.49998.49998.49998.49998.49998
4.2.49999.49999.49999.49999.49999.49999.49999.49999.49999.49999
4.3.49999.49999.49999.49999.49999.49999.49999.49999.49999.49999
4.4.49999.49999.49999.49999.49999.49999.49999.49999.49999.49999
4.5.49997.49997.49997.49997.49997.49997.49997.49997.49997.49997
4.6.49998.49998.49998.49998.49998.49998.49998.49998.49998.49998
4.7.49999.49999.49999.49999.49999.49999.49999.49999.49999.49999
4.8.49999.49999.49999.49999.49999.49999.49999.49999.49999.49999
4.9.499995.499995.499995.499995.499995.499995.499995.499995.499995.499995
5.0.499997

脚注

出典

参考文献

関連項目

外部リンク