Tom, thanks for spotting a bug in glmmML; internally, glmmML sorts data by cluster. The bug is that I missed to sort the weights accordingly. Weights are produced when the response is a two-column matrix, as in your case. So, glmmML calls glm.fit with wrong weights, and the warning comes from that. Unfortunately, this also affects the rest of the fitting procedure. I will fix this asap.
Meanwhile, a possible workaround: Sort your data.frame by cluster before calling glmmML. Göran On Wed, Jun 22, 2011 at 5:13 AM, Tom Kraft <thomas.s.kr...@dartmouth.edu> wrote: > Dear all, > > This question is basic but I am stumped. After running the below, I receive > the message: "non-integer #successes in a binomial glm!" > > model1 <- > glmmML(y~Brood.Size*Density+Date.Placed+Species+Placed.Emerging+Year+rate.of.parperplot, > data = data, cluster= data$Patch, family=binomial(link="logit")) > > My response variable is sex ratio, and I have learned quickly not to use > proportion values as data. Instead, my response variable y is a 2 column > matrix that looks like the following but with more rows. Column 1 would be > number of males and column 2 would be number of females in a brood. > > [,1] [,2] > [1,] 18 19 > [2,] 7 37 > [3,] 5 26 > [4,] 4 16 > [5,] 6 19 > [6,] 4 15 > [7,] 15 14 > [8,] 15 29 > > All the numbers in both columns are integers. I am able to use this same > format for my response variable when using the glmer() function with no > problem. In the documentation for glmmML, it says "For the binomial > families, the response can be a two-column matrix". Isn't that what I am > using? Is it immediately obvious to anyone what is wrong with my input? Or > perhaps I am doing something else wrong? Thanks so much in advance, > > Tom > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help@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. > -- Göran Broström ______________________________________________ R-help@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.
