【发布时间】:2015-03-21 00:49:19
【问题描述】:
我试图了解 JAGS 如何计算偏差和 Deviance Information Criterion (DIC)。
如果您只有一个链,R 中 rjags 中的 dic.samples 函数会引发错误。
stop("2 or more parallel chains required")
计算DIC的公式涉及计算模型中参数的预期偏差和预期值处的偏差。我想这样的估计可以从一条链中获得(尽管有时可能会出现收敛问题,只能通过一条以上的链来识别)。
手册指出:
pD 监控器估计对有效数量的贡献 通过比较观察到的随机节点的参数 (pD) [3] 跨多个链的偏差 [1]。它是通过使用创建的 选项类型(pD)。如果模型只有一个链,那么一个 pD 监视器 无法定义。
所以我的问题是:
- 为什么 JAGS 需要至少两条链来计算 DIC?
- 或者更具体地说,为什么 pD 监视器需要两条链?
【问题讨论】: