如何在共享ggplot2配置的特定图表中移除scale_fill_*相关的NA分组？

背景 #

数据集的分组变量存在NA值，ggplot2默认会将NA视为一个分组，但当前图表不需要这个效果。所有图表使用统一的标准化图层配置（保存为空白gg对象），如果在全局配置中移除NA会丢失其他图表的有效数据，不想单独复制数据集，寻求ggplot2原生的最优方案。

可复现示例 #

library(tidyverse)

gg <- airquality %>%
  slice_head(n = 150) %>%
  mutate(group = as.character(sample(c(1:2, NA), 150, replace = TRUE))) %>%
  ggplot() +
  theme_bw() +
  scale_fill_brewer(palette = "Paired")

# 生成包含NA分组的密度图
gg + geom_density(aes(Solar.R, fill = group), alpha = .7)

已尝试的无效方法 #

• 在scale层添加na.translate = FALSE：仅移除图例中的NA分组，但NA对应的密度曲线仍存在
• 为scale_fill设置na.value = NA：无明显效果，且该scale不支持NULL或"blank"值
• 在美学映射中使用na.rm = TRUE：geom_density()不支持此参数

解决方案 #

以下方法均不会影响全局的gg配置，仅针对当前特定图表生效：

方法1：在geom层的data参数中过滤NA #

利用tidyverse的公式语法，直接在当前geom的数据中剔除分组为NA的观测：

gg + geom_density(aes(Solar.R, fill = group), alpha = .7,
                  data = ~ filter(.x, !is.na(group)))

方法2：使用stat_density的subset参数 #

geom_density底层依赖stat_density，而stat支持subset参数，直接指定保留非NA的分组：

gg + geom_density(aes(Solar.R, fill = group), alpha = .7,
                  subset = !is.na(group))

方法3：在美学映射中过滤分组变量 #

在fill的映射里直接将NA值转为无效的映射值，配合scale的na.translate = FALSE彻底移除NA相关元素：

gg + geom_density(aes(Solar.R, fill = ifelse(is.na(group), NA, group)), alpha = .7) +
  scale_fill_brewer(palette = "Paired", na.translate = FALSE)

这三种方法都不需要单独复制或修改原始数据集，仅对当前调用的图层生效，完美适配共享全局配置的场景。

内容的提问来源于stack exchange，提问作者philsf