test_cosinor function

Test for differences in a cosinor model