See my answer here: Using all independent variables with @formula in a multiple linear model - #2 by nilshg
Basically what the preceding posts say - you want to construct Term
objects from symbols, in the other thread I’m doing this for all columns names(df)
but that of course easily generalizes to a vector holding whatever subset of columns you’re interested in.