Question

以下代码完美无缺，但最后一行除外，似乎被忽略了。代码创建一个由三列组成的数据帧，然后继续绘制2x2网格上的列。最后一行应该是整个多面板情节的标题。但是，它对结果没有影响。怎么了？我该如何纠正？

p <- c(3, 5, 10, 20, 50, 100)
n <- c(100, 5000, 100000)
f <- function(x, y){ return ((1 - 2^(-1/y))^(1/x))}

d <- as.data.frame(outer(p, n, FUN = f))
dimnames(d) <- list(p, n)

par(mfrow = c(2,2))    
lapply(colnames(d),
       function(x) plot(p, d[,x], type = "b",
                        main = paste("#points = ", x),
                        xlab = "Dim",
                        ylab = "Med Dist"))
mtext("Densities", outer = TRUE, cex = 1.5)

Answer 1

尝试

par(oma = c(0, 0, 2, 0))
par(mfrow = c(2,2))    
lapply(colnames(d),
      function(x) plot(p, d[,x], type = "b",
                    main = paste("#points = ", x),
                    xlab = "Dim",
                    ylab = "Med Dist"))
title("Densities", outer=TRUE)

多面板图上的标题

1 个答案: