r - 如何使用 ggplot2 将点和线叠加到等值线图上？

Question

我想用我想要突出显示的特定点(这些点存储在不同的数据集中)来注释等高线图。当我尝试时，出现错误:

Error: Aesthetics must either be length one, or the same length as the dataProblems:z

但是，当我尝试制作一个可重现的示例时，我收到了不同的错误:

Error in eval(expr, envir, enclos) : object 'z' not found

可重现示例的代码如下:

library(mnormt)
library(dplyr)
library(ggplot2)

f <- function(x, y) {
    dmnorm(x = c(x, y),
            mean = c(0, 0),
            varcov = diag(2))
}
f <- Vectorize(f)


xmesh <- seq(from = -3, to = 3, length.out = 100)
ymesh <- seq(from = -3, to = 3, length.out = 100)
dummy <- expand.grid(x = xmesh, y = ymesh)
dummy$z <- f(dummy$x, dummy$y)

stuff <- data_frame(x = c(0, 0, 1),
                    y = c(0, -1, -1),
                    point = c("O", "P", "Q"))

dummy %>%
    ggplot(aes(x = x, y = y, z = z)) +
    stat_contour(aes(color = ..level..)) +
    labs(color = "density") + 
    geom_point(data = stuff, mapping = aes(x = x, y = y, color = point))

Answer 1

ggplot 将第一次 ggplot 调用中的 aes 传递给其余的 geoms，除非另有说明。因此，错误告诉您它无法在 stuff 中找到 z，并且它仍然认为 z 应该是初始调用中的 z。

有多种方法可以解决这个问题，我认为最简单的解决方法是分别为每个几何对象提供其数据:

ggplot() +
  stat_contour(data = dummy, aes(x = x, y = y, z = z, color = ..level..)) +
  labs(color = "density") + 
  geom_point(data = stuff, aes(x = x, y = y, fill = factor(point)), pch = 21)

注意。你也有一个问题 colour cannot be mapped in two different geoms ，所以我使用 pch 和 fill 修复了它。