Ad: Re: R crashes for large formulas in lm() (PR#8180)

And some more informastion I forgot.
R does not crash if I write out the formula:

set.seed(123)
x1 <- runif(1000)
x2 <- runif(1000)
x3 <- runif(1000)
x4 <- runif(1000)
x5 <- runif(1000)
x6 <- runif(1000)
x7 <- runif(1000)
x8 <- runif(1000)
y <- rnorm(1000)
fit <- lm(y~(x1*x2*x3*x4*x5*x6*x7*x8)^2)
-> R crashes

fit <- lm(y~x1+x2+x3+x4+x5+x6+x7+x8
               +x1:x2+x1:x3+x1:x4+x1:x5+x1:x6+x1:x7+x1:x8
               +x2:x3++x2:x4+x2:x5+x2:x6+x2:x7+x2:x8
               +x3:x4+x3:x5+x3:x6+x3:x7+x3:x8
               +x4:x5+x4:x6+x4:x7+x4:x8
               +x5:x6+x5:x7+x5:x8
               +x6:x7+x6:x8
               +x7:x8)
-> R does not crash
This is the same formula, at least it should be.

terms(y~(x1*x2*x3*x4*x5*x6*x7*x8)^2, simplify=TRUE)

Uwe Ligges <ligges at statistik.uni-dortmund.de>
05.10.2005 12:13

       Til:    Prof Brian Ripley <ripley at stats.ox.ac.uk>
       cc:     hallgeir.grinde at elkem.no, R-bugs at biostat.ku.dk
       Emne:   Re: [Rd] R crashes for large formulas in lm() (PR#8180)


Prof Brian Ripley wrote:

On Wed, 5 Oct 2005 hallgeir.grinde at elkem.no wrote:

Full_Name: Hallgeir Grinde
Version: 2.1.1
OS: Windows XP
Submission from: (NULL) (144.127.1.1)


While using lm(y~(x*z*c*...*v)^2) R crashes/closes if the numbers of

variables

are at least 8.

OK, let's try to reproduce that:

x1 <- runif(1000)
x2 <- runif(1000)
x3 <- runif(1000)
x4 <- runif(1000)
x5 <- runif(1000)
x6 <- runif(1000)
x7 <- runif(1000)
x8 <- runif(1000)
y <- rnorm(1000)
fit <- lm(y~(x1*x2*x3*x4*x5*x6*x7*x8)^2)

No crash, a quite reasonable fit.

Can we please have a reproducible example, as we do ask?

Hmm, crashes for me as well with R-2.1.1 and R-2.2.0 beta (2005-09-27
r35682M) on WinNT 4.0, SP6.


Let's make it reproducible:

set.seed(123)
x1 <- runif(1000)
x2 <- runif(1000)
x3 <- runif(1000)
x4 <- runif(1000)
x5 <- runif(1000)
x6 <- runif(1000)
x7 <- runif(1000)
x8 <- runif(1000)
y <- rnorm(1000)
fit <- lm(y~(x1*x2*x3*x4*x5*x6*x7*x8)^2)


Uwe Ligges



                               NOTICE
              Please immediately e-mail back to sender
              if you are not the intended recipient.

              Thereafter delete the e-mail along with
              any attachments without making copies.

              Elkem reserves all rights of privilege,
              confidentiality and copyright.