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

Question

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

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产生以下方差表。我想知道我是否做错了，或者是否有这种方法的替代方法，用于将值提取到向量中？

标签： rvectoranova

Answer 1

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

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

-检查OP的方法

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

---

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

str(vartest)

这将给出一个关于如何提取组件的想法