我想创建一个 ggplot,其中统计参数根据美学映射因素而变化。具体来说,我想使用 stat_density2d() 创建一个等高线图,我想将离散因子映射到颜色,并且我想为每个因子级别指定不同的中断值。

这是一个最小的工作示例:

d <- data.frame(x=c(rnorm(500), rnorm(500, 2)), y=rnorm(1000), z=c(rep("a", 500), rep("b", 500)))

ggplot(d, aes(x, y, fill=z)) +
      stat_density2d(breaks=.05, geom="polygon", alpha=.5)

这就是我要做的,除了不同因素的休息时间是相同的。我正在寻找一种方法来为每个因素指定不同的中断值。一种方法是为每个因素创建一个单独的层:
ggplot() +
      stat_density2d(data=d[d$z=="a",], aes(x, y), breaks=.05, geom="polygon", alpha=.5, fill="red") +
      stat_density2d(data=d[d$z=="b",], aes(x, y), breaks=.1, geom="polygon", alpha=.5, fill="blue")

但这行不通,因为我丢失了图例,而且对于具有两个以上因子水平的情况来说太麻烦了。

我想结合以上两种情况的优点。我几乎肯定我已经看到了一种完成我正在尝试做的事情的方法,但现在似乎无法找到它,因为它是相关的。有人对可能的解决方案有想法吗?

最佳答案

几年前就记得这个 - 所以直接从 Kohske's answer 获取解决方案

# Data
set.seed(1)
d <- data.frame(x=c(rnorm(500), rnorm(500, 2)), y=rnorm(1000),
                                   z=c(rep("a", 500), rep("b", 500)))

library(plyr)
library(ggplot2)

hls <- mapply(function(x, b) stat_density2d(data = x, breaks = b,
                    geom="polygon", alpha=.5), dlply(d, .(z)), c(0.05, 0.1))

ggplot(d, aes(x, y, fill=z)) + hls + scale_fill_manual(values=c("red", "blue"))
# Compare
ggplot() +
  stat_density2d(data=d[d$z=="a",], aes(x, y), breaks=.05, geom="polygon",
                                                      alpha=.5, fill="red") +
  stat_density2d(data=d[d$z=="b",], aes(x, y), breaks=.1, geom="polygon",
                                                     alpha=.5, fill="blue")

关于r - ggplot2:非美学参数会因因素而异吗?,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/27327198/

10-12 17:39