全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件
3499 5
2013-06-12
运行程序时候提示undefined real result,修改过初值了,但是还是跑步起来,请问如何解决?


model{
c <- 1000
for (i in 1:N){
zeros[i] <- 0
zeros[i] ~ dpois(zeros.mean[i])
zeros.mean[i] <- -l[i]+c
l[i] <- log(lambda*beta*gama*pow((x[i]/alpha),(beta-1))*exp(pow((x[i]/alpha),beta))*exp(lambda*alpha*(1-exp(pow((x[i]/alpha),beta))))*pow((1-exp(lambda*alpha*(1-exp(pow((x[i]/alpha),beta))))),(gama-1)));
}
# prior
alpha ~ dgamma(400,10);
beta ~ dunif(0,10);
lambda ~dunif(0,1);
gama ~dunif(0,1);
}


list(x=c(0.0006,0.1281,0.1487,0.2499,0.5107,1.2060,2.9581,3.5664,4.3006, 5.5184),N=10)

list(alpha = 40, beta = 1,lambda= 0.5,gama=0.5)

二维码

扫码加我 拉你入群

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

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

全部回复
2013-6-13 20:47:11
程序没有问题啊, WinBUGS跑不起来, 用 OpenBUGS. 下图是运转1000的截图.
Untitled.png
二维码

扫码加我 拉你入群

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

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

2013-6-14 06:50:21
wwwujl 发表于 2013-6-13 20:47
程序没有问题啊, WinBUGS跑不起来, 用 OpenBUGS. 下图是运转1000的截图.
我下载一个试试……
为什么winbugs提示不行呢?
二维码

扫码加我 拉你入群

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

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

2013-11-18 16:36:16
没有zero[i]的值当然跑不起来
二维码

扫码加我 拉你入群

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

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

2015-4-12 18:13:58
运行到2100次就自动停止了
二维码

扫码加我 拉你入群

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

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

2015-4-23 00:22:12
这个程序应该是可以运行的,可以试试将C变大,换成10000,还有l[i]后应该不需要分号,最后要保证l[i]的表达式没有推错。
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群