在R中找到最佳线性模型的命令

Lin*_*ter 4 model r linear-regression lm

有没有办法让R运行所有可能的模型(数据集中变量的所有组合)以生成最佳/最精确的线性模型,然后输出该模型?

我觉得有办法做到这一点,但是我很难找到信息。

Kon*_*rad 6

有许多方法可以实现,但是对于执行此操作的简单方法,我建议您看一下该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)