如何在R中保存线性模型plot(model)生成的所有绘图？

方法1：逐张保存为单独的PNG文件 #

利用plot.lm()的which参数指定要绘制的子图，循环生成并保存每张图：

save_model_plots <- function(model, prefix = "model_plot") {
  # 遍历4张子图
  for (i in 1:4) {
    # 设置带序号的输出文件名
    png(paste0(prefix, "_", i, ".png"))
    # 绘制第i张残差图
    plot(model, which = i)
    # 关闭绘图设备
    dev.off()
  }
}

# 测试代码
x <- 1:10
y <- 1:10
model <- lm(x ~ y)
save_model_plots(model)

运行后会生成model_plot_1.png到model_plot_4.png四个文件，分别对应四张残差图。

方法2：保存所有图到单个PDF文件 #

PDF支持多页绘图，直接用pdf()包裹plot(model)即可一次性保存所有图：

save_model_plots_pdf <- function(model, file = "model_plots.pdf") {
  pdf(file)
  plot(model)
  dev.off()
}

# 测试代码
save_model_plots_pdf(model)

打开生成的PDF文件，里面会包含4页，每页对应一张残差图。

方法3：将所有图合并为单张图片 #

如果想把4张图放到同一张图片里，可以用par()设置绘图布局：

save_combined_plot <- function(model, file = "combined_plot.png") {
  png(file, width = 800, height = 800)
  # 设置2x2的绘图布局
  par(mfrow = c(2, 2))
  plot(model)
  par(mfrow = c(1, 1)) # 恢复默认布局
  dev.off()
}

# 测试代码
save_combined_plot(model)

生成的图片会以2行2列的形式展示所有4张残差图。

内容的提问来源于stack exchange，提问作者stats_b