robustbase adjbox segfault - memory not mapped

"B" == Baan  <baanbapat at gmail.com>
     on Mon, 4 Mar 2013 22:47:10 +0530 writes:

     B> Thank you Martin. Look forward to the fix.

Committed to the R-forge version of robustbase.

It was a simple integer overflow, indeed,
"necessarily" happening when the sample size was >= 2^16.5.

I'm planning to submit  robustbase_0.9-7  to CRAN today.
Martin

     B> Regards
     B> Baan


     B> On Monday 04 March 2013 10:19 PM, Martin Maechler wrote:

     >>>>>>> "B" == Baan  <baanbapat at gmail.com>
     >>>>>>> on Mon, 4 Mar 2013 15:02:02 +0530 writes:

     B> Hi, I encountered a segfault, memory not mapped error
     B> when using adjbox in robustbase. In trying to recreate
     B> the issue I found that the error occurs only for large
     B> sample size. Here is the code.

     >>

     >> >> require(robustbase)

     B> Loading required package: robustbase

     >> >> x <- rnorm(10)
     >> >> y <- rep(1, 10)
     >> >> adjbox(x ~ y) ## gives a plot
     >> >> x <- rnorm(10000)
     >> >> y <- rep(1, 10000)
     >> >> adjbox(x ~ y) ## gives a plot
     >> >> x <- rnorm(100000)
     >> >> y <- rep(1, 100000)
     >> >> adjbox(x ~ y)

     >>

     B> *** caught segfault ***
     B> address 0xfffffffcc47af530, cause 'memory not mapped'

     >>
     >>

     B> Traceback:
     B> 1: .C(mc_C, x, n, eps = eps, iter = c.iter, medc = double(1))
     B> 2: mcComp(x, doReflect, eps1 = eps1, eps2 = eps2, maxit = maxit,
     B> trace.lev = trace.lev)
     B> 3: mc.default(x, ..., na.rm = TRUE)
     B> 4: mc(x, ..., na.rm = TRUE)
     B> 5: adjboxStats(unclass(groups[[i]]), coef = range, doReflect = doReflect)
     B> 6: adjbox.default(split(mf[[response]], mf[-response]), ...)
     B> 7: adjbox(split(mf[[response]], mf[-response]), ...)
     B> 8: adjbox.formula(x ~ y)
     B> 9: adjbox(x ~ y)

     >>
     >> Indeed, I (as maintainer of robustbase) can reproduce the
     >> segfault.... *even* though you did not specify the random seed...
     >>
     >> So this should be fixed ... hopefully within a week or so,
     >> but I am not promising anything, given my busy schedule!
     >>
     >> Martin Maechler,
     >> ETH Zurich
     >>
     >> [................]
     >>

     B> My setup details:

     >>

     B> R --version
     B> R version 2.15.2 (2012-10-26) -- "Trick or Treat"

     >>

     B> Package:            robustbase
     B> Version:            0.9-5
     B> Date:               2012-03-01
     B> Packaged:           2013-03-01 16:34:03 UTC; maechler
     B> NeedsCompilation:   yes
     B> Repository:         CRAN
     B> Date/Publication:   2013-03-01 18:31:33
     B> Built:              R 2.15.2; x86_64-pc-linux-gnu; 2013-03-04 05:54:20
     B> UTC; unix

     >>
     >>

     B> Platform: x86_64-pc-linux-gnu (64-bit)
     B> uname -a
     B> Linux R 2.6.32-5-amd64 #1 SMP Mon Feb 25 00:26:11 UTC 2013 x86_64 GNU/Linux
     B> Debian squeeze

     >>

     B> Could someone pls help.

     >>

     B> Regards
     B> Baan