Bir beta dağılımının iki kantili parametrelerini belirler mi?


9

Açık aralıkta iki kantil ve karşılık gelen konumlarını (her biri) konumlarda bu olan bir beta dağıtımının parametrelerini her zaman bulabilir miyim?(q1,q2)(l1,l2)(0,1)


1
Hayır, temel karşı örnek (q1, q2) = (0,1) ve (l1, l2) = (0,1) parametreleri ne olursa olsun.
Tim

1
@Zaman Sanırım, noktanızı görüyorum, ancak karşı örneğiniz belirttiğim koşulları karşılamıyor (örneğin, konumlar açık aralıkta ). (0,1)
Bota

1
Sanırım bunu sayısal olarak yapabilirsiniz (ve benzersiz bir çözüm olacak), ancak biraz çaba gerektirecektir.
Glen_b-Monica'yı yeniden başlat

1
Ben de düşünüyorum - sayısal çözme zor değil, ama benzersizlik için bir argüman bulmak kolay değil.
Elvis

1
@Elvis aslında, her iki değişkenin (OP ve ) günlüklerine bakarak bunu yapmanın bir yolu olabileceğinden şüpheleniyorum . lq
Glen_b -Reinstate Monica

Yanıtlar:


9

Verilerin bariz tutarlılık gereksinimlerini karşılaması koşuluyla cevap evettir . Argüman basit, basit bir yapıya dayanır, ancak bazı kurulumlar gerektirir. Bir derece göze gerçeğine gelir: parametresi artan bir Beta dağılımı daha geniş için yoğunluğu (PDF) değerini arttırır daha küçük ; ve artırmak tam tersini yapar: ne kadar küçük olursa PDF'nin değeri o kadar artar.a(a,b)xxbx

Ayrıntılar aşağıdadır.


Let, istenen miktarsal be ve istenen miktarsal olarak ile ve (bu yüzden) . Daha sonra Beta dağılımının bu miktarlara sahip olduğu benzersiz ve vardır.q1x1q2x21>q2>q1>01>x2>x1>0ab(a,b)

Bunu göstermedeki zorluk, Beta dağılımının bir inatçı normalleştirici sabit içermesidir. Tanımlamayı hatırlayın: ve için Beta dağılımının yoğunluk işlevi vardır (PDF)a>0b>0(a,b)

f(x;a,b)=1B(a,b)xa1(1x)b1.

Normalleştirme sabiti Beta işlevidir

B(a,b)=01xa1(1x)b1dx=Γ(a)Γ(b)Γ(a+b).

yi ve göre doğrudan ayırt etmeye çalışırsak, her şey dağınık olur bir gösteri girişiminde bulunmanın kaba kuvvet yolu olacaktır.f(x;a,b)ab

Beta işlevini analiz etmekten kaçınmanın bir yolu, niceliklerin göreceli alanlar olduğunu belirtmektir . Yani,

qi=F(xi;a,b)=0xif(x;a,b)dx01f(x;a,b)dx

için . Örneğin, bir Beta dağılımının PDF ve kümülatif dağılım fonksiyonu (CDF) , bunun için ve .i=1,2F(1.15,0.57)x1=1/3q1=1/6

Şekil 1

Yoğunluk fonksiyonu solda çizilmiştir. olan alan sol eğri altındaki kırmızı ile gösterilen, göreceli eğrinin altındaki toplam alan. solundaki alandır , yine kırmızı ve mavi bölgeler toplamına eşit toplam alana göre . Sağdaki CDF ve üzerinde iki ayrı noktanın nasıl işaretlendiğini gösterir.xf(x;a,b)q1x1q2x2(x1,q1)(x2,q2)

Bu şekilde, olarak sabitlendi , olarak seçildi ve daha sonra Beta üzerinde bulunan değeri bulundu CDF.(x1,q1)(1/3,1/6)a1.15b(x1,q1)(a,b)

Lemma : Böyle bir her zaman bulunabilir.b

