Использование функции summary.lm в Rapache - PullRequest
Новая
       26

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

0 голосов
/ 22 октября 2009

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

Кто-нибудь сталкивался с подобной проблемой с R и rapache? Я относительно новичок в 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))

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

Спасибо Bharani

Ответы [ 4 ]

2 голосов
/ 22 октября 2009

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

1 голос
/ 22 октября 2009

Я проверил следующий пример и 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
0 голосов
/ 23 октября 2009

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

0 голосов
/ 23 октября 2009

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

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

-Bharani

...