R-方差分析提取 p_value 的一种方式 [英] R- One way anova extracting p_value

问题描述

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

I'm trying to do a one-way anova on several row of a dataset and extract the p_value to use it afterwards.

这是我所做的:

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

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

With this formula i'm able to extract the pvalue but it comes with other elements:

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

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

What I would like to have as a result is only this in a list. How could I extract it ?

推荐答案

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

Consider using broom. With tidy(), you can extract only the p.value field:

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

这篇关于R-方差分析提取 p_value 的一种方式的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持IT屋！