ansaurus

Question

Answer 1

+1 A:

Did you consider contacting the rapache Google Group as the rapache home page suggests? You may find more experienced reader there than here.

Dirk Eddelbuettel 2009-10-22 14:03:38

Thanks - just posted in the rapache group

Bharani 2009-10-22 14:12:00

Answer 2

+1 A:

I have tested the following example and cat(summary(fit)$adj.r.squared) works in my (default) setup (latest rapache 1.1.8 and R 2.9.2 under Ubuntu 9.04)

ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14)
trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69)
group <- gl(2,10,20, labels=c("Ctl","Trt"))
weight <- c(ctl, trt)
fit <- lm(weight ~ group - 1) # omitting intercept
setContentType("text/html")
cat('<HTML><BODY>')
cat(summary(fit)$adj.r.squared)
cat('</BODY></HTML>\n')
DONE

rcs 2009-10-22 14:41:23

Did you have to configure R with lapack? bcos with out that it failed the first time and then i did configure with lapack and it started working - just that it is failing when i call summary

Bharani 2009-10-23 04:53:00

I have the same version of both R and rapache but on RHE 5.1

Bharani 2009-10-23 05:01:54

I didn't compile R, just rapache. I'm using the R Ubuntu packages from CRAN (cran.r-project.org/bin/linux/…). BTW, there are also Red Hat packages.

rcs 2009-10-23 06:43:35

Answer 3

A:

Just found out that it is not with rapache. It is failing in R itself

 *** caught segfault ***
 address (nil), cause 'memory not mapped'

 Traceback:
  1: .Call("La_chol2inv", x, size, PACKAGE = "base")
  2: chol2inv(Qr$qr[p1, p1, drop = FALSE])
  3: summary.lm(fit)
  4: summary(fit)
  5: cat(summary(fit)$adj.r.squared)

 Possible actions:
 1: abort (with core dump, if enabled)
 2: normal R exit
 3: exit R without saving workspace
 4: exit R saving workspace

not sure what that means though

-Bharani

Bharani 2009-10-23 05:12:38

Hmm, curious. The "memory not mapped" error means that your code is reading from or writing to memory that it doesn't own.

rcs 2009-10-23 06:55:20

Answer 4

A:

I finally figured out the problem. Reading the discussion i wrongly libRlapck.so to lapack.so. Looks like that was causing problems. Did a clean install of R again and then modified apache to explicity load the libraries then it all worked Thanks - Bharani

Bharani 2009-10-23 09:32:27

ansaurus

tags:

views:

answers:

Using summary.lm function in rapache

related questions