用 R ggplot 绘制数百个图

Question

我有 141 个数值变量的季度预测值和实际值，这些变量是针对 3 个不同国家/地区测量的。我想找出我的预测对于哪些变量“足够好”。当我进行一些早期探索时，我试图为 3 * 141 = 423 对（国家，变量）中的每一对绘制实际值和预测。我不知道如何用 ggplot2 制作那么多图。

我已经用所有变量融化了数据框，所以它现在有列：

 $ Date    : 'yearqtr' num  2017 Q1 2017 Q2 2017 Q3 2017 Q4 ...
 $ Country : Factor w/ 3 levels "CN","DE","US": 1 1 1 1 1 1 1 1 1 1 ...
 $ variable: Factor w/ 141 levels "BCU_D","BUSCONFM",..: 1 1 1 1 1 1 1 1 1 1 
 $ value   : num  21382 56639 45900 71196 -34100 ...
 $ isActual: logi  TRUE TRUE TRUE TRUE TRUE TRUE ...

然后，我尝试了一个简单的 ggplot 调用：

ggplot(meltedData, aes(x = Date, y = value, color = isActual)) + geom_line() + facet_wrap(Country ~ variable, ncol = 3) + scale_x_yearqtr(format = "%YQ%q", n = 5)

我得到的是一些无用的灰色框。

我想得到如上所述的 141 * 3 图表。

我想制作一个 141 页的 pdf 文档，为每个变量逐页绘制一组图，如果这样更容易的话。

谢谢！

============

@Gregor 建议的解决方案（循环变量）很令人满意：

for (this_variable in unique(meltedData$variable)) {
  g = ggplot(meltedData[meltedData$variable == this_variable, ],
             aes(x = Date, y = value, color = isActual)) +
               geom_line() +
               facet_wrap(. ~ Country, nrow = 3) + 
               scale_x_yearqtr(format = "%YQ%q", n = 5) + 
               labs(title = this_variable)
  ggsave(filename = paste0(this_variable, ".png"), plot = g, path = "output_graphs")
}

Answer 1

facet_wrap（或facet_grid）将为您提供一个绘图网格 - 您将无法在行之间进行分页。对于 141 行，您的选择是创建一个 非常 高的图像文件，可能太笨重而无法使用。

更好的方法是将每一行（或一组行）保存到不同的文件中。由于我没有您的数据，我无法测试解决方案，但这里有一个使用内置数据的示例：

library(ggplot2)
for (i in unique(mtcars$cyl)) {
  g = ggplot(mtcars[mtcars$cyl == i, ], aes(x = wt, y = mpg)) + 
    geom_point() +
    facet_wrap(~ am) +
    labs(title = paste("Cylinder:", i))
  ggsave(filename = paste0("cylinder_", i, ".png"), plot = g)
}

根据您的数据，我认为它看起来像这样（未经测试）：

for(this_var in unique(meltedData$variable)) {
  g = ggplot(meltedData[meltedData$variable == this_var, ],
             aes(x = Date, y = value, color = isActual)) + 
    geom_line() + 
    facet_wrap(~ Country, ncol = 3) +
    scale_x_yearqtr(format = "%YQ%q", n = 5) +
    labs(title = this_var)
  ggsave(filename = paste0(this_var, ".png"), plot = g)
}