R中绘制双Y轴根系与地上新梢生长速率对比图及报错修复

错误原因 #

• 括号配对错误：第一行ggplot()的参数未闭合，缺少aes()对应的右括号
• 美学映射位置错误：color = Treatment、group = Treatment属于需要映射到数据列的参数，必须放在aes()内部，放在外部时R会从全局环境找Treatment对象，自然找不到就报错
• 双Y轴映射逻辑错误：次轴仅为刻度标注，第二个折线（侧梢生长）的y值需要先转换为和左轴（根系生长，毫米/天）同单位才能正确渲染，同时需要提前定义转换系数coeff
• geom_point()未指定y映射，无法确定对应的数据列

修正后代码 #

# 定义单位转换系数：1毫米=0.1厘米，左轴毫米转右轴厘米乘以0.1
coeff <- 0.1

ggplot(Growth_rateSeptC18, aes(x = Date)) +
  # 根系生长折线：左轴单位毫米/天，直接用原数值
  geom_line(aes(y = rootlength_root, color = Treatment, group = Treatment)) +
  # 侧梢生长折线：原单位是厘米/天，乘10转成毫米匹配左轴刻度，加线型区分两类指标
  geom_line(aes(y = AvgLoverall * 10, color = Treatment, group = Treatment, linetype = "侧梢生长")) +
  # 次轴设置必须包裹在scale_y_continuous内
  scale_y_continuous(sec.axis = sec_axis(~.*coeff, name="侧梢生长 - 厘米/天")) +
  # 分别为两类指标添加点，用不同形状区分
  geom_point(aes(y = rootlength_root, color = Treatment), position = position_dodge(0.2), size = 1) +
  geom_point(aes(y = AvgLoverall * 10, color = Treatment), position = position_dodge(0.2), size = 1, shape = 2) +
  labs(x="日期", y = "根系生长 - 毫米/天") +
  theme(
    panel.grid.major = element_blank(), 
    panel.grid.minor = element_blank(), 
    panel.background = element_blank(),
    axis.line = element_line(colour = "black") # 可选添加坐标轴边框，更符合科研绘图规范
  )

补充说明 #

如果需要展示3次重复的均值波动，可将折线和点替换为stat_summary()函数自动汇总重复值绘图，避免多条重复线干扰展示效果。

内容的提问来源于stack exchange，提问作者Anne Kearney Logan