为什么模型跑不了，是数据的问题吗？求助，急！谢谢

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

model{
for(t in 1:T){
   for(j in 1:N){
     x[j,t]~dnorm(mux[j,t],psi[j])
     ephat[j,t]<-x[j,t]-mux[j,t]
}
     mux[1,t]<-lam01*E[t]+lam11
     mux[2,t]<-lam22*E[t]+lam21
     mux[3,t]<-lam32*E[t]+lam31
     mux[4,t]<-lam42*E[t]+lam41
     E[t]~dnorm(0,psd)
}
#priors on observation model
    lam01~dnorm(0.0,1.0)
    lam11~dnorm(0.0,1.0)
    lam21~dnorm(0.0,1.0)
    lam31~dnorm(0.0,1.0)
    lam41~dnorm(0.0,1.0)
    lam22~dnorm(0.0,1.0)
    lam32~dnorm(0.0,1.0)
    lam42~dnorm(0.0,1.0)
    psd~dgamma(1.0,1.0)
   for(j in 1:N){
    psi[j]~dnorm(1.0,1.0)
}
}#end of model

Data
list(T=31,N=4,  
  x=structure(
      .Data=c(
20.224 , 0.215 , 1.240 , 2.252 , 4.256 , 0.266 , 0.309 , 0.351 , 0.377 , 0.403 , 0.432 , 0.439 , 0.454 ,
0.438 , 0.485 , 0.531 , 0.627 , 0.660 , 0.637 , 0.684 , 0.504 , 0.511 , 0.527 , 0.556 , 0.592 , 0.625 ,
0.648 , 0.667 , 0.688 , 0.705 , 4.716 ,0.097 , 0.112 , 0.125 , 0.150 , 0.145 , 0.144 , 0.167 , 0.229 ,
0.251 , 0.256 , 0.254 , 6.245 , 0.298 , 0.332 , 0.339 , 0.319 , 0.423 , 0.387 , 0.339 , 0.341 , 0.318 ,
0.333 , 0.396 , 8.385 , 0.427 , 0.519 , 0.598 , 0.638 , 0.665 , 0.648 , 0.598 ,0.311 , 0.282 , 0.255 ,
0.240 , 9.228 , 0.229 , 0.228 , 0.222 , 0.207 , 0.182 , 0.157 , 0.157 , 0.157 , 0.145 , 0.129 , 0.123 ,
0.108 , 0.103 , 5.104 , 0.110 , 0.117 , 0.128 , 0.135 , 0.149 , 0.157 , 0.160 , 0.165 , 0.173 , 0.183 ,
0.206 , 0.204 ,10.521 , 0.533 , 0.554 , 0.592 , 0.630 , 0.640 , 0.629 , 0.639 , 0.643 , 0.639 , 0.682 ,
0.744 , 20.717 , 20.697 , 0.683 , 0.728 , 0.780 , 0.784 , 0.793 , 0.793 , 0.774 , 0.743 , 0.737 , 0.758 ,
0.759 , 40.790 , 0.825 , 0.835 , 0.861 , 0.888 , 0.859  ),
       .Dim=c(4,31)))

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：model dnorm mode Norm For 模型

Hello

You should correct the model following the error message from BUGS. There is a website summarizing possible solutions. Good luck!
http://statisticalestimation.blogspot.com.au/2014/06/common-error-messages-in-bugs-and-their.html

scure 发表于 2015-1-15 12:36
Hello

You should correct the model following the error message from BUGS. There is a website summ ...

thank you！