R语言重复测量ANOVA的TukeyHSD事后检验报错求助

核心问题分析 #

你的脚本存在三个关键问题：

• 因变量类型错误：vo2列被定义为字符型（数值加了引号），R会将其当作分类变量处理，导致ANOVA完全失效。
• 重复测量模型结构错误：原脚本的Error(factor(id))未正确建模被试内变异，重复测量ANOVA需要用Error(id/treat)来区分被试间和被试内的变异来源。
• TukeyHSD函数不兼容：直接对带Error层的aov对象调用TukeyHSD会报错，这类对象需要专门的事后检验工具处理。

修正后的完整代码 #

第一步：修正数据并构建正确的重复测量ANOVA模型

my_data <- data.frame(
  stringsAsFactors = FALSE,
  id = c(1L,2L,3L,4L, 5L,1L,2L,3L,4L,5L,1L,2L,3L,4L,5L,1L,2L,3L,4L,5L),
  treat = c("o","o","o","o","o","j","j","j","j","j","z","z","z","z","z","w","w","w","w","w"),
  vo2 = c("47.48","42.74","45.23","51.65","49.11","51.00","43.82","49.88","54.61","52.20","51.31",
          "47.56","50.69","54.88","55.01","51.89","46.10","50.98","53.62","52.77"))

# 把vo2转为数值型（关键修正）
my_data$vo2 <- as.numeric(my_data$vo2)

# 构建正确的重复测量ANOVA模型
rm_result <- aov(vo2 ~ treat + Error(id/treat), data = my_data)
summary(rm_result)

第二步：执行Tukey HSD事后检验（推荐用emmeans包）

emmeans包能更好地处理重复测量模型的事后比较，自动考虑数据的相关性：

# 首次使用需安装包
# install.packages("emmeans")

library(emmeans)

# 生成边际均值并做Tukey校正的两两比较
emmeans(rm_result, pairwise ~ treat, adjust = "tukey")

备选方案：提取模型内效应后用TukeyHSD

如果坚持使用基础包的TukeyHSD，可以提取模型中被试内的效应部分：

# 提取被试内的aov结果
within_aov <- rm_result$`Within`
TukeyHSD(within_aov, "treat", ordered = TRUE)

结果说明 #

• emmeans的输出会包含各处理组的边际均值、标准误，以及所有两两比较的差异值、标准误、t值和校正后的p值。
• 修正后的模型能正确分解重复测量设计中的变异，确保ANOVA和事后检验结果的可靠性。

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