2

我试图了解 JAGS 如何计算偏差和偏差信息标准(DIC)。

如果您只有一个链，R中的dic.samples函数会引发错误。rjags

stop("2 or more parallel chains required")

计算DIC的公式涉及计算模型中参数的预期偏差和预期值处的偏差。我想这样的估计可以从单个链中获得（尽管有时可能存在只能通过多个链来识别的收敛问题）。

该手册指出：

pD 监视器通过比较跨多个链 [1] 的偏差来估计观察到的随机节点对参数有效数量 (pD) [3] 的贡献。它是通过使用选项 type(pD) 创建的。如果模型只有一个链，则无法定义 pD 监视器。

所以我的问题是：

  • 为什么 JAGS 需要至少两条链来计算 DIC？
  • 或者更具体地说，为什么 pD 监视器需要两条链？
4

1 回答 1

1

有许多不同的方法来计算 pD - Martyn Plummer 在“模型复杂性和拟合的贝叶斯测量”论文的讨论中描述了 JAGS 使用的方法：http: //onlinelibrary.wiley.com/doi/10.1111 /1467-9868.00353/摘要

获取 PDF，然后搜索“Plummer” - 第一个结果给出了 JAGS 中使用的计算的详细信息。

于 2015-04-07T06:59:54.697 回答