将多个对齐的绘图放置在一页上时避免浪费空间

我想将四个图放在一个页面上。轴标签应仅打印在最边缘，即x仅用于底部图表的轴标签，以及y仅左图的轴标签。这既适用于整个轴的名称，也适用于各个刻度线。我可以使用以下代码生成一些内容：

pdf(file = "ExampleOutput.pdf",
    width = 6.61,
    height = 6.61,
    pointsize = 10
    )
set.seed(42)
catA <- factor(c("m100", "m500", "m1000", "m2000", "m3000", "m5000"))
catB <- factor(20:28)
samples <- 100
rsample <- function(v) v[ceiling(runif(samples, max=length(v)))]
Tab <- data.frame(catA = rsample(catA),
                  catB = rsample(catB),
                  valA = rnorm(samples, 150, 8),
                  valB = pmin(1,pmax(0,rnorm(samples, 0.5, 0.3))))
par(mfrow = c(2,2))
for (i in 0:3) {
  x <- Tab[[1 + i %% 2]]
  plot(x, Tab[[3 + i %/% 2]],
       xlab = if (i %/% 2 == 1) "Some Categories" else NULL,
       ylab = if (i %% 2 == 0) "Some Values" else NULL,
       axes = FALSE
       )
  axis(side = 1,
       at=1:nlevels(x),
       labels = if (i %/% 2 == 1) levels(x) else FALSE)
  axis(side = 2, labels = (i %% 2 == 0))
  box(which = "plot", bty = "l")
}
par(mfrow = c(1,1))
dev.off()

我欢迎有关如何改进绘图命令的建议，也许可以避免手动耗尽左下角的轴和 L 。但这只是一个除此之外的事情。

该序列的结果如下所示：

这里的问题是huge浪费的空白量。我的印象是，R 为轴和刻度标签保留了空间，即使它们没有被使用。由于浪费了空间，对于左下图，只有每秒x蜱虫实际上被标记了，这在这里真的很糟糕。

我想生成一个类似的图，但没有那么多空白。实际的图应该具有相同的大小，因此它们可以正确排列，但标签的空间应该仅位于外部。我想象这样的布局（在 GIMP 中创建的模型）：

我怎样才能实现这样的布局？

这是对您显示的一般图的轻微修改，假设 y 和 x 轴标签属于所有图。它使用外边距来包含轴标签，我们将其添加为title()使用论证outer = TRUE。效果有点像中的标签ggplot2 or lattice plots.

这里的关键行是：

op <- par(mfrow = c(2,2),
          oma = c(5,4,0,0) + 0.1,
          mar = c(0,0,1,1) + 0.1)

它设置绘图参数（调用之前的值存储在op）。我们用5 and 4第 1 条和第 2 条边上的线作为外边距，这是通常的数量mar范围。绘制区域边距 (mar）各 1 行添加到顶部和右侧，以便在图之间留出一点空间。

添加轴标签after the for()循环与

title(xlab = "Some Categories",
      ylab = "Some Values",
      outer = TRUE, line = 3)

整个脚本是：

set.seed(42)
catA <- factor(c("m100", "m500", "m1000", "m2000", "m3000", "m5000"))
catB <- factor(20:28)
samples <- 100
rsample <- function(v) v[ceiling(runif(samples, max=length(v)))]
Tab <- data.frame(catA = rsample(catA),
                  catB = rsample(catB),
                  valA = rnorm(samples, 150, 8),
                  valB = pmin(1,pmax(0,rnorm(samples, 0.5, 0.3))))
op <- par(mfrow = c(2,2),
          oma = c(5,4,0,0) + 0.1,
          mar = c(0,0,1,1) + 0.1)
for (i in 0:3) {
  x <- Tab[[1 + i %% 2]]
  plot(x, Tab[[3 + i %/% 2]], axes = FALSE)
  axis(side = 1,
       at=1:nlevels(x),
       labels = if (i %/% 2 == 1) levels(x) else FALSE)
  axis(side = 2, labels = (i %% 2 == 0))
  box(which = "plot", bty = "l")
}
title(xlab = "Some Categories",
      ylab = "Some Values",
      outer = TRUE, line = 3)
par(op)

产生