有没有办法让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
为了进一步简化,您可以进行更多的手动模型比较:
fit1 <- lm(y ~ x1 + x2 + x3 + x4, data=mydata)
fit2 <- lm(y ~ x1 + x2, data=mydata)
anova(fit1, fit2)
关于r - 在R中找到最佳线性模型的命令,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/33740012/