cosinor.lm function

Fit cosinor model