Использование функции summary.lm в rapache

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

  •  05-07-2019
  •  | 
  •  

Вопрос

я установил рапаш и я пытаюсь поместить линейную модель в файл сценария R.Я настроил RFileHandler в http.conf.Когда я пытаюсь вызвать сводку (модель), она выдает мне ошибку ошибки сегмента (я вижу это в файле журнала Apache).Я предполагаю, что он пытается напечатать на консоли и поэтому терпит неудачу.

Кто-нибудь сталкивался с подобной проблемой с R и rapache?Я относительно новичок в R, и summary делает много вещей, которые не представлены напрямую как функции, поэтому я надеюсь, что смогу заставить его работать.

Вот мой 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

если я заменю 

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

с этим

    cat(coef(fit))

это работает!

Спасибо Бхарани

Это было полезно?

Решение 4

Я наконец понял проблему.Чтение обсуждение i ошибка libRlapck.so в lapack.so.Похоже, это вызвало проблемы.Сделал чистую установку R снова, а затем модифицировал apache для явности загрузите библиотеки, то все это работало Спасибо - Бхарани

Другие советы

Вы рассматривали возможность связаться с Рапаш Группа Google как Домашняя страница рапача предполагает?Возможно, там вы найдете более опытного читателя, чем здесь.

Я протестировал следующий пример и cat(summary(fit)$adj.r.squared)работает в моей настройке (по умолчанию) (последние версии Rapache 1.1.8 и R 2.9.2 под 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

Только что узнал, что это не с рапачем.Он терпит неудачу в самом 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

хотя не уверен, что это значит

-Бхарани

Лицензировано под: CC-BY-SA с атрибуция
Не связан с StackOverflow
scroll top