从步骤函数中提取被拒绝变量的名称

Question

我在glm对象上使用step函数

stepglm<-step(glmmodel)

并且步骤拒绝来自glmmodel的2个变量。有没有办法从stepglm对象中提取被拒绝的变量的名称？似乎stepglm只是一个glm对象所以也许没有简单的方法可以做到这一点？

Answer 1

在MASS包中尝试stepAIC。这是一个例子：

library(MASS)

example(birthwt) # construct bwt

birthwt.glm <- glm(low ~ ., family = binomial, data = bwt)
birthwt.step <- stepAIC(birthwt.glm, trace = FALSE)
as.character(birthwt.step$anova$Step)[-1]
## [1] "- ftv" "- age"