>>>>> "B" == Baan <baanba...@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 <baanba...@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 ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.