mcmc - 为先验定义约束时，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)模型中。