R绘图添加y轴最大值20padding报错及ggplot2实现咨询

错误原因说明 #

Error in which.max(density(y3)) : 'list' object cannot be coerced to type 'double'报错触发原因是：density()函数返回结果为列表结构，包含x轴取值、y轴密度值、带宽等多个属性，直接传入which.max()会触发类型匹配错误，需要显式提取密度值字段$y后再计算最大值。
后续出现的y轴上限过高/最大值低估问题，是因为未统一计算5条密度曲线的全局最大密度值，仅基于单条曲线设置范围导致的。

原生R绘图修复代码 #

# 1. 统一计算5组数据密度曲线的全局最大y值
max_density_y <- max(
  density(y1)$y,
  density(y2)$y,
  density(y3)$y,
  density(y4)$y,
  density(y5)$y
)
# 2. 上边界增加20的padding
ylim_upper <- max_density_y + 20

# 3. 绘图
plot(density(y1), ylim = c(0, ylim_upper), main = "多组分布密度曲线", col = 1)
lines(density(y2), col = 2)
lines(density(y3), col = 3)
lines(density(y4), col = 4)
lines(density(y5), col = 5)
# 可按需补充图例、坐标轴标签等配置
legend("topright", legend = paste0("分组",1:5), col = 1:5, lty = 1)

ggplot2版本实现代码 #

ggplot2基于长表数据绘图时会自动计算所有分组的密度值范围，无需手动逐个计算最大值：

library(ggplot2)
library(tidyr)

# 假设原始宽表数据框名为raw_df，5列分别对应5组数据y1-y5
# 转换为ggplot2适配的长表格式
df_long <- pivot_longer(raw_df, cols = everything(), names_to = "分组", values_to = "观测值")

ggplot(df_long, aes(x = 观测值, color = 分组)) +
  geom_density(linewidth = 1) +
  # 仅在y轴上边界添加20的padding，下边界不额外扩展
  scale_y_continuous(expand = expansion(add = c(0, 20))) +
  labs(title = "多组分布密度曲线", x = "观测值", y = "密度") +
  theme_bw()

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