从模型摘要列表中保存单个模型输出

Question

我在回归模型中循环了一个变量列表，并将输出保存为列表，如何将列表中的模型输出到单个文本文件中？

  chr  start     end     
1 chr7 123268311 123268334
2 chr7 123272067 123272207 
3 chr7 123272905 123273039

谢谢！

Answer 1

也许更好的方法是使用：

for(i in seq_along(mlist_summ)){
  stargazer(mlist_summ[[i]]$coefficients, type = "text", out = paste0("summary_", i, ".txt"))
}

这样您就可以控制文件名。

包含lm和mtcars数据的完整示例：

library(stargazer)
data(mtcars)

varlist <- list("cyl", "disp", "hp")

mlist <- lapply(varlist, function(x) {
  lm(substitute(mpg ~ i, 
                  list(i = as.name(x))), data = mtcars)
})

mlist_summ <- lapply(mlist, summary)

for(i in seq_along(mlist_summ)){
  stargazer(mlist_summ[[i]]$coefficients, type = "text", out = paste0("summary_", i, ".txt"))
}

另一种方法是首先提取系数：

mlist_summ <- lapply(mlist, function(x){
  summary(x)$coefficients})

并将该列表提供给观星者：

stargazer(mlist_summ, type = "text", out = "summary.txt")

这样所有模型都在一个文件中。