winbugs 做SV-N的 初始值设定的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

求助，求助，哪位winbugs的高人可以给予指导？做SV-N模型，模型代码什么的如下：
model{
for(i in 1:n)
{p<-1/exp(theta)
y~dnorm(0,p)
}
mu~dnorm(0,0.01)
phi1~dbeta(20,1.5)
itau2~dgamma(2.5,0.025)
phi<-2*phi1-1
tau<-sqrt(1/itau2)
theta0~dnorm(mu,itau2)
thmean[1]<-mu+phi*(theta0-mu)
theta[1]~dnorm(thmean[1],itau2)
for(j in 2:n)
{thmean[j]<-mu+phi*(theta[j-1]-mu)
theta[j]~dnorm(thmean[j],itau2)
}
}
list(n=20,y=c(-0.009644843,0.003862035,0.009047287,0.052496355,0.028357443,-0.016457656,0,0.000985286,0.03131365,0,
-0.024182912,0.002311462,-0.005053068,-0.011504959,0.03262237,0.06839109,-0.079104272,-0.00740191,0.001627922,0.011514828)
)
list(
mu=0,phi1=0.975,itau2=0.02)

现在关键就在这个初始值的设定，如果我这个初始值的list设定去掉 直接用随机生成初始值的话，我是可以做出结果的，但是如果我如上写的话（因为我只想要设置这几个参数的初始值，其他的参数初始值都随机）我点击load init按钮出现的是this chain contains uninitialized variables ，我再点击gen inits ，然后继续操作，迭代后结果在sample monitor tool 对话框中density，stats等选项是灰的，只有trace 和history 选项可以查看，这是这么一回事？那么初始值应该怎么设定呢？这种情况下LIST语句要怎么写才好？请求好心高手指点呀？

扫码加我 拉你入群

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

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

分享0 收藏4 回帖

关键词：winbugs WINBUG BUGS bug Win 模型

还是么的人回答呀 现在的问题是我可以做出来结果 但是点DIC的时候显示 不能够计算DIC的值 是怎么回事？是我初始值错误码？还是迭代次数不够呢？

梦苏＾.＾ 发表于 2013-4-29 12:26
还是么的人回答呀 现在的问题是我可以做出来结果 但是点DIC的时候显示 不能够计算DIC的值 是怎么回事？是我 ...

您好！那个初值设定怎么解决的呢？我也遇到了这样的问题，郁闷好多天了没解决

601361436 发表于 2013-5-10 11:42
您好！那个初值设定怎么解决的呢？我也遇到了这样的问题，郁闷好多天了没解决

目前还没有解决，同郁闷中

我也碰过类似的问题，迭代次数不够，迭代4000次左右，应该就可以了

我也碰过类似的问题，迭代次数不够，迭代4000次左右，应该就可以了

yoyoyowt 发表于 2013-6-3 15:09
我也碰过类似的问题，迭代次数不够，迭代4000次左右，应该就可以了

谢谢

我是刚刚用winbugs的新手，请教下下面这几句是什么意思，能不能注释下
tau<-sqrt(1/itau2)
theta0~dnorm(mu,itau2)
thmean[1]<-mu+phi*(theta0-mu)
theta[1]~dnorm(thmean[1],itau2)
for(j in 2:n)
{thmean[j]<-mu+phi*(theta[j-1]-mu)
theta[j]~dnorm(thmean[j],itau2)

您好！设置这几个参数的初始值，点击load init按钮出现的是this chain contains uninitialized variables ，这个问题你解决了吗？我也遇到了，怎么办呢？好急啊！

Jczoe 发表于 2015-3-6 22:59
您好！设置这几个参数的初始值，点击load init按钮出现的是this chain contains uninitialized variables ， ...

你再按gen init就可以了如果你是一条chain的话~这是因为你模型中含有通过初定值推导出的参数。但是如果你是多条chain，我忘了，你百度一下，之前看到有一个博客有写。