为先验定义约束时，JAGS 中出现错误消息尝试定义节点

Question

我正在尝试从先前的分布中进行采样，并将一些总和和差异限制为大于 0，但是，当我运行通常在 OpenBugs/WinBugs 中运行的此代码时，我收到以下错误“尝试重新定义节点'。我想知道是否有任何可能的方法来定义我的约束以防止发生此错误。任何帮助将不胜感激。

谢谢

S[1,1] ~ dbeta(8.10, 32.81)
d[1,1] ~ dbeta(4.36, 32.61)
d[1,2] ~ dnorm(0.097, 250)
d[2,2] ~ dbeta(2.23, 20.10)

Z1 <- 1
Z1 ~ dbern(constraint) 

#Constraints

constraint <- step(S[1,1] + d[1,2])*
                       step(1 - S[1,1] - d[1,2])*
                       step(S[1,1] - d[1,1])*
                       step(S[1,1] + d[1,2] - d[2,2])

Answer 1

一种方法是在数据中将Z1定义为 1 并将其从模型中取出。然后，您可以在模型中留下Z1~dbern(constraint)。