r - 如何绘制对数后验的轨迹（直至成比例）？

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)，因为对象有自己的绘图方法。

环顾四周，我发现这个小插曲也可以回答您可能遇到的其他问题。