R中基于cut分箱的直方图为何不显示x轴区间标签

复现代码 #

set.seed(50)
data <- sample(50)
res <- cut(data, breaks = seq(0,50, by= 12), labels = c("<12","12-24","24-35","36-50"))
hist(as.numeric(res))

问题成因 #

• 从代码逻辑拆解，标签不显示的核心是两个环节的属性丢失：
  1. cut()函数输出的res是带标签的因子类型，每个分类值对应预设的区间语义。很多人对as.numeric()处理因子的行为存在误解，这个操作不会把因子代表的区间边界转成实际数值，只会返回每个值在因子水平中的位置序号，转完的结果就是1、2、3、4这类纯序号值，和原本的区间标签完全脱钩
  2. hist()的设计目标是绘制连续数值的分布，它只会接收纯数值输入自动计算刻度，不会关联之前丢失的因子标签，最终x轴只会展示1-4附近自动生成的数值刻度，不会出现你预设的区间文本

可落地的修正方案 #

• 方案1：针对cut生成的分组数据，本质是离散分类频数统计，改用适配因子的barplot()绘制，直接保留原因子的标签属性：

set.seed(50)
data <- sample(50)
res <- cut(data, breaks = seq(0,50, by= 12), labels = c("<12","12-24","24-35","36-50"))
barplot(table(res),
        xlab = "数据区间",
        ylab = "频数",
        main = "分组频数分布")

• 方案2：如果坚持使用hist()绘制原生连续直方图，不需要提前用cut做分组，直接在hist的参数中指定匹配的区间断点，手动绑定自定义标签即可：

set.seed(50)
data <- sample(50)
hist(data,
     breaks = seq(0,50, by = 12),
     xaxt = "n", # 临时隐藏默认x轴刻度
     main = "原始数据分布直方图")
# 在每个区间的中点位置绑定自定义文本标签
axis(1, at = c(6, 18, 30, 43), labels = c("<12","12-24","24-35","36-50"))

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