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

[Andrew Robinson]

Thanks, Martin.  It was vanilla.


> traceback()
4: .Call(mer2_getPars, mer)
3: as.double(start)
2: nlminb(.Call(mer2_getPars, mer), function(x) .Call(mer2_deviance, 
       .Call(mer2_setPars, mer, x), as.integer(0)), lower = ifelse(const, 
       0, -Inf), control = list(trace = cv$msVerbose, iter.max = cv$msMaxIter, 
       rel.tol = abs(0.001/.Call(mer2_deviance, mer, 0))))
1: lmer2(Reaction ~ Days + (Days | Subject), sleepstudy)


(I will add it to my list of New Year's Resolutions!)



Andrew



On Sat, Jan 27, 2007 at 04:20:30PM +0100, Martin Maechler wrote:
> >>>>> "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

-- 
Andrew Robinson  
Department of Mathematics and Statistics            Tel: +61-3-8344-9763
University of Melbourne, VIC 3010 Australia         Fax: +61-3-8344-4599
http://www.ms.unimelb.edu.au/~andrewpr
http://blogs.mbs.edu/fishing-in-the-bay/