As described previously,
9 the PCI tracings were processed by a semiautomated algorithm to calculate axial length, corresponding to the distance from the corneal surface to a reflective surface located at the region of the interface of Bruch’s membrane and the retinal pigment epithelium interface.
9 10 11 This protocol provides an axial length that differs from that of ultrasonography, because ultrasound measures to the retinal surface.
12
Average daily axial length and IOP were calculated for each data set by using the mean of all measurement series taken on each day for each eye. High–low differences during the day for each data set for both axial length and IOP were calculated as the difference between the mean value at the time of largest parameter value and that of the smallest. A one-way analysis of variance (ANOVA) with replicate measures using a generalized linear model (SAS 8.2; SAS Institute, Inc., Cary, NC) was fit to each data set to determine whether the axial length measured at any of the time points differed significantly from the others. We used a criterion of P < 0.05 from the ANOVA to identify data sets that showed significant daily high–low differences in either parameter. Unless otherwise noted, all mean data are presented as the mean ± SD.
To obtain estimates of the period and phase of the daily fluctuations, the data sets with a statistically significant high–low difference in axial length and/or IOP were modeled with a sine curve
(Table 1) . Sine and cosine curves are traditionally used to model diurnal rhythms.
6 13 14 To achieve this modeling, an adjusted value was calculated by subtracting the mean daily axial length or IOP from the measured value at each time point, and sine curve functions were fit to the adjusted value versus time of day by a nonlinear model (PROC NLIN; SAS 8.2; SAS Institute, Inc.). The following equation was used to curve-fit the data for both axial length and IOP:
\[y\ {=}\ (a/2)\ {\cdot}\ \mathrm{sin}(2{\pi}\ {\cdot}\ \mathrm{time}/b{+}c),\]
where
y represents the adjusted parameter (either axial length or IOP),
a represents the peak-to-trough difference,
b represents the period, and
c represents the phase of the sine curve. The period was constrained in the model to be 24 ± 12 hours. The model yielded estimates with 95% confidence intervals for the amplitude of the peak-to-trough difference, the period and phase for each individual and, as indicators of goodness-of-fit of the model, the correlation coefficient (
R 2) and the probability of the model fit. In addition, the time of maximum axial length or IOP was estimated by solving the equation sin (2π · time/
b̂ + ĉ) = 1 for time, with the constraint of time between 0 and 24 hours and where
b̂ and ĉ were the period and phase estimated from the sin-curve fitting, respectively.