请教winbugs关于杠杆SV的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

在论坛里下的代码如下
model
{                                                               
for (i in 1:n) {
Ymean[i]<-rho/tau*exp(0.5*theta[i])*(theta[i+1]-mu-phi*(theta[i]-mu));
Yisigma2[i] <- 1/(exp(theta[i])*(1-rho*rho));
y[i]~ dnorm(Ymean[i],Yisigma2[i]);   
    }     
isigma2<-(1-phi*phi)*itau2;     ***这里为什么要出现这个呢？                                                   
theta0 ~ dnorm(mu,isigma2);  
thetamean[1] <- mu + phi*(theta0-mu);   
theta[1] ~ dnorm(thetamean[1],itau2);                                 
for (i in 2:(n+1)) {                                                                    
thetamean[i] <- mu + phi*(theta[i-1]-mu);                                   
theta[i]~dnorm(thetamean[i],itau2);}                                 
phi1 ~ dbeta(20,1.5);     
phi <- 2*phi1-1;     
mu ~ dnorm(0,0.04);                                                              
itau2 ~ dgamma(2.5,0.025);                                                      
tau <- sqrt(1/itau2) ;
rho ~ dunif(-1,1)                                                        
}


update时出现cannot sample from slice for node itau2？？？是怎么回事，而且运行反应特别慢

还是为什么要出来isigma2这个呢

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：winbugs WINBUG BUGS bug Win

程序我用wingus跑過了，沒有問題，lz可以去下個最新的winbugs試試

我也出现了和楼主同样的问题  请问楼楼怎么解决的？

ice486250 发表于 2013-3-24 04:31
程序我用wingus跑過了，沒有問題，lz可以去下個最新的winbugs試試

我已经更新到1.4.3了 但是还是不行，请问你是用的这个程序吗  ？itau2可以运行吗？

陈小小uu 发表于 2015-6-6 11:09
我也出现了和楼主同样的问题  请问楼楼怎么解决的？

忘记了