Несколько определений узла W [1]

Question

Я пытаюсь оценить дельта [j, k] на сумме условия (delta [j, 1: 5]) = 0 для каждого «j». но когда я компилирую код, программным выходом является «несколько определений узла W [1]». Может кто-нибудь мне помочь?

model { 
for (j in 1:p){ 
for (k in 1:5){ 
Z[j, k]<- sum(delta[j,1:k]) 
} 
for (i in 1:n){ 
Y[i , j] ~ dcat (prob [i , j , 1: 5]) 
}} 
for (i in 1:n){ 
theta [i] ~ dnorm (0.0 , 1.0) 
} 
for (i in 1:n){ 
for (j in 1:p){ 
for (k in 1:5){ 
eta[i , j , k] <- alpha [j] * (k*theta [i] - k*beta [j]+Z[j, k]) 
psum [i , j , k] <- sum(eta[i , j , 1: k]) 
exp.psum[i , j , k]<- exp(psum [i , j , k]) 
prob [i , j , k] <- exp.psum[i , j , k]/sum(exp.psum [i , j , 1:5]) 
} } } 
for (j in 1:p){ 
W[j] <- sum(delta [j, 1:5]) 
W[j]<- 0 
alpha [j] ~ dlnorm (0.83 , pr.alpha) 
beta [j] ~ dnorm (-1.73 , pr.beta) 
delta[j,1] <- 0.0 
for (k in 2:5){ 
delta [j , k] ~ dnorm (0.02 , pr.delta) 
} } 
pr.alpha <- pow(1.2 , -2) 
pr.beta <- pow(0.7, -2) 
pr.delta <- pow(1.3, -2) 
} 

благодаря

Answer 1

ОШИБКИ не позволяют перезаписывать детерминированных узлы, ... у вас есть есть W[j] <- дважды в последнем for цикле.

Я думаю, что существует много способов написать код для удовлетворения вашего состояния. Например, вы могли бы использовать другое распределение для delta или установить delta[1] быть остаток требуется, чтобы получить все delta просуммировать 0 после моделирования delta[2] к delta[5]