物件
平均は確率質量中心、つまり1次モーメント です。 中央値は、逆像 F −1 (1/2) です。 レートパラメータλ を持つ指数分布に従う確率変数X の平均値または期待値は 、次式で与えられる。 E [ X ] = 1 λ 。 {\displaystyle \operatorname {E} [X]={\frac {1}{\lambda }}.}
以下の 例を考慮すると、これは理にかなっています。1時間に平均2回の電話を受ける人は、連続する電話の間隔が0.5時間、つまり30分になると予想できます。
X の分散は 次のように与えられる。 バラ [ X ] = 1 λ 2 、 {\displaystyle \operatorname {Var} [X]={\frac {1}{\lambda ^{2}}},} つまり、標準偏差 は平均値と等しい。
X のモーメントは 、n ∈ N {\displaystyle n\in \mathbb {N} } は E [ X n ] = n ! λ n 。 {\displaystyle \operatorname {E} \left[X^{n}\right]={\frac {n!}{\lambda ^{n}}}.}
X の中心モーメント は、n ∈ N {\displaystyle n\in \mathbb {N} } は μ n = ! n λ n = n ! λ n ∑ k = 0 n ( − 1 ) k k ! 。 {\displaystyle \mu _{n}={\frac {!n}{\lambda ^{n}}}={\frac {n!}{\lambda ^{n}}}\sum _{k=0}^{n}{\frac {(-1)^{k}}{k!}}.} ここで 、! nは n の階乗 です。
X の中央値 は次のように与えられる。 m [ X ] = ln ( 2 ) λ < E [ X ] 、 {\displaystyle \operatorname {m} [X]={\frac {\ln(2)}{\lambda }}<\operatorname {E} [X],} ここでln は 自然対数 を表します。したがって、平均値と中央値の 絶対差は | E [ X ] − m [ X ] | = 1 − ln ( 2 ) λ < 1 λ = σ [ X ] 、 {\displaystyle \left|\operatorname {E} \left[X\right]-\operatorname {m} \left[X\right]\right|={\frac {1-\ln(2)}{\lambda }}<{\frac {1}{\lambda }}=\operatorname {\sigma } [X],}
中央値と平均値の不平等 に従って。
指数型確率変数の無記憶性 指数分布に従う確率変数T は次の関係式に従う。 教授 ( T > s + t ∣ T > s ) = 教授 ( T > t ) 、 ∀ s 、 t ≥ 0. {\displaystyle \Pr \left(T>s+t\mid T>s\right)=\Pr(T>t),\qquad \forall s,t\geq 0.}
これは、相補累積分布関数 を考慮することで理解できる。 教授 ( T > s + t ∣ T > s ) = 教授 ( T > s + t ∩ T > s ) 教授 ( T > s ) = 教授 ( T > s + t ) 教授 ( T > s ) = e − λ ( s + t ) e − λ s = e − λ t = 教授 ( T > t ) 。 {\displaystyle {\begin{aligned}\Pr \left(T>s+t\mid T>s\right)&={\frac {\Pr \left(T>s+t\cap T>s\right)}{\Pr \left(T>s\right)}}\\[4pt]&={\frac {\Pr \left(T>s+t\right)}{\Pr \left(T>s\right)}}\\[4pt]&={\frac {e^{-\lambda (s+t)}}{e^{-\lambda s}}}\\[4pt]&=e^{-\lambda t}\\[4pt]&=\Pr(T>t).\end{aligned}}}
Tを ある初期時刻からの事象発生までの待ち時間として解釈する場合、この関係は、 Tが 初期期間s 内に事象が観測されないことを条件としている場合、残りの待ち時間の分布は元の無条件分布と同じであることを意味します。例えば、30秒経過しても事象が発生しない場合、発生までに少なくとも10秒以上かかるという条件付き確率 は、初期時刻から10秒以上経過後に事象が観測されるという無条件確率と等しくなります。
指数分布と幾何分布は、 記憶を持たない唯一の確率分布 である。
したがって、指数分布は必然的に、一定の故障率 を持つ唯一の連続確率分布でもある。
分位数 異常値判定のためのテューキー基準。 Exp( λ )の 分位関数 (逆累積分布関数)は F − 1 ( p ; λ ) = − ln ( 1 − p ) λ 、 0 ≤ p < 1 {\displaystyle F^{-1}(p;\lambda )={\frac {-\ln(1-p)}{\lambda }},\qquad 0\leq p<1}
したがって、四分位数 は次のとおりです。
第1四分位数:ln(4/3)/ λ 中央値 :ln(2)/ λ 第3四分位数:ln(4)/ λ その結果、四分位範囲 は ln(3)/ λ となります。
条件付きバリュー・アット・リスク(予想不足額)条件付きリスク値(CVaR)は、Exp( λ )の 期待ショートフォール またはスーパークオンタイルとしても知られており、次のように導出されます。[ 4 ]
q ¯ α ( X ) = 1 1 − α ∫ α 1 q p ( X ) d p = 1 ( 1 − α ) ∫ α 1 − ln ( 1 − p ) λ d p = − 1 λ ( 1 − α ) ∫ 1 − α 0 − ln ( y ) d y = − 1 λ ( 1 − α ) ∫ 0 1 − α ln ( y ) d y = − 1 λ ( 1 − α ) [ ( 1 − α ) ln ( 1 − α ) − ( 1 − α ) ] = − ln ( 1 − α ) + 1 λ {\displaystyle {\begin{aligned}{\bar {q}}_{\alpha }(X)&={\frac {1}{1-\alpha }}\int _{\alpha }^{1}q_{p}(X)dp\\&={\frac {1}{(1-\alpha )}}\int _{\alpha }^{1}{\frac {-\ln(1-p)}{\lambda }}dp\\&={\frac {-1}{\lambda (1-\alpha )}}\int _{1-\alpha }^{0}-\ln(y)dy\\&={\frac {-1}{\lambda (1-\alpha )}}\int _{0}^{1-\alpha }\ln(y)dy\\&={\frac {-1}{\lambda (1-\alpha )}}[(1-\alpha )\ln(1-\alpha )-(1-\alpha )]\\&={\frac {-\ln(1-\alpha )+1}{\lambda }}\\\end{aligned}}}
緩衝型超過確率(bPOE)超過確率の緩衝値は、CVaRが閾値と等しくなる確率レベルから1を引いた値である。x {\displaystyle x} 以下のように導出される:[ 4 ]
p ¯ x ( X ) = { 1 − α | q ¯ α ( X ) = x } = { 1 − α | − ln ( 1 − α ) + 1 λ = x } = { 1 − α | ln ( 1 − α ) = 1 − λ x } = { 1 − α | e ln ( 1 − α ) = e 1 − λ x } = { 1 − α | 1 − α = e 1 − λ x } = e 1 − λ x {\displaystyle {\begin{aligned}{\bar {p}}_{x}(X)&=\{1-\alpha |{\bar {q}}_{\alpha }(X)=x\}\\&=\{1-\alpha |{\frac {-\ln(1-\alpha )+1}{\lambda }}=x\}\\&=\{1-\alpha |\ln(1-\alpha )=1-\lambda x\}\\&=\{1-\alpha |e^{\ln(1-\alpha )}=e^{1-\lambda x}\}=\{1-\alpha |1-\alpha =e^{1-\lambda x}\}=e^{1-\lambda x}\end{aligned}}}
カルバック・ライブラー情報量ナッツ の方向性のあるカルバック・ライブラーダイバージェンス e λ {\displaystyle e^{\lambda }} (分布を「近似」する)e λ 0 {\displaystyle e^{\lambda _{0}}} (「真の」分布)は次のように与えられる。 Δ ( λ 0 ∥ λ ) = E λ 0 ( ログ p λ 0 ( x ) p λ ( x ) ) = E λ 0 ( ログ λ 0 e λ 0 x λ e λ x ) = ログ ( λ 0 ) − ログ ( λ ) − ( λ 0 − λ ) E λ 0 ( x ) = ログ ( λ 0 ) − ログ ( λ ) + λ λ 0 − 1. {\displaystyle {\begin{aligned}\Delta (\lambda _{0}\parallel \lambda )&=\mathbb {E} _{\lambda _{0}}\left(\log {\frac {p_{\lambda _{0}}(x)}{p_{\lambda }(x)}}\right)\\&=\mathbb {E} _{\lambda _{0}}\left(\log {\frac {\lambda _{0}e^{\lambda _{0}x}}{\lambda e^{\lambda x}}}\right)\\&=\log(\lambda _{0})-\log(\lambda )-(\lambda _{0}-\lambda )E_{\lambda _{0}}(x)\\&=\log(\lambda _{0})-\log(\lambda )+{\frac {\lambda }{\lambda _{0}}}-1.\end{aligned}}}
指数分布に従う確率変数の最小値 X 1 、 ...、X n をレートパラメータλ 1 、 ...、λ n を持つ独立な 指数分布確率変数とする。すると[ 6 ] ミニ { X 1 、 … 、 X n } {\displaystyle \min \left\{X_{1},\dotsc ,X_{n}\right\}} は指数分布に従い、パラメータは λ = λ 1 + ⋯ + λ n 。 {\displaystyle \lambda =\lambda _{1}+\dotsb +\lambda _{n}.}
これは、相補累積分布関数 を考慮することで理解できる。 教授 ( ミニ { X 1 、 … 、 X n } > x ) = 教授 ( X 1 > x 、 … 、 X n > x ) = ∏ 私 = 1 n 教授 ( X 私 > x ) = ∏ 私 = 1 n exp ( − x λ 私 ) = exp ( − x ∑ 私 = 1 n λ 私 ) 。 {\displaystyle {\begin{aligned}&\Pr \left(\min\{X_{1},\dotsc ,X_{n}\}>x\right)\\={}&\Pr \left(X_{1}>x,\dotsc ,X_{n}>x\right)\\={}&\prod _{i=1}^{n}\Pr \left(X_{i}>x\right)\\={}&\prod _{i=1}^{n}\exp \left(-x\lambda _{i}\right)=\exp \left(-x\sum _{i=1}^{n}\lambda _{i}\right).\end{aligned}}}
最小値を達成する変数の指標は、カテゴリ分布に従って分布する。 教授 ( X k = ミニ { X 1 、 … 、 X n } ) = λ k λ 1 + ⋯ + λ n 。 {\displaystyle \Pr \left(X_{k}=\min\{X_{1},\dotsc ,X_{n}\}\right)={\frac {\lambda _{k}}{\lambda _{1}+\dotsb +\lambda _{n}}}.}
証明は、私 = アルグミン 私 ∈ { 1 、 ⋯ 、 n } { X 1 、 … 、 X n } {\displaystyle I=\operatorname {argmin} _{i\in \{1,\dotsb ,n\}}\{X_{1},\dotsc ,X_{n}\}} 。 それから、 教授 ( 私 = k ) = ∫ 0 ∞ 教授 ( X k = x ) 教授 ( ∀ 私 ≠ k X 私 > x ) d x = ∫ 0 ∞ λ k e − λ k x ( ∏ 私 = 1 、 私 ≠ k n e − λ 私 x ) d x = λ k ∫ 0 ∞ e − ( λ 1 + ⋯ + λ n ) x d x = λ k λ 1 + ⋯ + λ n 。 {\displaystyle {\begin{aligned}\Pr(I=k)&=\int _{0}^{\infty }\Pr(X_{k}=x)\Pr(\forall _{i\neq k}X_{i}>x)\,dx\\&=\int _{0}^{\infty }\lambda _{k}e^{-\lambda _{k}x}\left(\prod _{i=1,i\neq k}^{n}e^{-\lambda _{i}x}\right)dx\\&=\lambda _{k}\int _{0}^{\infty }e^{-\left(\lambda _{1}+\dotsb +\lambda _{n}\right)x}dx\\&={\frac {\lambda _{k}}{\lambda _{1}+\dotsb +\lambda _{n}}}.\end{aligned}}}
対照的に、最大 { X 1 、 … 、 X n } {\displaystyle \max\{X_{1},\dotsc ,X_{n}\}} X 1 、 ...、X n が すべてパラメータ 0 を持たない場合、指数分布には従わない。[ 7 ]
統計的推論 以下では、確率変数X がレート パラメータ λ の指数分布に従うと仮定し、x 1 、 … 、 x n {\displaystyle x_{1},\dotsc ,x_{n}} はXからの n 個の独立したサンプルであり、サンプル平均はx ¯ {\displaystyle {\bar {x}}} 。
パラメータ推定
最尤推定法 λの最尤推定量は次のように構築される。
変数から抽出された 独立同分布の 標本x = ( x 1 , ..., x n )が与えられた場合の λ の尤度関数は次のようになります。 L ( λ ) = ∏ 私 = 1 n λ exp ( − λ x 私 ) = λ n exp ( − λ ∑ 私 = 1 n x 私 ) = λ n exp ( − λ n x ¯ ) 、 {\displaystyle L(\lambda )=\prod _{i=1}^{n}\lambda \exp(-\lambda x_{i})=\lambda ^{n}\exp \left(-\lambda \sum _{i=1}^{n}x_{i}\right)=\lambda ^{n}\exp \left(-\lambda n{\overline {x}}\right),}
どこ: x ¯ = 1 n ∑ 私 = 1 n x 私 {\displaystyle {\overline {x}}={\frac {1}{n}}\sum _{i=1}^{n}x_{i}} は標本平均です。
尤度関数の対数の導関数は次のとおりです。 d d λ ln L ( λ ) = d d λ ( n ln λ − λ n x ¯ ) = n λ − n x ¯ { > 0 、 0 < λ < 1 x ¯ 、 = 0 、 λ = 1 x ¯ 、 < 0 、 λ > 1 x ¯ 。 {\displaystyle {\frac {d}{d\lambda }}\ln L(\lambda )={\frac {d}{d\lambda }}\left(n\ln \lambda -\lambda n{\overline {x}}\right)={\frac {n}{\lambda }}-n{\overline {x}}\ {\begin{cases}>0,&0<\lambda <{\frac {1}{\overline {x}}},\\[8pt]=0,&\lambda ={\frac {1}{\overline {x}}},\\[8pt]<0,&\lambda >{\frac {1}{\overline {x}}}.\end{cases}}}
したがって、レートパラメータの 最尤推定値は次のようになる。 λ ^ mle = 1 x ¯ = n ∑ 私 = 1 n x 私 {\displaystyle {\widehat {\lambda }}_{\text{mle}}={\frac {1}{\overline {x}}}={\frac {n}{\sum _{i=1}^{n}x_{i}}}}
これは不偏推定量 ではありません λ 、 {\displaystyle \lambda ,} それでもx ¯ {\displaystyle {\overline {x}}} は、不偏な[ 13 ] MLE [ 14 ] である。 1 / λ {\displaystyle 1/\lambda } そして分布の平均。バイアスを減らす方法があります。λ ^ mle {\displaystyle {\widehat {\lambda }}_{\text{mle}}} バイアス補正された最尤推定量 が得られる。
サンプルサイズが2より大きいと仮定すると、最尤推定値に対する補正係数を用いて、 平均二乗誤差 の近似最小値(バイアス・バリアンスのトレードオフ も参照)を見つけることができます。λ ^ = ( n − 2 n ) ( 1 x ¯ ) = n − 2 ∑ 私 = 1 n x 私 {\displaystyle {\widehat {\lambda }}=\left({\frac {n-2}{n}}\right)\left({\frac {1}{\bar {x}}}\right)={\frac {n-2}{\sum _{i=1}^{n}x_{i}}}} これは逆ガンマ分布 の平均と分散から導き出され、インビガンマ ( n 、 λ ) {\textstyle {\mbox{Inv-Gamma}}(n,\lambda )} [ 15 ]
UMVUE 一様最小分散不偏推定 量λ {\displaystyle \lambda } は
λ ^ u m v u e = n − 1 ∑ 私 = 1 n x 私 、 {\displaystyle {\widehat {\lambda }}_{\mathrm {umvue} }={\frac {n-1}{\sum _{i=1}^{n}x_{i}}},}
これは、最尤推定量とは対照的に、定義上不偏である。
フィッシャー情報 量私 ( λ ) {\displaystyle {\mathcal {I}}(\lambda )} レートパラメータの推定値の場合λ {\displaystyle \lambda } 次のように与えられます。 私 ( λ ) = E [ ( ∂ ∂ λ ログ f ( x ; λ ) ) 2 | λ ] = ∫ ( ∂ ∂ λ ログ f ( x ; λ ) ) 2 f ( x ; λ ) d x {\displaystyle {\mathcal {I}}(\lambda )=\operatorname {E} \left[\left.\left({\frac {\partial }{\partial \lambda }}\log f(x;\lambda )\right)^{2}\right|\lambda \right]=\int \left({\frac {\partial }{\partial \lambda }}\log f(x;\lambda )\right)^{2}f(x;\lambda )\,dx}
分布を代入して解くと次のようになる。 私 ( λ ) = ∫ 0 ∞ ( ∂ ∂ λ ログ λ e − λ x ) 2 λ e − λ x d x = ∫ 0 ∞ ( 1 λ − x ) 2 λ e − λ x d x = λ − 2 。 {\displaystyle {\mathcal {I}}(\lambda )=\int _{0}^{\infty }\left({\frac {\partial }{\partial \lambda }}\log \lambda e^{-\lambda x}\right)^{2}\lambda e^{-\lambda x}\,dx=\int _{0}^{\infty }\left({\frac {1}{\lambda }}-x\right)^{2}\lambda e^{-\lambda x}\,dx=\lambda ^{-2}.}
これは、指数分布の各独立サンプルが未知のレートパラメータについて持つ情報量を決定する。λ {\displaystyle \lambda } 。
信頼区間 指数分布のレートパラメータの正確な 100(1 − α)% 信頼区間は次のように与えられます。[ 16 ] 2 n λ ^ mle χ α 2 、 2 n 2 < 1 λ < 2 n λ ^ mle χ 1 − α 2 、 2 n 2 、 {\displaystyle {\frac {2n}{{\widehat {\lambda }}_{\textrm {mle}}\chi _{{\frac {\alpha }{2}},2n}^{2}}}<{\frac {1}{\lambda }}<{\frac {2n}{{\widehat {\lambda }}_{\textrm {mle}}\chi _{1-{\frac {\alpha }{2}},2n}^{2}}}\,,} これはまた、 2 n x ¯ χ α 2 、 2 n 2 < 1 λ < 2 n x ¯ χ 1 − α 2 、 2 n 2 、 {\displaystyle {\frac {2n{\overline {x}}}{\chi _{{\frac {\alpha }{2}},2n}^{2}}}<{\frac {1}{\lambda }}<{\frac {2n{\overline {x}}}{\chi _{1-{\frac {\alpha }{2}},2n}^{2}}}\,,} ここで、χ 2 p , v は自由度 v のカイ二乗分布 の 100( p ) パーセンタイル、n は観測数、x バーは標本平均です。χ 2 p , v 分布の正規近似を用いることで、正確な区間端点の簡単な近似値を導出できます。この近似値を用いると、 95 % 信頼区間に対して以下の値が得られます。 λ より低い = λ ^ ( 1 − 1.96 n ) λ アッパー = λ ^ ( 1 + 1.96 n ) {\displaystyle {\begin{aligned}\lambda _{\text{lower}}&={\widehat {\lambda }}\left(1-{\frac {1.96}{\sqrt {n}}}\right)\\\lambda _{\text{upper}}&={\widehat {\lambda }}\left(1+{\frac {1.96}{\sqrt {n}}}\right)\end{aligned}}}
この近似は、少なくとも15〜20個の元素を含むサンプルに対しては許容できる可能性がある。[ 17 ]
共役事前分布を用いたベイズ推論 指数分布の共役事前 分布はガンマ分布 です(指数分布はガンマ分布の特殊な場合です)。ガンマ確率密度関数の以下のパラメータ化が有用です。
ガンマ ( λ ; α 、 β ) = β α Γ ( α ) λ α − 1 exp ( − λ β ) 。 {\displaystyle \operatorname {Gamma} (\lambda ;\alpha ,\beta )={\frac {\beta ^{\alpha }}{\Gamma (\alpha )}}\lambda ^{\alpha -1}\exp(-\lambda \beta ).}
事後分布 pは 、上記で定義した尤度関数とガンマ事前分布を用いて次のように表すことができる。
p ( λ ) ∝ L ( λ ) Γ ( λ ; α 、 β ) = λ n exp ( − λ n x ¯ ) β α Γ ( α ) λ α − 1 exp ( − λ β ) ∝ λ ( α + n ) − 1 exp ( − λ ( β + n x ¯ ) ) 。 {\displaystyle {\begin{aligned}p(\lambda )&\propto L(\lambda )\Gamma (\lambda ;\alpha ,\beta )\\&=\lambda ^{n}\exp \left(-\lambda n{\overline {x}}\right){\frac {\beta ^{\alpha }}{\Gamma (\alpha )}}\lambda ^{\alpha -1}\exp(-\lambda \beta )\\&\propto \lambda ^{(\alpha +n)-1}\exp(-\lambda \left(\beta +n{\overline {x}}\right)).\end{aligned}}}
これで事後確率密度pは 、不足している正規化定数を除いて特定されました。ガンマ分布の形をしているため、これは簡単に補完でき、次の式が得られます。
p ( λ ) = ガンマ ( λ ; α + n 、 β + n x ¯ ) 。 {\displaystyle p(\lambda )=\operatorname {Gamma} (\lambda ;\alpha +n,\beta +n{\overline {x}}).}
ここで、ハイパーパラメータ α は事前観測の数、β は事前観測の合計と解釈できます。事後平均は次のようになります。 α + n β + n x ¯ 。 {\displaystyle {\frac {\alpha +n}{\beta +n{\overline {x}}}}.}
較正事前分布を用いたベイズ推論 指数分布は、グループ 構造を持つ統計分布の1つです。グループ構造の結果として、指数分布には関連するハール測度 があり、それは1 / λ 。 {\displaystyle 1/\lambda .} ベイズ予測において、ハール測度を 事前分布 (ハール事前分布として知られる)として使用すると、基礎となる真のパラメータ値に関わらず、完全に較正された確率が得られます。 [ 18 ] [ 19 ] [ 20 ] 完全に較正された確率は、予測された確率がサンプル外イベントの頻度と完全に一致するという性質を持ちます。指数分布の場合、ハール事前分布を使用して生成されたベイズ予測の正確な式は、次のように与えられます。
p H 1 1 r − p r 私 o r ( x n + 1 ∣ x 1 、 … 、 x n ) = n n + 1 ( x ¯ ) n ( n x ¯ + x n + 1 ) n + 1 。 {\displaystyle p_{\rm {Haar-prior}}(x_{n+1}\mid x_{1},\ldots ,x_{n})={\frac {n^{n+1}\left({\overline {x}}\right)^{n}}{\left(n{\overline {x}}+x_{n+1}\right)^{n+1}}}.}
これは、事前予測の較正の一例であり、較正精度を向上させるため(そしてこの場合は較正を完璧にするため)に事前分布が選択されます。指数分布に対するHaar事前分布を用いた事前予測の較正は、R ソフトウェアパッケージfitdistcpに実装されています。
同様の予測は、以下の予測の セクションで説明するように、他の多くの視点からも導き出すことができる。
発生と応用
出来事の発生 指数分布は、均質なポアソン過程 における到着間隔時間の長さを記述する際に自然に現れる。
指数分布は、離散的な プロセスが状態変化を起こすのに必要なベルヌーイ試行 の回数を表す幾何分布 の連続版と見なすことができる。一方、指数分布は、連続的なプロセスが状態変化を起こすのにかかる時間を表す。
現実世界では、一定のレート(または単位時間あたりの確率)という仮定はめったに満たされません。たとえば、着信電話のレートは時間帯によって異なります。しかし、平日の午後2時から4時までのように、レートがほぼ一定の時間帯に注目すれば、指数分布は次の電話がかかってくるまでの時間の近似モデルとして使用できます。指数分布に近い変数が得られる以下の例にも、同様の注意点があります。
放射性粒子が崩壊する までの時間、またはガイガーカウンター のクリック音の間隔。電話を1回受けてから次の電話を受けるまでの時間 簡略化された信用リスクモデルにおける、債務不履行(企業債権者への支払い)までの期間 指数変数は、 DNA 鎖上の突然変異 間の距離や、特定の道路における轢死事故の間隔など、単位長さあたりの一定の確率で特定の事象が発生する状況をモデル化するためにも使用できます。
待ち行列理論 では、システム内のエージェントのサービス時間(例えば、銀行の窓口係などが顧客にサービスを提供するのにかかる時間など)は、指数分布に従う変数としてモデル化されることが多い。(例えば、顧客の到着は、到着が独立で同一の分布に従う場合、ポアソン分布 によってモデル化される。)複数の独立したタスクのシーケンスと考えることができるプロセスの長さは、アーラン分布 (複数の独立した指数分布に従う変数の合計の分布)に従う。
信頼性理論 および信頼性工学 においても、指数分布は広く利用されています。この分布は無記憶性を持つため、信頼性理論で使用されるバスタブ曲線の一定 ハザード率 部分をモデル化するのに適しています。また、信頼性モデルに故障率を 簡単に追加できるため、非常に便利です。しかしながら、指数分布は生物や技術機器の全体的な寿命をモデル化するのに適していません。なぜなら、これらの「故障率」は一定ではなく、非常に若いシステムと非常に古いシステムでは故障が多く発生するからです。
年間最大1日降雨量に累積指数分布を当てはめた。 物理学 において、一定の温度 と圧力 で均一な重力場に ある気体を 観測すると、様々な分子の高さも近似的に指数分布に従うことが知られています。これは気圧公式 として知られています。これは、後述するエントロピーの性質の結果です。
水文学 では、指数分布は、日降雨量や河川流量の月間および年間最大値などの変数の極値を分析するために使用されます。[ 21 ]
青色の図は、年間最大1日降水量を順位付けしたデータに指数分布を当てはめた例を示しており、二項分布 に基づく90%信頼区間も表示されています。降水量データは、 累積頻度分析 の一環として、位置をプロットする ことで表されています。 手術室管理においては、典型的な作業内容を持たない 手術のカテゴリー(例えば、救急室におけるあらゆる種類の手術を含む)における手術時間の分布。
予測 未知の指数分布からn個 のデータ点のサンプルを観測した後、同じソースからの将来のデータについて予測を行うのが一般的なタスクです。将来のサンプルに対する一般的な予測分布は、レートパラメータλ の適切な推定値を指数密度関数に代入することによって形成される、いわゆるプラグイン分布です。一般的な推定値は最尤法の原理によって提供されるものであり、これを使用すると、観測されたサンプルx = ( x 1 , ..., x n )を条件とした 将来のサンプルx n +1 に対する予測密度が得られます。p M L ( x n + 1 ∣ x 1 、 … 、 x n ) = ( 1 x ¯ ) exp ( − x n + 1 x ¯ ) 。 {\displaystyle p_{\rm {ML}}(x_{n+1}\mid x_{1},\ldots ,x_{n})=\left({\frac {1}{\overline {x}}}\right)\exp \left(-{\frac {x_{n+1}}{\overline {x}}}\right).}
ベイズアプローチは、推定パラメータの不確実性を考慮した予測分布を提供するが、これは事前分布の選択に大きく依存する可能性がある。
主観的ベイズアプローチで生じる事前分布の選択の問題から解放された予測分布は
p C N M L ( x n + 1 ∣ x 1 、 … 、 x n ) = n n + 1 ( x ¯ ) n ( n x ¯ + x n + 1 ) n + 1 、 {\displaystyle p_{\rm {CNML}}(x_{n+1}\mid x_{1},\ldots ,x_{n})={\frac {n^{n+1}\left({\overline {x}}\right)^{n}}{\left(n{\overline {x}}+x_{n+1}\right)^{n+1}}},}
これは次のように考えることができる
中心量の分布から得られる頻度論的信頼分布 x n + 1 / x ¯ {\displaystyle {x_{n+1}}/{\overline {x}}} ; [ 22 ] プロファイル予測尤度は、 x n +1 とλ の同時尤度からパラメータλ を 最大化によって除去することによって得られます。[ 23 ] 客観的なベイズ予測事後分布は、非情報的なジェフリーズ事前分布 1/ λ を使用して得られ、 この場合、右ハール事前分布に等しくなります。右ハール事前分布を使用して生成された予測は、完全に較正された確率を与えることが保証されています。[ 24 ] [ 25 ] 情報理論的考察に基づく条件付き正規化最大尤度(CNML)予測分布。[ 26 ] 予測分布の精度は、レートパラメータλ₀を持つ真の指数分布と、サンプルxに基づく予測分布との間の距離または乖離を用いて測定できます。カルバック・ライブラー情報量は、2つの分布間の差を測る、パラメータ化を必要としない一般的な尺度です。レートパラメータλ₀ を持つ 指数分布 と予測分布pとの間のカルバック・ライブラー情報量をΔ( λ₀ || p ) とすると 、 次の ことが示されます。
E λ 0 [ Δ ( λ 0 ∥ p M L ) ] = ψ ( n ) + 1 n − 1 − ログ ( n ) E λ 0 [ Δ ( λ 0 ∥ p C N M L ) ] = ψ ( n ) + 1 n − ログ ( n ) {\displaystyle {\begin{aligned}\operatorname {E} _{\lambda _{0}}\left[\Delta (\lambda _{0}\parallel p_{\rm {ML}})\right]&=\psi (n)+{\frac {1}{n-1}}-\log(n)\\\operatorname {E} _{\lambda _{0}}\left[\Delta (\lambda _{0}\parallel p_{\rm {CNML}})\right]&=\psi (n)+{\frac {1}{n}}-\log(n)\end{aligned}}}
ここで、期待値はレートパラメータλ 0 ∈ (0, ∞) を持つ指数分布に関して取られ、ψ( · )はディガンマ関数です。CNML予測分布は、すべてのサンプルサイズ n > 0 について、平均カルバック・ライブラー情報量の観点から、最尤プラグイン分布よりも明らかに優れています。
乱数生成 指数分布に従う乱数 を生成する概念的に非常に単純な方法は、逆変換サンプリング に基づいています。単位区間(0, 1)上の 一様分布 から抽出された乱数U が与えられた場合、乱数は
T = F − 1 ( U ) {\displaystyle T=F^{-1}(U)}
は指数分布を持ち、ここでF −1 は分位関数 であり、次のように定義される。
F − 1 ( p ) = − ln ( 1 − p ) λ 。 {\displaystyle F^{-1}(p)={\frac {-\ln(1-p)}{\lambda }}.}
さらに、U が (0, 1) 上で一様分布であれば、1 − U も一様分布になります。これは、次のように指数分布に従う変量を生成できることを意味します。
T = − ln ( U ) λ 。 {\displaystyle T={\frac {-\ln(U)}{\lambda }}.}
指数分布に従う変数を生成する他の方法については、Knuth [ 27 ] および Devroye [ 28 ]によって議論されている。
ソートルーチンを使用せずに、順序付けされた指数変量のセットを生成する高速な方法も利用可能です。[ 28 ]