R语言执行TukeyHSD检验出现undefined columns selected错误求解

报错的核心触发原因有两个：

• 变量名包含空格，TukeyHSD函数解析带空格的列名时无法正确匹配数据框的列，触发列不存在的报错
• 拟合aov模型时没有显式指定data参数，函数后续调用时无法找到原始数据的对应列

解决步骤 #

1. 先对原始数据集的列名重命名，替换所有空格和特殊符号

# 请将下方的df替换为你自己的原始数据集名称
colnames(df) <- c("freq_ingest", "hyd_reg", "depth_zone", "st_loc")

2. 显式指定data参数重新拟合方差分析模型
   从你的ANOVA结果来看所有交互项p值均大于0.05，无统计显著性，可简化模型减少Tukey检验的无效输出：

# 全三因素交互模型写法
data.aov <- aov(freq_ingest ~ hyd_reg * depth_zone * st_loc, data = df)

# 简化主效应模型写法（移除无意义的高阶交互）
data.aov_simple <- aov(freq_ingest ~ hyd_reg + depth_zone + st_loc, data = df)

3. 重新执行Tukey检验即可正常输出结果

# 执行检验
tukey_result <- TukeyHSD(data.aov_simple)
# 输出检验结果
print(tukey_result)

# 如需生成多重比较字母标记
library(multcompView)
multcompLetters(extract_p(tukey_result))

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