如何在ggplot2的facet_wrap()中控制各分面的坐标轴刻度数量？

我在使用facet_wrap绘制多分组核密度估计图时，遇到了不同分组方差差异过大的问题：方差极小的分组对应的分面板x轴标签会堆叠在一起，无法正常阅读。最小复现代码如下：

library(tidyverse)
x1 <- rnorm(1e4)
x2 <- rnorm(1e4,mean=2,sd=0.00001)

data.frame(x=c(x1,x2),group=c(rep("1",length(x1)),rep("2",length(x2)))) %>%
  ggplot(.) + geom_density(aes(x=x)) + facet_wrap(~group,scales="free")

最直接的解决思路是增大画布尺寸让所有标签完全显示，但我的分面数量很多，这个方案不实用。我希望能控制每个分面x轴的刻度数量，比如所有分面的x轴只展示2个刻度，请问有没有办法实现这个需求？

已尝试的方案 #

1. 直接添加+ scale_x_continuous(n.breaks = 2)，看似符合需求，但实际运行后没有达到预期效果：
   
2. 参考自定义breaks函数的方案修改代码后，确实得到了2个刻度，但刻度值的小数位数过多，不符合使用要求：

equal_breaks <- function(n = 3, s = 0.5, ...){
  function(x){
    # 重缩放
    d <- s * diff(range(x)) / (1+2*s)
    seq(min(x)+d, max(x)-d, length=n)
  }
}

data.frame(x=c(x1,x2),group=c(rep("1",length(x1)),rep("2",length(x2)))) %>%
  ggplot(.) + geom_density(aes(x=x)) + facet_wrap(~group,scales="free")  + scale_x_continuous(breaks=equal_breaks(n=3, s=0.05), expand = c(0.05, 0))

你可以在自定义breaks函数的基础上，增加有效数字保留的逻辑，即可同时满足「固定刻度数量」和「刻度值格式整洁」的需求。
以下是修改后的可运行代码：

library(tidyverse)
# 自定义breaks函数，新增有效数字控制参数digits
custom_breaks <- function(n = 2, s = 0.05, digits = 3, ...){
  function(x){
    d <- s * diff(range(x)) / (1 + 2 * s)
    breaks <- seq(min(x) + d, max(x) - d, length = n)
    # 保留指定有效数字，避免小数位过长
    signif(breaks, digits = digits)
  }
}

# 绘图代码
x1 <- rnorm(1e4)
x2 <- rnorm(1e4,mean=2,sd=0.00001)
data.frame(x=c(x1,x2),group=c(rep("1",length(x1)),rep("2",length(x2)))) %>%
  ggplot(.) + 
  geom_density(aes(x=x)) + 
  facet_wrap(~group,scales="free") + 
  scale_x_continuous(
    breaks = custom_breaks(n = 2, s = 0.05, digits = 3),
    expand = c(0.05, 0)
  )

如果需要更灵活的刻度格式控制，还可以搭配scales包的标签函数使用，比如用scales::label_number(accuracy = 0.001)统一控制小数位数，或者用scales::label_auto()自动适配不同分面的数值范围。

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