Это, вероятно, то, что вы хотите, но может потребоваться некоторые изменения, чтобы получить в полной мере то, что вам нужно.
df <- data.frame(WD = c(1,1,0,0,0,1,1,1,0,0),
E1 = rnorm(10,0,1),
E2 = rnorm(10,0,1))
mod <- function(x){
lm(WD ~ x, data = df)
}
sapply(df[setdiff(names(df),"WD")],mod)