在WinBUGS中定义模型的一些问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

刚开始学R和WinBUGS，在WinBUGS中定义模型时遇到了困难。首先我在R里面模拟了一批数据，代码如下：
pnum<-100
inum<-20
hnum<-2
pdf<-pnum-1
idf<-inum-1
hdf<-hnum-1
pvar<-16
ivar<-4
hvar<-1
pivar<-64
phvar<-2
ihvar<-3
pihvar<-144
pz<-rnorm(pnum,0,(pvar)^0.5)
iz<-rnorm(inum,0,(ivar)^0.5)
hz<-rnorm(hnum,0,(hvar)^0.5)
piz<-array(rnorm((pnum*inum),0,(pivar)^0.5),dim=c(pnum,inum))
phz<-array(rnorm((pnum*hnum),0,(phvar)^0.5),dim=c(pnum,hnum))
ihz<-array(rnorm((inum*hnum),0,(ihvar)^0.5),dim=c(inum,hnum))
pihz<-array(rnorm((pnum*inum*hnum),0,(pihvar)^0.5),dim=c(pnum,inum,hnum))
x<-array(0,dim=c(pnum,inum,hnum))
for(p in 1:pnum){
for(i in 1:inum){
for(h in 1:hnum){
x[p,i,h]<-0+pz[p]+iz[i]+hz[h]+piz[p,i]+phz[p,h]+ihz[i,h]+pihz[p,i,h]
}
}
}


然后我想通过这批模拟数据在WinBUGS中估计pvar,ivar,hvar,pivar,phvar和pihvar，模型应该怎么写呢？

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：winbugs WINBUG BUGS bug Win WinBUGS 模型