R中gglikert绘制偶数级Likert图的排序与总计显示问题

1. 修正满意度占比计算逻辑 #

你怀疑的sum(x %in% 4:5)/length(x[!is.na(x)])确实是问题核心——当前代码错误将第4级（Neutral）纳入了满意范畴，需要把统计范围调整为仅第5、6级（Satisfied和Very Satisfied）：

data_fun <- function(x) {
  # 仅统计Satisfied(5)和Very Satisfied(6)的有效占比
  valid_count <- sum(!is.na(x))
  satisfied_count <- sum(x %in% 5:6, na.rm = TRUE)
  return(satisfied_count / valid_count)
}

2. 修复分面排序 #

偶数级Likert排序失效的根源是排序依据的统计值错误，修正计算逻辑后，需基于正确的满意度占比对分面变量重新排序。假设你的分面分组变量是group，处理方式如下：

library(dplyr)

# 计算每个分组的满意度占比，作为排序依据
sorting_table <- your_data %>%
  group_by(group) %>%
  summarise(satisfaction_rate = data_fun(likert_score)) %>%
  arrange(desc(satisfaction_rate)) # 按满意度从高到低排序

# 将分组变量转换为按满意度排序的因子
your_data$group <- factor(your_data$group, levels = sorting_table$group)

3. 调整总计显示（以ggplot可视化为例） #

如果用ggplot绘制Likert图，需要过滤掉Neutral及以下的类别，仅保留目标两级，并计算其合计占比：

library(ggplot2)

your_data %>%
  # 仅保留Satisfied和Very Satisfied两级
  filter(likert_score %in% 5:6) %>%
  mutate(likert_label = case_when(
    likert_score == 5 ~ "Satisfied",
    likert_score == 6 ~ "Very Satisfied"
  )) %>%
  ggplot(aes(x = group, fill = likert_label)) +
  geom_bar(position = "stack") +
  # 添加满意两级的总计百分比标签
  stat_summary(
    aes(label = paste0(round(..y../sum(..y..)*100, 1), "%")),
    fun = sum, geom = "text", position = position_stack(vjust = 0.5)
  ) +
  facet_wrap(~group, scales = "free_x") +
  coord_flip()

关键注意事项 #

• 确认Likert编码对应关系：如果Satisfied和Very Satisfied对应其他层级（比如非5、6级），只需调整代码中的5:6为对应数字即可。
• 排序逻辑可按需调整：如果需要从低到高排序，将arrange(desc(satisfaction_rate))改为arrange(satisfaction_rate)即可。
• 过滤数据时严格锁定目标层级，确保总计不会混入Neutral或不满意类别的数据。

内容的提问来源于stack exchange，提问作者Homer Jay Simpson