R scheint von der zur Ausgabe nette Zusammenfassung Plots der Lage zu sein , bugs
und jags
Objekte , die von den Funktionen erzeugt R2WinBUGS :: Bugs und R2jags: jags .
Ich benutze jedoch das rjags
Paket. Wenn ich versuche, die Ergebnisse der Funktion rjags::coda.samples
anhand R2WinBUGS::plot.mcmc.list
der Ergebnisse zu zeichnen, handelt es sich um Diagnosediagramme (Parameterdichte, Kettenzeitreihen, Autokorrelation) für jeden Parameter.
Unten ist die Art der Handlung aufgeführt, die ich aus Andrew Gelmans Tutorial "Ausführen von WinBuugs und OpenBugs von R" erstellen möchte . Diese wurden unter Verwendung der hergestellt plot.pugs
.
Das Problem ist, dass plot.bugs
ein bugs
Objekt als Argument verwendet wird, während plot.mcmc.list
die Ausgabe von verwendet wird coda.samples
.
Hier ist ein Beispiel (aus dem coda.samples
):
library(rjags)
data(LINE)
LINE$recompile()
LINE.out <- coda.samples(LINE, c("alpha","beta","sigma"), n.iter=1000)
plot(LINE.out)
Was ich brauche ist auch nicht
- Eine Möglichkeit, ein ähnliches, informationsreiches, einseitiges Zusammenfassungsdiagramm zu erstellen, das dem von ähnelt
plot.bugs
- eine Funktion, die
LINE.out
in ein Bugs-Objekt konvertiert wird oder
quelle
mcmc.list
(soweit ich das beurteilen kann).