我有几个因变量(我们称它们为dep1,dep2和dep3),我想在它们上回归几个自变量。我有自变量a,b和c,我总是希望包含在回归中。但是,只有当数字与因变量的数量匹配时,才应使用自变量d1,d2,d3。这是:
lm(dep1 ~ a + b + c + d1)
lm(dep2 ~ a + b + c + d2)
lm(dep3 ~ a + b + c + d3)
更确切地说,我的因变量和自变量都是月度时间序列(xts对象)。我合并了因变量以获得矩阵“dep”。 a,b,c,d1,d2和d3存储在矩阵“indep”中。
是否可以创建一个自动执行此操作的循环?
答案 0 :(得分:0)
我们可以使用paste
v1 <- paste0("dep", 1:3)
v2 <- paste0("d", 1:3)
v3 <- paste(v1, paste(paste(letters[1:3], collapse=" + "), v2, sep=" + "), sep="~")
lapply(v3, function(x) lm(x))