r-通过单个模型为从循环中提取的多个迭代创建摘要



我创建了一个函数,该函数在循环中使用不同的列作为因变量来运行lm回归的多次迭代。我正在提取每次迭代的摘要和关系图,但我无法创建所有迭代结果的单个摘要表。由于我只有8个专栏,我认为这是可以做到的。这是我的函数,数据低于

quantmodel<-function(a){
i<-1
a <- janitor::clean_names(a)
colnames1 <- colnames(a)
lm_model <- linear_reg() %>% 
set_engine('lm') %>%
set_mode('regression')

out_lst <- vector('list', ncol(a))

for (i in seq_along(a)) {
lm_fit <- lm_model %>% 
fit(as.formula(paste(colnames1[i], "~ .")), data = a)

#Saving relevance plot of each parameter
temp_plot = vip(lm_fit ,geom = "col", aesthetics = list(color = "black", fill = "black"))
ggsave(temp_plot, file=paste0("plot_", i,".png"), width = 14, height = 10, units = "cm")
#Saving pdf of individual summaries
pdf(paste0("summary_", colnames1[i],".pdf"), width = 10,height = 3) 
grid.table(coef(summary(lm_fit$fit)))
dev.off()
paste0("m",i)<-lm_fit$fit
}

}
quantmodel(set1)

数据:第1组(前3列(

家庭的影响朋友的影响休闲的影响
211
121

由于您没有提供最小工作示例,我们无法正确诊断您的问题。

也就是说,一种选择是确保您的函数返回一个模型列表,然后将其提供给modelsummary函数。在这个示例代码中,请注意末尾的return()调用和modelsummary()调用:

quantmodel<-function(a){
i<-1
a <- janitor::clean_names(a)
colnames1 <- colnames(a)
lm_model <- linear_reg() %>% 
set_engine('lm') %>%
set_mode('regression')

out_lst <- list()

for (i in seq_along(a)) {
lm_fit <- lm_model %>% 
fit(as.formula(paste(colnames1[i], "~ .")), data = a)
out_lst[[i]] <- lm_fit

#Saving relevance plot of each parameter
temp_plot = vip(lm_fit ,geom = "col", aesthetics = list(color = "black", fill = "black"))
ggsave(temp_plot, file=paste0("plot_", i,".png"), width = 14, height = 10, units = "cm")
#Saving pdf of individual summaries
pdf(paste0("summary_", colnames1[i],".pdf"), width = 10,height = 3) 
grid.table(coef(summary(lm_fit$fit)))
dev.off()
paste0("m",i)<-lm_fit$fit
}
return(out_lst)

}
library(modelsummary)
models <- quantmodel(set1)
modelsummary(dvnames(models))

最新更新