Üstel dağılım

Olasılık kuramı ve istatistik bilim dallarında üstel dağılımı bir sürekli olasılık dağılımları grubudur. Sabit ortalama değişme haddinde ortaya çıkan bağımsız olaylar arasındaki zaman aralığını modelleştirirken bir üstel dağılım doğal olarak ortaya çıkar.

Üstel

Olasılık yoğunluk fonksiyonu
Yığmalı dağılım fonksiyonu
Parametreler	${\displaystyle \lambda >0\,}$ oran veya ters ölçek (reel)
Destek	${\displaystyle [0,\infty )\!}$
Olasılık yoğunluk fonksiyonu (OYF)	${\displaystyle \lambda e^{-\lambda x}}$
Birikimli dağılım fonksiyonu (YDF)	${\displaystyle 1-e^{-\lambda x}}$
Ortalama	${\displaystyle {\frac {1}{\lambda }}\,}$
Medyan	${\displaystyle {\frac {\ln(2)}{\lambda }}\,}$
Mod	${\displaystyle 0\,}$
Varyans	${\displaystyle \lambda ^{-2}\,}$
Çarpıklık	${\displaystyle 2\,}$
Fazladan basıklık	${\displaystyle 6\,}$
Entropi	${\displaystyle 1-\ln(\lambda )\,}$
Moment üreten fonksiyon (mf)	${\displaystyle \left(1-{\frac {t}{\lambda }}\right)^{-1}\,}$
Karakteristik fonksiyon	${\displaystyle \left(1-{\frac {it}{\lambda }}\right)^{-1}\,}$

Tipik karakteristikler

Olasılık yoğunluk fonksiyonu

Bir üstel dağılım için olasılık yoğunluk fonksiyonu şu şekli alır:

