如何绘制对数后验轨迹（最大比例）？

Question

我正在努力从JAGS模型中获取一些后验轨迹图。获得对数后验轨迹图（达到比例）的最佳方法是什么？

mydata = read.table('aaa', header = F)
jmodel = jags.model(file = '111.model', data = list(x = mydata, N = 50, k=2), n.chains = 2)
samps_f = coda.samples(jmodel, variable.names=c('a'), n.iter=1e3)

Answer 1

我假设您正在使用rjags。 您可以简单地执行plot(samps_f)，因为该对象具有其自己的plot方法。

环顾四周，我发现this vignette可能还会回答您可能遇到的其他问题。