A sample of the objective dynamic measurement data obtained from the refractor is shown in
Figure 1 . Note the asymmetries between accommodation and disaccommodation in some cycles and the usual fluctuations in accommodation response.
The records were analyzed using fast Fourier analysis and also with the peak-fitting module in a software program (Origin Pro 7; Microcal Software Inc., Northampton, MA). The peak-fitting module was used to detect the peak locations automatically, by means of an irregular sinusoid function, once the baseline refractive error was specified. No averaging was performed, and each response was analyzed individually. The time constants for accommodation and disaccommodation were calculated by assuming that the accommodation response during the changes could be described by the exponential equations
\[\mathrm{A}\ {=}\ A_{H}\ {-}\ \mathrm{a}\ {\times}\ \mathrm{e}^{{-}t/{\tau}}\]
for accommodation and
\[\mathrm{A}\ {=}\ A_{L}\ {+}\ \mathrm{a}\ {\times}\ \mathrm{e}^{{-}t/{\tau}}\]
for disaccommodation where,
A is the accommodative response;
A H and
A L are the accommodative responses at the near and distance peaks, respectively;
t is time in seconds after commencement of the response; and τ is the time constant in seconds. The amplitude of accommodative response (
a =
A H −
A L) for each accommodative half-cycle was calculated from the data obtained from the peak locations. The amplitude data were then used to calculate the response levels when 10% and 90% of the accommodative response amplitude was reached, and the refractor records were analyzed to find the corresponding times:
t 10 and
t 90. These values were then used to calculate the time constants using a derivation from the above equation where,
\[{\tau}\ {=}\ \frac{\mathrm{t}_{90}\ {-}\ \mathrm{t}_{10}}{\mathrm{ln}9}.\]
On the above assumption of an exponential response change, the maximum velocities (
V) of accommodation and disaccommodation were calculated as
V = (
dA/
dt)
t =0 =
a/τ.
A fast Fourier transform (FFT) was performed for data obtained over approximately 41 seconds, which produced 1024 data points. The amplitude of the Fourier transform at each frequency was determined by calculating the square root of the sum of squares of the imaginary and the real components of the FFT. The fundamental frequency peak obtained for this data set was then converted to the corresponding rate in cycles per minute.