${\displaystyle f(x;\lambda )=\left\{{\begin{matrix}\lambda e^{-\lambda x}&,\;x\geq 0,\\0&,\;x<0.\end{matrix}}\right.}$

Burada λ > 0 dağılım için tek parametredir ve çok zaman oran parametresi olarak anılır. Dağılım için destek [0,∞) aralığında verilir. Eğer X rassal değişkeni bu üstel dağılım gösteriyorsa bu şöyle yazılır:

X ~ Üstel(λ).

Ancak bir diğer şekilde değişik parametreleme ile ise üstel dağılım için olasılık yoğunluk fonksiyonu şöyle ifade edilir:

${\displaystyle f(x;\beta )=\left\{{\begin{matrix}{\frac {1}{\beta }}e^{-x/\beta }&,\;x\geq 0,\\0&,\;x<0.\end{matrix}}\right.}$

Burada β > 0 bir ölçek parametresidir ve yukarıda tanımlanan oran parametresi olan λ'nın bir üstü değeri çarpım tersi, yani β=1/λ; dır. Bu çeşit tanımlamada β kalım parametresi çünkü eğer bir rassal değişken X bir biyolojik veya mekanik sistem M için ömür geçirme zaman uzunluğu ise ve X ~ Üstel(β) ise

${\displaystyle \mathbb {E} [X]=\beta }$

yani M için beklenen hayatta kalım süresi zaman birimleri ile β olur.

Bu ikinci şekilde tanımlama bazen birinci tanımlamadan daha kullanışlı olur ve bazı istatistikçiler bu ikinci tanımı üstel dağılım için standard tanım kabul etmektedirler.

Bu gerçek dikkat çekilmesi gereken bir konu olarak burada işaret edilmektedir. Çünkü iki değişik tanım bazen bir kavram karmaşaklığına neden olmaktadır. Genel olarak üstel dağılımı kullanan istatistikçi birinci tanım kullanırsa

X ~ Üstel(λ)

ve ikinci tanımı kullanırsa

X ~ Üstel(β)

yazılır ve β=1/λ olur.

Yığmalı dağılım fonksiyonu

Genel olarak kullanılan bir yönteme göre yığmalı dağılım fonksiyonu şu ifade ile verilir:

${\displaystyle F(x;\lambda )=\left\{{\begin{matrix}1-e^{-\lambda x}&,\;x\geq 0,\\0&,\;x<0.\end{matrix}}\right.}$

Ortaya çıkma ve uygulanma

Bir homojen Poisson süreçde varışlar arasındaki zaman dönemlerini tanımlarken üstel dağılım doğal olarak ortaya çıkar.

Üstel dağılım geometrik dağılımin sürekli dağılımlara uzantısı olarak görülebilir. Geometrik dağılım durumu değiştirmek gereken Bernoulli süreçlerinin sayısını tanımlar ve bu yüzden bir ayrık süreçtir. Buna karşılık, durumu değiştirmek için sürekli bir süreç için geçen zamanı tanımlar.

Pratik gerçek hayatta bir değişme oranının (veya her zaman birimi içinde olasılığın) gerçekleşmesi çok nadirdir. Örneğin, bir mobil telefona gelen çağrılar birim saatin gün içindeki yerine değişir. Fakat araştırmamızı günün öyle bir zaman aralığına odaklayabiliriz ki (diyelim öğleden sonra 2 ile 4), bu zaman aralığından gelen telefon çağrı ortalamları kabaca sabit olabilir. Üstel dağılım o halde iyi bir yaklaşık model olarak kullanılabilir ve en son çağrıdan sonra ne zaman aralığından sonra bir yeni çağrının geleceği hakkında üstel dağılım kullanarak tahmin yapabiliriz.

Benzer şekilde uzun ve karmaşık varsayımlar ve açıklamalar pratikte yaklaşık olarak üstel dağılım gösteren değişkenlere da uygulan şu olaylar için de uygulanabilir:

• bir radyoaktif parçacığın bozunmasına kadar geçen zaman veya bir geiger sayacının birbirini takibedecek düdük seslerinin arasında geçen zamanın tahmini;
• gelecek telefon çağrısını en son yaptığınız çağrıdan ne kadar zaman sonra yapacağınız;
• indirgenmis şekilde olan kredi rizikosu modelinde bir firmanın borçluları ile ilgili olarak en son borcunu ödeyemiyecegini bildiren borçludan ne zaman sonra bir başka daha borç ödeyemiyecek borçlu çıkacağını tahmin etmek.

Özellikler

Ortalama ve varyans

Bir λ oran parametresi ile üstel dağılım gösteren bir X rassal değişkeni için ortalama veya beklenen değer şöyle verilir:

${\displaystyle \mathrm {E} [X]={\frac {1}{\lambda }}.\!}$

Bu verilen pratik örneklerden sağduyu ile çıkarılabilir. Örneğin eğer telefon çağrı ortalama oranı saatte 3 ise (λ), her telefon çağrısı için ortalama 1/3 saat veya 20 dakika (β) beklemek gerekmektedir

X için varyans şöyle verilir

${\displaystyle \mathrm {Var} [X]={\frac {1}{\lambda ^{2}}}.\!}$

Belleksizlik

Üstel dağılımın bir önemli niteliği de belleksiz olmasıdır. Bu demektir ki eğer bir rassal değişken T üstel dağılım gösteriyorsa, onun koşullu olasılığı

${\displaystyle P(T>s+t\;|\;T>s)=P(T>t)\;\;{\hbox{butun}}\ s,t\geq 0.}$

ifadesine uygunluk gösterir. Buna göre, bir hizmet noktasındaki hizmet ve bekleme kuyruğu problemi örneği için bir koşullu olasılık olan ilk varışın 30 saniye geçtikten sonra ortaya çıkmadığını bilerek ilk varıştan sonra 10 saniyeden daha fazla beklemek gereğinin olasılığının, birinci varıştan sonra 10 saniyeden daha fazla bekleme gereğinin koşulsuz başlangıç olasılığı arasında bir fark yoktur. Bu çok kere olasılık hesaplarını ilk gören kişiler tarafından yanlış anlaşılmaktadır:

P(T > 40 | T > 30) = P(T > 10)

gerçeği

T>40 ve T>30

olayları birbirinden bağımsızdır anlamına gelmez. İlk varışa kadar T bekleme zamanının olasılık dağılımının belleksizlik karakteri olduğunu bildirmek

${\displaystyle \mathrm {(Dogru)} \ P(T>40\mid T>30)=P(T>10).}$

olur demektir; yoksa

${\displaystyle \mathrm {(Yanlis)} \ P(T>40\mid T>30)=P(T>40).}$

demek değildir çünkü bu ikinci ifade bağımsızlık kavramını açıklar ve burada olaylar bağımsız değildir.

Bütün mevcut dağılımlar arasında sadece üstel dağılımlar ve geometrik dağılımlar belleksizlik özelliği taşırlar.

Üstel dağılımının ayrıca sabit bir tehlike fonksiyonu bulunmaktadır.

Dörtebirlikler

Bir λ parametreli üstel dağılım için (ters yığmalı dağılım fonksiyonu) şudur:

${\displaystyle F^{-1}(p;\lambda )={\frac {-\ln(1-p)}{\lambda }},\!}$

burada 0 ≤ p < 1.

Onun için şu ifadeler dörttebirlikler verir:

birinci dörttebirlik : ${\displaystyle \ln(4/3)/\lambda \,}$

medyan : ${\displaystyle \ln(2)/\lambda \,}$

üçüncü dörttebirlik : ${\displaystyle \ln(4)/\lambda \,}$

Kullback-Leibler ayrılımı

'Gerçek' üstel dağılım olan Exp(λ₀) ile ('yaklaşık' dağılım) olan Exp(λ) arasında yönlendirilmiş Kullback-Leibler ayrılımı şöyle verilir:

${\displaystyle \Delta (\lambda ||\lambda _{0})=\log(\lambda )-\log(\lambda _{0})+{\frac {\lambda _{0}}{\lambda }}-1.}$

Maksimum entropi dağılımı

[0,∞) and mean μ, de destekli bulunan bütün sürekli olasılık dağılımları arasında sadece λ = 1/μ parametresi ile üstel dağılımın en yüksek entropisi bulunmaktadır.

Üstel rassal değişirlerin minimumu için dağılım

X₁, ..., X_n bağımsız oran parametreleri λ₁, ..., λ_n olan üstel olarak dağılım gösteren rassal değişkenler olsun. Bu halde

${\displaystyle \min\{\,X_{1},\dots ,X_{n}\,\}}$

ifadesi de üstel dağılımdır ve bu dağılımın parametresi

${\displaystyle \lambda =\lambda _{1}+\cdots +\lambda _{n}.\,}$

olur.

Fakat,

${\displaystyle \max\{\,X_{1},\dots ,X_{n}\,\}}$

üstel dağılım göstermez.

Parametre tahmin edilmesi

Verilmiş bir değişkenin üstel dağılım gösterdiği bilinmiş olsun ve oran parametresi olan λnın değerinin tahmin edilmesi gerekmektedir.

Maksimum olabilirlilik

İlgi gösterilen değişkenden bir bağımsız aynen dağılma gösteren örneklem x = (x₁, ..., x_n) olarak seçilsin; o halde λ için olabilirlilik fonksiyonu şöyle verilir:

${\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),}$

burada

${\displaystyle {\overline {x}}={1 \over n}\sum _{i=1}^{n}x_{i}}$

örnek ortalamasıdır.

Olabilirlik fonksiyonunun logaritmasının türevi şudur:

${\displaystyle {\frac {\mathrm {d} }{\mathrm {d} \lambda }}\ln L(\lambda )={\frac {\mathrm {d} }{\mathrm {d} \lambda }}\left(n\ln(\lambda )-\lambda n{\overline {x}}\right)={n \over \lambda }-n{\overline {x}}\ \left\{{\begin{matrix}>0&{\mbox{if}}\ 0<\lambda <1/{\overline {x}},\\\\=0&{\mbox{if}}\ \lambda =1/{\overline {x}},\\\\<0&{\mbox{if}}\ \lambda >1/{\overline {x}}.\end{matrix}}\right.}$

Bu nedenle oran parametresinin maksimum olabilirlilik tahmini şöyle verilir:

${\displaystyle {\widehat {\lambda }}={\frac {1}{\overline {x}}}.}$

Bayes tipi çıkarımsal analiz

Bir üstel dağılımın eşlenik önseli bir gamma dağılımı olur (çünkü üstel dağılım bir özel hal gamma dağılımıdır). Gamma olasılık dağılım fonksiyonunun şu çeşit parametrik tanımı analizde kullanılacaktır:

${\displaystyle \mathrm {Gamma} (\lambda \,;\,\alpha ,\beta )={\frac {\beta ^{\alpha }}{\Gamma (\alpha )}}\,\lambda ^{\alpha -1}\,\exp(-\lambda \,\beta ).\!}$

Bu halde p için sonsal dağılım yukarıda tanımlanan olabilirlilik fonksiyonu ve bir gamma önsel ile şöyle ifade edilebilir:

${\displaystyle p(\lambda )\propto L(\lambda )\times \mathrm {Gamma} (\lambda \,;\,\alpha ,\beta )}$

${\displaystyle =\lambda ^{n}\,\exp(-\lambda \,n{\overline {x}})\times {\frac {\beta ^{\alpha }}{\Gamma (\alpha )}}\,\lambda ^{\alpha -1}\,\exp(-\lambda \,\beta )}$

${\displaystyle \propto \lambda ^{(\alpha +n)-1}\,\exp(-\lambda \,(\beta +n{\overline {x}})).}$

Şimdi p için sonsal yoğunluk bir kayıp olmuş normalizasyon sabiti değerine kadar tanımlanmıştır.

Bunun dağılımı gamma olduğu için bu eksiklik hemen tamamlanabilir ve şu ifade elde edilir:

${\displaystyle p(\lambda )=\mathrm {Gamma} (\lambda \,;\,\alpha +n,\beta +n{\overline {x}}).}$

Burada parametre α önsel gözlemlerin sayısı olarak yorumlanabilir ve β önsel gözlemlerin toplamıdır.

Üstel değişebilirleri üretme

Üstel değişebilirler için üstel dağılım üreten kavramsal olarak bir basit yöntem ters dönüşüm örnekleme dayanır: Verilmiş olan bir birim aralıkta, yani [0,1] arasında, bulunan bir tekdüze dağılımdan çekilmiş U rassal değişebiliri verilmiş olsun,

${\displaystyle T=F^{-1}(U)\!}$

değişebiliri bir üstel dağılım gösterir ve ${\displaystyle F^{-1}}$ ifadesi

${\displaystyle F^{-1}(p)={\frac {-\ln(1-p)}{\lambda }}.\!}$

ile tanımlanmış bir kuantil fonksiyonu olur.

Bunun yanında, eğer U ${\displaystyle (0;1)}$ aralığında bir tekdüze dağılım gösterirse, ${\displaystyle 1-U}$ için de aynı özellik gerçektir. Bu demektir ki şu şekilde üstel değişebilirler üretilebilir:

${\displaystyle T={\frac {-\ln U}{\lambda }}.\!}$

Üstel değişebilirlerin diğer yöntemlerle üretilebilmesi Knuth (1998)de [1] ve Luc Devroye (1986) da [2] görülebilir.

Üstel değişebilirleri üretmek için bir hızlı yöntem zigurat algoritması iledir.

İlişkili dağılımlar

• Bir üstel dağılım bir gamma dağılımının bir özel halidir ve kullanılan parametre setine göre

${\displaystyle \alpha =1}$ veya ${\displaystyle k=1}$ olur.

• Hem bir üstel dağılım ve hem de bir gamma dağılım, faz-tipi dağılımın ozel halleridir.
• Eğer

${\displaystyle Y=X^{1/\gamma }\,}$ ve ${\displaystyle X\sim \operatorname {Ustel} (\lambda ^{-\gamma })}$

ise

${\displaystyle Y\sim \operatorname {Weibull} (\gamma ,\lambda )}$,

olur yani Y Weibull dağılım gösterir. Özellikle, her üstel dağılım da bir Weibull dağılımıdır.

• Eğer

${\displaystyle Y={\sqrt {2X/\lambda }}}$ ve ${\displaystyle X\sim \operatorname {Ustel} (\lambda )}$.

ise

${\displaystyle Y\sim \operatorname {Rayleigh} (1/\lambda )}$

olur; yani Y bir Rayleigh dağılımı gösterir.

• Eğer

${\displaystyle Y=\mu -\beta \log(X/\lambda )\,}$ ve ${\displaystyle X\sim \operatorname {Ustel} (\lambda )}$.

ise

${\displaystyle Y\sim \operatorname {Gumbel} (\mu ,\beta )}$,

olur yani Y Gumbel dağılımı gösterir.

• Eğer iki bağımsız üstel dağılımı olan ${\displaystyle X_{1}\,}$ ve ${\displaystyle X_{2}\,}$ için ${\displaystyle Y=X_{1}-X_{2}\,}$ ise

${\displaystyle Y\sim \operatorname {Laplace} }$

olur yani Y Laplace dağılımı gösterir.

• Bağımsız üstel dağılımlar olan ${\displaystyle X_{i}\,}$ için

${\displaystyle Y=\min(X_{1},\dots ,X_{N})}$

ise

${\displaystyle Y\sim \operatorname {Ustel} }$

olur; yani Y bir üstel dağılım gösterir.

• Eger

${\displaystyle Y=\exp(-X\lambda )\,}$ and ${\displaystyle X\sim \operatorname {Ustel} (\lambda )}$

ise

${\displaystyle Y\sim \operatorname {Tekduze} (0,1)}$,

olur yani Y tekdüze dağılım gösterir.

• Eğer

${\displaystyle X\sim \operatorname {Ustel} (\lambda =1/2)\,\;}$.

ise ${\displaystyle X\sim \chi _{2}^{2}}$ olur yani X için 2 serbestlik derecesi olan ki-kare dağılımı geçerlidir.

• ${\displaystyle X_{1}\dots X_{n}\sim \operatorname {Ustel} (\lambda )\,}$ üstel dağılımlı ve bağımsız olsun ve ${\displaystyle Y=\sum _{i=1}^{n}X_{i}\,}$ olsun; o halde ${\displaystyle Y\sim \operatorname {Gamma} (n,1/\lambda )\,}$
• ${\displaystyle X\sim \operatorname {CarpikLogistik} (\theta )\,}$ ise ${\displaystyle \operatorname {log} (1+e^{-X})\sim \operatorname {Ustel} (\theta )\,}$ olur

Kaynakça

1. Donald E. Knuth (1998). The Art of Computer Programming, Cilt 2: Seminumerical Algorithms, 3. ed. Boston: Addison-Wesley. ISBN 0-201-89684-2. Bak bölüm 3.4.1, say. 133.
2. Luc Devroye (1986). Tekdüze olmayan rassal değişebilir üretimi 5 Mayıs 2009 tarihinde Wayback Machine sitesinde arşivlendi.. New York: Springer-Verlag. ISBN 0-387-96305-7. Bak Bölüm IX 27 Eylül 2007 tarihinde Wayback Machine sitesinde arşivlendi., kısım 2, say. 392–401.