R , R2WinBUGS :: bugs ve R2jags: jags işlevleri tarafından oluşturulan nesnelerden bugs
ve jags
nesnelerden güzel özet grafikler çıkarabiliyor gibi görünüyor .
Ancak, rjags
paketi kullanıyorum . Sonuçları rjags::coda.samples
kullanarak fonksiyonun sonuçlarını çizmeye çalıştığımda R2WinBUGS::plot.mcmc.list
, her parametre için teşhis grafikleri (parametre yoğunluğu, zincir zaman serisi, otokorelasyon) vardır.
Andrew Gelman'ın "WinBuugs ve OpenBugs'ı R'den Çalıştırma" dersinden üretmek istediğim arsa türü aşağıdadır . Bunlar plot.pugs
.
Sorun, plot.bugs
bir bugs
nesneyi bağımsız değişken olarak alırken plot.mcmc.list
çıktısını alır olmasıdır coda.samples
.
İşte bir örnek (from coda.samples
):
library(rjags)
data(LINE)
LINE$recompile()
LINE.out <- coda.samples(LINE, c("alpha","beta","sigma"), n.iter=1000)
plot(LINE.out)
İhtiyacım olan şey de
- tarafından üretilene benzer, bilgi açısından zengin, tek sayfalık bir özet grafik oluşturmanın bir yolu
plot.bugs
- bir
LINE.out
hata nesnesine dönüştürülecek bir işlev veya
mcmc.list
(anlayabildiğim kadarıyla) ile başlama problemini çözmez .