如何对数据框多变量执行Levene's F检验并提取p值？

你的问题出在对lapply返回结果的处理上：leveneTest会返回一个包含检验结果的对象，lapply处理多个变量后得到的是列表（每个元素对应一个变量的检验结果），不能直接在整个列表后提取'Pr(>F)'，需要对列表中的每个元素单独提取p值。

正确实现方式 #

首先确保加载car包（leveneTest属于该包），然后用以下两种方法之一实现：

方法1：用sapply（和你之前做t检验的逻辑一致）

library(car)

df <- data.frame(EMO_ETA.3 = c(1, 4, 5, 6), EMO_ETA.6 = c(4, 1, 8, 9), zyg = c("1", "1", "2", "2"))

# 批量提取Levene检验的p值
sapply(df[c("EMO_ETA.3", "EMO_ETA.6")], function(x) {
  leveneTest(x ~ df$zyg)$'Pr(>F)'[1]
})

方法2：用lapply分步处理

先存储所有检验结果，再逐个提取p值：

library(car)

df <- data.frame(EMO_ETA.3 = c(1, 4, 5, 6), EMO_ETA.6 = c(4, 1, 8, 9), zyg = c("1", "1", "2", "2"))

# 先得到所有变量的Levene检验结果列表
levene_results <- lapply(df[c("EMO_ETA.3", "EMO_ETA.6")], leveneTest, group = df$zyg)

# 从列表中提取每个结果的p值
unlist(lapply(levene_results, function(res) res$'Pr(>F)'[1]))

原代码返回NULL的原因 #

lapply返回的是一个列表，列表的每个元素是单个变量的Levene检验结果（包含'Pr(>F)'字段），但列表本身没有'Pr(>F)'这个元素，所以直接写lapply(...)$'Pr(>F)'会返回NULL。

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