1

我正在制作一个表格，通过 texreg 包比较使用 r 和乳胶的不同线性模型（这是我第一次使用这个包）。我发现这个包非常有用，但我真的不明白为什么 RMSE 不包含在 extract.lm 方法中。

我以这种方式修改了 extract.lm 方法，将其作为选项包含在内：

extract.lm<-  function (model, include.rsquared = TRUE, include.adjrs = TRUE, include.rmse = TRUE, 
                include.nobs = TRUE, ...) 
      {
        s <- summary(model, ...)
        names <- rownames(s$coef)
        co <- s$coef[, 1]
        se <- s$coef[, 2]
        pval <- s$coef[, 4]
        rs <- s$r.squared
        adj <- s$adj.r.squared
        sig<-s$sigma ##added it
        n <- nobs(model)
        gof <- numeric()
        gof.names <- character()
        gof.decimal <- logical()
        if (include.rsquared == TRUE) {
          gof <- c(gof, rs)
          gof.names <- c(gof.names, "R$^2$")
          gof.decimal <- c(gof.decimal, TRUE)
        }
        if (include.adjrs == TRUE) {
          gof <- c(gof, adj)
          gof.names <- c(gof.names, "Adj. R$^2$")
          gof.decimal <- c(gof.decimal, TRUE)
        }
        if (include.rmse == TRUE) { ##added it
          gof <- c(gof, sig)
          gof.names <- c(gof.names, "RMSE")
          gof.decimal <- c(gof.decimal, TRUE)
        }
        if (include.nobs == TRUE) {
          gof <- c(gof, n)
          gof.names <- c(gof.names, "Num. obs.")
          gof.decimal <- c(gof.decimal, FALSE)
        }
        tr <- createTexreg(coef.names = names, coef = co, se = se, 
                           pvalues = pval, gof.names = gof.names, gof = gof, gof.decimal = gof.decimal)
        return(tr)
      }

然后我在控制台中运行它并搜索以明确覆盖它使用

setMethod("extract", signature=className("lm","stats"), definition=extract.lm, where="package:texreg")

这也是作者在rforge上建议的方式，但我有这个错误

Error in setMethod ...
  the environment "texreg" is locked; cannot assign methods for function

我也尝试使用 fixInNamespace()，但未显示 body 函数，即使在终端中插入函数的名称或使用 body() 函数可以访问函数的源（这是我访问的方式和修改它）。

我的错在哪里？？我怎样才能做到这一点？

注意：extract 是 S4 通用的。

请，我想在包环境中明确地覆盖该功能。

谢谢！

4

1 回答 1

2

该错误是由您尝试覆盖包环境中的对象引起的。texreg我建议在您的本地环境中创建该方法，当您与该方法结合使用时，R 将优先考虑该extract.lm方法：

setMethod("extract", signature = className("lm", "stats"), 
    definition = extract.lm)

请参阅《统计软件杂志》中文章的第 6 节：http: //www.jstatsoft.org/v55/i08/

于 2014-05-20T09:53:27.963 回答