如何在使用 lattice 制作的 levelplot 中添加线条（abline 不知何故不起作用）？

Question

我想在我的水平图上绘制水平线和垂直线，对应于从 74 到 76 的 x 值和从 28 到 32 的 y 值。下面是我的 R 代码。但是当我运行以下命令时，我得到了水平图但没有线条。我也没有收到来自 R 的错误。我安装的默认主题是将值映射到粉红色和青色的东西。我也尝试过使用面板功能，但也没有运气。

levelplot(d_fire_count_nom ~ longitude + latitude | factor(day)+factor(year), 
          data = asia,
          subset = (month == 10),  aspect="iso", contour = FALSE, layout=c(1,1), 
          main="If a fire occured in a region (low confidence) in October during 2001-2008", 
          scales=list(x=list(at=seq(from=60,to=98, by=1)), 
                      y=list(at=seq(from=5,to=38,by=1)),cex=.7, alternating=3), 
          xlim=c(60, 98), ylim=c(5, 38),
          abline=list(h=74:76, v=28:32, col="grey"))

Answer 1

lattice 图形不是这样工作的。事实上，如果您阅读 ?levelplot，您会发现名为 abline 的函数没有参数，所以我不确定您从哪里得到该语法。

您可以通过更改 panel 函数向 lattice 图形添加内容。有许多panel.* 函数可以做各种事情，比如绘制点、线、散点图平滑器等。在这种情况下，我们想使用panel.abline。所以我们定义了自己的panel函数。

这使用了?levelplot的第一个例子：

x <- seq(pi/4, 5 * pi, length.out = 100)
y <- seq(pi/4, 5 * pi, length.out = 100)
r <- as.vector(sqrt(outer(x^2, y^2, "+")))
grid <- expand.grid(x=x, y=y)
grid$z <- cos(r^2) * exp(-r/(pi^3))
levelplot(z~x*y, grid,
        panel = function(...){
            panel.levelplot(...)
            panel.abline(h = 2.5)
            panel.abline(v = 2.5)
        }, 
        cuts = 50, scales=list(log="e"), xlab="",
        ylab="", main="Weird Function", sub="with log scales",
        colorkey = FALSE, region = TRUE)

我们的新面板函数需要先绘制水平图，所以我们先调用panel.levelplot。然后我们想添加一些行，所以我们为此添加了panel.abline。