多个变量的 leveneTest

leveneTest for several variables

如何分别对每个因变量应用 leveneTest,然后如何提取每个分析的 P 值并将所有 p 值放入矩阵中。

responses <- as.matrix(mtcars[,-c(4,8,9)])
leveneTest(responses~as.factor(am)*as.factor(vs),data=mtcars)

你可以试试这个:

responses <- as.matrix(mtcars[,-c(4,8,9)])
results <- data.frame(var = colnames(responses), p = rep(NA,dim(responses)[2]))
results$p <- apply(responses,2,function(x) {leveneTest(x ~ as.factor(mtcars$am)*as.factor(mtcars$vs))[1,3]})