Chris Fonnesbeck
2010-Apr-26 04:31 UTC
[R] Unexpected warnings from summary() on mcmc.list objects
I am trying to get summary statistics from WinBUGS/JAGS output in the form of mcmc.list objects, using the summary() function. However, I get odd warning messages: Warning messages: 1: In glm.fit(x = X, y = Y, weights = weights, start = start, etastart = etastart, : algorithm did not converge 2: In glm.fit(x = X, y = Y, weights = weights, start = start, etastart = etastart, : algorithm did not converge I assume this is a bug, since glm's should not be fit when simply trying to summarize the elements of a list. Any idea what is going on here? Incidentally, is there additional documentation on the mcmc.list type, other than in the coda docs? The information there is pretty thin, and lacks illustrative examples. Thanks, Chris
Uwe Ligges
2010-Apr-26 12:37 UTC
[R] Unexpected warnings from summary() on mcmc.list objects
Please specify reproducible code that show the problem. When I do it, ewverything is fine on my data with the most recently released versions of R, R2WinBUGS and coda. But perhaps you have another package loaded with another summary.mcmc.list method? On 26.04.2010 06:31, Chris Fonnesbeck wrote:> I am trying to get summary statistics from WinBUGS/JAGS output in the > form of mcmc.list objects, using the summary() function. However, I > get odd warning messages: > > Warning messages: > 1: In glm.fit(x = X, y = Y, weights = weights, start = start, etastart > = etastart, : > algorithm did not converge > 2: In glm.fit(x = X, y = Y, weights = weights, start = start, etastart > = etastart, : > algorithm did not converge > > I assume this is a bug, since glm's should not be fit when simply > trying to summarize the elements of a list. Any idea what is going on > here? > > Incidentally, is there additional documentation on the mcmc.list type, > other than in the coda docs? The information there is pretty thin, and > lacks illustrative examples.Not that I know. Best, Uwe Ligges> Thanks, > Chris > > ______________________________________________ > R-help at r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code.
Ben Bolker
2010-Apr-26 14:03 UTC
[R] Unexpected warnings from summary() on mcmc.list objects
Chris Fonnesbeck <fonnesbeck <at> gmail.com> writes:> > I am trying to get summary statistics from WinBUGS/JAGS output in the > form of mcmc.list objects, using the summary() function. However, I > get odd warning messages: > > Warning messages: > 1: In glm.fit(x = X, y = Y, weights = weights, start = start, etastart > = etastart, : > algorithm did not converge[snip]> I assume this is a bug, since glm's should not be fit when simply > trying to summarize the elements of a list. Any idea what is going on > here?Are you using JAGS? Are you sure that these warning messages are really from summary() and not left over from the JAGS run? I think I have seen similar messages from JAGS ...