Tengo una lista de fórmulas que quiero ajustar a los datos, en lugar de ejecutar un bucle, me gustaría hacer esto de inmediato, por el bien del rendimiento. Las estimaciones aún deben estar separadas, no estoy tratando de estimar un SUR ni nada. El siguiente código hace lo que quiero

x <- matrix(rnorm(300),ncol=3)
y <- x %*% c(1,2,3)+rnorm(100)
formulae <-list(y~x[,1],
                y~x[,2],
                y~x[,1] + x[,2])
lapply(formulae,lm)

Desafortunadamente, esto se vuelve algo lento a medida que la duración deformulae aumentos ¿hay alguna manera de realmente vectorizar esto?

Si es de alguna ayuda, los únicos resultados delm Me importan los coeficientes y algunos errores estándar.