ggplot2绘制多基因表达曲线及添加均值线报错问题求助

咱们先拆解你遇到的两个核心问题：颜色映射报错和均值线位置偏离，一步步来搞定：

1. 修复"Discrete value supplied to continuous scale"错误 #

你之前把color="red"放在aes()内部了，这会让ggplot把"red"当成一个离散分类变量，去匹配你原本的基因颜色尺度（离散型），但咱们是想直接指定颜色，不是映射变量，所以要把color移到aes()外面。另外，均值线是单条线，需要用group=1告诉ggplot把所有时间点的均值归为同一组来连线：

# 修正后的均值线图层写法
plotC + geom_line(data=means, aes(x=variable, y=value, group=1), color="red", size=0.75)

2. 解决均值线位置偏离的问题 #

看你提供的means数据，数值明显不对：比如时间点5的三个基因表达量是-0.61、-6.26、-3.77，计算均值应该是(-0.61-6.26-3.77)/3 ≈ -3.54，但你的means里对应的value是-1.5，这说明你手动计算均值时出错了。

推荐直接用dplyr从原始数据按时间点自动计算均值，避免手动失误：

library(dplyr)
# 从df2按时间点分组计算均值
means_correct <- df2 %>%
  group_by(variable) %>%
  summarise(value = mean(value), .groups = "drop")

得到的正确均值数据应该是：

用这个数据绘制的均值线，就会出现在所有基因曲线的中间位置了。

完整可运行代码 #

把所有步骤整合起来，直接就能跑出正确的图：

library(ggplot2)
library(dplyr)

# 你的原始数据
df2 <- data.frame(
  variable = c(5,5,5,6,6,6,7,7,7),
  value = c(-0.610384193, -6.25967087, -3.773389731,
            -0.358879035, -6.066341017, -4.202998579,
            -0.103885903, -6.648844687, -5.041554127),
  gene = factor(c(1,2,3,1,2,3,1,2,3))
)

# 计算正确的时间点均值
means_correct <- df2 %>%
  group_by(variable) %>%
  summarise(value = mean(value), .groups = "drop")

# 绘制基因表达曲线 + 红色均值线
plotC <- ggplot(df2, aes(x=variable, y=value, group=gene, color=gene)) +
  geom_line(size=0.5, alpha=0.4) +
  # 添加均值线：color在aes外指定红色，group=1确保连线
  geom_line(data=means_correct, aes(x=variable, y=value, group=1), 
            color="red", size=0.75) +
  labs(x="Time Point", y="Normalized Expression", color="Gene") +
  theme_bw()

print(plotC)

额外小贴士 #

• 如果想让均值线更突出，可以加linetype="dashed"改成虚线，或者调整size/alpha参数
• 你原始means数据里的gene列（值为50）完全多余，均值线不需要按基因分组，计算时直接按variable分组即可

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