在ggplot日期轴图表中按网格坐标绘制箭头的技术问询

问题描述 #

绘制患者血液检测结果的时间序列图时，已完成基础绘图和参考范围阴影区域添加，但使用annotate()或geom_segment()时，因x轴为日期（YYYY-MM-DD）类型，直接用数字指定坐标会报错，希望能忽略轴单位像网格一样指定箭头坐标，最终目标是用箭头指向特定数据点标记治疗起始时间。

测试代码运行报错：

Error: Invalid input: date_trans works with objects of class Date only

报错原因 #

x轴使用了scale_x_date()进行日期类型转换，而annotate("segment")中直接传入x=1这类数字，并非Date类对象，导致类型不匹配触发报错。

解决方案 #

方案1：使用相对坐标（NPC单位） #

直接基于绘图区域的相对位置指定坐标，无需关注轴的实际单位。用unit()函数指定"npc"单位，取值范围为0到1：

• 0代表绘图区域最左侧/最底部
• 1代表绘图区域最右侧/最顶部

示例代码：

annotate("segment", 
         x = unit(0.1, "npc"), y = unit(0.2, "npc"), 
         xend = unit(0.3, "npc"), yend = unit(0.8, "npc"), 
         arrow = arrow(length = unit(0.15, "cm")))

方案2：直接指向特定数据点（更贴合业务需求） #

如果要标记治疗起始，直接提取目标数据点的日期和检测值作为坐标，无需手动估算位置。

示例步骤：

1. 提取目标数据点（以第3个数据点为例）：

target_point <- q[3, ]

2. 添加指向该点的箭头（可调整起始位置的偏移，比如提前15天）：

annotate("segment", 
         x = target_point$date - days(15), y = 3.5, 
         xend = target_point$date, yend = target_point$result, 
         arrow = arrow(length = unit(0.15, "cm")))

修改后的完整测试代码 #

library(ggplot2)
library(data.table)
library(lubridate)
library(ggprism)

result <- runif(25, min = 2.0, max = 3.5)
start_date <- ymd("2021-08-16")
end_date <- ymd("2022-10-29")
date <- sample(seq(start_date, end_date, by = "days"), 25, replace = TRUE)
q <- data.table(result, date) # 修正原代码中未定义的numbers

# 提取目标数据点（示例为第3个）
target_point <- q[3, ]

ggplot(q, aes(x = date, y = result)) + 
  geom_line() + 
  geom_point(aes(x = date, y = result), shape = 21, size = 3) +
  scale_x_date(limits = c(min(q$date), max(q$date)), 
               breaks = date_breaks("1 month"),
               labels = date_format("%b %Y")) +
  ylab("Corrected calcium (mmol/L)")+ # 修正原代码未闭合的括号
  xlab("Date of blood test")  +
  ylim(1,4)+ 
  geom_ribbon(aes(ymin=2.1, ymax=2.6), fill="grey", alpha=0.2, colour="grey")+ 
  geom_vline(xintercept=as.numeric(q$date[c(3, 2)]),
             linetype=4, colour="black") +
  theme(axis.text.x = element_text(angle = 45)) + theme_prism(base_size = 10) +
  # 使用方案2：指向目标数据点的箭头
  annotate("segment", 
           x = target_point$date - days(15), y = 3.5, 
           xend = target_point$date, yend = target_point$result, 
           arrow = arrow(length = unit(0.15, "cm")))

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