4

我在winbugs中对这段代码有一些问题。该模型在本质上是正确的并且数据已加载，但是当我编译时，软件输出是“节点 Z 的多个定义”。我不知道如何解决这个问题。

这是模型：

#BUGS Model
model {
for (i in 1:n){
  for (j in 1:p){
    Y[i , j] ~ dcat ( prob [i , j , 1: M[j]])
    B <- sum(alpha[j])  
  } 
  theta [i] ~ dnorm (0.0 , 1.0)
}
for (i in 1:n){
  for (j in 1:p){
    for (k in 1:M[j]){
      Z <- sum(delta [k ])
      eta [i , j , k] <- 1.7* alpha [j] * (B * (theta [i] - beta [j] ) + Z)
      exp.eta[i , j , k] <- exp( eta[i , j , k])
      psum[ i , j , k] <- sum(eta[i , j , 1:k])
      prob[i , j , k] <- exp.eta[i , j , k] / psum[i , j , 1:M[j]]
    }
  }
}
for (j in 1:p){
  alpha [j] ~ dnorm (0 , pr.alpha) I(0 , )
  for (k in 2:M[j]){
    delta [k] ~ dnorm (0.0 , 1.0)
  }
  for (k in 1:M[j]){
    beta [j] ~ dnorm (0 , pr.beta )
  }
}
delta [1] <- 0.0
pr.alpha <- pow(1 , -2)
pr.beta <- pow(1, -2)
}

#data
list(n=10, p=8)

M[] M[] M[] M[] M[] M[] M[] M[]
2 2 4 2 2 3 4 2 
2 1 1 2 1 2 2 3
1 2 1 3 1 1 4 4
2 1 1 2 1 1 2 4 
3 4 4 3 3 3 1 1 
4 3 4 4 4 4 4 4 
1 1 2 2 1 2 4 4 
2 1 1 3 1 4 2 4 
3 4 1 1 1 2 2 2 
2 2 2 1 4 4 4 4 
END

感谢所有会回答的人。

4

1 回答 1

3

您的问题在于在 BUGS 循环中多次定义某些节点。例如B在第一个i和j循环中定义 np 次。BUGS 不允许这样做。您不能覆盖节点值。你需要要么

1) 添加一些索引到B、和 ，Z以使 BUGS 能够在循环期间将模拟值存储在节点的元素内。例如用和用替换delta[k]beta[j]BB[i,j]ZZ[i,j,k]

或者

2) 移动B, Z,delta[k]和beta[j]to 仅覆盖它们已有索引的循环。即B，Z不在循环中，因为它们没有索引，delta[k]仅在for(k in 1:...)循环中。

该决定取决于您对模型的想法以及您想要存储的参数。

于 2014-05-28T09:03:22.377 回答