如何让ggplot2的geom_density接受数据框列作为y轴频率变量？

我明白你现在的困境：想用类似密度图的平滑曲线展示基因组区间的强度，但数据是非等宽bin，每个区间的高度由独立的强度值（V7）决定，不是计数，之前的代码要么报错要么效果不对。让我一步步帮你解决：

为什么你之前的代码出问题？ #

• 报错原因：geom_density是统计图层，它的y轴是内部计算的密度/统计量（比如..density..），不能直接在aes里指定原始数据的y=V7，这就是你遇到object 'y' not found的原因。
• 效果不佳原因：你之前的代码只用了区间起始位置V5作为x，没考虑区间宽度，导致bin位置偏移；同时没有把V7作为权重，所以峰值高度和实际强度无关。

解决方案：加权密度估计（无需重构等宽bin） #

我们的核心思路是：把每个区间的强度V7作为权重，让密度曲线的峰值高度和强度成正比，同时考虑区间的位置范围。

步骤1：预处理数据（仅计算辅助列，不改变原始结构） #

先算出每个区间的中点和结束位置，方便后续更准确的定位：

library(ggplot2)
library(dplyr) # 用于数据处理

# 计算区间结束位置和中点
locus_df_trim <- locus_df_trim %>%
  mutate(
    x_end = V5 + annot_width,
    x_mid = V5 + annot_width / 2
  )

步骤2：绘制加权平滑密度曲线 #

使用stat_density并指定weight=V7，让密度估计时优先考虑高强度的区间，同时用geom="area"实现填充效果，区分不同细胞类型：

ggplot(locus_df_trim, aes(x = x_mid, weight = V7, fill = V2, colour = V2)) +
  # 加权密度曲线，adjust调整平滑程度，alpha控制透明度
  stat_density(alpha = 0.3, adjust = 0.4, kernel = "gaussian", 
               geom = "area", position = "identity") +
  # 添加你需要的垂直参考线
  geom_vline(xintercept = 204738919, size = 1, colour = "#FF3721", linetype = "dashed") +
  theme_classic()

进阶：更准确的区间覆盖（可选） #

如果想让密度曲线更精准地反映整个区间的覆盖范围，而不是只用中点，可以把每个区间拆分成多个采样点（比如每个区间取10个点），再做加权密度估计：

# 拆分每个区间为多个采样点，保证总权重不变
expanded_df <- lapply(1:nrow(locus_df_trim), function(i) {
  row <- locus_df_trim[i,]
  # 生成区间内的10个均匀分布点
  x_vals <- seq(row$V5, row$x_end, length.out = 10)
  data.frame(
    V2 = row$V2,
    x = x_vals,
    weight = row$V7 / 10  # 拆分权重，每个点的权重是原强度的1/10
  )
}) %>% bind_rows()

# 绘制更精准的加权密度曲线
ggplot(expanded_df, aes(x = x, weight = weight, fill = V2, colour = V2)) +
  stat_density(alpha = 0.3, adjust = 0.4, kernel = "gaussian", 
               geom = "area", position = "identity") +
  geom_vline(xintercept = 204738919, size = 1, colour = "#FF3721", linetype = "dashed") +
  theme_classic()

额外：先看原始区间块（可选） #

如果你想先确认每个区间的位置和高度是否正确，可以先用geom_rect绘制原始区间：

ggplot(locus_df_trim, aes(xmin = V5, xmax = x_end, ymin = 0, ymax = V7, fill = V2)) +
  geom_rect(alpha = 0.3) +
  geom_vline(xintercept = 204738919, size = 1, colour = "#FF3721", linetype = "dashed") +
  theme_classic()

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