winbugs求助，自己建立的模型update不了，内贴程序希望能得到解决！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

运行程序时候提示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)

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：winbugs WINBUG update BUGS date update 程序 模型 winbugs

程序没有问题啊, WinBUGS跑不起来, 用 OpenBUGS. 下图是运转1000的截图.

wwwujl 发表于 2013-6-13 20:47
程序没有问题啊, WinBUGS跑不起来, 用 OpenBUGS. 下图是运转1000的截图.

我下载一个试试……
为什么winbugs提示不行呢？

没有zero[i]的值当然跑不起来

运行到2100次就自动停止了

这个程序应该是可以运行的，可以试试将C变大，换成10000，还有l[i]后应该不需要分号，最后要保证l[i]的表达式没有推错。