R sembra essere in grado di produrre simpatici grafici di riepilogo dagli oggetti bugs
e jags
generati dalle funzioni R2WinBUGS :: bugs e R2jags: jags .
Tuttavia, sto usando il rjags
pacchetto. Quando provo a tracciare i risultati della funzione rjags::coda.samples
usando R2WinBUGS::plot.mcmc.list
i risultati sono grafici diagnostici (densità dei parametri, serie temporali della catena, autocorrelazione) per ciascun parametro.
Di seguito è riportato il tipo di trama che vorrei produrre, dal tutorial di Andrew Gelman "Esecuzione di WinBuugs e OpenBugs da R" . Questi sono stati prodotti utilizzando il plot.pugs
.
Il problema è che plot.bugs
accetta un bugs
oggetto come argomento, mentre plot.mcmc.list
prende l'output di coda.samples
.
Ecco un esempio (dal coda.samples
):
library(rjags)
data(LINE)
LINE$recompile()
LINE.out <- coda.samples(LINE, c("alpha","beta","sigma"), n.iter=1000)
plot(LINE.out)
Neanche quello di cui ho bisogno
- un modo per generare un grafico riassuntivo di una pagina simile, ricco di informazioni, simile a quello prodotto da
plot.bugs
- una funzione che verrà convertita
LINE.out
in un oggetto bug o
mcmc.list
(per quanto posso dire).