使用R,我使用以下脚本w3
您可以使用sapply
遍历每个方差分析并提取p值。
sapply(w3, function(x) summary(x)[[1]][["Pr(>F)"]][[1]])
与mtcars
一起使用可复制的示例
w3 <- lapply(split(mtcars, mtcars$cyl), aov, formula=mpg ~ am)
sapply(w3, function(x) summary(x)[[1]][["Pr(>F)"]][[1]])
# 4 6 8
#0.0892 0.2209 0.8662