I am running a regression in a loop of the following form:
regression[[i]] <- glm(as.formula(paste(depvars[i],"~",indepvars[j],"+ var1 + var2
