geom_histogram特定binwidth下cut_short_scale()报错求助

问题原因 #

这个错误是scales包1.2.1之后版本的已知bug：当ggplot根据特定binwidth计算出的x轴刻度值，恰好落在cut_short_scale()内部的刻度分割边界附近，且涉及非10次幂的特殊刻度点时，函数内部的下标赋值逻辑会出现NA值，进而触发报错。

比如使用binwidth=90时，生成的x轴刻度组合刚好触发了这个逻辑漏洞；而binwidth=80或100时，刻度值的分布避开了这个边界场景，因此不会报错。

解决办法 #

1. 回滚到scales 1.2.1版本 #

这是已验证有效的临时解决方案，可通过devtools安装指定版本：

# 若未安装devtools先执行：install.packages("devtools")
devtools::install_version("scales", version = "1.2.1")

2. 使用自定义刻度标签函数替代cut_short_scale() #

手动实现类似的短刻度格式化逻辑，绕过官方函数的bug：

library(tidyverse)
library(AmesHousing)

# 自定义短刻度标签函数
my_short_label <- function(x) {
  case_when(
    x >= 1000 ~ str_c(round(x/1000, 1), "k"),
    x >= 1 ~ as.character(round(x)),
    TRUE ~ as.character(x)
  )
}

# 绘图调用自定义函数
ames_raw |> 
    mutate(`Bsmt Unf SF` = replace_na(`Bsmt Unf SF`, 0)) |>
    ggplot(aes(x = `Bsmt Unf SF`)) +
    geom_histogram(binwidth = 90) +
    scale_x_continuous(label = my_short_label)

3. 等待scales包官方修复 #

该bug属于包内部逻辑漏洞，后续版本大概率会修复，可关注scales包的更新日志，待修复后升级到最新版本即可。

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