如何在ggplot分面图中强制全零数据的Y轴从0开始？

问题背景 #

使用ggplot的facet_wrap绘制多站点时序图时，希望所有分面的Y轴从0开始，但当某站点数据全为0时，scale_y_continuous(limits = c(0,NA))会导致Y轴出现负数；expand_limits(y=0)无效，而设置scales='free_x'虽解决全零问题，但会挤压其他站点的X轴显示。

解决方案 #

方法1：使用ggh4x包为每个分面单独设置Y轴范围 #

通过手动计算每个站点的Y轴上下限，再利用ggh4x的facetted_pos_scales为每个分面配置独立的Y轴，既保证全零站点Y轴从0开始且无负数，又让其他站点Y轴自由缩放。

library(tidyverse)
library(ggh4x)

# 计算每个站点的Y轴范围：全零站点设极小上限，其他站点用数据最大值
y_limits <- dat %>%
  group_by(shortname) %>%
  summarise(
    y_min = 0,
    y_max = ifelse(all(result == 0), 0.001, max(result, na.rm = TRUE))
  ) %>%
  split(.$shortname) %>%
  map(~c(.$y_min, .$y_max))

# 绘制图形
ggplot(dat) +
  geom_point(aes(x=date, y=result, color=param_name), pch=19) +
  theme_bw() +
  facet_wrap(facets = ~shortname, scales = 'free') +
  scale_x_date(date_labels = '%y') +
  # 为每个分面设置专属Y轴范围
  facetted_pos_scales(
    y = lapply(names(y_limits), function(nm) {
      scale_y_continuous(limits = y_limits[[nm]], expand = c(0, 0))
    })
  ) +
  theme(legend.position = 'none') +
  xlab('Year') +
  ylab(expression(paste('Total Phosphorus (',mu,'g/l)')))

方法2：修改数据添加极小虚拟值（无需额外包） #

给全零站点的数据集添加一个极小的正值（仅用于确定Y轴范围，不实际显示），配合expand = c(0,0)避免Y轴自动扩展到负数。

library(tidyverse)

# 修改数据：为全零站点添加极小值用于Y轴范围计算
dat_mod <- dat %>%
  group_by(shortname) %>%
  mutate(
    result_mod = ifelse(all(result == 0), 0.001, result)
  ) %>%
  ungroup()

# 绘图：用修改后的字段确定Y轴，绘制原始数据
ggplot(dat_mod) +
  geom_point(aes(x=date, y=result, color=param_name), pch=19) +
  scale_y_continuous(limits = c(0, NA), expand = c(0, 0)) +
  theme_bw() +
  facet_wrap(facets = ~shortname, scales = 'free') +
  scale_x_date(date_labels = '%y') +
  theme(legend.position = 'none') +
  xlab('Year') +
  ylab(expression(paste('Total Phosphorus (',mu,'g/l)')))

原理说明 #

• 原方法中scale_y_continuous(limits = c(0,NA))在全零分面时，因数据范围为0，ggplot默认的expand参数会自动扩展Y轴范围，导致出现负数；
• 两种方法均通过为全零站点设置明确的Y轴上限（极小值），配合expand = c(0,0)阻止自动扩展，同时保留其他站点Y轴的自由缩放能力。

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