我正在使用 runjags 并试图为我的模型获取 DIC。
我知道如果我打印 runjags 对象,它将返回包含 DIC 的模型摘要,但我的模型有很多变量,因此计算摘要统计需要很长时间。当我使用 summarise=FALSE
时,不会返回 DIC。
我试过使用:extract(runJagsObj, what='dic')
但它返回此错误消息:"Error: Please supply column name"
有没有办法在摘要信息中查看 DIC 而不必计算每个变量的摘要统计数据?或者为什么提取物不起作用?
如果这是有用的信息,我正在使用 method=parallel
。
最佳答案
我终于弄清楚了 "Error: Please supply column name"
错误。
我正在使用 ggmcmc 进行绘图,这需要 tidyr。 Tidyr 与 extract
的命名重叠。我能够通过在解决错误的 ggmcmc 之后加载 runjags 来解决命名冲突。
关于r - 无法提取 DIC,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/31838807/