[R-sig-ME] glmer: Downdated X'X is not positive definite

[Ben Bolker]

Ross Boylan <ross at ...> writes:

 [snip]

> Ben Bolker wrote
> > It is indeed a little surprising (to me) that you're getting 
> > downdating problems with an intercept-only model.
> >    Hard to say without a closer look at the data.
> >    I know it won't solve your problem, but does fitting the model
> > _without_ the singleton clusters work?

  I didn't mean without the clusters, I meant dropping the singleton
clusters from the data ...

> glmer complained that there were no random terms, as I recall (not at the
system now).
> glm fit it without complaints.
> 
> SAS GLIMMIX fit a poisson with random intercepts to the data without incident.

  You could try glmmPQL if you like ... (it's not clear whether you
used PQL or Laplace approximation within GLIMMIX ...)

> OTOH, glmer did manage to fit poisson random intercept model to toy
> simulated data I constructed, even though the simulated data had lots of
> singleton "clusters".
> 
> >    Do you have a very small number of values >0?
> No.  N > 2,000 and c 1/3 are 1's.
> 

  [snip]

  Did you try with development lme4?  Does it make a difference?