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

Question

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

如何只打印回归模型的（调整后）R平方？在这个帖子中，是只打印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，你不能直接从模型的P值中得到 summary - 它的计算是在飞行过程中由 print.summary.lm 方法。你可以用

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