Hi Paul,
We could be of more assistance if you included some source code...
Are you doing the looping in R or in python?
-Greg
On 6/23/08 10:32AM , "Jung, Paul" <[EMAIL PROTECTED]> wrote:
> I'm having the following issue with rpy. I'm calling r.lm to do a bunch
> of regressions (several hundred thousand over a loop) but get a segfault
> as listed below. I think it is a memory issue because the same
> regression will work if I start the loop closer to where the segfault
> occurs. Also, I think the issue concerns data.frames since I don't get
> any problems when I use lsfit (which does the regression using vectors
> instead of a data frame). Any experience with such issues or
> suggestions as to where to look for a solution would be much
> appreciated.
>
> Thanks,
> Paul
>
> RHOME= /usr/local/R-2.6.0/lib/R
> RVERSION= 2.6.0
> RVER= 2060
> RUSER= /home/jung
> Loading Rpy version 2060 .. Done.
> Creating the R object 'r' .. Done
>
> *** caught segfault ***
> address 0xc5ac, cause 'memory not mapped'
>
> Traceback:
> 1: match.call()
> 2: function (formula, data, subset, weights, na.action, method = "qr",
> model = TRUE, x = FALSE, y = FALSE, qr = TRUE, singular.ok = TRUE,
> contrasts = NULL, offset, ...) { ret.x <- x ret.y <- y cl <-
> match.call() mf <- match.call(expand.dots = FALSE) m <-
> match(c("formula", "data", "subset", "weights", "na.action",
> "offset"), names(mf), 0) mf <- mf[c(1, m)] mf$drop.unused.levels
> <- TRUE mf[[1]] <- as.name("model.frame") mf <- eval(mf,
> parent.frame()) if (method == "model.frame") return(mf)
> else if (method != "qr") warning(gettextf("method = '%s' is not
> supported. Using 'qr'", method), domain = NA) mt <-
> attr(mf, "terms") y <- model.response(mf, "numeric") w <-
> as.vector(model.weights(mf)) if (!is.null(w) && !is.numeric(w))
> stop("'weights' must be a numeric vector") offset <-
> as.vector(model.offset(mf)) if (!is.null(offset)) { if
> (length(offset) == 1) offs!
> et <- rep(offset, NROW(y)) else if (length(offset) != NROW(y))
> stop(gettextf("number of offsets is %d, should equal %d (number of
> observations)", length(offset), NROW(y)), domain = NA)
> } if (is.empty.model(mt)) { x <- NULL z <-
> list(coefficients = if (is.matrix(y)) matrix(, 0, 3) else
> numeric(0), residuals = y, fitted.values = 0 * y, weights =
> w, rank = 0, df.residual = if (is.matrix(y)) nrow(y) else length(y))
> if (!is.null(offset)) { z$fitted.values <- offset
> z$residuals <- y - offset } } else { x <-
> model.matrix(mt, mf, contrasts) z <- if (is.null(w))
> lm.fit(x, y, offset = offset, singular.ok = singular.ok,
> ...) else lm.wfit(x, y, w, offset = offset, singular.ok =
> singular.ok, ...) } class(z) <- c(if (is.matrix(y))
> "mlm", "lm") z$na.action <- attr(mf, "na.action") z$offset <-
> offse!
> t z$contrasts <- attr(x, "contrasts") z$xlevels <- .getX!
> levels(m
> t, mf) z$call <- cl z$terms <- mt if (model) z$model <-
> mf if (ret.x) z$x <- x if (ret.y) z$y <- y if
> (!qr) z$qr <- NULL z}(y ~ x + x2 + x3, data = list(x2 = c(0a,
> d7, a3, 70, 3d, 0a, b7, 3f, ec, 51, b8, 1e, 85, eb, c9, 3f, 5d, 8f, c2,
> f5, 28, 5c, d3, 3f, 0a, d7, a3, 70, 3d, 0a, d7, 3f, c6, dc, b5, 84, 7c,
> d0, d7, 3f, 0b, 46, 25, 75, 02, 9a, d8, 3f, d8, 12, f2, 41, cf, 66, d9,
> 3f, 2d, 43, 1c, eb, e2, 36, da, 3f, 0b, d7, a3, 70, 3d, 0a, db, 3f, 71,
> ce, 88, d2, de, e0, db, 3f, 60, 29, cb, 10, c7, ba, dc, 3f, d7, e7, 6a,
> 2b, f6, 97, dd, 3f, d6, 09, 68, 22, 6c, 78, de, 3f), x3 = NULL, y =
> NULL, x = c(33, 33, 33, 33, 33, 33, d3, 3f, cd, cc, cc, cc, cc, cc, dc,
> 3f, 9a, 99, 99, 99, 99, 99, e1, 3f, 33, 33, 33, 33, 33, 33, e3, 3f, 85,
> eb, 51, b8, 1e, 85, e3, 3f, d7, a3, 70, 3d, 0a, d7, e3, 3f, 29, 5c, 8f,
> c2, f5, 28, e4, 3f, 7b, 14, ae, 47, e1, 7a, e4, 3f, cd, cc, cc, cc, cc,
> cc, e4, 3f, 1f, 85, eb, 51, b8, 1e, e5, 3f, 71, 3d, 0a, !
> d7, a3, 70, e5, 3f, c3, f5, 28, 5c, 8f, c2, e5, 3f, 15, ae, 47, e1, 7a,
> 14, e6, 3f)), weights = c(5L, 4L, 15L, 4L, 6L, 5L, 6L, 8L, 14L, 7L, 23L,
> 11L, 17L, 6L, 21L, 8L, 8L, 8L, 8L, 10L, 14L, 13L, 18L, 12L, 38L, 14L,
> 33L, 26L, 37L, 35L, 26L, 28L, 56L, 46L, 27L, 33L, 51L, 22L, 33L, 25L,
> 49L, 28L, 12L, 22L, 22L, 19L, 20L, 16L, 20L, 8L, 15L, 10L, 29L, 10L,
> 16L, 12L, 17L, 32L, 30L, 28L, 22L, 14L, 30L, 19L, 17L, 21L, 15L, 36L,
> 23L, 24L, 37L, 51L, 32L, 17L, 36L, 7L, 34L, 9L, 20L, 4L, 40L, 9L, 34L,
> 31L, 36L, 55L, 50L, 35L, 45L, 35L, 16L, 13L, 14L, 9L, 23L, 5L, 6L, 9L,
> 9L, 25L, 8L, 4L, 12L, 9L))
>
> Possible actions:
> 1: abort (with core dump, if enabled)
> 2: normal R exit
> 3: exit R without saving workspace
> 4: exit R saving workspace
>
> IMPORTANT: The information contained in this email and/or its attachments is
> confidential. If you are not the intended recipient, please notify the sender
> immediately by reply and immediately delete this message and all its
> attachments. Any review, use, reproduction, disclosure or dissemination of
> this message or any attachment by an unintended recipient is strictly
> prohibited. Neither this message nor any attachment is intended as or should
> be construed as an offer, solicitation or recommendation to buy or sell any
> security or other financial instrument. Neither the sender, his or her
> employer nor any of their respective affiliates makes any warranties as to the
> completeness or accuracy of any of the information contained herein or that
> this message or any of its attachments is free of viruses.
>
> -------------------------------------------------------------------------
> Check out the new SourceForge.net Marketplace.
> It's the best place to buy or sell services for
> just about anything Open Source.
> http://sourceforge.net/services/buy/index.php
> _______________________________________________
> rpy-list mailing list
> rpy-list@lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/rpy-list
>
--
Gregory R. Warnes, Ph.D
Program Director
Center for Computational Arts, Sciences, and Engineering
University of Rochester
Tel: 585-273-2794
Fax: 585-276-2097
Email: [EMAIL PROTECTED]
-------------------------------------------------------------------------
Check out the new SourceForge.net Marketplace.
It's the best place to buy or sell services for
just about anything Open Source.
http://sourceforge.net/services/buy/index.php
_______________________________________________
rpy-list mailing list
rpy-list@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/rpy-list
