Para modelos longos e repetidos, quero criar uma "macro" (chamada Stata e realizada com oglobal var1 var2 ...) que contém oregressores da fórmula do modelo.

Por exemplo, de

library(car)
lm(income ~ education + prestige, data = Duncan)

Eu quero algo como:

regressors <- c("education", "prestige")
lm(income ~ @regressors, data = Duncan)  

Eu poderia encontrar éessa abordagem. Mas minha aplicação nos regressores não funcionará:

reg = lm(income ~ bquote(y ~ .(regressors)), data = Duncan)

como me joga:

Error in model.frame.default(formula = y ~ bquote(.y ~ (regressors)), data =
Duncan,  :  invalid type (language) for variable 'bquote(.y ~ (regressors))'

Até a resposta aceita da mesma pergunta:

lm(formula(paste('var ~ ', regressors)), data = Duncan)

ataca e me mostra:

Error in model.frame.default(formula = formula(paste("var ~ ", regressors)),
: object is not a matrix`. 

E é claro que eu tenteias.matrix(regressors) :)

Então, o que mais posso fazer?