Analysis Method and Experimental Conditions Affect Computed Circadian Phase from Melatonin Data

Accurate determination of circadian phase is necessary for research and clinical purposes because of the influence of the master circadian pacemaker on multiple physiologic functions. Melatonin is presently the most accurate marker of the activity of the human circadian pacemaker. Current methods of analyzing the plasma melatonin rhythm can be grouped into three categories: curve-fitting, threshold-based and physiologically-based linear differential equations. To determine which method provides the most accurate assessment of circadian phase, we compared the ability to fit the data and the variability of phase estimates for seventeen different markers of melatonin phase derived from these methodological categories. We used data from three experimental conditions under which circadian rhythms - and therefore calculated melatonin phase - were expected to remain constant or progress uniformly. Melatonin profiles from older subjects and subjects with lower melatonin amplitude were less likely to be fit by all analysis methods. When circadian drift over multiple study days was algebraically removed, there were no significant differences between analysis methods of melatonin onsets (P = 0.57), but there were significant differences between those of melatonin offsets (P