R使用dplyr/rstatix汇总数据时自定义因子排序丢失问题求助

核心原因 #

你当前遇到的排序丢失问题由两个常见问题导致：

1. mutate()语法使用错误：转换因子时未使用=完成赋值，实际上转换后的因子并未覆盖原有列，你观察到的汇总前排序正常大概率为临时效果
2. ggpubr包的get_summary_stats()默认会将分组变量转为字符型，主动丢弃了因子的层级属性

具体解决办法 #

方案1：修正语法+汇总后重新指定因子层级（适配get_summary_stats） #

先修正mutate的赋值逻辑，汇总完成后再对结果的分组列重新指定因子层级即可：

# 提前定义好层级，避免重复写
date_levels <- c("4/8/2020","4/9/2020","4/10/2020")
plot_levels <- as.character(1:11)

df %>% 
  # 修正mutate赋值语法，用=覆盖原列
  mutate(
    Date.Coll = factor(Date.Coll, levels = date_levels), 
    Shoot.Plot = factor(Shoot.Plot, levels = plot_levels)
  ) %>%
  group_by(Date.Coll, Site, Eelgrass, Oyster, Shoot.Plot) %>%
  filter(
    is.na(BGB),
    # 简化筛选逻辑，用%in%替代多个|
    Date.Coll %in% date_levels
  ) %>% 
  select(AGB) %>% 
  get_summary_stats(type="mean_se") %>%
  # 汇总后重新指定因子层级，保留排序
  mutate(
    Date.Coll = factor(Date.Coll, levels = date_levels), 
    Shoot.Plot = factor(Shoot.Plot, levels = plot_levels)
  )

方案2：使用dplyr原生summarize替代get_summary_stats（自动保留因子属性） #

dplyr原生的分组汇总逻辑不会修改因子的层级属性，不需要二次转换：

date_levels <- c("4/8/2020","4/9/2020","4/10/2020")
plot_levels <- as.character(1:11)

df %>% 
  mutate(
    Date.Coll = factor(Date.Coll, levels = date_levels), 
    Shoot.Plot = factor(Shoot.Plot, levels = plot_levels)
  ) %>%
  group_by(Date.Coll, Site, Eelgrass, Oyster, Shoot.Plot) %>%
  filter(
    is.na(BGB),
    Date.Coll %in% date_levels
  ) %>% 
  # 原生计算均值和标准误，结果自动保留因子排序
  summarize(
    mean = mean(AGB, na.rm = T),
    se = sd(AGB, na.rm = T)/sqrt(n()),
    .groups = "drop"
  )

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