Cosinor Analysis of Biorhythm Data

This Demonstration models the 24-hour (circadian) variation of heart rate as a sinusoid with superimposed random noise. The methods of cosinor biorhythm analysis are used to estimate the magnitude, phase, and statistical significance of the modeled rhythm.
The cosinor method fits data to a cosine function, , where and are the mean and amplitude of the sinusoidal rhythm, is the period, and is the time of peak heart rate. In cosinor phraseology, the mean is called the mesor [1, 2, 3].
Plots are shown for the cosinor best-fit solution of and the bivariate error ellipse with amplitude and acrophase of 0.95 confidence intervals, which are dynamically computed and plotted on the polar graph.


  • [Snapshot]
  • [Snapshot]
  • [Snapshot]


The cosinor 0.95 confidence error ellipse is centered at the tip of the amplitude-acrophase vector on the polar plot. The eccentricity and orientation (i.e., tilt) of the ellipse depend on the degree of nonuniformity in data collection, as occurs, for example, when data is missing or otherwise irregularly sampled. When the error ellipse overlaps the pole the amplitude estimate fails the test of significance (). This indicates the absence of a sinusoidal rhythm. In the case of pole overlap, amplitude and acrophase confidence intervals do not have calculable bounds [1].
Solving the nonlinear estimating function is straightforward after linearizing the parameters and applying least-squares techniques. Computing the amplitude and acrophase confidence intervals is more difficult, since efficient close-formed solutions do not exist. We therefore derive the relevant nonlinear implicit equations and apply Mathematica numerical root-finding techniques. This approach is efficient and permits dynamic updating of parameter values in the Demonstration.
The implicit equations defining amplitude and acrophase confidence limits are relevant to detection of circle-ellipse collisions and locating line-ellipse tangency points.
[1] W. Nelson, Y. L. Tong, J. K. Lee, and F. Halberg, "Methods for Cosinor-Rhythmometry," Chronobiologia, 6(4), 1979 pp. 305–323.
[2] R. Ferrara and A. Antonellini, "Using the SAS System for COSINOR Analysis of 24 Hour Ambulatory Blood Pressure Monitoring (ABPM) Data," in SAS European Users Group International Proceedings (SEUGI 1993), Jersey, Channel Islands: SAS Institute, 1993. www.sascommunity.org/seugi/SEUGI1993/Using%20 the %20 SAS %20 System %20 for %20 COSINOR %20 Analysis %20 of %2024 %20 Hour %20 Ambulatory %20 Blood %20 Pressure %20 Monitoring %20 %28 ABPM %29 %20 Data.pdf.
[3] 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(3), 2005 pp. 160–64. advan.physiology.org/content/29/3/160.full.pdf+html.
    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.

Mathematica »
The #1 tool for creating Demonstrations
and anything technical.
Wolfram|Alpha »
Explore anything with the first
computational knowledge engine.
MathWorld »
The web's most extensive
mathematics resource.
Course Assistant Apps »
An app for every course—
right in the palm of your hand.
Wolfram Blog »
Read our views on math,
science, and technology.
Computable Document Format »
The format that makes Demonstrations
(and any information) easy to share and
interact with.
STEM Initiative »
Programs & resources for
educators, schools & students.
Computerbasedmath.org »
Join the initiative for modernizing
math education.
Step-by-step Solutions »
Walk through homework problems one step at a time, with hints to help along the way.
Wolfram Problem Generator »
Unlimited random practice problems and answers with built-in Step-by-step solutions. Practice online or make a printable study sheet.
Wolfram Language »
Knowledge-based programming for everyone.
Powered by Wolfram Mathematica © 2014 Wolfram Demonstrations Project & Contributors  |  Terms of Use  |  Privacy Policy  |  RSS Give us your feedback
Note: To run this Demonstration you need Mathematica 7+ or the free Mathematica Player 7EX
Download or upgrade to Mathematica Player 7EX
I already have Mathematica Player or Mathematica 7+