R tampaknya dapat menampilkan plot ringkasan yang bagus dari bugs
dan jags
objek yang dihasilkan oleh fungsi R2WinBUGS :: bugs dan R2jags: jags .
Namun, saya menggunakan rjags
paket itu. Ketika saya mencoba memplot hasil fungsi rjags::coda.samples
menggunakan R2WinBUGS::plot.mcmc.list
hasilnya adalah plot diagnostik (kepadatan parameter, rangkaian waktu rantai, autokorelasi) untuk setiap parameter.
Di bawah ini adalah jenis plot yang ingin saya hasilkan, dari tutorial Andrew Gelman "Menjalankan WinBuugs dan OpenBugs dari R" . Ini diproduksi dengan menggunakan plot.pugs
.
Masalahnya adalah bahwa plot.bugs
mengambil bugs
objek sebagai argumen, sementara plot.mcmc.list
mengambil output dari coda.samples
.
Berikut ini sebuah contoh (dari coda.samples
):
library(rjags)
data(LINE)
LINE$recompile()
LINE.out <- coda.samples(LINE, c("alpha","beta","sigma"), n.iter=1000)
plot(LINE.out)
Yang saya butuhkan adalah baik
- cara untuk menghasilkan plot ringkasan satu halaman yang mirip, kaya informasi, mirip dengan yang dihasilkan oleh
plot.bugs
- sebuah fungsi yang akan dikonversi
LINE.out
ke objek bug atau
sumber
mcmc.list
(sejauh yang saya tahu).