Usando a função de summary.lm em rapache

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

  •  05-07-2019
  •  | 
  •  

Pergunta

Eu instalei rapache e eu estou tentando ajustar um modelo linear dentro do R Arquivo de script. Eu tenha configurado o RFileHandler no http.conf. Quando eu estou tentando invocar o resumo (modelo) ele está me dando um erro de falha de segmento (i ver isso no arquivo de log apache). Eu estou supondo que ele está tentando imprimir para o console e é por isso que está a falhar.

Já alguém encontrou um problema semelhante com R e rapache? Eu sou relativamente novo para R e resumo está fazendo um monte de coisas que não estão expostos diretamente como funções, por isso estou esperando que eu poderia obtê-lo para o trabalho

Aqui está o meu script r 

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

Se eu substituir 

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

com este 

    cat(coef(fit))

que está trabalhando!

Graças Bharani

Foi útil?

Solução 4

Eu finalmente descobri o problema. Lendo o discussão i erroneamente libRlapck.so para lapack.so. Looks como que estava causando problemas. Será que uma instalação limpa do R novamente e então modificado apache para explicity carregar as bibliotecas, em seguida, tudo funcionou obrigado - Bharani

Outras dicas

Será que você considerar contactar o rapache Google Grupo como o home page rapache sugere? Você pode achar leitor mais experiente lá do que aqui.

Eu testei o exemplo e cat(summary(fit)$adj.r.squared) seguinte obras no meu (default) de configuração (última rapache 1.1.8 e R 2.9.2 sob 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

Só descobri que não é com rapache. Ele está falhando em si 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

não sei o que isso significa embora

-Bharani

Licenciado em: CC-BY-SA com atribuição
Não afiliado a StackOverflow
scroll top