Lin*_*ter 4 model r linear-regression lm
有没有办法让R运行所有可能的模型(数据集中变量的所有组合)以生成最佳/最精确的线性模型,然后输出该模型?
我觉得有办法做到这一点,但是我很难找到信息。
有许多方法可以实现,但是对于执行此操作的简单方法,我建议您看一下该glmulti软件包,本文对此进行了详细描述:
另外,可以在Quick-R网站上找到非常简单的模型选择示例:
# Stepwise Regression
library(MASS)
fit <- lm(y~x1+x2+x3,data=mydata)
step <- stepAIC(fit, direction="both")
step$anova # display results
Run Code Online (Sandbox Code Playgroud)
为了进一步简化,您可以进行更多的手动模型比较:
fit1 <- lm(y ~ x1 + x2 + x3 + x4, data=mydata)
fit2 <- lm(y ~ x1 + x2, data=mydata)
anova(fit1, fit2)
Run Code Online (Sandbox Code Playgroud)