`fit_model` - fitting the model

[gareth-j]

Create a model fitting function (for internal use) that can be used easily as it has sensible defaults.