I am often fitting models for more than one dependent variable in this way. It lets me get the coefficients quickly for plotting. But then I have a hard time passing these results to modelsummary() to print the full regression tables.