[嗨,我开始使用r,并且一直在分析我的数据。我有一个具有80列的数据框。第1列是因变量,从2到80列是自变量。我要执行78个多元线性回归,使模型的第一个自变量固定不变(第2列),并创建一个列表,在其中可以保存所有回归,以便以后可以使用AIC得分比较模型。我该怎么办?
这是我的循环
data.frame
for(i in 2:80)
{
Regressions <- lm(data.frame$column1 ~ data.frame$column2 + data.frame [,i])
}
以iris
数据集为例,您可以这样做:
lapply(seq_along(iris)[-c(1:2)], function(x) lm(data = iris[,c(1:2, x)]))
[[1]]
Call:
lm(data = iris[, c(1:2, x)])
Coefficients:
(Intercept) Sepal.Width Petal.Length
2.2491 0.5955 0.4719
[[2]]
Call:
lm(data = iris[, c(1:2, x)])
Coefficients:
(Intercept) Sepal.Width Petal.Width
3.4573 0.3991 0.9721
[[3]]
Call:
lm(data = iris[, c(1:2, x)])
Coefficients:
(Intercept) Sepal.Width Speciesversicolor Speciesvirginica
2.2514 0.8036 1.4587 1.9468