全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件
6283 14
2013-04-28
求助,求助,哪位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语句要怎么写才好?请求好心高手指点呀?



二维码

扫码加我 拉你入群

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

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

全部回复
2013-4-29 12:26:15
还是么的人回答呀 现在的问题是我可以做出来结果 但是点DIC的时候显示 不能够计算DIC的值 是怎么回事?是我初始值错误码?还是迭代次数不够呢?
二维码

扫码加我 拉你入群

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

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

2013-5-10 11:42:02
梦苏^.^ 发表于 2013-4-29 12:26
还是么的人回答呀 现在的问题是我可以做出来结果 但是点DIC的时候显示 不能够计算DIC的值 是怎么回事?是我 ...
您好!那个初值设定怎么解决的呢?我也遇到了这样的问题,郁闷好多天了没解决
二维码

扫码加我 拉你入群

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

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

2013-5-14 15:05:44
601361436 发表于 2013-5-10 11:42
您好!那个初值设定怎么解决的呢?我也遇到了这样的问题,郁闷好多天了没解决
目前还没有解决,同郁闷中
二维码

扫码加我 拉你入群

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

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

2013-6-3 15:08:35
我也碰过类似的问题,迭代次数不够,迭代4000次左右,应该就可以了
二维码

扫码加我 拉你入群

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

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

2013-6-3 15:09:31
我也碰过类似的问题,迭代次数不够,迭代4000次左右,应该就可以了
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

分享

扫码加好友,拉您进群