SV-T模型OPENBUGS程序不能运行，求助

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

这是我的程序，check model和load data之后load inits还是暗的

1. #model specification
   
2. model{
   
3. for(i in 1:N)
   
4. {y[i]~dt(0,p[i],omega)
   
5. p[i]<-exp(-theta[i])
   
6. }
   
7. theta[1]~dnorm(mu,itau2)
   
8. for(j in 2:N)
   
9. {theta[j]~dnorm(theta2[j],itau2)
   
10. theta2[j]<-mu+phi*(theta[j-1]-mu)}
    
11. phi<-2*phi1-1
    
12. tau<-sqrt(1/itau2)
    
13. mu~dnorm(0,0.01)
    
14. itau2~dgamma(2.5,0.025)
    
15. phi1~dbeta(20,1.5)
    
16. omega~dchisqr(8)
    
17. }
    
18. 
    
19. #Data
    
20. 
    
21. list(N=387,T=1,y = structure(.Data = c(0.000137324,

复制代码

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：openbugs SV-T Open BUGS bug 程序 模型

显示说我的N is not defined

请各位大虾帮个忙！！

.....

我想到个办法，直接把N=387带到model里去就可以了，不过不知道有哪些参数要设初始值

看一下重新修改过的程序，初始值设置还是有点问题

1. for(i in 1:387)
   
2. {y[i]~dt(0,p[i],omega)
   
3. p[i]<-exp(-theta[i])
   
4. }
   
5. theta[1]~dnorm(mu,itau2)
   
6. for(j in 2:387)
   
7. {theta[j]~dnorm(theta2[j],itau2)
   
8. theta2[j]<-mu+phi*(theta[j-1]-mu)}
   
9. phi<-2*phi1-1
   
10. tau<-sqrt(1/itau2)
    
11. mu~dnorm(0,0.01)
    
12. itau2~dgamma(2.5,0.025)
    
13. phi1~dbeta(20,1.5)
    
14. omega~dchisqr(8)
    
15. }
    
16. 
    
17. #Initial values
    
18. list(mu=0.01,itau2=0.1,phi1=2,omega=8)
    
19. 
    
20. #Data
    
21. 
    
22. list(y = c(0.000137324,

复制代码

按了load inits出现initial values loaded but chain contains uninitialized variables

没有人知道吗。。。自己顶一下

设置了初始值list(mu=0,itau2=0.5,theta0=0.97) 但是出现initial values loaded but chain contains uninitialized variables

我再来人工顶个贴。。。论文急着用啊。。。