[精彩WinBUGS答问]Variables that do not converge in winbugs

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

I conducted a Bayesian analysis in winbugs and then checked the convergence from the history plots. The regression coefficients look stabilized, but the variance parameters don't. I got the number of iteration up to 5000, but it still does not look stabilized. I attached the history plots of a regression coefficient and variances. Can you please tell me what I might be doing wrong? Or is is converged?

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Variables Variable converge winbugs WINBUG conducted history checked please number

The second and third picture don't seem to converge. You can look at Rhat value. If it is 1.0, then it is converged. If it is far away from that, say 1.5, probably it doesn't converge.

In papers it is usual to use like 200,000. Variance parameters are always the worse for estimation. You have a huge autocorrelation in the chains - you should increase your thin parameter at least 16 times. What is your thin parameter now?

Anyway, I think the variables would converge quite well if you let the chains run longer. So set the thin and do at least 50,000 iterations.

The model is likely to take hours to run, but even with that I suspect that phiA and tauA chains won't converge. They seem to exhibit poor mixing and autocorrelation. Sometimes long burn-ins and lots of thinning can solve this, but often doesn't in my experience. Centering or scaling to N(0, 1) can help. I'm actually looking for advice on what else to do to improve mixing and convergence in hierarchical Bayesian models, so someone else might be of more help.