从glm模型摘要中按p值排序xtable()输出

Question

我正在为不同的公司建模大量数据，对于每个公司，我需要快速确定最重要的模型参数。我希望看到一个xtable()输出的拟合模型，按p值递增的顺序排序所有系数（即，最重要的参数首先出现）。

x <- data.frame(a=rnorm(100), b=runif(100), c=rnorm(100), e=rnorm(100))
fit <- glm(a ~ ., data=x)
xtable(fit)

我猜我可以通过更改fit对象的结构来实现这样的事情。但我对结构不太熟悉，无法自信地进行任何更改。

有什么建议吗？

- Brandon Bertelsen

---

当您想要检查R对象时，str()是您的好朋友。 :-) - chl

Answer 1

这并不一定是最优雅的解决方案，但应该能够完成工作:

data(birthwt, package="MASS")
glm.res <- glm(low ~ ., data=birthwt[,-10])
idx <- order(coef(summary(glm.res))[,4])  # sort out the p-values
out <- coef(summary(glm.res))[idx,]       # reorder coef, SE, etc. by increasing p
library(xtable)
xtable(out)

- chl

---

看起来相当直接。正确使用提取器功能。有什么不喜欢的呢？ - IRTFM

完全符合我的要求。谢谢 chl。 - Brandon Bertelsen