Dear Dr Skaug and R users,
just discovered glmm.admb in R, and it seems a very useful tool.
However, I ran into a problem when I compare two models:
m1<-glmm.admb(survival~light*species*damage, random=~1, group="table",
data=bm, family="binomial", link="logit")
m1.1<-glmm.admb(survival~(light+species+damage)^2, random=~1,
group="table", data=bm, family="binomial", link="logit")
anova(m1, m1.1)
I get the following output with the warning
Analysis of Variance Table
Model 1: survival ~ light * species * damage
Model 2: survival ~ (light + damage + species)^2
NoPar LogLik Df -2logQ P.value
1 9.000 -103.307
2 9.000 -103.781 0 -0.948
Warning message:
NaNs produced in: pchisq(q, df, lower.tail, log.p)