R ดูเหมือนว่าจะสามารถที่จะมีความสุขเอาท์พุทแปลงสรุปจากbugs
และjags
วัตถุที่สร้างขึ้นโดยฟังก์ชั่นR2WinBUGS :: ข้อบกพร่องและR2jags: Jags
อย่างไรก็ตามฉันใช้rjags
แพ็คเกจ เมื่อฉันพยายามพล็อตผลลัพธ์ของฟังก์ชันrjags::coda.samples
โดยใช้R2WinBUGS::plot.mcmc.list
ผลลัพธ์คือพล็อตการวินิจฉัย (ความหนาแน่นของพารามิเตอร์, อนุกรมเวลาของโซ่, ความสัมพันธ์อัตโนมัติ) สำหรับแต่ละพารามิเตอร์
ด้านล่างนี้คือประเภทของพล็อตที่ผมอยากจะผลิตจากกวดวิชาแอนดรูว์เกลแมนของ"วิ่ง WinBuugs และ OpenBugs จาก R" plot.pugs
เหล่านี้ถูกผลิตโดยใช้
ปัญหาคือว่าplot.bugs
ใช้bugs
วัตถุเป็นอาร์กิวเมนต์ในขณะที่ใช้เวลาการส่งออกของplot.mcmc.list
coda.samples
นี่คือตัวอย่าง (จากcoda.samples
):
library(rjags)
data(LINE)
LINE$recompile()
LINE.out <- coda.samples(LINE, c("alpha","beta","sigma"), n.iter=1000)
plot(LINE.out)
สิ่งที่ฉันต้องการก็คือ
- วิธีสร้างพล็อตสรุปสรุปแบบหน้าเดียวที่คล้ายกับข้อมูลที่คล้ายกับที่สร้างโดย
plot.bugs
- ฟังก์ชั่นที่จะแปลง
LINE.out
เป็นวัตถุบั๊กหรือ
mcmc.list
(เท่าที่ฉันสามารถบอกได้)