R dường như có thể xuất ra các sơ đồ tóm tắt đẹp từ bugs
và jags
các đối tượng được tạo bởi các hàm R2WinBUGS :: bug và R2jags: jags .
Tuy nhiên, tôi đang sử dụng rjags
gói. Khi tôi cố gắng vẽ kết quả của hàm rjags::coda.samples
bằng R2WinBUGS::plot.mcmc.list
các kết quả là các ô chẩn đoán (mật độ tham số, chuỗi thời gian chuỗi, tự động tương quan) cho từng tham số.
Dưới đây là loại cốt truyện mà tôi muốn sản xuất, từ hướng dẫn "Chạy WinBuugs và OpenBugs từ R" của Andrew Gelman . Chúng được sản xuất bằng cách sử dụng plot.pugs
.
Vấn đề là plot.bugs
lấy một bugs
đối tượng làm đối số, trong khi plot.mcmc.list
lấy đầu ra của coda.samples
.
Đây là một ví dụ (từ coda.samples
):
library(rjags)
data(LINE)
LINE$recompile()
LINE.out <- coda.samples(LINE, c("alpha","beta","sigma"), n.iter=1000)
plot(LINE.out)
Cái tôi cần là
- một cách để tạo ra một cốt truyện tóm tắt tương tự, giàu thông tin, tương tự như một trang được tạo bởi
plot.bugs
- một chức năng sẽ chuyển đổi
LINE.out
thành một đối tượng lỗi hoặc
mcmc.list
(theo như tôi có thể nói).