[R-sig-ME] new version of lme4 - Log-Likelihood of model returned as "-Inf"

juwb08 at hampshire.edu juwb08 at hampshire.edu
Thu Feb 13 00:10:35 CET 2014


Hi All,

I was working with Poisson glmms in lme4 last year in 2013 using  
Windows XP, lme4 version 0.999375-16 (2008-06-23) and R version 2.15.1  
(2012-06-22), evaluating a set of candidate models with AIC to select  
the best ones.

This worked well until I updated lme4 and R last month.

I now updated to R version 3.0.2 (2013-09-25) -- "Frisbee Sailing" and  
lme4 version 1.0-5 (2013-10-24), still Windows XP, and am running the  
same models again on the same data, and lme4 now returns "-Inf" for  
the log likelihood, so it can't calculate AIC.

Using the previous versions, I got stable values for Log-Likelihoods  
and calculated AICs with them.

Attached is a sample of the dataset and the code I am using.

How can I get a Log-Likelihood estimate with the current versions of  
lme4 and R?

Thanks in advance.
Justin Baldwin

Hampshire College, MA, USA
-------------- next part --------------
setwd("C:/Documents and Settings/Justin/Desktop/Research/REU2")

library(lme4)

shits<-read.csv("PretEx-EXAMPLE.csv", header=T)
shits2<-subset(shits, GRT_min!="NA")
shits2$BatID<-as.factor(shits2$BatID)
shits2$PlantID_factor<-as.factor(shits2$PlantID_factor)

summary(glmer(GRT_min~Treatment_Factor+Mass_Consumed+FedByHand+(1|BatID/PlantID_factor), data=shits2, family="poisson", REML=F))



More information about the R-sig-mixed-models mailing list