Сообщение об ошибке в JAGS при определении ограничений для априорных точек Попытка переопределить узел

Я пытаюсь взять образцы из предыдущих дистрибутивов и ограничить некоторые суммы и различия, чтобы они были больше 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])

mcmc jags rjags winbugs openbugs
person Martin    schedule 23.03.2020    source источник

Ответы (1)


arrow_upward
0
arrow_downward

Один из способов - определить Z1as 1 в данных и извлечь его из модели. Затем вы можете оставить Z1~dbern(constraint) в модели.

person DaveArmstrong    schedule 01.10.2020