Bayes Tahmin Tahminlerini Anlama


9

Bayes'e Giriş kursuna katılıyorum ve tahmine dayalı dağılımları anlamakta zorluk çekiyorum. Neden faydalı olduklarını anlıyorum ve tanıma aşinayım, ama tam olarak anlamadığım bazı şeyler var.

1) Yeni gözlemlerin bir vektörü için doğru tahmini dağıtım nasıl elde edilir

bir örnekleme modeli ve önceki bir varsayalım . Gözlemler varsayın verilen şartlı bağımsızdır .p(yi|θ)p(θ)yiθ

verilerini gözlemledik ve önceki posterior .D={y1,y2,...,yk}p(θ)p(θ|D)

Yeni gözlemlerin bir vektörünü tahmin etmek istersek , ben bu formülü eşit değildir ki yani öngörülen gözlemler bağımsız değil, değil mi?N={y~1,y~2,...,y~n}

p(N|D)=p(θ|D)p(N|θ)dθ=p(θ|D)i=1np(y~i|θ)dθ,
i=1np(θ|D)p(y~i|θ)dθ,

Bunu söyle Sabit bir için Beta ( ) ve Binom ( ) . Bu durumda, 6 yeni benzetmek istersem, bunu doğru anlarsam, tek bir gözlem için posterior öngörüye karşılık gelen Beta-Binom dağılımından bağımsız olarak 6 çekilişi simüle etmek yanlış olur. Bu doğru mu? Gözlemlerin marjinal olarak bağımsız olmadığını nasıl yorumlayacağımı bilmiyorum ve bunu doğru anladığımdan emin değilim.θ|Da,bp(yi|θ)n,θny~

Posterior öngörülerden simülasyon

Çoğu zaman posterior öngörücüden veri simülasyonu yaptığımızda bu şemayı takip ederiz:

İçin 1 ila :bB

1) den .θ(b)p(θ|D)

2) Daha sonra den yeni verileri .N(b)p(N|θ(b))

Sezgisel görünse de, bu şemanın nasıl çalıştığını kanıtlamayı pek bilmiyorum. Ayrıca, bunun bir adı var mı? Bir gerekçe aramaya çalıştım ve farklı isimler denedim ama şansım yoktu.

Teşekkürler!


Stats.stackexchange.com/questions/72570/… ' da benzer bir soru sordum ama sizinki şu ana kadar daha fazla oy aldı gibi görünüyor.
John

Yanıtlar:


4

Diyelim ki , olduğu göz önüne alındığında koşullu olarak bağımsızdır . Ardından, ki burada ilk eşitlik toplam olasılık yasasından, ikincisi ürün kuralından ve üçüncüsü varsayılan koşullu bağımsızlıktan gelir:X1,,Xn,Xn+1Θ=θ

fXn+1X1,,Xn(xn+1x1,,xn)=fXn+1,ΘX1,,Xn(xn+1,θx1,,xn)dθ
=fXn+1Θ,X1,,Xn(xn+1θ,x1,,xn)fΘX1,,Xn(θx1,,xn)dθ
=fXn+1Θ(xn+1θ)fΘX1,,Xn(θx1,,xn)dθ,
Θ, dağılımını belirlemek için değerlerine ihtiyacımız yok .X1,,XnXn+1

İçin: simülasyon programı doğru , çizmek dağılımından , o zaman çekme dağıtımından . Bu size dağılımından size örneği verir. .i=1,,Nθ(i)ΘX1=x1,,Xn=xnxn+1(i)Xn+1Θ=θ(i){xn+1(i)}i=1NXn+1X1=x1,,Xn=xn


Eğer posterior kestirimi birden fazla dönemde alıyorsanız ne olacak? Ben kullanıyorum her biri için , ama yeni bir teta yeniden çizmek için mantıklı olabilir neden görebilirsiniz. θ(i)xn+j
John

2

Adım adım öngörülü dağılım üretmenin arkasındaki sezgiyi gözden geçirmeye çalışacağım.

Let bir olasılık dağılımı gelen gözlenen verilerin bir vektör olmak ve let geleceğin bir vektör olmak (veya dışı numunesi) biz tahmin etmek istiyoruz değer verir. Biz varsayalım aynı dağıtım gelen . Bu dağıtım hakkında bilgi almak için --- gibi MLE veya MAP tahmini --- en iyi tahminimizi kullanmak cazip olabilir . Ancak, bunu yapmak kaçınılmaz olarak hakkındaki belirsizliğimizi görmezden gelecektir . Bu nedenle, PROCEDE uygun bir şekilde arka dağılımı üzerindeki ortalama etmektir , yani . Dikkat edinyp(y|θ)y~y~yθθθp(θ|y)y~bağımsız verilen aynı dağıtım çekilen bağımsız bir numune olduğu varsayılmaktadır, . Böylece,yθy

p(y~|θ,y)=p(y~,y|θ)p(θ)p(θ,y)=p(y~|θ)p(y|θ)p(θ)p(y|θ)p(θ)=p(y~|θ).

arka kestirimci dağılımı ,y~

p(y~|y)=Θp(y~|θ,y)p(θ|y)dθ=Θp(y~|θ)p(θ|y)dθ

nerede desteğidir .Θθ

Şimdi, örnekleri den nasıl elde edebiliriz ? Açıkladığınız yönteme bazen aşağıdaki gibi çalışan kompozisyon yöntemi denir :p(y~|y)


s = 1,2, ..., S için

'yi dan çizinθ(s)p(θ|y)

çizmek deny~(s)p(y~|θ(s))


Burada, çoğu durumda, zaten den çekimlerimiz vardır , böylece sadece ikinci adım gereklidir.p(θ|y)

Bunun çalışmasının nedeni oldukça basit: İlk önce olduğuna dikkat edin . Bu nedenle, örnekleme bir parametre vektörü den ve daha sonra, örnek, bu vektör kullanılarak gelen eklem dağılımından örnekler verir . Bundan sonra, örneklenen değerler , marjinal dağılımdan örneklerdir .p(y~,θ|y)=p(y~|θ,y)p(θ|y)θ(s)p(θ|y)y~(s)p(y~|θ(s))=p(y~|θ(s),y)p(y~,θ|y)y~(s),s=1,2,...,Sp(y~|y)


1

İlk sorunuza cevap vermek için: evet, değerini bilmiyorsanız gözlemler bağımsız değildir . Diyelim ki, oldukça aşırı bir değeri olduğunu gözlemlediniz . bilinmeyen değerinin aşırı olduğunun bir göstergesi olabilir ve bu nedenle diğer gözlemlerin de aşırı olmasını beklemelisiniz.θy~1θ

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.