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