Cosinor Analysis Using Rotating Ellipse Tangents and Circle Collisions

Cosinor analysis uses a least squares method to fit a cosine wave to biorhythm time series data. The method can be applied to nonuniformly sampled studies with missing data points. The fitted cosine function is given by , where is the amplitude and is the acrophase of the time series [1], [2]. The acrophase is the time period in which the cycle peaks. The Demonstration plots the cosinor 95% bivariate error ellipse of the cosinor parameter estimates and dynamically computes associated confidence interval (CI) limits.

In the bookmarks, data is shown from cosinor analysis of 24-hour heart rate (HR) recordings in humans. HR is reported by its reciprocal parameter, RR interval, the time between consecutive heartbeats measured in milliseconds. Illnesses such as heart disease and depression disturb the amplitude and acrophase of 24-hour HR periodicity.

Using coordinate geometry analysis to characterize tangency conditions and circle-ellipse collisions, Mathematica root-solving methods are applied dynamically to determine acrophase and amplitude confidence interval (CI) limits. The cosinor error ellipse parametric curve, and , lies on the axis of the acrophase angle () offset by a tilt (). Cosinor studies with uneven and/or missing data produce a rotated, eccentric ellipse, which influences CI limits.

References

[1] M. A. Rol de Lama, J. P. Lozano, V. Ortiz, F. J. Sánchez-Vázquez, and J. A. Madrid, "How to Engage Medical Students in Chronobiology: An Example on Autorhythmometry," Advances in Physiology Education, 29, 2005 pp. 160–164. advan.physiology.org/content/29/3/160.full.pdf+html.