[R] Bug in lmer?

[Renaud Lancelot]

I got a crash too with my Win XP config running:

platform i386-pc-mingw32
arch     i386
os       mingw32
system   i386, mingw32
status   beta
major    2
minor    2.0
year     2005
month    09
day      24
svn rev  35666
language R

Package:       Matrix
Version:       0.98-7
Date:          2005-09-09

Package:       lme4
Version:       0.98-1
Date:          2005-07-27

The monitoring of the crash with Dr.MingW provided the following information:

Rgui.exe caused an Access Violation at location 1001ed89 in module
R.dll Writing to location 00000008.

Registers:
eax=00000000 ebx=00000040 ecx=011c8a18 edx=011c8058 esi=011c8a40 edi=00000018
eip=1001ed89 esp=0022dc20 ebp=0022dc48 iopl=0         nv up ei pl nz ac pe nc
cs=001b  ss=0023  ds=0023  es=0023  fs=003b  gs=0000             efl=00000212

Call stack:
1001ED89  R.dll:1001ED89  malloc_usable_size
1001FC2A  R.dll:1001FC2A  Rm_malloc
1009EF69  R.dll:1009EF69  Rf_allocVector
100389D8  R.dll:100389D8  do_makevector
100A6141  R.dll:100A6141  do_internal
10081FBF  R.dll:10081FBF  Rf_eval
10084C2B  R.dll:10084C2B  Rf_applyClosure
10081D50  R.dll:10081D50  Rf_eval
10084C2B  R.dll:10084C2B  Rf_applyClosure
10081D50  R.dll:10081D50  Rf_eval
10083985  R.dll:10083985  do_set
10081FBF  R.dll:10081FBF  Rf_eval
10083A55  R.dll:10083A55  do_begin
10081FBF  R.dll:10081FBF  Rf_eval
10084C2B  R.dll:10084C2B  Rf_applyClosure
10081D50  R.dll:10081D50  Rf_eval
10081FBF  R.dll:10081FBF  Rf_eval
10083985  R.dll:10083985  do_set
10081FBF  R.dll:10081FBF  Rf_eval
10083A55  R.dll:10083A55  do_begin
10081FBF  R.dll:10081FBF  Rf_eval
1008457D  R.dll:1008457D  do_if
1008489B  R.dll:1008489B  R_execMethod
016A2907  methods.dll:016A2907  R_standardGeneric
100A9254  R.dll:100A9254  do_standardGeneric
10082094  R.dll:10082094  Rf_eval
10084C2B  R.dll:10084C2B  Rf_applyClosure
10081D50  R.dll:10081D50  Rf_eval
10083985  R.dll:10083985  do_set
10081FBF  R.dll:10081FBF  Rf_eval
10083716  R.dll:10083716  Rf_DispatchOrEval
10081FBF  R.dll:10081FBF  Rf_eval
10083A55  R.dll:10083A55  do_begin
10081FBF  R.dll:10081FBF  Rf_eval
10081FBF  R.dll:10081FBF  Rf_eval
10083A55  R.dll:10083A55  do_begin
10081FBF  R.dll:10081FBF  Rf_eval
1008457D  R.dll:1008457D  do_if
1008489B  R.dll:1008489B  R_execMethod
016A2907  methods.dll:016A2907  R_standardGeneric
100A9254  R.dll:100A9254  do_standardGeneric
10082094  R.dll:10082094  Rf_eval
10084C2B  R.dll:10084C2B  Rf_applyClosure
10081D50  R.dll:10081D50  Rf_eval
10083985  R.dll:10083985  do_set
10081FBF  R.dll:10081FBF  Rf_eval
1009AAF5  R.dll:1009AAF5  Rf_ReplIteration
1009AC27  R.dll:1009AC27  Rf_ReplIteration
1009AEB7  R.dll:1009AEB7  run_Rmainloop
004013AA  Rgui.exe:004013AA
00401315  Rgui.exe:00401315
00401568  Rgui.exe:00401568
00401236  Rgui.exe:00401236
00401288  Rgui.exe:00401288
7C816D4F  kernel32.dll:7C816D4F  RegisterWaitForInputIdle

Best,

Renaud

2005/9/30, Martin Maechler <maechler at stat.math.ethz.ch>:
> >>>>> "Mark" == Mark Lyman <mlyman at byu.edu>
> >>>>>     on Thu, 29 Sep 2005 14:44:38 -0600 writes:
>
>     Mark> I am relatively new to R so I am not confident enough in what I am doing
>     Mark> to be certain this is a bug.
>
>
>     Mark> I am running R 2.1.1 on a Windows XP
>     Mark> machine and the lme4 package version 0.98-1.
>
> lme4 nowadays is heavily based on "Matrix" which version are you
> using there?
>
>     Mark> The following code fits the model I want using the
>     Mark> nlme package version 3.1-60.
>
>       < .............. >  {see a script at the end}
>
>     Mark> The problem is that when I try fitting the model using
>     Mark> the lmer function with the following code:
>
>     Mark> lmer(adg~trt+(1|loc)+(1|block:loc)+(1|loc:trt),mltloc)
>
>     Mark> I get this message from Windows and R closes.
>
>     Mark> >> R for Windows GUI front-end has encountered a problem and needs to
>     Mark> >> close.  We are sorry for the inconvenience.
>
> That definitely means there is a bug.
> The question is *where* the bug is:  "lme4", "Matrix", "R", "Windows".
>
> One first thing coming to mind is a mismatch of "lme4" and "Matri
>
>     Mark> This same code works on a Macintosh. So it doesn't
>     Mark> seem that I have made an error in my code.
>
> correct; I can also fit the model nice and quickly on Linux,
> and summary() confirms the same fit {with the "usual problem" of
> different estimates for the degrees of freedoms 'df'}.
>
> So currently the bug only shows on the Windows platform.
> Could it be that you have a mismatching package "Matrix" version
> there, but not on the Mac?
>
>     Mark> Also if anyone of the random effect terms is removed there is
>     Mark> no problem. Is this something that is being looked at?
>
> not yet, AFAIK.
>
>     Mark> Or I have I made a mistake somewhere? I have included
>     Mark> the data that I am using below.
>
> I'm putting the data and an R script up for FTP,
> so that you or others can run this ``from anywhere'' via
>
>  source("ftp://stat.ethz.ch/U/maechler/R/mltloc-ex.R", echo = TRUE)
>
> Maybe this helps diagnosis,
> Martin Maechler
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
>