model{
c<-10
for(i in 1:n){zeros<-0
zeros~dpois(zeros.mean)
zeros.mean<--log(l)+c
l<-alpha*beta*(pow(y,(beta-1)))*(exp(-alpha*((pow(y,beta)-(pow(t,beta))))))}
for(j in 2:n){log(beta[j])<-log(beta[j-1])+(eta[j])
eta[j]~dnorm(0,lamda)}
s<-sqrt(1/lamda)
lamda~dgamma(1,0.5)
alpha~dgamma(1,0.5)
beta[1]~dlnorm(0,0.1)
}
#Data list(n=19,y=c(0.3945,0.7435,0.9055,1.1075,1.4631,1.8113,1.9136,2.1390,2.3514,2.7121,2.7612,2.7807,2.8761,2.9606,2.9802,3.0674,3.1310,3.1719,3.2531))
#Initial values
list(alpha=44,mu=0.01,tau=0.01,beta=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1))
错在什么地方？