Spesifik olmak gerekirse, kez ve herkes için düzeltilsin. (Onlar izleyin resimlerde aynı kalır: Her üç durumda da, nispi alanı sola eşittir .) Herhangi biri için , Lemma eşsiz bir değer yoktur iddia , yazılı için , Beta dağılımının .(x1,q1)x1q1a>0bb(a),x1q1(a,b(a))

İlk olarak bu yüzden, not görmek için sıfır, bir yakın değerlere kadar olan tüm olasılık yığınları yaklaşımlar , nereden yaklaşımları . De sonsuz, bir yakın değerlere kadar olan tüm olasılık yığınları yaklaşımlar , nereden yaklaşımları . Arasında olarak işlev kesin artmaktadır .b0F(x1;a,b)1b1F(x1;a,b)0bF(x1;a,b)b

Bu iddia geometrik açıktır: bu demek anlamına gelir biz eğri altındaki sol alanında bakarsak toplam alanı altında göre eğriyi ve bunu için eğrisinin altındaki göreli alanla karşılaştırın , sonra ikinci alan nispeten daha büyük. Bu iki işlevin oranı . Bu işlev için eşit olan istikrarlı bir şekilde bırakarak Fonksiyon nedenle yükseklikleri olan göreceli olarak daha büyük birxxa1(1x)b1 xxa1(1x)b1b>b(1x)bb1x=0,0x=1.xf(x;a,b)yükseklikleri daha için soluna onlar için daha sağındaki Sonuç olarak, bölge solundaki eski olmalıdır nispeten sağındaki alandan daha büyük (Örneğin, bir Riemann toplamı kullanarak titiz bir argümana tercüme etmek kolaydır.)xf(x;a,b)xx1xx1.x1x1.

fonksiyonunun sırasıyla monoton olarak ve sınır değerleri ve olarak arttığını gördük . Aynı zamanda (açıkça) süreklidir. Sonuç olarak çok sayıda vardır ve bu sayı lemma kanıtlayan benzersizdir.bf(x1;a,b)01b0b,b(a)f(x1;a,b(a))=q1

Aynı argüman, arttıkça solundaki alanın arttığını gösterir. bx2 Dolayısıyla değerleri gibi bir numara bazı aralığı boyunca aralığı hemen hemen gelen gelişmeler hemen hemen hiç Sınırı olarak olanf(x2;a,b(a))a0.f(x2;a,b(a))a0q1.

Bir örnek yakın (bu eşittir ). İle ve (önceki şekildeki gibi), ve arasında neredeyse hiç alan yoka00.1x1=1/3q1=1/6b(a)0.02.x1x2:

şekil 2

CDF pratik olarak ve arasında düzdür nedenle pratik olarak sınır olarak ,x1x2,q2q1.a0q2q1.

Diğer uçta ise, yeterince büyük değerler kurşun isteğe bağlı olarak yakın Burada bir örnek daha önce olduğu gibi.aF(x2;a,b(a))1.(x1,q1)

Figür 3

Burada ve yaklaşık Şimdi esasen sağında neredeyse hiç alan yoktura=8b(a)10.F(x2;a,b(a))1:x2.

Sonuç olarak, ile arasında herhangi bir ve olana kadar ayarını Daha önce olduğu gibi, bu benzersiz olmalı, QED .q2q11aF(x2;a,a(b))=q2.a


RÇözüm bulmak için çalışma kodu , iki rasgele noktadan (kantil) ve beta dağıtım parametrelerini belirlemeαβ bölümünde yayınlanmıştır .


Bu cevap, sabit veya seçersek , benzersiz bir karşılık gelen değer bulacağımızı gösterir. , ve sabit bir alana sahip fonksiyonlar oluşturmak mümkün olacaktır . Bunun neden ve setinin benzersiz olduğunu garanti ettiğini hemen anlayamıyorum . Beni ayrıntılandırmaya ve aydınlatmaya istekli misiniz? ab[0,x1][x1,x2][x2,1]αβ
Ocak

@Jan " ve kümesi" ile ne demek istediğinizi açıklayabilir misiniz ? Bu semboller bu iş parçacığında hiçbir yerde görünmez. αβ
whuber
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.