用于提取值的 R 循环答案

【问题标题】：R Loop For Extracting Values用于提取值的 R 循环
【发布时间】：2020-02-06 23:50:56
【问题描述】：

hsb2 <- read.csv("https://stats.idre.ucla.edu/stat/data/hsb2.csv")
names(hsb2)
varlist <- names(hsb2)[8:11]

models <- lapply(varlist, function(x) {
  lm(substitute(read ~ i, list(i = as.name(x))), data = hsb2)
})

## look at the first element of the list, model 1
models[[1]]

上面的代码为不同的自变量生成了一系列简单的回归模型。我的首要任务是提取 varlist 中列出的每个变量的系数和标准误差。我的尝试如下所示。

ATTEMPT = lapply(1:length(models), function(x) {
                   cbind(cov, coef(summary(models[[x]]))[2,1:2])})

我希望的输出将显示三列——变量、系数、标准。错误：

【问题讨论】：

将数据重新整形为长格式并使用包 nlme 中的 lmList 会更容易。
看看broom 包。这些可能有用stackoverflow.com/a/55194281/786542 & stackoverflow.com/a/49484854/786542 & stackoverflow.com/a/54624775/786542

标签： r dplyr data.table lapply mapply

【解决方案1】：

怎么样：

ATTEMPT2 = lapply(1:length(models), function(x) {
                    cf <- coef(summary(models[[x]]))
                    data.frame(Variable=rownames(cf)[2], 
                               Estimate=cf[2,1],
                               Std.Error=cf[2,2])})
(df2 <- do.call("rbind", ATTEMPT2))
#   Variable  Estimate  Std.Error
# 1    write 0.6455300 0.06168323
# 2     math 0.7248070 0.05827449
# 3  science 0.6525644 0.05714318
# 4    socst 0.5935322 0.05317162

【讨论】：