我有一组因变量y1, y2, ...，一组独立变量x1,x2,...和一组控件d1,d2,...。这些都在一个里面data.table，让我们称之为它data。
我需要做一些事情
out1 <- lm(y1 ~ x1, data=data)
out2 <- lm(y1 ~ x1 + d1 + d2, data=data)
....
当然这不是很好，所以我正在考虑编写一个包含所有这些回归的列表，而不仅仅是遍历该列表。遵循以下原则
myRegressions <- list('out1' = y1 ~ x1, 'out2' = y1 ~ x1 + d1 + d2)
output <- NULL
for (reg in myRegressions)
{
output[reg] <- lm(myRegressions[[reg]])
}
这当然是行不通的：我无法构造列表，因为语法在之外无效lm()。这里有什么好的方法？
公式可以引用为：
myReg <- list('out1' = "mpg ~ cyl")
lm(myReg[[1]],data=mtcars)
Call:
lm(formula = myReg[[1]], data = mtcars)
Coefficients:
(Intercept) cyl
37.885 -2.876
本文收集自互联网，转载请注明来源。
如有侵权，请联系[email protected] 删除。
我来说两句