Michael dan Fraijo menyarankan bahwa hanya memeriksa apakah nilai parameter dari minat terkandung dalam beberapa wilayah yang kredibel adalah setara Bayesian dari interval kepercayaan pembalik. Saya agak skeptis tentang ini pada awalnya, karena tidak jelas bagi saya bahwa prosedur ini benar-benar menghasilkan tes Bayesian (dalam arti biasa).
Ternyata, itu terjadi - setidaknya jika Anda bersedia menerima jenis fungsi kerugian tertentu. Terima kasih banyak kepada Zen , yang menyediakan referensi ke dua makalah yang membangun koneksi antara wilayah HPD dan pengujian hipotesis:
Saya akan mencoba merangkumnya di sini, untuk referensi di masa mendatang. Dalam analog dengan contoh dalam pertanyaan asli, saya akan memperlakukan kasus khusus di mana hipotesisnya adalah mana adalah ruang parameter.Θ
H0:θ∈Θ0={θ0}andH1:θ∈Θ1=Θ∖Θ0,
Θ
Pereira & Stern mengusulkan metode untuk menguji hipotesis tersebut tanpa harus meletakkan probabilitas sebelumnya pada danΘ 1Θ0Θ1 .
Biarkan menunjukkan fungsi kerapatan dan tentukanθ T ( x ) = { θ : π ( θ | x ) > π ( θ 0 | x ) } .π(⋅)θ
T(x)={θ:π(θ|x)>π(θ0|x)}.
Ini berarti bahwa adalah wilayah HPD , dengan kredibilitas .T(x)P(θ∈T(x)|x)
Tes Pereira-Stern menolak ketika adalah "kecil" ( , katakanlah). Untuk posterior unimodal, ini berarti bahwa jauh di ekor posterior, membuat kriteria ini agak mirip dengan menggunakan nilai-p. Dengan kata lain, ditolak di level jika dan hanya jika tidak terdapat di wilayah HPD . P ( θ ∉ T ( x ) | x ) < 0,05 θ 0 Θ 0 5 % 95 %Θ0P(θ∉T(x)|x)<0.05θ0Θ05 %95 %
Biarkan fungsi tes menjadi jika diterima dan jika ditolak. Madruga et al. mengusulkan fungsi kerugian
dengan .1 Θ 0 0 Θ 0 L ( θ , φ , x ) = { a ( 1 - I ( θ ∈ T ( x ) ) , jika φ ( x ) = 0 b + c I ( θ ∈ ( T ( x )) ) , jika φ ( x ) = 1 , aφ1Θ00Θ0
L(θ,φ,x)={a(1−I(θ∈T(x)),b+cI(θ∈(T(x)),if φ(x)=0if φ(x)=1,
a,b,c>0
Meminimalkan kerugian yang diharapkan mengarah ke tes Pereira-Stern di mana ditolak jikaΘ0P(θ∉T(x)|x)<(b+c)/(a+c).
Sejauh ini, semuanya baik-baik saja. Tes Pereira-Stern sama dengan memeriksa apakah berada di wilayah HPD dan ada fungsi kerugian yang menghasilkan tes ini, yang berarti bahwa tes ini didasarkan pada teori keputusan.θ0
Bagian yang kontroversial adalah bahwa fungsi kerugian tergantung padax . Sementara fungsi kerugian seperti itu telah muncul dalam literatur beberapa kali, mereka tampaknya tidak diterima secara umum sebagai sangat masuk akal.
Untuk membaca lebih lanjut tentang topik ini, lihat daftar makalah yang mengutip Madruga et al. artikel .
Pembaruan Oktober 2012:
Saya tidak sepenuhnya puas dengan fungsi kerugian di atas, karena ketergantungannya pada membuat pengambilan keputusan lebih subjektif daripada yang saya inginkan. Saya menghabiskan lebih banyak waktu untuk memikirkan masalah ini dan akhirnya menulis catatan pendek tentang hal itu, diposting di arXiv sebelumnya hari ini .x
Misalkan menunjukkan fungsi kuantil posterior dari , sehingga . Alih-alih set HPD kami mempertimbangkan interval pusat (sama-tailed) . Untuk menguji menggunakan interval ini dapat dibenarkan dalam kerangka kerja decision-theoretic tanpa fungsi kerugian yang tergantung pada .qα(θ|x)θP(θ≤qα(θ|x))=α(qα/2(θ|x),q1−α/2(θ|x))Θ0x
Triknya adalah merumuskan kembali masalah pengujian hipotesis titik-nol sebagai masalah tiga keputusan dengan kesimpulan terarah. kemudian diuji terhadap dan .Θ0={θ0}Θ0Θ−1={θ:θ<θ0}Θ1={θ:θ>θ0}
Biarkan fungsi tes jika kita menerima (perhatikan bahwa notasi ini adalah kebalikan dari yang digunakan di atas!). Ternyata di bawah fungsi kerugian
tertimbang
the Bayes Tes adalah untuk menolak jika tidak dalam interval pusat.φ=iΘi0−1
L2(θ,φ)=⎧⎩⎨0,α/2,1,if θ∈Θi and φ=i,i∈{−1,0,1},if θ∉Θ0 and φ=0,if θ∈Θi∪Θ0 and φ=−i,i∈{−1,1},
Θ0θ0
Ini sepertinya fungsi kerugian yang cukup masuk akal bagi saya. Saya membahas kerugian ini, kehilangan dan pengujian Madruga-Esteves-Wechsler menggunakan set yang kredibel lebih lanjut dalam naskah di arXiv.