从 R 中的方差分析表中提取 p 值

【问题标题】：Extracting p-value from anova table in R从 R 中的方差分析表中提取 p 值
【发布时间】：2020-05-01 19:29:30
【问题描述】：

** 前面的步骤是这样的：

a <- rnorm(100, mean=2, sd=3)
response <- a * 1.5 + rnorm(100, mean=0, sd=1)
model <- lm(response ~ a)
vartest <- anova(model)

我想将 p 值提取到与 a 系数相关联的向量中，该值将小于 2.2e-16。

我的代码是：

vartest[1,5]
[1] 1.002182e-63

其中vartest 生成以下方差表。我想知道我是否做错了，或者是否有这种方法的替代方法，用于将值提取到向量中？

【问题讨论】：

【解决方案1】：

我们可以使用[[或$的列名直接提取

out <- vartest[["Pr(>F)"]][1]
is.vector(out)
#[1] TRUE

-检查 OP 的方法

identical(out, vartest[1,5])
#[1] TRUE

我们可以用str检查对象的结构

str(vartest)

这将提供有关如何提取组件的想法

【讨论】：