使用StepAIC后从模型中提取公式

时间:2020-09-11 14:11:56

标签: r dplyr glm

我正在为R运行stepAIC。我正在尝试从返回的值中提取最佳模型的公式以及最佳模型本身。有什么想法吗?

示例:

library(tidyverse)
library(MASS)
df<-mtcars 
mdl <- glm.nb(mpg ~ cyl +hp +drat+wt+qsec+gear+am +carb , data = df)
mdl <- stepAIC(mdl, trace = FALSE)

编辑 我知道这段代码可以调整很多值,我只需要公式即可。

mdl$terms
mpg ~ hp + wt
attr(,"variables")
list(mpg, hp, wt)
attr(,"factors")
    hp wt
mpg  0  0
hp   1  0
wt   0  1
attr(,"term.labels")
[1] "hp" "wt"
attr(,"order")
[1] 1 1
attr(,"intercept")
[1] 1
attr(,"response")
[1] 1
attr(,".Environment")
<environment: R_GlobalEnv>
attr(,"predvars")
list(mpg, hp, wt)
attr(,"dataClasses")
      mpg        hp        wt 
"numeric" "numeric" "numeric" 

1 个答案:

答案 0 :(得分:1)

怎么样

formula(mdl)

? .....