[本週專題]WinBUGS: Undefined Real Result?

2701

收藏 2014-06-26

Dear lister
I have run a zero-inflated Poisson model in WinBUGS without problems, and now I am trying to run its equivalent negative binomial. However, I get an "undefined real result" trap message over and over. Here is the data set http://www.megafileupload.com/en/file/474809/data1.html. The trap arises from the Pr(Y > 0) part of the model since the Pr(Y = 0) runs without problems. I've tried different inits and priors to no avail.
Thanks for your help!

复制代码

扫码加我拉你入群

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

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

全部回复

Lisrelchen

2014-6-26 21:17:26

Ok, I found a missing zero in one of the priors:

b.descart ~ dnorm(.0,0.001)
Apparently:

BUGS doesn't take .0 as a 0.0, and
WinBUGS doesn't take .0 as a syntax error
Another point: the model neither runs with more than three covariates nor more than two categorical covariates. Any clue?

As an aside, the value I gave to the constant K makes the deviance and DIC to get wildly high values.

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群