我想对相同的结果进行线性回归,并且在每个模型中使用一些协变量减去一个协变量。我已经看过这个page上的例子,但这可能没有提供我想要的东西。
示例数据
a <- data.frame(y = c(30,12,18), x1 = c(7,6,9), x2 = c(6,8,5),
x3 = c(4,-2,-3), x4 = c(8,3,-3), x5 = c(4,-4,-2))
m1 <- lm(y ~ x1 + x4 + x5, data = a)
m2 <- lm(y ~ x2 + x4 + x5, data = a)
m3 <- lm(y ~ x3 + x4 + x5, data = a)
我怎样才能在短时间内运行这些模型,而且不会一次又一次地重复相同的协变量?
答案 0 :(得分:2)
关注this example你可以这样做:
lapply(1:3, function(i){
lm(as.formula(sprintf("y ~ x%i + x4 + x5", i)), a)
})