[R-sig-ME] Do glmer and glmmadmb calculate log likelihood on thesame scale ?
David Duffy
David.Duffy at qimr.edu.au
Fri Jun 22 04:13:35 CEST 2012
On Thu, 21 Jun 2012, Joshua Wiley wrote:
> Just a note that those are not the same models.
Huh.
> OR you could change glmmadmb() to
>
> fm2 <- glmmadmb(y~Base*trt+Age+Visit+(Visit|subject),
> data=epil2, family="poisson", corStruct = "full")
> however, on my machine, there are warnings about the matrix not being
> positive definite, and the model never converges.
I haven't updated glmmADMB recently (0.6.4 2011-08-30), so it
gives the warning, but continues merrily to give:
(Intercept) Visit
(Intercept) 0.24931 0.24931
Visit 0.54572 0.54572
and the same likelihood as before
Log-likelihood: -655.41
glmer was
Groups Name Variance Std.Dev. Corr
subject (Intercept) 0.24912 0.49912
Visit 0.53962 0.73459 0.011
'log Lik.' -262.6586 (df=9)
Your suggestion Visit + (1 | Visit) + (0 + Visit | subject)
gave me
Groups Name Variance Std.Dev.
subject Visit 0.77301 0.87921
Visit (Intercept) 5.1845e-16 2.2769e-08
'log Lik.' -421.1128 (df=8)
I have to admit, the direction of the change in loglikelihood here
confuses me a little ;)
My lme4 is also a little antiquated: 0.999375-42 (2011-10-02).
--
| David Duffy (MBBS PhD) ,-_|\
| email: davidD at qimr.edu.au ph: INT+61+7+3362-0217 fax: -0101 / *
| Epidemiology Unit, Queensland Institute of Medical Research \_,-._/
| 300 Herston Rd, Brisbane, Queensland 4029, Australia GPG 4D0B994A v
More information about the R-sig-mixed-models
mailing list