在循环中保存图R

Question

昨天我将R升级到版本3.0.0，将ggplot2升级到版本0.9.3.1（并对我的脚本做了一些小修改）。现在，当我尝试保存绘图时出现错误 - 不幸的是，错误不会以较小的数据框来重现，因此我已经包含了生成相同大小的代码。

library("ggplot2") 

# Create data frame 
# Time interval ID (x) 
bin.ts.avg <- as.data.frame(rep(1:18, 31)) 
names(bin.ts.avg) <- "x" 
# Time (sequence of 10 minuter intervals between 7am and 10am) 
tt.month.bins <- seq(from=as.POSIXct("2012-01-01 GMT"), to=as.POSIXct("2012-01-01 GMT") + 60*60*24*31, by="10 mins") 
tt.month.bins <- tt.month.bins[-length(tt.month.bins)] 
temp <- as.numeric(format(tt.month.bins, "%H")) 
ind <- which(temp >=7 & temp <= 9) 
tt.month.bins <- tt.month.bins[ind] 
bin.ts.avg$dep <- tt.month.bins 
# Value (with some NA) 
bin.ts.avg$tt <- runif(558, min=2.5, max=5) 
bin.ts.avg$tt[trunc(runif(200, min=1, max=558))] <- NA 
# Day 
bin.ts.avg$depday <- rep(1:31, each=18) 

for (i in 1:2){ 
    if (1){ 
    hist(rnorm(100)) 
    dev.print(file="MyHist.png",device=png, bg="white", width=640, height=352) 

    p <- ggplot(bin.ts.avg, aes(x, tt)) + geom_point() +geom_line() + facet_grid(.~depday) 
    p <- p + ggtitle("10 minute averages")+ xlab("Hour") + ylab("Values")  
    p <- p + scale_x_continuous(breaks=c(min(bin.ts.avg$x), max(bin.ts.avg$x)), labels=c("7", "10")) 
    print(p) 
    dev.print(file="MyGGPlot.png",device=png, bg="white", width=640, height=352) 
    } 
} 

在运行该脚本时，收到以下错误消息：

错误UseMethod（“深度”）：无适用方法为“深度” 应用于类“NULL的目的“

但是，如果我逐行运行脚本，一切运行正常（下图）。 现在，如果我改变for循环和使用dev.copy，而是ggsave dev.print，如下

for (i in 1:2){ 
    if (1){ 
    hist(rnorm(100)) 
    dev.copy(file="MyHist.png",device=png, bg="white", width=640, height=352) 
    dev.off() 

    p <- ggplot(bin.ts.avg, aes(x, tt)) + geom_point() +geom_line() + facet_grid(.~depday) 
    p <- p + ggtitle("10 minute averages")+ xlab("Hour") + ylab("Values")  
    p <- p + scale_x_continuous(breaks=c(min(bin.ts.avg$x), max(bin.ts.avg$x)), labels=c("7", "10")) 
    print(p) 
    ggsave(filename="MyGGPlot.png") 
    } 
} 

在试图打开“MyGGPlot.png”使用画图，我收到一条错误信息，说明

A sharing violation occurred while accessing <filename> 

我使用RStudio版本0.97.449运行脚本。 关于我需要更改以保存当前地块的任何想法？

Answer 1

dev.copy后，一对夫妇点

使用graphics.off()的。这将关闭所有图形设备。你也可以拨打dev.off()两次（但graphics.off()是一种包装，基本上都会叫dev.off()足够的时间以关闭所有图形设备

ggsave不需要print ED对象（这是不其中FAQ 7.22是相关的情况下） 。

默认是last_plot是创建的最后一个对象ggplot，改性的或印刷的，所以创造p值足以 ggsave('filname.png')保存该对象。

for (i in 1:2){ 
    if (1){ 
    hist(rnorm(100)) 
    dev.copy(file="MyHist.png",device=png, bg="white", width=640, height=352) 
    graphics.off() 

    p <- ggplot(bin.ts.avg, aes(x, tt)) + geom_point() +geom_line() + facet_grid(.~depday) 
    p <- p + ggtitle("10 minute averages")+ xlab("Hour") + ylab("Values")  
    p <- p + scale_x_continuous(breaks=c(min(bin.ts.avg$x), max(bin.ts.avg$x)), labels=c("7", "10")) 
    # no need to print p 
    ggsave(filename="MyGGPlot.png") 
    # note specifying p is redundant but explicit. 
    # ggsave(filename = 'MyGGplot.png', plot = p) 
    } 
}