从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)

我想提取与a系数相关联的p值进入向量，该系数是一个<2.2e-16的值。

我的代码是：

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

其中vartest产生以下方差表。我想知道我是否做错了，还是有其他方法可以将值提取到向量中？

Answer 1

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

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

-检查OP的方法

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

---

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

str(vartest)
这将提供有关如何提取组件的想法