從階梯函數

Question

提取拒絕變量的名稱我使用GLM對象

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"