我创建了一个函数,该函数在循环中使用不同的列作为因变量来运行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列(
家庭的影响 | 朋友的影响 | 休闲的影响 |
---|---|---|
2 | 1 | 1 |
1 | 2 | 1 |
由于您没有提供最小工作示例,我们无法正确诊断您的问题。
也就是说,一种选择是确保您的函数返回一个模型列表,然后将其提供给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))