R语言中如何提取单因素方差分析的p值并生成纯数值列表？

你的代码思路没问题，但提取p值的方式导致结果附带了多余的元数据标签。咱们可以调整自定义函数，精准提取bt.factor对应的纯p值数值：

修改后的代码 #

# 重写函数，只提取组间效应的p值数值
anova_pval <- function(x) {
  # 定位到方差分析表中bt.factor行的Pr(>F)数值
  summary(aov(x ~ bt.factor))[[1]]$`Pr(>F)`[1]
}

# 逐行计算并得到纯p值向量
anv.pval <- apply(golubALL, 1, anova_pval)

为什么原代码会带额外元素？ #

原代码里的summary(aov(x ~ bt.factor))[[1]]["Pr(>F)"]返回的是带行/列标签的子表，包含了bt.factor和Residuals两行的p值信息，所以输出会附带这些元数据。

修改后的代码做了两个关键优化：

• 用$Pr(>F)``直接定位到p值列
• 用[1]只取第一行（也就是bt.factor对应的组间效应p值）
• 最终返回的是纯数值，没有多余的标签干扰

可选：去掉结果的行名 #

如果连结果向量里的探针名（比如1414_at）也不需要，可以用unname()去除：

anv.pval <- unname(apply(golubALL, 1, anova_pval))

内容的提问来源于stack exchange，提问作者Solal