fit_model_and_process function

Fit the cosinor GLMM model using the output from update_formula_and_data() and a new formula