Question

我尝试估计每个“j”的条件和（delta [j，1：5]）= 0的delta [j，k]。但是当我编译代码时，软件输出是“节点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

BUGS不允许您覆盖确定性节点，...在最后一个W[j] <-循环中，您有for次。

我想有很多方法可以编写满足您条件的代码。例如，您可以使用delta的不同分布或将delta[1]设置为在模拟delta到{{1}后将所有delta[2]求和为0所需的余数}

节点W的多个定义[1]

1 个答案: