你如何在 R 中的多图形环境中画一条线？

Question

举个很简单的例子，mfrow=c(1,3)；每个图都是不同的直方图；我将如何绘制一条跨越所有 3 个数字的水平线（类似于abline(h=10)）？ （也就是说，甚至它们之间的边距。）显然，我可以为每个图形添加一个 abline，但这不是我想要的。我可以想到一个非常复杂的方法来做到这一点，真的只有一个数字，并使用polygon 等在其中绘制每个“数字”。这太荒谬了。有没有简单的方法可以做到这一点？

Answer 1

正如@joran 所说，grid 图形系统可以更灵活地控制单个设备上多个绘图的排列。

这里，我先用grconvertY()查询y轴上高度为50的位置以“标准化设备坐标”为单位。 （即作为绘图设备总高度的比例，0=底部，1=顶部）。然后我使用 grid 函数来：（1）推送一个填充设备的viewport； (2) 在grconvertY() 返回的高度绘制一条线。

## Create three example plots
par(mfrow=c(1,3))
barplot(VADeaths, border = "dark blue") 
barplot(VADeaths, border = "yellow") 
barplot(VADeaths, border = "green") 

## From third plot, get the "normalized device coordinates" of 
## a point at a height of 50 on the y-axis.
(Y <- grconvertY(50, "user", "ndc"))
# [1] 0.314248

## Add the horizontal line using grid
library(grid)
pushViewport(viewport())
grid.lines(x = c(0,1), y = Y, gp = gpar(col = "red"))
popViewport()

编辑：@joran 询问如何绘制一条从第一个图的 y 轴延伸到第三个图的最后一个条形边缘的线。这里有几个选择：

library(grid)
library(gridBase)
par(mfrow=c(1,3))

# barplot #1
barplot(VADeaths, border = "dark blue") 
X1 <- grconvertX(0, "user", "ndc")
# barplot #2
barplot(VADeaths, border = "yellow") 
# barplot #3
m <- barplot(VADeaths, border = "green") 
X2 <- grconvertX(tail(m, 1) + 0.5, "user", "ndc") # default width of bars = 1
Y <- grconvertY(50, "user", "ndc")

## Horizontal line
pushViewport(viewport())
grid.lines(x = c(X1, X2), y = Y, gp = gpar(col = "red"))
popViewport()

最后，这是一个几乎等效的、更普遍有用的方法。它使用了 @mdsumner 的答案中链接到的文章中 Paul Murrell 演示的函数 grid.move.to() 和 grid.line.to()：

library(grid)
library(gridBase)
par(mfrow=c(1,3))

barplot(VADeaths); vps1 <- do.call(vpStack, baseViewports())
barplot(VADeaths) 
barplot(VADeaths); vps3 <- do.call(vpStack, baseViewports())

pushViewport(vps1)
Y <- convertY(unit(50,"native"), "npc")
popViewport(3)

grid.move.to(x = unit(0, "npc"), y = Y, vp = vps1)
grid.line.to(x = unit(1, "npc"), y = Y, vp = vps3, 
             gp = gpar(col = "red"))

Answer 2

这是我能做的最好的事情，而不用更深思熟虑：

par(mfrow = c(1,3),xpd = NA)

for (i in 1:3){
    x <- rnorm(200,i)
    hist(x)
    if (i == 1) segments(par("usr")[1],10,30,10)
}

我不确定如何确保线路在正确的位置结束而不进行修补。在每个区域中绘制一个段可以解决这个问题，但会引入正确排列高度的问题。但这至少是一个很好的起点。

我猜这在 grid 图形中更容易，但我必须做一些研究来验证。

Answer 3

Paul Murrell 的这篇文章展示了使用 grid 图形在两个不同的坐标系之间绘制线条，在这种情况下，线条的端点在两个单独的子图的原始空间中指定：

保罗·穆雷尔。网格图形包。 R 新闻，2(2):14-19，2002 年 6 月

在 PDF 文章的第 17 页：

http://cran.r-project.org/doc/Rnews/Rnews_2002-2.pdf