R中ANOVA后执行Tukey检验提示无因子报错如何解决？

问题根源 #

你的sleep_category变量当前为数值型，R的aov()函数会默认将数值型自变量识别为连续协变量，而非多水平分组因子：

• 原ANOVA结果中sleep_category自由度为1，不符合3水平分类变量自由度应为2的预期，说明原模型本质是线性回归而非多组方差分析，结果不代表三组间的整体差异
• TukeyHSD要求自变量为因子类型才能开展多组两两比较，因此触发报错

解决步骤 #

1. 将sleep_category转换为因子类型

# 基础转换，自动识别取值为水平
sleep_analysis$sleep_category <- factor(sleep_analysis$sleep_category)

# 如需自定义水平对应标签可使用以下写法，标签内容可根据你的实际分组含义修改
sleep_analysis$sleep_category <- factor(sleep_analysis$sleep_category,
                                        levels = c(1,2,3),
                                        labels = c("轻度睡眠不足", "正常睡眠", "超长睡眠"))

2. 重新拟合ANOVA模型

res.aov_distance <- aov(distance ~ sleep_category, data = sleep_analysis)

此时运行summary(res.aov_distance)可看到sleep_category的自由度变为2，说明模型已按多组方差分析逻辑拟合。
3. 运行Tukey检验

TukeyHSD(res.aov_distance)

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