メインコンテンツまでスキップ

切断正規分布(Truncated normal distribution)

腕をバッサリ切り落とし切断正規分布​

上限と下限​

さて、突然ですが袋に30グラムと書いてある某メーカーヒット商品のスナック菓子がここにあります。
突撃訪問にて工場で次々と作られる製品の重さを量ってみるとピッタリ30グラムなわけもなく、29.98グラムのものもあれば、30.12グラムのものもあったりします。 もちろん30グラム周辺にデータは固まっていて、極端に多すぎるものや少なすぎるものは少なくなっています。そうです。ピンときましたか？これはまさに「正規分布」をなす典型的な場合だったのです！（テンション高めに左側の図をご覧ください。）

では、忘れていましたが、我々の手元にあるこの製品の重さはどうなっているのでしょうか。
早速大人買いをして、重さの分布を見てみることにしましょう。さすが30グラムと書いてあるだけのことはあって、大体どれも30グラムに近い値ですが、ピッタリ30グラムなわけもなく、29.98グラムのものもあれば、30.12グラムのものもありますね（ま、当たり前なんですが）。 こうして努力と不意の出費の結果出来あがった分布図は右側のものになりました。
ここで左と異なる点がはっきりします。我々の手元に届いたという条件の下では 30±330\pm 3グラムより誤差が大きいものが存在しないのです！すると分布は裾をバッサリ切り取られたような形になってしまいます。 しかし存在範囲では依然、正規分布であることは間違いないのです。
そう、この分布こそが切断正規分布なのです。

条件付確率​

さていきなりここで問題です。図を見ると右の分布は左右の両側を切り取っただけでなく、背が高くなっていることに気がつきましたか？これは何故でしょうか？
そうです、確率の基礎を熟知しているあなたなら（アナタですよ、アナタ！）もう明らかでしょう。そう、確率密度関数の曲線の下の面積は1でないといけないからです（ちゃんと覚えてましたか？）。
もともと面積が1だった正規分布の曲線（図の左）から単純にに左右を切り取ると、面積は1より小さくなってしまいますね。そこで、面積が1になるように縦に引き伸ばしているのです。すごいでしょう？
そしてこれは、条件付確率の事実を表現しているのです。つまり、この状況下では

「我々の手元に来た、という条件の下で重さが xx 以下になる確率は？」

という問いを表現したものが求める分布関数であったのです。まさに条件付確率ではないでしょうか。
「我々の手元に来る確率」は左の図で左右の裾を切り取った部分の面積に他ならないですよね（平均±3\pm 3グラムに収まる確率）。従って、分布関数は、単に左右を切り取ったものを「我々の手元に来る確率」で割ったものになるのです。

その他​

  • 実は上限と下限は必ずしも分布の山を挟む必要はないのです。下図のように山を挟まないで分布の裾の部分だけを残したような切断正規分布ももちろん考えられます（何に使われるかは現在調査中ですが）。

  • 上限と下限を定めた切断正規分布を「両側切断正規分布」、上限または下限のみを定めた切断正規分布を「片側切断正規分布」と呼ぶことがあります。

分布の形状​

基本情報​

  • 4つのパラメータ a,b,m,σa, b,m,\sigma が必要です (どうやって求めるの？).

    a<ba<b σ>0\sigma>0

  • 有限区間 a≤x≤ba\leq x \leq b で定義された連続分布です。

  • 平均対して対称にも非対称にもなり得ます。

確率​

  • 累積分布関数

    F(x)=1Δ[Φ(x−mσ)−Φ(A)]F(x)=\frac{1}{\Delta}\left[\Phi\left(\frac{x-m}{\sigma}\right)-\Phi(A)\right]

    ここで

    Δ=Φ(B)−Φ(A)\Delta=\Phi(B)-\Phi(A) A=a−mσ,B=b−mσA=\frac{a-m}{\sigma},B=\frac{b-m}{\sigma}

    、 Φ(⋅)\Phi(\cdot) は 標準正規分布の累積分布関数です。

  • 確率密度関数

    f(x)=1σΔϕ(x−mσ)f(x)=\frac{1}{\sigma\Delta}\phi\left(\frac{x-m}{\sigma}\right)

  • Excel での累積分布関数 (c.d.f.) と 確率密度関数 (p.d.f.)の求め方

