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

[mad][em39]

按gen inits 一直说unable to generate initial values for node是什么情况？怎么解决啊？？？

贴上最新的code 求问问题在哪。。。。

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

复制代码

楼主问题解决了么？

sheepyang1992 发表于 2015-5-30 20:45
按gen inits 一直说unable to generate initial values for node是什么情况？怎么解决啊？？？

楼主问题解决了么？能否分先一下经验

sheepyang1992 发表于 2015-5-29 12:32
我再来人工顶个贴。。。论文急着用啊。。。

楼主最后解决了吗 我正好也要写关于SV模型的论文

你需要增加一条链 WINBUGS中有一个 gen init 点击 如果显示 initial values generated：model initialized就ok了 代表预设数值已经一一对应了

为什么mu~dnorm(0,0.01)？按照Kim(1998)的文献不应该是mu服从方差为100的正态分布吗，那应该假设为mu~dnorm(0,100)呀