r - lm 中 "call"的提取函数

Question

是否有用于线性模型的通用提取器函数来从模型中提取调用？我正在寻找跨数据帧列表应用模型，然后提取对该模型的调用(我在许多数据帧和模型上做了很多次，所以寻找一种简单的方法来跟踪哪个模型是哪个)

示例数据:

set.seed(1234)
library(plyr)
    data <- data.frame(
        biz = sample(c("telco","shipping","tech"), 50, replace = TRUE),
        region = sample(c("mideast","americas"), 50, replace = TRUE),
        orders = sample(1:50, 50, replace=TRUE),
        revenue = sample(100:150, 50, replace=TRUE)
        )

现在，如果我想通过“区域”对这些数据应用一个简单的 lm:

modlist <- dlply(data, 'region', function(df){
    summary(lm(revenue ~ orders, data=df))
    })

对于具有标准提取函数的东西，比如“coefficients” = “coef”，这很简单:

coefs <- ldply(modlist, coef)

但我找不到任何方法来提取“调用”，例如“lm(formula = revenue ~ orders, data = df)”。有任何想法吗？

Answer 1

只需编写您自己的函数，以与 stats:::print.lm 相同的方式提取调用

coefs <- ldply(modlist, function(x) deparse(x$call))