R parece ser capaz de generar buenos gráficos de resumen de los objetos bugs
y jags
generados por las funciones R2WinBUGS :: bugs y R2jags: jags .
Sin embargo, estoy usando el rjags
paquete. Cuando trato de trazar los resultados de la función rjags::coda.samples
usando R2WinBUGS::plot.mcmc.list
los resultados, hay diagramas de diagnóstico (densidad de parámetros, series temporales de cadena, autocorrelación) para cada parámetro.
A continuación se muestra el tipo de trama que me gustaría producir, del tutorial de Andrew Gelman "Ejecutar WinBuugs y OpenBugs desde R" . Estos fueron producidos mediante el uso de plot.pugs
.
El problema es que plot.bugs
toma un bugs
objeto como argumento, mientras plot.mcmc.list
toma la salida de coda.samples
.
Aquí hay un ejemplo (del coda.samples
):
library(rjags)
data(LINE)
LINE$recompile()
LINE.out <- coda.samples(LINE, c("alpha","beta","sigma"), n.iter=1000)
plot(LINE.out)
Lo que necesito es
- una forma de generar un diagrama de resumen de una página similar, rico en información, similar al producido por
plot.bugs
- una función que se convertirá
LINE.out
en un objeto de errores o
fuente
mcmc.list
(por lo que puedo decir).