刚开始学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，模型应该怎么写呢？