Winbugs 代码编写求助

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

本人运行了依据自己编写的一段程序，在模型check和load data 都没有问题，但是到了compile这一步，一点击下去，WinBUGS就马上没有反应了，卡住了，不知道是什么问题，哪位高手能否帮忙看一下。还有我用的的数据d和e是两个19X101的矩阵，我的数据在data list 里面那样输入不知道是否正确？
或者有哪位对WinBUGS比较熟练可以帮我实现这个代码，可以支付报酬。我可以把条件分布，先验分布和数据都整理出来，这个代码应该不难，两个小时之内应该可以搞定，可以联系QQ:1027331092.

程序如下：

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)))     #中间数据省略掉

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：winbugs WINBUG BUGS bug Win 酒店英文 迈阿密 程序 故事 艺术

我读了一下（注意：我没去理解你的程序做什么，只是看语法），至少发现以下语法错误（我不理解为何你说CHECK没问题，check应该报错的）：
1. 第5行，即 lambda[i,j]<-e[i,j]*exp(a[j]+b[j]*k ，最后一个k不合语法, 不应该是单个变量，应该是个数组中元素，例如k[j]或k[i]之类
2. 第7行，a~dnorm(0,invsigma2.a)，同样第一个a不合语法，不应该是单个变量，应该是个数组中元素，例如a[i]之类
3.第8行，同理，同样第一个b不合语法.