WinBUGS パッケージを使用して差分項目機能モデルに取り組んでいます。より単純なモデルの構築に成功しましたが、「expected collection operator 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))