我有一个简单的问题,我在下面使用线性模型
model1 <- lm(output ~. , data=df)
我得到一些输出,我想直接提取线性方程,而不是看到原始形式的输出
我知道扫帚包装我最擅长的是
library(broom)
tidy(mod1)
我要以表单
输出output=1.1*var1 + 1.2*var2+....
答案 0 :(得分:4)
您应该看一下等价软件包(https://github.com/datalorax/equatiomatic)
下面的代码应该可以解决您的问题,我以虹膜为例
remotes::install_github("datalorax/equatiomatic")
library(equatiomatic)
mod1 <- lm(Petal.Width~.-Species, data= iris)
extract_eq(mod1, use_coefs = TRUE, wrap=TRUE,fix_signs = FALSE,ital_vars = TRUE,terms_per_line = 4)