【请教】winbugs初始值、DIC，收敛等问题设置问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

rt，刚才贴的程序不对，改了下：
model{
    yisigma2[1]<-1/exp(theta[1])
    ymean[1]<-beta0+d*exp(theta[1])
    y[1]~dnorm(ymean[1],yisigma2[1])
for(t in 2:n){
    yisigma2[t]<-1/exp(theta[t])
    ymean[t]<-beta0+beta1*(y[t-1]-beta0)+d*exp(theta[t])
    y[t]~dnorm(ymean[t],yisigma2[t])
   }
mu~dnorm(0,0.1)
phi~dbeta(20,1.5)
itau2 ~dgamma(2.5,0.025)
beta0~dnorm(0,1.0E-6)
beta1~dnorm(0,1.0E-6)
d~dnorm(0,1.0E-6)
#betal1<- 2*beta1star-1
tau<-sqrt(1/itau2)
theta0~dnorm(mu,itau2)
thmean[1]<-mu+phi*theta0
theta[1]~dnorm(thmean[1],itau2)
for(t in 2:n){
    thmean[t]<-mu+phi*(theta[t-1]-mu)
    theta[t]~dnorm(thmean[t],itau2)
    }
}
1.关于初始值：用winbugs做SV-M模型，总是提示“this chain contains unitialized varibles”。
程序如下：
我设置的初始值是：list(itau2=0.5, mu=-10, phi=0.975, d=6,beta0=0,beta1=0,theta0=0.1)

2.我的DIC值出不来，翻了之前的帖子，看来是程序没有运行掉进了trap，但程序验证说是语法没有问题。所以不知道了。
3.关于收敛性问题，因为我做的是单链，G_R统计量使用多链问题，那单链的收敛问题如何解决呢。

谢谢大家了。

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：winbugs WINBUG BUGS 初始值 dic 如何 chain 程序 模型