r - rjags 中的总和为零约束模型

Question

我正在尝试将 ANOVA 模型放入 rjags。模型是这样的

for (r in 1:nE){
  for ( j in 1:nP){
    for ( i in 1:nA){
      logit(p[i,j,r]) <- mu[r] + theta[i,r] + varphi[j,r] + psi[(nA-i)+j,r]
    }
  }
}

我需要适应这些限制

for (r in 1:nE){
  theta[nA,r] <- 0 - sum(theta[1:(nA-1), r])
  varphi[nP,r] <- 0 - sum(varphi[1:(nP-1), r])
  psi[nK,r] <- 0 - sum(psi[1:(nK-1), r])
}

这是该模型的零约束总和。但是，rjags 给了我信息

"Compilation error on line 14. Attempt to redefine node varphi[16,1]"

如果我删除约束部分，模型编译得很好但不会收敛。在 BUGS 中，模型被接受。

如何在 rjags 中实现这些约束？

Answer 1

您有两个选择:

创建新的零和变量(不要重新定义旧的)并监视它们。

重新参数化模型以使用角约束而不是总和为零。

1 是更好的选择，因为新变量可能收敛得更快。见 Ntzoufras' book, (ch.5) ，秒。 5.4.2 讨论以及相关的错误代码。这也应该适用于锯齿，虽然我没有检查过。