将多个 glm 图导出为 PNG？

Question

所以，

我正在尝试导出线性模型的图。当我将其制作为 PDF 时，PDF 有四页和四个不同的图表。当我导出为 PNG 时，我只得到第一个图表。如何导出以便将所有四个图表作为单独的 PNG 文件获取？

什么适用于 PDF：

lrfitX11SUMS <- glm(X11SUMS ~ POLITA + CIRCULATION + COMPAPER + TrafficRankUS, data=NewspTest)

    summary(lrfitOTONE)
    pdf("/Users/william/Desktop/output/lmfitOTONE1.pdf")
    plot(lrfitOTONE)
    dev.off()

什么不适用于 PNG（花了两个小时在互联网和绘图文档中进行挖掘，但无济于事）：

lrfitX11SUMS <- glm(X11SUMS ~ POLITA + CIRCULATION + COMPAPER + TrafficRankUS, data=NewspTest)

summary(lrfitOTONE)
png("/Users/william/Desktop/output/lmfitOTONE1.png", width=720, height=720, pointsize=16)
png("/Users/william/Desktop/output/lmfitOTONE2.png", width=720, height=720, pointsize=16)
png("/Users/william/Desktop/output/lmfitOTONE3.png", width=720, height=720, pointsize=16)
png("/Users/william/Desktop/output/lmfitOTONE4.png", width=720, height=720, pointsize=16)
plot(lrfitOTONE)
dev.off()

如何获取我的图像？

谢谢，

-Wm

Answer 1

PDF 允许多页文档。 PNG 图像从根本上与这个想法不兼容。阅读 ?png 并意识到需要查看 filename 参数，将会引导您前往 ?postscript 了解详细信息。

您想要类似的内容：

png("/Users/william/Desktop/output/lmfitOTONE%1d.png", width=720, 
    height=720, pointsize=16)
plot(lrfitOTONE)
dev.off()

其中文件名中的 %1d 是一个通配符，它​​扩展为 1 位数字值，这样您就可以得到带有您想要的名称的四个数字。您对 png() 的 4 次调用设置了四个单独的设备，只有后者被使用并随后关闭，其他三个保持打开状态。

Answer 2

像这样的事情：

setwd("/Users/william/Desktop/output/")
tmpf <- function(i) {
   png(paste("lmfitOTONE",i,".png",sep=""), width=720, height=720, pointsize=16)
}
wplot <- c(1,2,3,5) ## see ?plot.lm for definition of 'which'
for (i in seq_along(wplot)) {
   tmpf(i); plot(lrfitOTONE, which=wplot[i]); dev.off()
}

关键是实现 plot.lm （这是 plot 使用的方法应用于 glm 对象，因为 glm 是 lm 的子类，没有自己特定的绘图方法）根据 which 参数显示诊断图，默认值为which 与 wplot 相同 多于。所以基本上：弄清楚如何创建各个子图。