[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