[R] glmm.admb - bug and possible solution??

Hans Skaug Hans.Skaug at mi.uib.no
Tue Feb 21 13:31:50 CET 2006


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




More information about the R-help mailing list