请教关于WinBUGS软件的model程序问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

以下是一个WinBUGS软件的model程序，是参照范例Rats: a normal hierarchical model写成的，当WinBUGS软件在进行模型检测的时候，总是出现"expected a comma"的错误信息，请各位专家帮忙看看model程序的问题在哪里？导致出现"expected a comma"的错误信息，多谢！
model{
#specify the multilevel model, z and d are the number of first-level units and second-level units, respectively.
for(j in 1:z){
for(i in 1:d){
#specify the first level model: m[j,i]=beta3[j]+e2[j,i]
m[j,i] ~ dnorm(beta3[j], prec.m)
#specify the first level model: y[j,i]=beta1[j]+beta2[j]*m_1[j,i]+e1[j,i]  
y[j,i] ~ dnorm(mean.y[j,i], prec.y)
mean.y[j,i]<- beta1[j]+beta2[j]*m_1[j,i]
}
#specify the second level model: beta2[j]=beta2+e3[j]
beta2[j] ~ dnorm(beta2, prec.beta2)
#specify the second level model: beta3[j]=beta3+alpha*x[j]+e4[j]
beta3[j] ~ dnorm(mean.beta3[j], prec.beta3)
mean.beta3[j]<-beta3+alpha*x[j]
#specify the second level model: beta1[j]=beta1+tau.p*x[j]+beta*m_2[j]+e5[j]
beta1[j] ~ dnorm(mean.beta1[j], prec.beta1)
mean.beta1[j]<-beta1+tau.p*x[j]+beta*m_2[j]
}
#prior distribution of parameters.
beta2 ~ dnorm(0,1.0E-6)
beta1 ~ dnorm(0,1.0E-6)
tau.p ~ dnorm(0,1.0E-6)
beta ~ dnorm(0,1.0E-6)
beta3 ~ dnorm(0,1.0E-6)
alpha ~ dnorm(0,1.0E-6)
prec.y ~ dgamma(0.001,0.001)
prec.m ~ dgamma(0.001,0.001)
prec.beta1 ~ dgamma(0.001,0.001)
prec.beta2 ~ dgamma(0.001,0.001)
prec.beta3 ~ dgamma(0.001,0.001)
#define the effects as function of parameters
theta <- alpha*beta
}

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Winbugs软件 winbugs WINBUG model mode 请教 程序 软件 model winbugs

问题已经解决，谢谢各位专家！

你好！楼上的同学， 加一下我qq 327871324，很着急想请教一下软件的使用。非常感谢！