如何在R循环中批量执行Kruskal检验？原ANOVA方法报错

报错原因分析 #

你的代码存在两个核心问题：

1. 公式方向错误：reformulate(paste(i,'~Species'))生成的是Species ~ 变量名的反向公式，而我们需要的是变量名 ~ Species（因变量在前，分组变量在后），这会导致检验函数无法识别正确的变量关系。
2. 参数传递错误：base R的kruskal.test默认优先接收x和y参数，若你像aov那样把data放在前面，R会调用默认的非公式方法，从而报x缺失的错误。

解决方案代码 #

方案1：使用base R的kruskal.test

vars = c("Sepal.Length" , "Sepal.Width" ,"Petal.Length" ,  "Petal.Width")

for (i in vars) {
  # 生成正确的公式：因变量 ~ 分组变量
  f <- reformulate("Species", response = i)
  # 或用as.formula直接拼接：f <- as.formula(paste(i, "~ Species"))
  
  # 按公式方法的正确参数顺序调用
  test <- kruskal.test(formula = f, data = iris)
  
  cat("\n=== 变量:", i, "===\n")
  print(test)
}

方案2：使用rstatix::kruskal_test

library(rstatix)

vars = c("Sepal.Length" , "Sepal.Width" ,"Petal.Length" ,  "Petal.Width")

for (i in vars) {
  f <- reformulate("Species", response = i)
  # rstatix的参数顺序是data在前，formula在后
  test <- kruskal_test(data = iris, formula = f)
  
  cat("\n=== 变量:", i, "===\n")
  print(test)
}

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