WinBUGS中hierarchial bayesian model问题求助

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

#model
model{
for (i in 1:N){risk[i,1:2] ~ dmnorm(beta[i,1:2],inverse(structure(.Data=c(risk[i,3],risk[i,4],risk[i,5],risk[i,6]),.Dim=c(2,2))))}
for (i in 1:N){beta[i,1:2] ~ dmnorm(mu[1:2],pre.sigma[1:2,1:2])}
mu[1:2] ~ dmnorm(c(0,0),inverse(structure(.Data=c(1,0,0,1),.Dim=c(2,2))))
pre.sigma[1:2,1:2] ~ dwish(structure(.Data=c(1,0,0,1),.Dim=c(2,2)),2)
pre.sigma[1:2,1:2] <- inverse(sigma[1:2,1:2])
}

#data
list(N=3,risk=structure(.Data=c(0.0001164745,-0.0006360886,5.454619*10^(-8),-2.15988*10^(-6),-2.15988*10^(-6),0.02313066,
0.0002386487,-0.0001279721,1.419456*10^(-7),-3.196332*10^(-6),-3.196332*10^(-6),0.03024943,
-6.353768*10^(-5),-0.0001478020,2.418857*10^(-7),-4.84376*10^(-6),-4.84376*10^(-6),0.08127031),.Dim=c(N,6)))
)

#initial values
list(beta=structure(.Data=rep(c(1,0),N),.Dim=c(N,2)),mu=c(1,1),sigma2=structure(.Data=c(1,0,0,1),.Dim=c(2,2)))

运行代码check model时即报错expected variable name,小弟初次接触WinBUGS，请大家多多指教

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Bayesian winbugs WINBUG model Bayes WinBUGS Bayesian hierarchial model

求大侠出手解救啊！！！