¿Exporta múltiples diagramas de glm como PNG?

Question

SO,

Estoy tratando de exportar gráficos de mi modelo lineal. Cuando lo hago como PDF, el PDF tiene cuatro páginas y cuatro gráficos diferentes. Cuando exporto como PNG, solo obtengo el primer gráfico. ¿Cómo exporto para que obtenga los cuatro gráficos como archivos PNG separados?

Lo que funcionó con el PDF:

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

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

Lo que no funcionaba con PNG (y pasó dos horas cavando alrededor en el Internet y en la documentación gráfica en vano):

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() 

¿Cómo obtengo mis imágenes?

Gracias,

-Wm

Answer 1

algo como esto:

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() 
} 

La clave es darse cuenta de que plot.lm (que es el método utilizado por plot aplicado a un objeto glm, porque glm es una subclase de lm y no tiene su propio método de trazado particular) muestra gráficas de diagnóstico basadas en el argumento which, y el valor predeterminado de which es lo mismo que wplot anterior. Entonces, básicamente: descubra cómo crear subtramas individuales.

Answer 2

Un PDF permite documentos de varias páginas. Una imagen PNG es fundamentalmente incompatible con esta idea. Leyendo ?png y apreciando la necesidad de mirar el argumento filename le habría dirigido a ?postscript para más detalles.

quieres algo así como:

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

donde el %1d en el nombre del archivo es un comodín que se expande a un valor numérico de 1 dígito de tal manera que se obtiene cuatro cifras con los nombres que quería. Sus 4 llamadas al png() configuraron cuatro dispositivos separados, solo el último de los cuales se usó y posteriormente se cerró, los otros tres permanecieron abiertos.