Utilizzo della funzione summary.lm in rapache

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

  •  05-07-2019
  •  | 
  •  

Domanda

Ho installato rapache e sto cercando di adattare un modello lineare all'interno della R file di script. Ho configurato RFileHandler in http.conf. Quando sto cercando di invocare il riepilogo (modello), mi viene visualizzato un errore di errore del segmento (lo vedo nel file di registro di Apache). Immagino che stia cercando di stampare sulla console ed è per questo che non riesce.

Qualcuno ha riscontrato un problema simile con R e rapache ? Sono relativamente nuovo su R e il sommario sta facendo molte cose che non sono direttamente esposte come funzioni, quindi spero di riuscire a farlo funzionare

Ecco il mio 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 sostituisco 

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

con questo 

    cat(coef(fit))

funziona!

Grazie Bharani

È stato utile?

Soluzione 4

Alla fine ho capito il problema. Lettura della discussione i erroneamente libRlapck.so su lapack.so. Sembra che fosse la causa i problemi. Ha eseguito di nuovo un'installazione pulita di R e quindi ha modificato apache in modo esplicito caricare le librerie quindi tutto ha funzionato Grazie  - Bharani

Altri suggerimenti

Hai considerato di contattare il rapache Google Group come home page rapache suggerisce? Puoi trovare un lettore più esperto lì che qui.

Ho testato il seguente esempio e cat (summary (fit) $ adj.r.squared)  funziona nella mia configurazione (predefinita) (ultime rapache 1.1.8 e R 2.9.2 sotto 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

Ho appena scoperto che non è con rapache. Sta fallendo in R stesso 

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

non sono sicuro di cosa significhi però

-Bharani

Autorizzato sotto: CC-BY-SA insieme a attribuzione
Non affiliato a StackOverflow
scroll top