Usando la función summary.lm en rapache

StackOverflow https://stackoverflow.com/questions/1607413

  •  05-07-2019
  •  | 
  •  

Pregunta

He instalado rapache y estoy tratando de ajustar un modelo lineal dentro de la R archivo de comandos. He configurado el RFileHandler en http.conf. Cuando estoy intentando invocar el resumen (modelo), me está dando un error de falla de segmento (veo esto en el archivo de registro de apache). Supongo que está intentando imprimir en la consola y es por eso que está fallando.

¿Alguien ha encontrado un problema similar con R y rapache ? Soy relativamente nuevo en R y el resumen es hacer muchas cosas que no están directamente expuestas como funciones, por lo que espero poder hacerlo funcionar

Aquí está mi script rp 

mydata <- read.table("/home/user/test.csv", header = TRUE, sep = ",")
fit <- lm(y~x1+x2+x3, data = mydata)
setContentType("text/html")
cat('<HTML><BODY>')
cat(summary(fit)$adj.r.squared)
cat('</BODY></HTML>\n')
DONE

si sustituyo 

    cat(summary(fit)$adj.r.squared)

con esto 

    cat(coef(fit))

está funcionando!

Gracias Bharani

¿Fue útil?

Solución 4

Finalmente me di cuenta del problema. Leyendo la discusión i incorrectamente libRlapck.so a lapack.so. Parece que eso estaba causando problemas. Hice de nuevo una instalación limpia de R y luego modifiqué apache para explicitar cargar las bibliotecas entonces todo funcionó Gracias  - Bharani

Otros consejos

¿Consideró ponerse en contacto con el rapache Google Group como la página de inicio de rapache sugiere? Puede encontrar lectores más experimentados allí que aquí.

He probado el siguiente ejemplo y cat (resumen (ajuste) $ adj.r.squared)  funciona en mi configuración (predeterminada) (las últimas Rapache 1.1.8 y R 2.9.2 bajo 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

Acabo de descubrir que no es con rapache. Está fallando en la propia R 

 *** 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

no estoy seguro de qué significa eso

-Bharani

Licenciado bajo: CC-BY-SA con atribución
No afiliado a StackOverflow
scroll top