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

Extracting p-value from anova table in R

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

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 产生以下方差 table。 我想知道我是否做错了,或者是否有替代方法可以将值提取到向量中?

我们可以使用 [[$

直接提取列名
out <- vartest[["Pr(>F)"]][1]
is.vector(out)
#[1] TRUE

-检查OP的方法

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

我们可以用str

检查对象的结构
str(vartest)

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