请高手帮忙看看下面程序哪里出现了问题，公式和完整的程序见附件，非常感谢！！！
无法运行的程序：
另外，如果我想让WINBUGS指定初始值，是不是就不用写Inits?
附件在这里：

MODEL
{
# DINA MODEL
for (i in 1:N) {
theta~dnorm(0,0.0001)
# compute alpha
for (k in 1:M) {
logit(alpha.p[i,k])<-1.7*(theta-L0[k])
alpha[i,k]~dbern(alpha.p[i,j])
}
}
# compute eta
for (i in 1:N) {
for (j in 1:T) {
eta[i,j]<-1
for (k in 1:M) {
eta3[i,j]<-pow(alpha[i,k],q[j,k])
eta[i,j]<-eta[i,j]*eta3[i,j]
p[i,j]<-pow(g[j],(1-eta[i,j]))*pow((1-s[j]),eta[i,j])
response[i,j]~dbern(p[i,j])
}
}
}

#priors
for (j in 1:T) {
s[j]~dunif(0,0.6) #此处应该是4参数的beta分布，但是WINBUGS好像不能处理
g[j]~dunif(0,0.6) #此处应该是4参数的beta分布，但是WINBUGS好像不能处理
}
for (k in 1:M) {
L0[k]~dlnorm(0,0.0001)
}
}
DATA
list(N = 1000, T = 15, M=5,
response = structure(.Data = c(
1,0,1,1,0,1,0,1,0,1,1,0,1,0,0,
……
1,1,1,1,1,1,1,0,1,0,1,1,1,1,1,
1,1,0,1,1,0,1,1,0,0,0,0,0,0,0),.Dim=c(1000,15)),
q = structure(.Data = c(
1,0,0,0,0,
0,1,0,0,0,
0,0,1,0,0,
0,0,0,1,0,
0,0,0,0,1,
1,1,0,0,0,
0,1,1,0,0,
0,0,1,1,0,
0,0,0,1,1,
1,1,1,0,0,
0,1,1,1,0,
0,0,1,1,1,
1,1,1,1,0,
0,1,1,1,1,
1,1,1,1,1), .Dim=c(15,5))
)