Извлечение p-значения из таблицы anova в R

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. </p>

Код, который у меня есть:

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)

, и это даст представление о том, как извлечь компоненты