Membuat Bayesian sebelum dari hasil yang sering

13

Bagaimana seharusnya seseorang mengubah hasil yang sering menjadi Bayesian prior?

Pertimbangkan skenario umum yang cukup umum: Eksperimen dilakukan di masa lalu dan hasilnya pada beberapa parameter diukur. Analisis dilakukan dengan metodologi frequentist. Interval kepercayaan untuk diberikan dalam hasil.ϕϕ

Saya sekarang sedang melakukan beberapa percobaan baru di mana saya ingin mengukur beberapa parameter lain, misalnya dan . Eksperimen saya berbeda dari penelitian sebelumnya --- tidak dilakukan dengan metodologi yang sama. Saya ingin melakukan analisis Bayesian, jadi saya perlu menempatkan prior pada dan .θϕθϕ

Tidak ada pengukuran sebelumnya yang telah dilakukan, jadi saya menempatkan yang tidak informatif (katakan seragamnya) sebelumnya. θ

Seperti disebutkan, ada hasil sebelumnya untuk , diberikan sebagai interval kepercayaan. Untuk menggunakan hasil itu dalam analisis saya saat ini, saya perlu menerjemahkan hasil frequentist sebelumnya menjadi informatif sebelum analisis saya. ϕ

Salah satu opsi yang tidak tersedia dalam skenario ini adalah mengulangi analisis sebelumnya yang mengarah ke pengukuran dengan cara Bayesian. Jika saya bisa melakukan ini, akan memiliki posterior dari percobaan sebelumnya yang kemudian akan saya gunakan sebagai sebelumnya, dan tidak akan ada masalah.ϕ ϕ

Bagaimana saya harus menerjemahkan CI yang sering ke dalam distribusi Bayesian sebelumnya untuk analisis saya? Atau dengan kata lain, bagaimana saya bisa menerjemahkan hasil yang paling sering pada ke posterior di yang kemudian akan saya gunakan sebagai prioritas dalam analisis saya?ϕϕϕ

Wawasan atau referensi apa pun yang membahas masalah jenis ini dipersilakan.

bill_e
sumber
Distribusi sebelum, atau posterior?
Tim
diedit untuk kejelasan, lebih baik?
bill_e
Bisakah Anda memiliki seragam dari-infinity hingga + infinity
mdewey
Tidak yakin apa hubungannya dengan meta-analisis. Bisakah Anda mengklarifikasi
mdewey
3
Anda sedang mencari gaya priors, Welch and Peers yang cocok. Lihatlah ulasan ini: projecteuclid.org/euclid.lnms/1215091929
Zen

Jawaban:

3

Versi singkat: Ambil Gaussian yang berpusat pada perkiraan sebelumnya, dengan std. dev. sama dengan CI.

Versi lama: Let menjadi nilai sebenarnya dari parameter, dan membiarkan φ perkiraan yang Anda miliki. Asumsikan sebuah apriori sebelum seragam P ( φ ) = c t . Anda ingin tahu distribusi φ 0 mengingat bahwa perkiraan φ telah diperoleh:ϕ0ϕ^P(ϕ)=ctϕ0ϕ^

Sekarang satu-satunya ketergantungan padaφ0adalah dalam jangkaP( φ |φ0), sisanya adalah normalisasi konstan. Dengan asumsi φ adalah estimator maksimum likelihood (atau beberapa estimator yang konsisten lainnya), kita dapat menggunakan fakta-fakta berikut:

P(ϕ0|ϕ^)=P(ϕ^|ϕ0)P(ϕ0)P(ϕ^)=P(ϕ^|ϕ0)ctP(ϕ^)
ϕ0P(ϕ^|ϕ0)ϕ^
  1. Seiring dengan meningkatnya jumlah pengamatan, MLE secara Gaussian asimptotik,
  2. Secara asimptotik tidak bias (berpusat pada nilai sebenarnya ),ϕ0
  3. Ia berfluktuasi sekitar dengan varian sama dengan Informasi Fisher terbalik dari pengamatan sebelumnya, dan itulah yang akan saya gunakan sebagai CI (kuadrat).ϕ0

Cara lain untuk menggambarkannya: Posterior Bayesian dan distribusi estimator yang konsisten dan efisien menjadi sama tanpa gejala.

Alex Monras
sumber
Saya harus menambahkan bahwa solusi ini adalah untuk 68% CI, yaitu 1 sigma. Jika interval kepercayaan Anda adalah 95% Anda berada di dua sigma, jadi Anda harus membagi CI dengan 2, jika mereka di 99,7%, maka mereka adalah 3 sigmas, jadi Anda harus membaginya dengan 3. en.wikipedia.org/wiki/ 68% E2% 80% 9395% E2% 80% 9399.7_rule
Alex Monras
Saya berkomentar persis apa yang ada di komentar Anda :-) Mungkin Anda harus menambahkannya ke balasan Anda. Saya akan ...
Rolazaro Azeveires
1

Tergantung. Dalam beberapa kasus sederhana dengan data terdistribusi normal, ketika Anda memiliki frequentist interval kepercayaan didasarkan pada -Distribusi, yang sesuai posterior marginal dari analisis Bayesian akan menjadi bergeser, Rescaled mahasiswa t -Distribusi dengan quantiles cocok dengan batas kepercayaan frequentist, lihat https : //en.wikipedia.org/wiki/Student%27s_t-distribution#Bayesian_inference . Demikian pula, jika Anda memiliki interval kepercayaan berulang untuk beberapa parameter varians σ 2 yang diturunkan melalui distribusi chi-square dari beberapa kuantitas penting seperti S 2 ( n - p ) / σ 2ttσ2S2(np)/σ2, posterior marginal Bayesian yang sesuai akan menjadi "kuadrat rescaled chi-square" (distribusi Gamma terbalik), sekali lagi dengan kuantil yang cocok dengan batas kepercayaan sering (dengan ketentuan bahwa skala non-informatif sebelum digunakan). 1/σ2

Jarle Tufto
sumber