[R-sig-ME] lme4 broken ?

lancelot renaud.lancelot at cirad.fr
Fri Apr 13 08:55:24 CEST 2012


I re-installed the previous version of Matrix (Matrix_1.0-5.tar.gz 
available on CRAN) from source, and it now works (see below).

I confirm there is a problem - at least on Windows - between lme4 (CRAN 
and R-Forge present versions) and the present version of Matrix 1.0-6.

All the best,

Renaud

 > sessionInfo()
R version 2.15.0 (2012-03-30)
Platform: i386-pc-mingw32/i386 (32-bit)

locale:
[1] LC_COLLATE=French_France.1252  LC_CTYPE=French_France.1252 
LC_MONETARY=French_France.1252 LC_NUMERIC=C 
LC_TIME=French_France.1252

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base

other attached packages:
[1] lme4_0.999902344-0 Matrix_1.0-5       lattice_0.20-6

loaded via a namespace (and not attached):
[1] grid_2.15.0    MASS_7.3-17    minqa_1.2.0    nlme_3.1-103 
splines_2.15.0 tools_2.15.0


 > library(lme4)
Le chargement a nécessité le package : lattice
Le chargement a nécessité le package : Matrix
 > example(lmer)

lmer> ## linear mixed models - reference values from older code
lmer> (fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy))
Linear mixed model fit by REML ['lmerMod']
Formula: Reaction ~ Days + (Days | Subject)
    Data: sleepstudy

REML criterion at convergence: 1743.628

Random effects:
  Groups   Name        Variance Std.Dev. Corr
  Subject  (Intercept) 612.10   24.741
           Days         35.07    5.922   0.066
  Residual             654.94   25.592
Number of obs: 180, groups: Subject, 18

Fixed effects:
             Estimate Std. Error t value
(Intercept)  251.405      6.825   36.84
Days          10.467      1.546    6.77

Correlation of Fixed Effects:
      (Intr)
Days -0.138

lmer> (fm2 <- lmer(Reaction ~ Days + (1|Subject) + (0+Days|Subject), 
sleepstudy))
Linear mixed model fit by REML ['lmerMod']
Formula: Reaction ~ Days + (1 | Subject) + (0 + Days | Subject)
    Data: sleepstudy

REML criterion at convergence: 1743.669

Random effects:
  Groups   Name        Variance Std.Dev.
  Subject  (Intercept) 627.58   25.052
  Subject  Days         35.86    5.988
  Residual             653.58   25.565
Number of obs: 180, groups: Subject, 18

Fixed effects:
             Estimate Std. Error t value
(Intercept)  251.405      6.885   36.51
Days          10.467      1.560    6.71

Correlation of Fixed Effects:
      (Intr)
Days -0.184

lmer> anova(fm1, fm2)
Data: sleepstudy
Models:
fm2: Reaction ~ Days + (1 | Subject) + (0 + Days | Subject)
fm1: Reaction ~ Days + (Days | Subject)
     Df    AIC    BIC  logLik deviance  Chisq Chi Df Pr(>Chisq)
fm2  5 1762.0 1778.0 -876.00   1752.0
fm1  6 1763.9 1783.1 -875.97   1751.9 0.0639      1     0.8004



Le 12/04/2012 20:11, Joehanes, Roby (NIH/NHLBI) [F] a écrit :
> Hi Renaud,
>
> It works okay in my computer (Linux). Maybe you want to update all underlying libraries?
>
> Hope this helps,
> Roby
>
>   (fm1<- lmer(Reaction ~ Days + (Days|Subject), sleepstudy))
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>

-- 
Renaud Lancelot
EDENext Project, coordinator: http://www.edenext.eu/

CIRAD, UMR15, Campus International de Baillarguet TA A-DIR / B
F34398 Montpellier

Tel.  +33 4 67 59 37 17  -  Fax  +33 4 67 59 37 98
Secr. +33 4 67 59 37 37  - Cell. +33 6 77 52 08 69



More information about the R-sig-mixed-models mailing list