Utilisation de la fonction summary.lm dans rapache

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

  •  05-07-2019
  •  | 
  •  

Question

J'ai installé rapache et j'essaie d'intégrer un modèle linéaire à l'intérieur du R Fichier de script. J'ai configuré le RFileHandler dans le fichier http.conf. Lorsque j'essaie d'appeler le résumé (modèle), cela me donne une erreur de segment (je le vois dans le fichier journal Apache). Je devine qu'il essaie d'imprimer sur la console et c'est pourquoi il échoue.

Quelqu'un at-il rencontré un problème similaire avec R et rapache ? Je suis relativement nouveau en R et summary fait beaucoup de choses qui ne sont pas directement exposées en tant que fonctions, alors j’espère pouvoir le faire fonctionner

Voici mon 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

si je remplace 

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

avec ceci 

    cat(coef(fit))

ça marche!

Merci Bharani

Était-ce utile?

La solution 4

J'ai finalement compris le problème. Lecture de la discussion i à tort libRlapck.so à lapack.so. On dirait que ça causait problèmes. A fait une nouvelle installation de R puis modifié apache à explicity charger les bibliothèques puis tout a fonctionné Merci  - Bharani

Autres conseils

Avez-vous envisagé de contacter le groupe Google rapache en tant que page d'accueil de rapache suggère? Vous y trouverez peut-être un lecteur plus expérimenté qu’ici.

J'ai testé l'exemple suivant et cat (résumé (ajustement) $ adj.r.squared)  fonctionne dans ma configuration (par défaut) (les dernières versions de rapache 1.1.8 et R 2.9.2 sous 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

Je viens de découvrir que ce n’est pas avec le rapache. Il échoue dans R lui-même 

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

ne sais pas ce que cela signifie bien

-Bharani

Licencié sous: CC-BY-SA avec attribution
Non affilié à StackOverflow
scroll top