本人运行了依稀自己编写的一段程序，在模型check和load data 都没有问题，但是到了compile这一步，一点击下去，WinBUGS就马上没有反应了，卡住了，不知道是什么问题，哪位高手能否帮忙看一下。还有我用的的数据d和e是两个19X101的矩阵，我的数据在data list 里面那样输入不知道是否正确？

程序如下：

MODEL{

for(i in 1:19){
for(j in 1:101){
d[i,j]~dpois(lambda[i,j])
lambda[i,j]<-e[i,j]*exp(a[j]+b[j]*k[i])
}}

for(i in 2:101){
a[i]~dnorm(0,invsigma2.a)
b[i]~dnorm(0.0099,invsigma2.b)
}

a[1]~dnorm(0,invsigma2.a)
b[1]<-1-sum(b[2:101])


k[3]<-5*mu+k[1]+w[2]
r1<-invsigma2/5
w[2]~dnorm(0,r1)


k[2]<-3*mu+k[2]+w[1]
r2<-invsigma2/3
w[1]~dnorm(0,r2)


k[1]<-(0-sum(k[2:19]))


mu~dnorm(mu0,invsigma2.mu)
invsigma2~dgamma(alpha,beta)


for(j in 1:16){
k[j+3]<-mu+k[j+2]+w[j+2]
}

for(j in 3:18){
w[j]~dnorm(0,invsigma2)
}

}

data list(invsigma2=4,mu0=10,alpha=2.1,beta=4,invsigma2.mu=8,
invsigma2.a=5,invsigma2.b=6,  d=structure(.Data=c(353,......1406),.Dim=c(19, 101)),e= structure(.Data=c(167199.74,......411.01),.Dim=c(19, 101)))     #中间数据省略掉