model{
#Likelihood:
for(i in 1:n){
count[i]~dpois(mu[i])
mu[i]<-u[i]*lambda[i]
u[i]~dbern(p[i])
logit(p[i])<-alpha00+alpha[1,X[i,1]]+alpha[2,X[i,2]]+alpha[3,X[i,3]]+alpha[4,X[i,4]]
logit(lambda[i])<-beta00+beta[1,X[i,1]]+beta[2,X[i,2]]+beta[3,X[i,3]]+beta[4,X[i,4]]
zdp[i]<-1-p[i]+p[i]*exp(-lambda[i])}
mzdp<-mean(zdp[])
#Priors;
for(j in 1:4){
alpha[j,1]~dnorm(0,0.001)   alpha[j,2]<-0
alpha[j,3]~dnorm(0,0.001)
beta[j,1]~dnorm(0,0.001)   beta[j,2]<-0
beta[j,3]~dnorm(0,0.001) }
alpha00~dnorm(0,0.001)
beta00~dnorm(0,0.001)}
list(n=40,X=structure(.Data=c(
0.43,
0.61,
0.29,
0.16,
0.70,
0.35,
0.45,
0.05,
0.10,
0.14,
0.04,
0.62,
0.15,
0.72,
0.93,
0.58,
0.26,
0.72,
0.04,
0.10,
0.73,
0.23,
0.22,
0.61,
0.95,
0.24,
0.34,
0.05,
0.74,
0.12,
0.06,
0.19,
0.10,
0.68,
0.69,
0.87,
0.35,
0.50,
0.70,
0.21,
0.77,
0.07,
0.53,
0.66,
0.94,
0.32,
0.52,
0.77,
0.13,
0.56,
0.98,
0.33,
0.54,
0.51,
0.28,
0.35,
0.92,
0.41,
0.10,
0.33,
0.62,
0.64,
0.75,
0.20,
0.14,
0.45,
0.37,
0.49,
0.10,
0.73,
0.04,
0.17,
0.29,
0.67,
0.76,
0.80,
0.70,
0.05,
0.39,
0.23,
0.73,
0.88,
0.18,
0.28,
0.70,
0.90,
0.12,
0.51,
0.12,
0.92,
0.52,
0.27,
0.64,
0.64,
0.85,
0.07,
0.41,
0.57,
0.84,
0.19,
0.24,
0.19,
0.99,
0.27,
0.67,
0.67,
0.60,
0.93,
0.84,
0.95,
0.54,
0.87,
0.71,
0.38,
0.48,
0.01,
0.19,
0.58,
0.98,
0.66,
0.74,
0.29,
0.37,
0.66,
0.04,
0.81,
0.56,
0.60,
0.76,
0.11,
0.38,
0.75,
0.95,
0.44,
0.71,
0.58,
0.05,
0.53,
0.36,
0.48,
0.98,
0.68,
0.70,
0.03,
0.01,
0.78,
0.89,
0.53,
0.20,
0.68,
0.44,
0.41,
0.50,
0.27,
0.06,
0.56,
0.76,
0.06,
0.72,
0.16),
.Dim=c(40,4)))
list(
alpha=structure(.Data=c(1.000,0,1.000,1.000,0,1.000,1.000,0,1.000,1.000,0,1.000),
.Dim=c(4,3)),
beta=structure(.Data=c(0.500,0,0.500,0.500,0,0.500,0.500,0,0.500,0.500,0,0.500),
.Dim=c(4,3)),
alpha00=1.000,beta00=0.500)
上面模型应该是对的，是在文献上的代码，但是数据是我自己弄得，想让X~dunif(0,1)
,初始部分也是自己弄得，compile时老说“array index is not an integer”
请大家帮忙解答下哈。。。。