R语言for循环执行Kruskal-Wallis检验遇model.frame.default错误

你循环代码报错的核心原因是：直接在公式中使用字符串变量i时，kruskal_test会将i当作字面量的变量名，而非循环中索引对应的列名。比如第一次循环时，公式被解析为"Observed" ~ Month，R会尝试寻找名为"Observed"的变量（而非列名为Observed的列），这个变量不存在，就会触发"变量长度不匹配"的错误。

下面提供几种修正方案：

方案1：修改for循环，用reformulate构造公式 #

Indices <- c("Observed", "Shannon", "InvSimpson", "Evenness")
result.kruskal <- data.frame()

for (i in Indices) {
  # 用reformulate生成正确的公式
  formula <- reformulate(response = i, termlabels = "Month")
  kruskal <- df %>% kruskal_test(formula)
  result.kruskal <- rbind(result.kruskal, kruskal)
}

方案2：用as.formula拼接公式字符串 #

Indices <- c("Observed", "Shannon", "InvSimpson", "Evenness")
result.kruskal <- data.frame()

for (i in Indices) {
  # 拼接公式字符串并转为公式对象
  formula <- as.formula(paste(i, "~ Month"))
  kruskal <- df %>% kruskal_test(formula)
  result.kruskal <- rbind(result.kruskal, kruskal)
}

方案3：用purrr替代for循环（更简洁的tidyverse风格） #

library(purrr)
library(rstatix)

Indices <- c("Observed", "Shannon", "InvSimpson", "Evenness")

result.kruskal <- map_dfr(Indices, function(index) {
  df %>% kruskal_test(reformulate(response = index, termlabels = "Month"))
})

这三种方法都能正确生成每个指数对应的X ~ Month公式，从而顺利执行Kruskal-Wallis检验并合并结果。

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