[R-sig-ME] New version of lme4 - memory error

[Martin Maechler]

>>>>> "Andrew" == Andrew Robinson <A.Robinson at ms.unimelb.edu.au>
>>>>>     on Sat, 27 Jan 2007 10:27:14 +1100 writes:

    Andrew> Hi Doug, sorry, it's me having problems again :(

    Andrew> I can install and load the new lme4 package with no
    Andrew> trouble, but when I try to run the examples, I get:

    >> require(lme4)
    Andrew> Loading required package: lme4
    Andrew> Loading required package: Matrix
    Andrew> Loading required package: lattice
    Andrew> [1] TRUE
    >> sessionInfo()
    Andrew> R version 2.4.1 Patched (2007-01-25 r40572) 
    Andrew> i386-unknown-freebsd6.1 

    Andrew> locale:
    Andrew> C

    Andrew> attached base packages:
    Andrew> [1] "stats"     "graphics"  "grDevices" "utils"     "datasets"
    Andrew> "methods"  
    Andrew> [7] "base"     

    Andrew> other attached packages:
    Andrew> lme4      Matrix     lattice 
    Andrew> "0.9975-11"  "0.9975-8"   "0.14-16" 
    >> fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy)
    >> fm1 <- lmer2(Reaction ~ Days + (Days|Subject), sleepstudy)
    Andrew> Error in as.double(start) : Calloc could not allocate (169499040 of 4)
    Andrew> memory

Hmm, I can't replicate your problem.  I get

 > system.time( fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy))
 [1] 0.112 0.001 0.124 0.000 0.000
 > system.time( fm1.2 <- lmer2(Reaction ~ Days + (Days|Subject), sleepstudy))
 [1] 0.056 0.000 0.056 0.000 0.000
 > summary(fm1.2)
 Linear mixed-effects model fit by REML 
   AIC  BIC logLik MLdeviance REMLdeviance
  1754 1770 -871.8       1752         1744
 Random effects:
  Groups   Name        Variance Std.Dev. Corr  
  Subject  (Intercept) 612.114  24.7409        
	   Days         35.072   5.9222  0.066 
  Residual             654.937  25.5917        
 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
 > 



    Andrew> Does anyone else find this?  Please let me know what else I can do to
    Andrew> help.

You can type  
        traceback()

after the error (which you should get a habit of doing ;-)
which might be revealing though I doubt it a bit in this case.

Are you sure that you started R as "R --vanilla", i.e.
that you have *not* loaded an .RData of a previous session
accidentally?

Regards,
Martin

    Andrew> Cheers,
    Andrew> Andrew