在rapache中使用summary.lm函数

时间:2009-10-22 13:57:32

标签: r rapache

我已安装rapache,我正在尝试在R脚本文件中插入线性模型。我在http.conf中配置了RFileHandler。当我试图调用摘要(模型)时,它给我一个段错误错误(我在apache日志文件中看到这一点)。我猜它正在尝试打印到控制台,这就是它失败的原因。

是否有人遇到与R和rapache类似的问题?我对R比较新,总结是做了很多没有直接作为函数公开的东西所以我希望我可以让它工作

这是我的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 个答案:

答案 0 :(得分:2)

您是否考虑过rapache Google Grouprapache home page建议的联系?你可能会找到比这里更有经验的读者。

答案 1 :(得分:1)

我测试了以下示例和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

答案 2 :(得分:0)

刚刚发现它不是rapache。 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

答案 3 :(得分:0)

我终于找到了问题所在。阅读discussion i 错误的libRlapck.so到lapack.so。看起来这是造成的 问题。 是否彻底安装了R然后修改了apache以明确 加载库然后它都工作 谢谢   - 巴拉尼