[R-sig-ME] lmer() with 'na.action=na.exclude'; error with summary()

John Maindonald john.maindonald at anu.edu.au
Mon Sep 8 06:56:49 CEST 2014


The following demonstrates the issue:

> library(DAAG)
> science.lmer <- lmer(like ~ sex + PrivPub + (1 | school) + 
+                      (1 | school:class), data = science,  
+                      na.action=na.exclude) 
> summary(science.lmer)
Linear mixed model fit by REML ['lmerMod']
Formula: like ~ sex + PrivPub + (1 | school) + (1 | school:class)
   Data: science

REML criterion at convergence: 5546.5

Scaled residuals: 
Error in quantile.default(resids) : 
  missing values and NaN's not allowed if 'na.rm' is FALSE
> ## Suppress details of residuals
> summary(science.lmer, show.resids=FALSE)
Linear mixed model fit by REML ['lmerMod']
Formula: like ~ sex + PrivPub + (1 | school) + (1 | school:class)
   Data: science

REML criterion at convergence: 5546.5

Scaled residuals: 
Error in quantile.default(resids) : 
  missing values and NaN's not allowed if 'na.rm' is FALSE
> 

With na.action = na.omit, there is no problem.

> sessionInfo()
R version 3.1.1 Patched (2014-08-22 r66458)
Platform: x86_64-apple-darwin10.8.0 (64-bit)

locale:
[1] en_AU.UTF-8/en_AU.UTF-8/en_AU.UTF-8/C/en_AU.UTF-8/en_AU.UTF-8

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

other attached packages:
[1] lme4_1.1-7      Rcpp_0.11.2     Matrix_1.1-4    DAAG_1.20       lattice_0.20-29

loaded via a namespace (and not attached):
[1] grid_3.1.1          latticeExtra_0.6-26 MASS_7.3-34         minqa_1.2.3         nlme_3.1-117       
[6] nloptr_1.0.4        RColorBrewer_1.0-5  splines_3.1.1       tools_3.1.1    

I doubt that there is an intention for summary.merMod() to throw an error
lmer() has been called with 'na.action=na.exclude’.  It should certainly not 
throw an error with the argument 'show.resids=FALSE’.

John Maindonald             email: john.maindonald at anu.edu.au
phone : +61 2 (6125)3473    fax  : +61 2(6125)5549
Centre for Mathematics & Its Applications, Room 1194,
John Dedman Mathematical Sciences Building (Building 27)
Australian National University, Canberra ACT 0200.



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