如何用ggplot2绘制分物种的花瓣长度分箱频率分面图？

你遇到的核心问题是：cut()函数把连续的花瓣长度转换成了因子型的分箱标签，所以ggplot默认将X轴识别为离散变量，没法呈现连续的数值刻度。这里有两种简洁的解决方案，帮你实现连续X轴的分面分箱频率图：

方法1：直接用ggplot原生直方图函数（推荐） #

不需要手动计算频数，geom_histogram可以直接指定分箱区间，自动生成连续X轴的直方图，同时结合分面实现按物种拆分：

library(ggplot2)

# 定义分箱区间
br <- seq(1, 6, by = 0.4)

ggplot(iris, aes(x = Petal.Length)) +
  # 指定分箱断点，设置柱子样式
  geom_histogram(breaks = br, fill = "#2c3e50", color = "white") +
  # 按物种分面
  facet_wrap(~Species, ncol = 1) +
  # 自定义轴标签
  labs(x = "花瓣长度 (cm)", y = "频数") +
  # 让X轴刻度和分箱起点对齐
  scale_x_continuous(breaks = br)

这个方法的优势是代码更简洁，不需要额外处理数据，X轴天然是连续型的，完全符合你的需求。

方法2：改造现有数据流程（适合需要手动控制频数的场景） #

如果一定要保留原有的table计算频数的逻辑，你需要把因子型的分箱标签转换成连续的分箱中点，作为X轴的数值：

library(ggplot2)
library(dplyr)
library(stringr)

br <- seq(1, 6, by = 0.4)

# 处理数据：分箱、计算频数、提取分箱中点
df <- iris %>%
  # 生成分箱标签
  mutate(pl_bin = cut(Petal.Length, br, include.lowest = TRUE)) %>%
  # 按物种和分箱统计频数
  count(Species, pl_bin, name = "Freq") %>%
  # 从分箱标签中提取上下限，计算中点
  mutate(
    # 拆分区间字符串，比如"(1,1.4]"拆成c("1", "1.4")
    bin_limits = str_extract_all(pl_bin, "\\d+\\.*\\d*", simplify = TRUE),
    pl_mid = (as.numeric(bin_limits[,1]) + as.numeric(bin_limits[,2])) / 2
  )

# 绘图：用分箱中点作为连续X轴
ggplot(df, aes(x = pl_mid, y = Freq)) +
  geom_col(width = 0.4, fill = "#2c3e50") + # width匹配分箱步长
  facet_wrap(~Species) +
  labs(x = "花瓣长度 (cm)", y = "频数") +
  scale_x_continuous(breaks = br)

关键说明： #

• 两种方法都能实现连续型X轴，方法1更高效，适合大多数场景；
• 如果需要调整分箱的样式（比如柱子宽度、颜色），两种方法都可以通过geom_histogram或geom_col的参数实现；
• scale_x_continuous(breaks = br)可以让X轴刻度和你定义的分箱起点完全对齐，更直观。

内容的提问来源于stack exchange，提问作者P.Chakytei