From murdoch at stats.uwo.ca  Mon Jun 26 15:14:29 2006
From: murdoch at stats.uwo.ca (Duncan Murdoch)
Date: Mon, 26 Jun 2006 09:14:29 -0400
Subject: [Rd] a new way to crash R? (PR#8981)
In-Reply-To: <20060614204445.3908C3F033@slim.kubism.ku.dk>
References: <20060614204445.3908C3F033@slim.kubism.ku.dk>
Message-ID: <449FDDB5.6060103@stats.uwo.ca>

On 6/14/2006 4:44 PM, jritter at umn.edu wrote:
> This is a multi-part message in MIME format.
> --------------010500030900050501010804
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
> Content-Transfer-Encoding: 7bit
> 
> Dear R Team,
> 
> First, thank you for incredibly useful software!
> 
> Now the bad news:  The attached script (or the original version with 
> real data) will reliably crash R on my machine.  I am using:
> 
> R version: either 2.2.1 or 2.3.1
> Windows 2000 Professional, Service Pack 4
> 512 MB of RAM

I don't have any Win2K machines, but on an XP laptop with only 256M, I 
see the crash, just after failing a couple of times with out-of-memory 
errors.  I'll see if I can track it down.

Duncan Murdoch

> 
> On my machine the script will crash R on line 42 [  probits21 <- 
> lapply(...  ].  
> 
> In both this script and the "real" one, memory.size() indicates that R 
> is about 512 MB just before executing the command that crashes it.  In 
> other contexts, R has been able to allocate virtual memory on this 
> computer without any trouble.
> 
> Thanks again for all your work on R.
> 
> Regards,
> 
> Joe Ritter
> 
> --------------------------------------------------------------------------
> Joseph A. Ritter
> Associate Professor
> Humphrey Institute of Public Affairs
> University of Minnesota
> Minneapolis, MN  55455
> U.S.A.
> 612-625-0442
> 
> --------------010500030900050501010804
> Content-Type: text/plain;
>  name="crash.r"
> Content-Transfer-Encoding: 7bit
> Content-Disposition: inline;
>  filename="crash.r"
> 
> x1 <- rnorm(10000)
> x2 <- rnorm(10000)
> x3 <- rnorm(10000)
> x4 <- rnorm(10000)
> x5 <- rnorm(10000)
> x6 <- rnorm(10000)
> x7 <- rnorm(10000)
> x8 <- rnorm(10000)
> y <- round(runif(10000))
> 
> specs <- list()
> specs$A <- formula(y~x1+x2+x3+x4+x5+x6+x7+x8) 
> specs$B <- specs$A
> specs$C <- specs$A
> specs$D <- specs$A
> specs$E <- specs$A
> specs$F <- specs$A
> specs$H <- specs$A
> specs$I <- specs$A
> specs$J <- specs$A
> 
> probits1 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits2 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits3 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits4 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits5 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits6 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits7 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits8 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits9 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits10 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits11 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits12 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits13 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits14 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits15 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits16 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits17 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits18 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits19 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits20 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits21 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits22 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits23 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits24 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits25 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits26 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits27 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits28 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits29 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits30 <- lapply(specs, glm, family=binomial(link="probit"))      
> probits31 <- lapply(specs, glm, family=binomial(link="probit"))      
> 
> --------------010500030900050501010804--
> 
> ______________________________________________
> R-devel at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel



[Rd] a new way to crash R? (PR#8981)
Duncan Murdoch