ggplot2：绘制多个页面

Question

我正在使用ggplot2绘制图形，我有facet_wrap函数来制作多个图形（~51），但它们都出现在一个页面上。现在搜索后，我发现ggplot2无法在多个页面上放置图形。

有办法做到这一点吗？我查看了这个问题（Multiple graphs over multiple pages using ggplot）并尝试了代码，但收效甚微。

这是我的图形代码，它在一个页面上生成~51个图形，使它们非常小而且难以看到，如果我可以在pdf中每页打印1个图形，那就太棒了：

ggplot(indbill, aes(x = prey, y = weight), tab) +
geom_polygon(aes(group = load, color = capture), fill = NA, size = 0.75) +
facet_wrap(~ individual) +
theme(axis.ticks.x = element_blank(),
    axis.text.x = element_text(size=rel(0.5)),
    axis.ticks.y = element_blank(),
    axis.text.y = element_blank()) +
xlab("") + ylab("") +
guides(color = guide_legend(ncol=2)) +
coord_radar()

如果有人可以写一些代码并向我解释，那就太棒了！

谢谢！

Answer 1

一种选择是使用您现在使用的相同代码，一次只绘制六个individual级别。您只需要对数据的每个子集进行一次迭代几次。您尚未提供样本数据，因此以下是使用Baseball数据框的示例：

library(ggplot2)
library(vcd) # For the Baseball data
data(Baseball)

pdf("baseball.pdf", 7, 5)
for (i in seq(1, length(unique(Baseball$team87)), 6)) {
   print(ggplot(Baseball[Baseball$team87 %in% levels(Baseball$team87)[i:(i+5)], ], 
                  aes(hits86, sal87)) + 
    geom_point() +
    facet_wrap(~ team87) +
    scale_y_continuous(limits=c(0, max(Baseball$sal87, na.rm=TRUE))) +
    scale_x_continuous(limits=c(0, max(Baseball$hits86))) +
    theme_bw())
}
dev.off()

上面的代码将生成一个包含四页图表的PDF文件，每个图表有一个页面的六个面。您还可以创建四个单独的PDF文件，每组六个方面一个：

for (i in seq(1, length(unique(Baseball$team87)), 6)) {
pdf(paste0("baseball_",i,".pdf"), 7, 5)
  ...ggplot code...
dev.off()
}

如果您需要更多灵活性，另一个选项是为构面变量的每个级别（即每个唯一值）创建单独的图，并将所有单个图保存在列表中。然后，您可以在每页上布置任意数量的图。这可能有点过头了，但here's an example灵活性派上用场。

首先，让我们创建所有的情节。我们将team87用作我们的分面列。所以我们想为team87的每个级别制作一个图。我们通过将数据分割为team87并为数据的每个子集制作单独的图来做到这一点。

在下面的代码中，split将数据拆分为每个级别team87的单独数据框。 lapply包装器将每个数据子集依次提供到ggplot中，以便为每个团队创建一个绘图。我们将输出保存在plist中，这是一个（在这种情况下）24个图表的列表。

plist = lapply(split(Baseball, Baseball$team87), function(d) {
  ggplot(d, aes(hits86, sal87)) + 
    geom_point() +
    facet_wrap(~ team87) +
    scale_y_continuous(limits=c(0, max(Baseball$sal87, na.rm=TRUE))) +
    scale_x_continuous(limits=c(0, max(Baseball$hits86))) +
    theme_bw() +
    theme(plot.margin=unit(rep(0.4,4),"lines"),
          axis.title=element_blank())
})

现在我们将在PDF文件中列出六个图表。下面是两个选项，一个包含四个单独的PDF文件，每个包含六个图，另一个包含一个四页的PDF文件。我还粘贴在底部的一块地块上。我们使用grid.arrange来绘制图表，包括使用left和bottom参数来添加轴标题。

library(gridExtra)

# Four separate single-page PDF files, each with six plots
for (i in seq(1, length(plist), 6)) {
  pdf(paste0("baseball_",i,".pdf"), 7, 5)
  grid.arrange(grobs=plist[i:(i+5)], 
               ncol=3, left="Salary 1987", bottom="Hits 1986")
  dev.off()
}

# Four pages of plots in one PDF file
pdf("baseball.pdf", 7, 5)
for (i in seq(1, length(plist), 6)) {
  grid.arrange(grobs=plist[i:(i+5)], 
               ncol=3, left="Salary 1987", bottom="Hits 1986")
}
dev.off()

Answer 2

有多种分页方法：ggforce或gridExtra::marrangeGrob。另请参见此answer。

ggforce：

library(ggplot2)
library(ggforce)

# Standard facetting: too many small plots
ggplot(diamonds) +
  geom_point(aes(carat, price), alpha = 0.1) +
  facet_wrap(~cut:clarity, ncol = 3)

# Pagination: page 1
ggplot(diamonds) +
  geom_point(aes(carat, price), alpha = 0.1) +
  facet_wrap_paginate(~cut:clarity, ncol = 3, nrow = 3, page = 1)

# Pagination: page 2
ggplot(diamonds) +
  geom_point(aes(carat, price), alpha = 0.1) +
  facet_wrap_paginate(~cut:clarity, ncol = 3, nrow = 3, page = 2)

# Works with grid as well
ggplot(diamonds) +
  geom_point(aes(carat, price), alpha = 0.1) +
  facet_grid_paginate(color~cut:clarity, ncol = 3, nrow = 3, page = 4)

gridExtra：

library(gridExtra)

set.seed(123)
pl <- lapply(1:11, function(.x) 
  qplot(1:10, rnorm(10), main=paste("plot", .x)))

ml <- marrangeGrob(pl, nrow=2, ncol=2)

## non-interactive use, multipage pdf
## ggsave("multipage.pdf", ml)
## interactive use; calling `dev.new` multiple times
ml

由reprex package（v0.2.0.9000）创建于2018-08-09。

Answer 3

类似的东西：

by(indbill, indbill$individual, function (x){
    ggplot(x, aes(x = prey, y = weight), tab) +
    geom_polygon(aes(group = load, color = capture), fill = NA, size = 0.75) +
    theme(axis.ticks.x = element_blank(),
        axis.text.x = element_text(size=rel(0.5)),
        axis.ticks.y = element_blank(),
        axis.text.y = element_blank()) +
    xlab("") + ylab("") +
    guides(color = guide_legend(ncol=2)) +
    coord_radar()
}