"Matt Loveland" <[EMAIL PROTECTED]> writes: > Hi, > > I wrote a few days ago about an error message I'm getting when I use GLMM > from lme4 to do random effects modelling. > When I add random effects, I get the following error message: Error in > "EMsteps<-"(`*tmp*`, value = control) : invalid source matrix. > > (I wanted to note that I've only just started to learn about random effects > modelling and R, so am a novice at both, but R/lme4 have been great up to > this hiccup) > > I was asked to provide a reproducable example, so here is my try, using the > guImmun data frame included with the package: > > library(lme4) > data(guImmun) > > fm1 = GLMM(immun ~ kid2p + mom25p + momWork, > data = guImmun, family = binomial, > random = ~1|mom) > > fm2 = GLMM(immun ~ kid2p + mom25p + momWork, > data = guImmun, family = binomial, > random = ~kid2p|mom) > ** > fm1 produces output, and fm2 gives me the following error message: > > Error in "EMsteps<-"(`*tmp*`, value = control) : > invalid source matrix > > Is my model statement incorrect? Are there possible data problems I should > direct my attention to? > Any help will be appreciated.
> I'm using R 1.8.1 I get the same thing with R 1.9.0. Looks like it might be a bug, but I notice that quite few moms have more than a single observation: > with(guImmun, table(table(mom))) 1 2 3 1063 500 32 and digging a bit deeper, I see that only 287 moms have observations with different values of kid2p. This kind of data pattern has been known to confuse lme-type algorithms before. We probably need to leave further investigation to the package author. -- O__ ---- Peter Dalgaard Blegdamsvej 3 c/ /'_ --- Dept. of Biostatistics 2200 Cph. N (*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918 ~~~~~~~~~~ - ([EMAIL PROTECTED]) FAX: (+45) 35327907 ______________________________________________ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html