Печатайте только p-значения линейной модели в r

Question

Есть ли способ распечатать только p-значения линейных моделей в R? Мне нужно только p-значение глобальной модели, а НЕ коэффициентов.

Как только распечатать (скорректированный) R-квадрат регрессионной модели? В этом посте есть решение для только печать прил. r².

Я пробовал summary(model)$p.value, но он не работает.

Answer 1

Самый простой способ - использовать пакет broom:

library(broom)
# coefficient p values
tidy(model)$p.value

# model p value
glance(model)$p.value

Оставаясь на базе, вы можете использовать

# coefficient p values
summary(model)$coefficients[, "Pr(>|t|)"]

В base не получить p-значение модели прямо из summary - оно рассчитывается на лету методом print.summary.lm. Вы можете повторить этот расчет с помощью

with(summary(model), 
  pf(
    fstatistic[1L], 
    fstatistic[2L],
    fstatistic[3L],
    lower.tail = FALSE
  )
)