R , кажется , чтобы иметь возможность вывод хороших краткие заговоры от bugs
и jags
объектов , порожденные функциями R2WinBUGS :: жучков и R2jags: откосы .
Тем не менее, я использую 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
(насколько я могу судить).