[R] robustbase adjbox segfault - memory not mapped

[Baan]

Glad to know. Thanks.

Regards
Baan

On Wednesday 06 March 2013 02:15 PM, Martin Maechler wrote:
>>>>>> "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