是否将多个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个调用设置了4个独立的设备，只有后者被使用并随后关闭，其他3个保持打开状态。

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的子类，没有自己特定的plot方法)根据which参数显示诊断图，which的默认值与上面的wplot相同。所以基本上:弄清楚如何创建单独的子图。