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