如何编程实现ggplot2的geom_density触达X轴（含自由刻度分面）

ggplot2默认的geom_density不会像base R的plot(density())那样将密度曲线延伸到密度为0的X轴边界，即使使用outline.type参数也无法解决。以下是两种可靠的实现方法，支持单图和分面（包括scales="free"）场景：

方法一：手动计算密度范围并设置X轴限制 #

适用于快速调整现有ggplot代码，无需重构数据：

单图场景 #

先计算目标变量的密度对象，提取其X轴全范围（包含密度为0的边界），再传递给xlim()：

library(ggplot2)

# 计算密度并获取X轴范围
dens <- density(diamonds$carat, adjust = 5)
x_full_range <- range(dens$x)

# 绘制密度图并设置X轴范围
ggplot(diamonds, aes(carat)) +
  geom_density(adjust = 5, outline.type = "full") +
  xlim(x_full_range)

分面（scales="free"）场景 #

为每个分面分组单独计算密度范围，结合ggh4x包的facetted_pos_scales为每个分面设置独立X轴限制：

library(ggplot2)
library(dplyr)
library(ggh4x)

# 计算每个分面分组的密度X范围
density_ranges <- diamonds %>%
  group_by(cut) %>%
  summarise(
    x_min = min(density(carat, adjust = 5)$x),
    x_max = max(density(carat, adjust = 5)$x)
  )

# 绘制分面密度图
ggplot(diamonds, aes(carat)) +
  geom_density(adjust = 5) +
  facet_wrap(~cut, scales = "free_x") +
  # 为每个分面设置对应X轴范围
  facetted_pos_scales(
    x = lapply(split(density_ranges, density_ranges$cut), function(group) {
      scale_x_continuous(limits = c(group$x_min, group$x_max))
    })
  )

方法二：基于base R密度数据重构绘图数据 #

完全自定义密度曲线的起点和终点，确保100%匹配base R的效果，适合需要精细控制的场景（如Shiny应用）：

单图场景 #

直接使用base R生成的密度数据，手动添加X轴两端密度为0的点：

library(ggplot2)

# 生成含边界点的密度数据
dens <- density(diamonds$carat, adjust = 5)
dens_data <- data.frame(x = dens$x, y = dens$y) %>%
  rbind(
    data.frame(x = min(dens$x), y = 0),
    data.frame(x = max(dens$x), y = 0)
  ) %>%
  arrange(x)

# 绘制密度曲线和填充区域
ggplot(dens_data, aes(x, y)) +
  geom_line(color = "black") +
  geom_area(fill = "gray70", alpha = 0.5)

分面（scales="free"）场景 #

为每个分面分组生成完整的密度数据（含边界点），再绘图：

library(ggplot2)
library(dplyr)

# 生成每个分面分组的完整密度数据
dens_facet_data <- diamonds %>%
  group_by(cut) %>%
  group_modify(~{
    dens <- density(.x$carat, adjust = 5)
    data.frame(x = dens$x, y = dens$y) %>%
      rbind(
        data.frame(x = min(dens$x), y = 0),
        data.frame(x = max(dens$x), y = 0)
      ) %>%
      arrange(x)
  })

# 绘制分面密度图
ggplot(dens_facet_data, aes(x, y)) +
  geom_line(color = "black") +
  geom_area(fill = "gray70", alpha = 0.5) +
  facet_wrap(~cut, scales = "free_x")

关键说明 #

• 两种方法都确保X轴范围覆盖密度曲线降至0的边界，解决了outline.type无效的问题；
• 方法二更适合Shiny应用，因为可以提前计算好数据，避免绘图时的动态计算延迟；
• 分面场景下scales="free_x"会自动适配每个分组的X轴范围，无需额外调整。

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