ggplot2绘制时间序列曲线：0-50毫秒区间阴影添加报错求助

核心问题点 #

• 语法拼接错误：scale_y_continuous(limits = c(0, max(data$force)))这一行末尾缺失+，导致后续的geom_point、theme_classic()等图层无法和前面的ggplot对象正确拼接，这是触发「无法合并ggproto对象」错误的直接原因。
• 废弃参数误用：geom_area里的geom_params是ggplot2旧版本参数，现已废弃，直接将填充、透明度等参数放在geom_area函数括号内即可。
• 区间筛选逻辑错误：原代码用ifelse处理x轴的方式会导致图形异常，应该直接针对0-50毫秒的数据集子集绘制阴影区域。

修正后的完整代码 #

# 原始数据
milliseconds <- c(0,50,100,150,200,250,300)
force <- c(20,75,120,260,400,500,600)
data <- data.frame(miliseconds = milliseconds, force)  # 统一列名拼写（与原代码保持一致）

# 绘制图形
RFD <- ggplot(data = data, mapping = aes(x = miliseconds, y = force)) +
  geom_line(color="black", size=2) +
  # 针对0-50毫秒的子集绘制阴影区域
  geom_area(data = subset(data, miliseconds >= 0 & miliseconds <= 50),
            fill = "red", alpha = 0.5) +
  scale_y_continuous(limits = c(0, max(data$force))) +  # 补上缺失的+号
  geom_point(size=4, color="black") +
  theme_classic() +
  ggtitle("Rate of Force Development")  # 修正标题拼写错误（Developement→Development）

# 补充x轴和y轴设置
RFD + scale_x_continuous(name = "Milliseconds",
                         breaks = seq(0, 300, 50),
                         limits=c(0, 300)) +
  scale_y_continuous(name = "Force")

额外说明 #

• 修正后的代码通过subset筛选0-50毫秒的数据子集，geom_area会自动连接该区间的所有数据点并填充曲线下方区域，阴影贴合度更准确。
• 标题拼写错误已修正（不影响代码运行，但更符合英文规范）。

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