r - 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ます。しかし、私は自信を持って何かを変更できるほど構造に精通していません。

提案？

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)