如何为ggstats::gglikert分面李克特图添加频数标签

报错原因 #

annotate()是全局注释函数，不会识别分面的分组逻辑。当你用facet_rows分面时，ggplot会为每个分面重复加载所有注释数据——比如你有2个分面、30行数据，ggplot就会要求注释数据的长度是30×2=60，但你只传入了30行，自然触发美学属性长度不匹配的报错。

可行方案：用geom_text替代annotate #

不要用annotate，改用geom_text，结合gglikert的底层计算数据，同时同步exclude_fill_values参数，保证标签位置和条形完全对齐。

具体操作步骤 #

• 提取绘图底层数据：先运行gglikert生成基础图，再用ggplot_build()拿到已经计算好条形位置的数据集，里面包含每个条形的x/y坐标、频数等信息。
• 同步position参数：在geom_text的position中传入和gglikert一致的position_likert(exclude_fill_values = 你的排除值)，避免标签位置错乱。
• 映射频数标签：把频数字段映射到label美学，调整vjust/size让标签显示更美观。

示例代码 #

假设你的数据是df，分组变量group，李克特题项item，响应response，排除的层级是c("NA", "不适用")：

library(ggstats)
library(ggplot2)

# 基础likert分面图
base_plot <- gglikert(df, item ~ response,
                      facet_rows = vars(group),
                      exclude_fill_values = c("NA", "不适用"))

# 提取绘图用的底层数据（包含已计算的位置）
plot_data <- ggplot_build(base_plot)$data[[1]]

# 关联频数（这里plot_data的count列就是实际频数，若需自定义可关联原始数据）
plot_data$freq_label <- plot_data$count

# 添加频数标签
base_plot + 
  geom_text(aes(label = freq_label),
            position = position_likert(exclude_fill_values = c("NA", "不适用")),
            vjust = 0.5, size = 3)

手动计算频数的情况 #

如果底层数据的count不是你需要的频数，可先用dplyr统计：

library(dplyr)
freq_df <- df %>%
  filter(!response %in% c("NA", "不适用")) %>%
  count(group, item, response, name = "freq")

# 把统计结果和plot_data关联（需匹配group、item、response字段）
plot_data <- plot_data %>%
  left_join(freq_df, by = c("group", "item", "response"))

之后再把freq映射到geom_text的label即可。

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