R-方差分析提取 p_value 的一种方式

Question

我正在尝试对数据集的几行进行单向方差分析，然后提取 p_value 以供使用。

这是我所做的：

anova <- function(x) {summary(aov(x ~ bt.factor))[[1]]["Pr(>F)"]}
anv.pval <- apply(golubALL, 1, anova)

使用这个公式，我可以提取 pvalue 但它带有其他元素：

$`1414_at`
            Pr(>F)
bt.factor   0.7871
Residuals

结果我想要的只是列表中的这个。我怎么能提取它？

Answer 1

考虑使用broom. 使用tidy()，您只能提取p.value字段：

require(broom)
a <- aov(mpg ~ wt, mtcars)

tidy(a)
#        term df    sumsq     meansq statistic      p.value
# 1        wt  1 847.7252 847.725250  91.37533 1.293959e-10
# 2 Residuals 30 278.3219   9.277398        NA           NA

tidy(a)$p.value
# [1] 1.293959e-10           NA