AB
1データ説明
22.5対象となる値
31分布のパラメータ Min の値
44分布のパラメータ Max の値
53分布のパラメータ M の値
60.9分布のパラメータ Sigma の値
7数式説明（計算結果）
8=NTTRUNCNORMDIST(A2,A3,A4,A5,A6,TRUE)上のデータに対する累積分布関数の値
9=NTTRUNCNORMDIST(A2,A3,A4,A5,A6,FALSE)上のデータに対する確率密度関数の値

Sample distribution

分位点​

  • 累積分布関数の逆関数

    F−1(P)=σΦ−1[ΔP+Φ(A)]+mF^{-1}(P)=\sigma\Phi^{-1}\left[\Delta P+\Phi(A)\right]+m

    ここで

    Δ=Φ(B)−Φ(A)\Delta=\Phi(B)-\Phi(A) A=a−mσ,B=a−mσA=\frac{a-m}{\sigma},B=\frac{a-m}{\sigma}

    で、 Φ(⋅)\Phi(\cdot) は標準正規分布の累積分布関数です。

  • Excel での分位点の求め方

AB
1データ説明
20.5上のデータに対する確率密度関数の値
31分布のパラメータ Min の値
44分布のパラメータ Max の値
53分布のパラメータ M の値
60.9分布のパラメータ Sigma の値
7数式説明（計算結果）
8=NTTRUNCNORMINV(A2,A3,A4,A5,A6)上のデータに対する累積分布関数の逆関数の値

分布の特徴​

平均 -- 分布の"中心"はどこ？ (定義)​

  • 分布の平均 は次式で与えられます。

    m+ϕ(A)−ϕ(B)Δσm+\frac{\phi(A)-\phi(B)}{\Delta}\sigma

    ここで

    Δ=Φ(B)−Φ(A)\Delta=\Phi(B)-\Phi(A) A=a−mσ,B=a−mσA=\frac{a-m}{\sigma},B=\frac{a-m}{\sigma}

    で、Φ(⋅)\Phi(\cdot) と ϕ(⋅)\phi(\cdot) はそれぞれ、標準正規分布の確率密度関数と累積分布関数です。

  • Excel での計算法

AB
1データ説明
21分布のパラメータ Min の値
34分布のパラメータ Max の値
43分布のパラメータ M の値
50.9分布のパラメータ Sigma の値
6数式説明（計算結果）
7=NTTRUNCNORMMEAN(A2,A3,A4,A5)上のデータに対す分布の平均

標準偏差 -- 分布はどのくらい広がっているか（定義）​

  • 分布の分散 は次式で与えられます。

    [1+Aϕ(A)−Bϕ(B)Δ−(ϕ(A)−ϕ(B)Δ)2]σ2\left[1+\frac{A\phi(A)-B\phi(B)}{\Delta}-\left(\frac{\phi(A)-\phi(B)}{\Delta}\right)^2\right]\sigma^2

    ここで

    Δ=Φ(B)−Φ(A)\Delta=\Phi(B)-\Phi(A) A=a−mσ,B=a−mσA=\frac{a-m}{\sigma},B=\frac{a-m}{\sigma}

    で、 Φ(⋅)\Phi(\cdot) と ϕ(⋅)\phi(\cdot) はそれぞれ、標準正規分布の確率密度関数と累積分布関数です。

    標準偏差 は 分散の正の平方根です。

  • Excel での計算法

AB
1データ説明
21分布のパラメータ Min の値
34分布のパラメータ Max の値
43分布のパラメータ M の値
50.9分布のパラメータ Sigma の値
6数式説明（計算結果）
7=NTTRUNCNORMSTDEV(A2,A3,A4,A5)上のデータに対する分布の標準偏差

