ANOVA分析后filter()函数报错的原因及解决方法

报错原因 #

你遇到的错误是因为rstatix::anova_test在分组操作后返回的是**grouped_anova_test类对象**——它虽然继承了data frame的结构，但属于rstatix自定义的特殊对象，dplyr的filter函数无法直接处理这类带自定义类的对象，底层调用vec_slice时会因对象类型不匹配抛出错误。

解决办法 #

这里提供两种实用的处理方案：

方案1：转换为普通数据框后再筛选 #

将grouped_anova_test对象转换为普通的tibble/data frame，就能正常使用dplyr的筛选函数了：

library(dplyr)
library(rstatix)

stress.effect <- performance %>%
  group_by(time, gender) %>%
  anova_test(dv = score, wid = id, between = stress)

# 转换为普通tibble后执行筛选
stress.effect %>% 
  as_tibble() %>% 
  filter(time == "t2")

方案2：提前过滤数据，只分析目标分组 #

既然你只需要t2时间点的结果，更高效的方式是在分组前就把数据过滤为t2子集，直接得到目标结果：

library(dplyr)
library(rstatix)

# 提前筛选t2数据，再分组做ANOVA
stress.effect_t2 <- performance %>%
  filter(time == "t2") %>%
  group_by(gender) %>%
  anova_test(dv = score, wid = id, between = stress)

# 直接查看t2的结果
stress.effect_t2

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