我正在为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"
答案 0 :(得分:1)
怎么样
formula(mdl)
? .....