应该将 xts 类数据绘制到 png 文件的 R 脚本不会生成 png 文件

Question

我有一个工作脚本，它使用循环生成许多 xts 类数据的 png 文件图。现在脚本抛出一个错误，如果我注释掉抛出错误的行（调用 abline() ）然后脚本执行但不生成 png 文件。该问题似乎涉及绘制 xts 类数据和/或循环或脚本。

在 stackoverflow 上搜索并没有提供对此问题的解决方案或参考。我在以下示例中重现了该问题。在实践中，脚本会在循环和重要数据中使用不同的文件名。

# put following code in 'myscript.R' and execute using source('myscript.R',print.eval=TRUE) or  source('myscript.R')
# xts class data
data <- xts(seq(1:10),order.by=as.Date(seq(1:10)))

# a non xts version of same data
#data <- seq(1:10)

for(i in 1:1) {
filename <- 'myfile.png'
png(filename)
plot(data)
lines( (data-1),col='red')
abline(h=1)
dev.off()
}
# The call to abline in above script with xts class data gives error 'plot.new has not been called yet'
# If comment out the call to abline it completes but doesn't produce a png file
# script works fine with abline for non xts data

Answer 1

使用xts >= 0.10.1，这会将您想要的内容保存到文件中

for(i in 1:1) {
    filename <- 'myfile.png'
    png(filename)
    plot(data)
    print(lines( (data-1),col='red', on = 1))
    print(lines(xts(x = rep(1, NROW(data)), order.by = index(data)),col='green', on=  1)    )
    dev.off()
}

使用print 调用额外的行。我也会使用lines 作为水平线，而不是abline，因为这与使用 xts 绘图更一致。

如果你这样做 print(abline(h=1)) 也可以避免你的错误