歪度 -- 分布はどちらに偏っているか(定義)​

  • 分布の歪度 は次式で与えられます。

    −1V3/2[2Δ0+(3Δ1−1)Δ0+Δ2]-\frac{1}{V^{3/2}}\left[2\Delta_0+(3\Delta_{1}-1)\Delta_0+\Delta_2\right]

    ここで

    z(x)=ϕ(x)Δz(x)=\frac{\phi(x)}{\Delta} Δk=Bkz(B)−Akz(A)\Delta_k=B^kz(B)-A^kz(A) V=1−Δ1−Δ02V=1-\Delta_1-\Delta_0^2 Δ=Φ(B)−Φ(A)\Delta=\Phi(B)-\Phi(A) A=a−mσ,B=a−mσA=\frac{a-m}{\sigma},B=\frac{a-m}{\sigma}

    で、 Φ(⋅)\Phi(\cdot) と ϕ(⋅)\phi(\cdot) はそれぞれ、標準正規分布の確率密度関数と累積分布関数です。

  • Excel での計算法

AB
1データ説明
21分布のパラメータ Min の値
34分布のパラメータ Max の値
43分布のパラメータ M の値
50.9分布のパラメータ Sigma の値
6数式説明（計算結果）
7=NTTRUNCNORMSKEW(A2,A3,A4,A5)上のデータに対する分布の歪度

尖度 -- 尖っているか丸まっているか (定義)​

  • 分布の尖度 は次式で与えられます。

    1V2[−3Δ04−2Δ02(3Δ1+1)−4Δ2Δ0−3Δ1−Δ3+3]−3\frac{1}{V^2}\left[-3\Delta_0^4-2\Delta_0^2(3\Delta_1+1)-4\Delta_2\Delta_0-3\Delta_1-\Delta_3+3\right]-3

    ここで

    z(x)=ϕ(x)Δz(x)=\frac{\phi(x)}{\Delta} Δk=Bkz(B)−Akz(A)\Delta_k=B^kz(B)-A^kz(A) V=1−Δ1−Δ02V=1-\Delta_1-\Delta_0^2 Δ=Φ(B)−Φ(A)\Delta=\Phi(B)-\Phi(A) A=a−mσ,B=a−mσA=\frac{a-m}{\sigma},B=\frac{a-m}{\sigma}

    で、Φ(⋅)\Phi(\cdot) と ϕ(⋅)\phi(\cdot) はそれぞれ、標準正規分布の確率密度関数と累積分布関数です。

  • Excel での計算法

AB
1データ説明
21分布のパラメータ Min の値
34分布のパラメータ Max の値
43分布のパラメータ M の値
50.9分布のパラメータ Sigma の値
6数式説明（計算結果）
7=NTTRUNCNORMKURT(A2,A3,A4,A5)上のデータに対する分布の尖度

乱数​

  • 乱数 x は一様乱数 U に対して次式で生成されます（逆関数法） :

    x=σΦ−1[ΔU+Φ(A)]+mx=\sigma\Phi^{-1}\left[\Delta U+\Phi(A)\right]+m

    ここで

    Δ=Φ(B)−Φ(A)\Delta=\Phi(B)-\Phi(A) A=a−mσ,B=a−mσA=\frac{a-m}{\sigma},B=\frac{a-m}{\sigma}

    で、Φ(⋅)\Phi(\cdot) と ϕ(⋅)\phi(\cdot) はそれぞれ、標準正規分布の確率密度関数と累積分布関数です。

  • Excel での乱数生成法

AB
1データ説明
21分布のパラメータ Min の値
34分布のパラメータ Max の値
43分布のパラメータ M の値
50.9分布のパラメータ Sigma の値
6数式説明（計算結果）
7=NTRANDTRUNCNORM(100,A2,A3,A4,A5,0) 切断正規乱数を Mersenne Twister アルゴリズムで生成します。

メモ： この使用例の数式は、配列数式として入力する必要があります。使用例を新規ワークシートにコピーした後、A7:A106 のセル範囲 (配列数式が入力されているセルが左上になる) を選択します。F2 キーを押し、Ctrl キーと Shift キーを押しながら Enter キーを押します。この数式が配列数式として入力されていない場合、単一の値 2 のみが計算結果として返されます。

関連 NtRand 関数​

参照​