R help - May 2005 - memory error message using MASS and GLMMGibbs

Hello,

I was just testing the MASS code examples for chapter 10 (Random and Mixed
Effects) and I have pasted the following code in an R session (2.1.0 in
windows 2000 professional; I have also Xemacs + ESS installed, but I was not
using them at that time; my machine has quite a lot of RAM):

library(MASS)
library(lattice)
library(nlme)
library(GLMMGibbs)
# declare a random intercept for each subject
epil$subject <- Ra(data = factor(epil$subject))
glmm(y ~ lbase*trt + lage + V4 + subject, family = poisson,
     data = epil, keep = 100000, thin = 100)

and then an Application Error appears: 

"The instruction at "0x1001edc9" referenced memory at
"0x00000008". The
memory could not be "written". "

It does not take long for this message to appear (less than 1s after I type
Enter).

Any help is welcome

Jordi




--------------------------------------------------------------------------------
The information contained herein is confidential and is inte...{{dropped}}

This works for me (after quite a while) using the pre-compiled binary from 
CRAN on Windows XP.

It is random, so try a different seed.

If you consult the rw-FAQ you will see how to get a more usable error 
message.

On Tue, 3 May 2005, Molins, Jordi wrote:
>
> Hello,
>
> I was just testing the MASS code examples for chapter 10 (Random and Mixed
> Effects) and I have pasted the following code in an R session (2.1.0 in
> windows 2000 professional; I have also Xemacs + ESS installed, but I was
not
> using them at that time; my machine has quite a lot of RAM):
>
> library(MASS)
> library(lattice)
> library(nlme)
> library(GLMMGibbs)
> # declare a random intercept for each subject
> epil$subject <- Ra(data = factor(epil$subject))
> glmm(y ~ lbase*trt + lage + V4 + subject, family = poisson,
>     data = epil, keep = 100000, thin = 100)
>
> and then an Application Error appears:
>
> "The instruction at "0x1001edc9" referenced memory at
"0x00000008". The
> memory could not be "written". "
>
> It does not take long for this message to appear (less than 1s after I type
> Enter).
>
> Any help is welcome
>
> Jordi
>
>
>
>
>
--------------------------------------------------------------------------------
> The information contained herein is confidential and is inte...{{dropped}}
>
> ______________________________________________
> 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
>
-- 
Brian D. Ripley,                  ripley at stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595