r - 估计 zoib beta 回归模型的 DIC

Question

我有下面的数据和代码，希望您能帮助估计DIC（或AIC）β回归模型zoib：

library(zoib)
data("GasolineYield", package = "zoib")
re.md <- zoib(yield ~ temp | 1 | 1, data=GasolineYield, 
joint = FALSE, random=1, EUID=GasolineYield$batch, 
zero.inflation = FALSE, one.inflation = FALSE, 
n.iter=3200, n.thin=15, n.burn=200)
sample2 <- re.md$coeff
summary(sample2)

按照包的小插图和作者发表的文章，我实现了下面的代码来获取，DIC但我得到一个错误：

dic.samples(sample2, n.iter=201, thin = 1, type="pD")
Error in nchain(model) : Invalid JAGS model object in nchain

我用文档中的一个对象运行了相同的代码coda.samples {rjags}并且它有效：

data(LINE)
LINE$recompile()
LINE.out <- coda.samples(LINE, c("alpha","beta","sigma"), n.iter=1000)
summary(LINE.out)
dic.samples(LINE, n.iter=201, thin = 1, type="pD")

然后我检查了上面每个对象的类LINE，sample2，re.md， 和LINE.out，输出是：

class(LINE)
"jags"
class(re.md)
"list"
class(sample2)
"mcmc.list"
class(LINE.out)
"mcmc.list"

这表明我的错误可能是因为该对象sample2不是 class jags。

因此，我将不胜感激有关如何获得sample2可以接受的形式的任何想法dic.samples ，以便能够获得DIC（或AIC）我的rd.md模型。