本人正在用WinBUGS做一个贝叶斯的问题，由于初次使用这个软件，不是很熟悉，导致目前在程序部分卡主了停滞不前。现在想请熟悉软件操作的大神帮忙做一下，过程很简单，一两个小时就能搞定。我会把样本分布，先验分布，以及样本数据都给出来，你需要做的是只要把这些步骤编成能运行的程序即可。
所需技术：WinBUGS以及贝叶斯的相关知识。
另外你可以给予现金报酬，QQ：1027331092
本人自己编写的程序如下，可以在此基础上进行改进编写，下面的程序大部分是没有错的，所用到的分布也都是简单常见的，但在compile这一步老是卡住，编译不了。
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)
}}

for(i in 2:101){
a~dnorm(0,invsigma2.a)
b~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)))     #中间数据省略掉