Pekerjaan rumah: Analisis Data Bayesian: Priors pada kedua parameter binomial

8

Berikut ini adalah masalah dari Bayesian Data Analysis 2nd ed , hal. 97. Andrew Gelman belum memasukkan solusinya dalam panduan di situs webnya dan itu membuatku gila sepanjang hari. Secara harfiah sepanjang hari.

Untuk beberapa data , dimodelkan sebagai distribusi binomial dengan parameter populasi dan probabilitas , yang keduanya tidak diketahui. Masalahnya mengatur pertanyaan dengan informasi ini: (1) Menentukan sebelum pada sulit, karena hanya mengambil bilangan alami positif, sehingga diperlakukan sebagai , di mana \ mu tidak dikenal. (2) Untuk mendefinisikan prior on (N, \ theta) , kita memiliki \ lambda = \ mu \ theta . (Logikanya di sini adalah bahwa mungkin lebih mudah untuk merumuskan suatu pertimbangan mengingat harapan pengamatan tanpa syarat, daripada rata-rata N yang tidak teramati.yNθNPr(N|μ)=Poisson(μ)μ(N,θ)λ=μθN.) (3) Potensial noninformatif sebelum adalah p(λ,θ)1/λ .

Bagian dari masalah yang saya tutupi adalah bagaimana mengubah variabel dan menentukan p(N,θ) .

Pendekatan yang saya coba adalah menulis p(N,θ|λ)p(λ,θ) , dan menghilangkan \ lambda yang tidak diinginkan λmelalui integrasi, yaitu p(N,θ)=0CμN/(exp(μ)λN!)dλ , dan mengganti μ dengan relasi μ=λ/θ . Pendekatan ini direduksi menjadi p(N,θ)=C/(N+1) , di mana C adalah konstanta proporsionalitas yang diperkenalkan dari (3).

Hasil ini menjadi perhatian saya, karena ini menyiratkan bahwa probabilitas gabungan dari beberapa nilai θ dan N hanya tergantung pada N , dan bukan pada θ . Selain itu, beberapa lonceng yang tidak jelas terdengar dari kalkulus multivariabel saya yang cukup tua, berusaha mengingatkan saya tentang Jacobian dan mengoordinasikan transformasi, tetapi saya tidak yakin bahwa pendekatan integrasi ini bahkan tepat.

Saya menghargai bantuan dan wawasan Anda.

Sycorax berkata Reinstate Monica
sumber
Dalam hal ini, mengapa tidak mengirim email ke Andrew? Dia mungkin ingin memperbaiki kelalaian.
Glen_b -Reinstate Monica

Jawaban:

10

Saya melakukan semua pertanyaan dari empat bab pertama enam tahun lalu. Inilah yang saya miliki:

p(μ,θ)|λμ|p(λ,θ)=μ1.

Begitu

p(N,θ)=0p(μ,N,θ)dμ=0p(μ,θ)Pr(N|μ)dμ0μ1(μNN!eμ)dμ=(N1)!N!=N1

Anda tidak perlu khawatir bahwa tidak bergantung pada . Ini hanya berarti bahwa prior untuk adalah seragam pada , yang keren untuk parameter Bernoulli.p(N,θ)θθ[0,1]

Cyan
sumber