如何让ggplot散点图坐标轴与geom_rect元素精准对齐？

错误原因 #

设置limits = c(0, 8000)报错的核心原因是：log2(0)的计算结果为负无穷，并非有限数值，ggplot无法基于该值生成坐标轴刻度，因此抛出'from' must be a finite number错误。

解决方案 #

1. 替换坐标轴下限为极小正数 #

将0替换为接近0的小数值（比如1，现有数据中最小的stroma为8.27、epi为9.72，1不会截断任何数据），此时log2(1)=0，属于有限值，可正常生成刻度。

2. 让矩形匹配数据实际最大范围 #

不要硬编码8000、4000这类固定值，而是从所有关联数据框中计算x、y轴的实际最大值，让矩形边界与数据范围完全对齐，消除空白。

3. 优化代码（可选） #

将矩形的参数整理为单独数据框，减少重复代码，提升可读性。

修改后的完整代码 #

# 定义截断值
epi_cutoff <- 200
stroma_cutoff <- 200

# 计算所有数据的x、y最大值
max_stroma <- max(c(final_df$stroma, crc_vals$stroma, crc_met_vals$stroma))
max_epi <- max(c(final_df$epi, crc_vals$epi, crc_met_vals$epi))

# 绘制图表
ggplot(final_df, aes(x=as.numeric(stroma), y=as.numeric(epi), 
                     label = ID)) +
  # 基于计算出的范围定义矩形，下限用1替代0
  geom_rect(aes(xmin=1, xmax=stroma_cutoff, ymin=1, ymax = epi_cutoff),
            alpha = 0.005, fill = "blue2") +
  geom_rect(aes(xmin=stroma_cutoff, xmax=max_stroma, ymin=1, ymax = epi_cutoff),
            alpha = 0.005, fill = "yellow1") +
  geom_rect(aes(xmin=1, xmax=max_stroma, ymin=epi_cutoff, ymax = max_epi),
            alpha = 0.005, fill = "red1") +
  geom_point(size = 2) +
  geom_point(data = crc_vals,
             aes(x=as.numeric(stroma), y=as.numeric(epi)),
             color = 'green',
             size = 3) +
  geom_point(data = crc_met_vals,
             aes(x=as.numeric(stroma), y=as.numeric(epi)),
             color = 'purple',
             size = 3) +
  # 设置坐标轴范围，expand参数消除边缘空白
  scale_x_continuous(trans='log2', limits = c(1, max_stroma), expand = c(0,0)) +
  scale_y_continuous(trans='log2', limits = c(1, max_epi), expand = c(0,0)) +
  labs(y='Feature1',x='Feature2',
       title = 'Chart\nblue2 0.005 - yellow1 0.005 - red1 0.005') +
  theme_bw()

关键细节说明 #

• expand = c(0,0)：移除坐标轴两端的默认空白，让矩形与图表边缘完全对齐。
• 动态计算最大值：确保矩形范围始终匹配数据的实际边界，后续数据更新无需手动调整硬编码数值。
• 下限用1替代0：既规避了log2(0)的无穷值问题，又不会截断现有有效数据。

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