Jeffreys Prior untuk distribusi normal dengan mean dan varian tidak diketahui

12

Saya membaca tentang distribusi sebelumnya dan saya menghitung Jeffrey sebelumnya untuk sampel variabel acak yang terdistribusi normal dengan rerata tidak diketahui dan tidak diketahui. Menurut perhitungan saya, berikut ini berlaku untuk Jeffreys sebelumnya: Di sini, adalah matriks informasi Fisher.

p(μ,σ2)=det(I)=det(1/σ2001/(2σ4))=12σ61σ3.
I

Namun, saya juga sudah membaca publikasi dan dokumen yang menyatakan

seperti Jeffreys sebelumnya untuk kasus distribusi normal dengan mean dan varians tidak dikenal. Apa yang sebenarnya dilakukan Jeffrey sebelumnya?

Nussig
sumber

Jawaban:

7

Saya pikir perbedaan dijelaskan oleh apakah penulis mempertimbangkan kepadatan lebih dari atau kepadatan lebih dari . Mendukung interpretasi ini, hal yang pasti yang ditulis oleh Kass dan Wassermann adalah sedangkan Yang dan Berger menulis σσ2

π(μ,σ)=1/σ2,
π(μ,σ2)=1/σ4.
A. Donda
sumber
2
Terima kasih, saya mengabaikan ini. Namun, ini masih tidak menjelaskan perbedaan antara dan . 1/σ31/σ4
Nussig
3
Sebenarnya, memiliki prior sama dengan memiliki prior , karena properti reparametrization dari Jeffrey sebelumnya: dengan matriks Jacobian dari , yaitu . π(μ,σ)=1/σ2π(μ,σ2)=1/σ3
π(μ,σ)=π(μ,σ2)det(Jf)1σ32σ1σ2
Jff:(μ,σ)(μ,σ2)
Jf=(1002σ)
Nussig
3
@Nussig, saya memeriksa perhitungan, dan saya pikir Anda benar tiba di . Anda juga benar bahwa reparametrization hanya sejumlah faktor . Mempertimbangkan hal ini, perhitunganmu sesuai dengan Kass dan Wassermann, dan aku hanya bisa menebak bahwa Yang dan Berger melakukan kesalahan. Ini masuk akal juga karena yang pertama adalah makalah jurnal yang ditinjau secara teratur dan yang terakhir adalah konsep semacam koleksi formula. 1/σ31/σ
A. Donda
3
Kass dan Wassermann juga mencatat bahwa Jeffreys memperkenalkan aturan yang dimodifikasi, sesuai dengan parameter lokasi dan skala mana yang harus diperlakukan secara terpisah. Ini mengarah ke dan karenanya , tetapi masih belum ke . π(μ,σ)=1/σπ(μ,σ2)=1/σ2π(μ,σ2)=1/σ4
A. Donda
2
Jim Berger masih seorang ilmuwan aktif, jadi untuk memastikan Anda dapat mengecek langsung dengannya: stat.duke.edu/~berger
A. Donda
4

Jawaban yang ada sudah dengan baik menjawab pertanyaan awal. Sebagai seorang fisikawan, saya hanya ingin menambahkan argumen dimensionalitas pada diskusi ini. Jika Anda mempertimbangkan dan untuk menggambarkan distribusi variabel acak dalam ruang 1D nyata dan diukur dalam meter, mereka memiliki dimensi dan . Untuk memiliki prior yang benar secara fisik, Anda perlu memiliki dimensi yang tepat, yaitu satu-satunya kekuatan dimungkinkan secara fisik dalam prior non-parametrik adalah: dan .μσ2[μ]m[σ2]m2σ

π(μ,σ)1/σ2
π(μ,σ2)1/σ3
Dr_Zaszuś
sumber
Mengapa ada dalam ekspresi kedua? σ3
cerebrou
3

1σ3 adalah Jeffreys sebelumnya. Namun dalam prakteknya cukup sering digunakan karena mengarah ke posterior yang relatif sederhana, "intuisi" dari prior ini adalah sesuai dengan flat prior on .1σ2log(σ)

Jorne Biccler
sumber
1
Terima kasih, @Noshgul. Saya mendapatkan poin tentang flat sebelumnya pada . Namun, dapatkah Anda menguraikan 'posterior yang relatif sederhana'? Jika saya tidak salah, hasil Jeffrey sebelumnya adalah posterior inversi normal , yaitu Sebelum harus menghasilkan normal-inverse- posterior, juga, hanya dengan parameter yang berbeda. χ 2 ( μ , σ 2 ) | D N χ - 1 ( ¯ X , n , n , 1log(σ)χ21/σ2χ2
(μ,σ2)|DNχ1(X¯,n,n,1n(XiX¯)2).
1/σ2χ2
Nussig
1
Ooh, ya itu mengarah ke inversi normal . Saya hanya merasa lebih alami bahwa marginal dari adalah inversi dengan n-1 bukannya n derajat kebebasan. Bagaimanapun, saya tentu tidak ingin menyiratkan bahwa prior lainnya akan menyebabkan distribusi yang mengganggu. Sejujurnya saya tidak tahu posterior karya Jeffry sebelumnya, saya juga tidak terlalu memikirkannya ketika saya menulis posting. σ 2 χ 2χ2(X¯,n,n1,s2)σ2χ2
Jorne Biccler