只打印 r 中线性模型的 p 值

Question

有没有办法在 R 中只打印线性模型的 p 值？我只需要全局模型的 p 值，而不是系数。

在这个 post 中是只打印 adj 的解决方案。 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 中，您无法直接从 summary 中获取模型 p 值 - 它是通过 print.summary.lm 方法动态计算的。您可以使用

复制此计算

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

Only print p-values of a linear model in r