leveneTest用于几个变量

Question

如何单独对每个因变量应用leveneTest，然后如何提取每个分析的P值并将所有p值都放在矩阵中。

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

Answer 1

你可以试试这个：

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]})