我在 R 中使用逐步回归,代码如下:
model_scicareer_all <- glm(scicareer ~ .,
family = binomial(link = "logit"), data = clean_data)
summary(model_scicareer_all)
intercept_only <- glm(scicareer ~ -1,
family = binomial(link = "logit"), data = clean_data)
forward <- step(intercept_only, direction='forward', scope=formula(model_scicareer_all), trace=1000)
代码运行没有问题。但是,我想获得每一步的系数。那可能吗? (我知道 SPSS 就是这样做的。)
我尝试了不同的输入,但无法打印每一步的系数
step
有一个 trace
参数,但似乎它不返回系数。请改用如图所示的 trace
函数。
trace(glm, exit = quote(print(coef(returnValue()))))
fm <- glm(mpg ~., data = mtcars)
step(fm)
untrace(flm)