fit_models function

Fit models in parallel across a list-column of datasets