list - WinBugs的“Expected collection operator c”错误

Question

我正在使用WinBUGS包来开发不同项目的功能模型。我已经成功构建了更简单的模型，但是我也遇到了错误“期望的集合运算符c”。
如果您能帮助我，我将非常高兴。

# Model
Model
{

for (j in 1:40){
for (k in 1:5){
y[j,k] ~ dbern(p[j,k])
logit(p[j,k])<- u2[j] - beta[k] + gamma[k]*grp[j]
}
}
# Random effect prior
for (j in 1:40){
u2[j] ~ dnorm(0,tau.u2)
}
# Fixed effect prior
for (k in 1:5){
beta[k] ~ dnorm(0, .0001)
gamma[k] ~ dnorm(0,.0001)
}

tau.u2 <- pow(sigma.u2, -2)
sigma.u2 ~ dunif (0, 100)
}

# Data
list(y=structure(
  .Data=c(1, 1, 1, 1, 1, 
            1, 0, 1, 0, 1, 
            0, 1, 1, 1, 1, 
            1, 0, 0, 1, 1, 
            1, 0, 0, 0, 1, 
            0, 1, 1, 1, 0, 
            0, 1, 1, 1, 1, 
            1, 0, 0, 0, 1, 
            1, 0, 1, 1, 1, 
            0, 1, 0, 0, 0, 
            1, 1, 1, 0, 1, 
            0, 1, 1, 0, 1, 
            1, 1, 1, 1, 1, 
            0, 1, 0, 1, 0, 
            1, 1, 1, 1, 1, 
            0, 0, 1, 1, 1, 
            0, 0, 0, 1, 0, 
            1, 1, 1, 0, 0, 
            1, 0, 0, 0, 0, 
            0, 0, 0, 0, 0, 
            0, 1, 0, 0, 0, 
            0, 1, 0, 0, 0, 
            0, 0, 0, 0, 1,
            0, 0, 0, 0, 0, 
            0, 1, 0, 0, 0, 
            0, 0, 0, 0, 0, 
            1, 0, 0, 1, 0, 
            1, 1, 1, 1, 0, 
            1, 1, 0, 0, 1, 
            0, 1, 1, 1, 1, 
            1, 1, 0, 0, 0, 
            0, 0, 0, 1, 0,
            0, 1, 1, 1, 0, 
            1, 0, 1, 1, 0, 
            1, 0, 1, 0, 0, 
            0, 0, 1, 1, 0, 
            0, 0, 1, 1, 0, 
            0, 1, 0, 1, 0, 
            1, 1, 1, 1, 0, 
            0, 0, 0, 1, 1), 
        .Dim=c(40,5)), grp=c(1, 0, 1, 0, 0, 1, 0, 1, 0, 0, 1, 1, 0, 1, 1, 1, 0, 1, 0, 0, 0, 1, 1,      1, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 1, 1, 1, 0))

#Initital values
list(beta=c(0, 0, 0, 0, 0), gamma=c(0, 0, 0, 0, 0), sigma.u2=1, u2=0))

Answer 1

在您的模型中，u2是40个值的 vector ，而您仅为其初始值提供了一个数字。该错误应该将光标置于u2初始值之后，这可以告诉您错误在哪里。 “期望的集合运算符c”意味着它期望看到u2 = c(0,0,0，...，即初始值的 vector )。

尽管由于您已经提供了sigma.u2的初始值，所以让WinBUGS自动生成u2的初始值可能是安全的。