Thanks for sharing the code and for the {tidycensus} pkg, it's quite informative and also a great opportunity to use {geofacet}

(2bis) Gamma regression

If lapply doesn't count as explicit loop, I'll do something like this

You can use .SDcols = patterns(....) to select the columns you want using regular expression, something like this should work.