Dear R list, A problem with the package glmmADMB was pointed out by R. Bagchi some time ago. We have now fixed the problem and a new version of the package can be downloaded from the usual place:
http://otter-rsch.com/admbre/examples/glmmadmb/glmmADMB.html The likelihood values that previously were off by a constant have now been corrected. Please let me know if there are other problems. hans >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) _____________________________ Hans Julius Skaug Department of Mathematics University of Bergen Johannes Brunsgate 12 5008 Bergen Norway ph. (+47) 55 58 48 61 ______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
