Eu tenho uma lista de fórmulas que quero ajustar aos dados, em vez de executar um loop que gostaria de fazer isso de uma só vez, por uma questão de desempenho. As estimativas ainda devem ser separadas, não estou tentando estimar um SUR ou algo assim. O código a seguir faz o que eu quero

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)

Infelizmente, isso fica um pouco lento, pois o comprimento daformulae aumenta existe uma maneira de realmente vetorizar isso?

Se for de alguma ajuda, os únicos resultados delm Eu me preocupo com os coeficientes e alguns erros padrão.