|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1From the Glaucoma Institute, SUNY State College of Optometry, New York, New York; and the 2School of Optometry, Indiana University, Bloomington, Indiana.
| Abstract |
|---|
|
|
|---|
METHODS. Twenty patients with glaucoma and 20 control subjects were tested with a CSP protocol implemented on a monitor-based testing station. In the protocol 26 locations were tested over the central visual field with Gabor patches with a peak spatial frequency of 0.4 cyc/deg and a two-dimensional spatial Gaussian envelope, with most of the energy concentrated within a 4° circular region. Threshold was estimated by a staircase method. Patients and 10 age-similar control subjects were also tested on conventional automated perimetry (CAP), with the 24-2 pattern with the SITA Standard testing strategy. The neuroretinal rim area of the patients was measured with a retinal tomograph (Retina Tomograph II [HRT]; Heidelberg Engineering, Heidelberg, Germany). A Bland-Altman analysis of agreement was used to assess test-retest variability, compare depth of defect shown by the two perimetric tests, and investigate the relations between contrast sensitivity and neuroretinal rim area.
RESULTS. Variability showed less dependence on defect depth for CSP than for CAP (z = 9.3, P < 0.001). Defect depth was similar for CAP and CSP when averaged by quadrant (r = 0.26, P > 0.13). The relation between defect depth and rim area was more consistent with CSP than with CAP (z = 9, P < 0.001).
CONCLUSIONS. The implementation of CSP was successful in reducing test-retest variability in glaucomatous defects. CSP was in general agreement with CAP in terms of depth of defect and was in better agreement than CAP with HRT-determined rim area.
Pearson et al.5 showed that large chromatic stimuli can provide low test-retest variability while maintaining sensitivity to glaucomatous loss. They hypothesized that the key aspect of the stimulus was the nature of cortical pooling and that a similar result would hold for any stimulus for which detection was mediated by cortical cells tuned to low spatial frequencies. This prediction was confirmed by Pan et al.,6 who used both chromatic stimuli and a Gabor stimulus with a peak spatial frequency of 0.5 cyc/deg.
Another form of sinusoidal stimuli with low spatial frequencies (0.25 and 0.50 cyc/deg) is used in frequency-doubling perimetry (FDP), and the results of studies on test-retest variability with FDP are consistent with the predictions of Pearson et al.5 FDP defects are generally similar in depth to defects with conventional perimetry, when expressed as log contrast sensitivity, and test-retest variability is independent of sensitivity.7 8 However, normal between-subject and test-retest variability are high for FDP. A recent study9 showed large, normal, between-subject variability on FDP: an MD of –5 dB was within 95% confidence limits of normal, as was a pattern standard deviation (PSD) of 5 dB (up to 8 dB in some cases). The 90% confidence interval for test-retest variability on FDP is ±4 dB even in areas with normal sensitivity.7 8 FDP uses high temporal frequencies with the goal of tapping spatially nonlinear retinal ganglion cells. However, a recent study demonstrated that spatially nonlinear cells cannot be the physiological substrate of the frequency-doubling illusion, so the use of high temporal frequency modulation may not be necessary.10 The use of high temporal frequencies makes FDP susceptible to the effects of reduced retinal illumination, due to factors such as pupillary miosis and increase in lens density.11
The goal of the present study was to design a new perimetric protocol, by using contrast sensitivity perimetry (CSP),12 which could improve clinical management of glaucoma by building on findings of these earlier studies. We used a stimulus with a peak spatial frequency in the range used by FDP, but without the sharp stimulus edges and high temporal frequencies used in FDP. The primary purpose of this design was to decrease test-retest variability in defective areas while maintaining good ability to detect glaucomatous loss. The success of the design was measured by comparing results of the new CSP test with results from CAP and with measures of neuroretinal rim area by using scanning laser tomography.
| Methods |
|---|
|
|
|---|
To compare the new CSP test with CAP in detecting depth of defect and variability in areas with low sensitivity, we recruited 20 patients with glaucoma; ages were 44 to 79 years (mean, 63 ± 7). Inclusion criteria were: diagnosis of glaucoma and undergoing treatment at UOC, clear ocular media, spherical correction within 6 D and cylindrical correction within 3 D, IOP below 22 mm Hg, and visual field defects consistent with a diagnosis of glaucoma.13 Subjects clinical characteristics including refractive error are listed in Table 1 . The spherical-cylindrical equivalent of refractive error ranged from –3.25 to +4.75 D across all subjects, which would not be expected to produce significant optical distortions of the stimuli.14 More than half of the patients had a superior or inferior nasal defect; half of the patients had a superior or inferior arcuate scotoma; six patients had combinations of defects including enlarged blind spot and nasal, arcuate, double arcuate, and/or paracentral defects. All but one of the patients had primary open-angle glaucoma; one patient had a diagnosis of normal-tension glaucoma. Diagnosis of glaucoma was established by a treating clinician and confirmed by the second author, based on medical and family history, slit lamp biomicroscopy (including gonioscopy), applanation tonometry, dilated funduscopy, stereoscopic ophthalmoscopy of the optic disc with a 78-D lens, stereo photographs of the optic nerve, optic nerve imaging, and visual field performance. The intraocular pressure before treatment initiation was >21 mm Hg for all but one patient and was successfully controlled by topical medications: two patients used three different drops, and the rest used only one or two different drops. Exclusion criteria were: eye disease (other than glaucoma), systemic disease, or medication known to affect visual function.
|
Apparatus
A field perimeter (Humphrey Field Analyzer II, HFA; Carl Zeiss Meditec, Dublin, CA) was used with the SITA 24-2 testing strategy, to provide a standard evaluation of the central visual field (±24°) of patients with glaucoma and age-similar control subjects. This is a self-calibrating apparatus, and every 6 months a company representative is scheduled for maintenance and evaluation.
A custom-built system based on the Visual Stimulus Generator (VSG2/5, Cambridge Research Systems Ltd.; Cambridge, UK) was used for specialized testing. Dual 8-bit video DACs provided 15-bit output resolution per phosphor. Stimuli from the VSG were displayed on a 21-in. monitor (F500 Trinitron CRT; Sony, Tokyo, Japan). The resolution of the monitor was set to 800 x 600 pixels, with a 150-Hz frame rate. The apparatus was calibrated with a photometer system (Opti-Cal; Cambridge Research Systems Ltd.), to measure luminance versus DAC values for each phosphor, to calculate the transfer functions, and to produce RGB gamma correction look-up tables. Mean luminance was measured periodically throughout the study with an LS-100 photometer (Minolta, Ramsey, NJ). A customized program was written in C++ (Borland, Cupertino CA).
Each participant was seated with forehead and chin on a head and chin rest at a fixed distance of 33 cm for the HFA and 40 cm for the VSG. One eye at a time was tested, while the other was patched. Closed-circuit video systems were used to monitor fixation stability on both devices. An appropriate corrective lens was placed at the fixed lens holder in front of the participant for each device. All subjects were refracted before testing, to determine refractive correction. For HFA and VSG testing of patients and age-similar controls, the near correction was determined based on age, distance correction, and viewing distance. All young subjects wore contact lenses with no other correction in place. For the VSG testing, a near visual acuity card was used to ensure the participants ability to read at least the 20/25 line monocularly at 40 cm with the correction in place.
A retinal tomograph (Retina Tomograph II [HRT II]; Heidelberg Engineering, Heidelberg, Germany) was used to acquire images of the optic nerve head and for quantitative analysis of the neuroretinal rim. A company representative performs annual calibration and maintenance. All images were obtained by the same researcher with a test distance of 1.0 to 1.5 cm from the eye, as the subject was viewing a fixation target. The spherical equivalent of the subjects prescription was dialed on the adjustable scanning lens, to focus the image. The quality of an image was assessed by the HRT software and the researcher. Images with a standard deviation of height measure that was more than 40 µm were excluded as unreliable, and the test was repeated. This method resulted in images for all patients that met the criterion of SD no more than 40 µm. The contour line on the optic disc edge was drawn by a single experienced clinician (MWD). The standard reference plane was used.
Stimulus
Two types of stimuli were used, Gabor patches and size III increments. Gabor patches were used only on the VSG apparatus. Size III stimuli, circular achromatic increments with a diameter of 0.43°, were used on both the HFA and VSG. On the HFA, size III stimuli were presented on a 10-cd/m2 uniform background, as a rectangular temporal pulse of 200 ms duration with Weber contrast between 31,500% and 0.3% in 0.1-log unit steps. On the VSG, size III stimuli were presented on a background with a uniform luminance of 12 cd/m2 and with Weber contrast between 900% and 0.1% in 0.15-log unit steps. The size III stimuli and the Gabor patches were presented on the VSG screen with a temporal Gaussian envelope of 600 ms (time constant of 100 ms); 68% of energy was within the central 200 ms.
The Gabor patches had a peak spatial frequency of 0.4 cyc/deg and were presented on a 55 cd/m2 background with Weber contrast between 100% and 0.1% in 0.15-log unit steps. The patches had a two-dimensional spatial Gaussian envelope with 83% of the energy within a circular region 4° in diameter. Spatial bandwidth was 1.5 octave.
![]() | (1) |
(1.5°) is the space constant,
(0.375 cyc/deg) is the spatial frequency, and C is the contrast. On the HFA, there were 54 locations separated by 6° in an offset grid over the central visual field: 27° nasally and 21° temporally, superiorly, and inferiorly (Fig. 1) .
|
Threshold Algorithm
For the HFA, stimulus presentation and threshold estimation were determined with the Swedish Interactive Thresholding Algorithm (SITA).15 For the VSG, stimulus presentation was determined by a staircase method. At each location in the visual field, the first Gabor patch was presented with the maximum contrast of 100%. The first size III stimulus was presented on the screen with 900% contrast. In the subsequent presentations at that location, the contrast of the stimulus was varied with an adaptive staircase16 procedure driven by the subjects responses and controlled by the computer. The contrast was decreased in steps of 0.3 log unit until the subject did not respond (first reversal); then the contrast of the stimulus was increased by 0.15-log-unit steps until the subject responded (second reversal); the rest of the staircase used steps of 0.15 log unit. Threshold for each location was calculated as the mean of the reversals at a 0.15-log unit step size.
On the VSG, at the locations where the patient did not respond to stimuli at the maximum contrast, the stimuli were presented at least seven more times at maximum contrast. If no response was obtained after eight stimulus presentations, then a psychometric function could not be fit, and the sensitivity was recorded as "not seen." In addition, sensitivities of 0.11 log unit or smaller were recorded as "not seen."
When the threshold is near the maximum contrast available, the estimates of fluctuation cannot be accurate, therefore cutoff limits were established for excluding data for study of variability. For the VSG data, locations scored as "not seen" were excluded. For the HFA data, if at least one of the two thresholds was in the range of 0 to 11 dB, then the location was excluded. For depth of defect comparison and analysis of structure-function relations, only unreliable points were excluded.
Limits for reliability of data from the HFA were set by the rate of fixation loss and the false-negative and -positive rates. In addition, the HFA data from locations above and below the blind spot were considered unreliable and were excluded from data analysis. The VSG allowed additional reliability indices.
Fixation stability for both VSG and HFA was monitored by observing the participant directly through a closed-circuit video system and by using the Heijl-Krakau method.17 The first author was present during all the testing, and there was continuous observation of the subjects fixation. If the subject started moving the eye under examination, he or she was instructed to maintain fixation as much as possible.
In the Heijl-Krakau method, stimuli are periodically presented throughout the experiment at the physiological blind spot. A positive response to a stimulus presented at the blind spot indicates either the patients poor fixation or the patients tendency to press the button without seeing a stimulus. During this study, there were two instances when a subjects responses to these stimuli were greater than 20%. These subjects were retested on a different day, and if reliable visual field results were obtained, then the subjects were included in the study. One young control subject was excluded based on a high rate of fixation losses, and one patient with glaucoma was retested; during the second test, reliable data were obtained and used in analysis.
On the HFA the location of the center of the blind spot was determined by projecting the stimulus at the standard location of the blind spot and if the subject did not respond to the stimulus, then this location was taken as an individuals blind spot. If an individual responded to at least two stimuli presented at the assumed blind spot, then the testing was suspended and the blind spot was mapped.
On the VSG, the center of the blind spot was estimated for each individual before testing. The blind spot algorithm started testing at 6° temporally and 1° below horizontal midline. A 0.5° diameter circular stimulus was presented at the maximum contrast. At first, the stimulus was presented in 1° increments moving temporally until the subject stopped responding and then began responding again at more temporal locations. The midpoint between the two edges of the "nonseeing" region was considered to be the horizontal center of the blind spot. The vertical parameters were then determined in a similar manner by presenting the stimulus along the line perpendicular to the horizontal center of the blind spot. On the VSG, to evaluate fixation stability during an experiment, Gabor patches or size III stimuli (during the testing of the young subjects with size III stimuli) were presented periodically at the individuals blind spot. Gabor patches projected at the blind spot had most of the energy concentrated over the central 2°, with peak spatial frequency of 0.5 cyc/deg and spatial bandwidth of 3.0 octaves.
The false-positive error score characterizes the tendency of the participant to press the button even when there is no stimulus presented. For the VSG system, the false-positive rate was calculated as the number of responses to blank stimuli divided by the number of blank stimuli presented. For SITA, the false-positive rate was determined based on the number of times the patient pressed a button when it was not expected.18 A false-positive rate higher than 20% was used to exclude data. There were no points excluded based on the false-positive criterion.
For the HFA, the false-negative rate represents the number of times that a subject fails to respond to a stimulus eight times brighter than the threshold that has been already determined at this test point location. A false-negative rate higher than 20% was used to exclude data obtained from normal subjects only. During this study, no visual fields were excluded based on this criterion.
Visual field results of patients with glaucoma who had high false-negative rates were not considered unreliable.13 19 20 21 22 A higher false-negative rate in eyes with glaucomatous field loss compared with unaffected eyes may be due to the increased variability in threshold values typically found in such eyes.19 Only two patients had false-negative rates greater than 20%, and when their data were removed, the results did not change.
Data from the VSG were analyzed using five parameters as indices of reliability: rate of fixation loss, false-positive rate, false-negative rate, slope of the psychometric function, and a distribution parameter. The rates of fixation loss and false-positive responses were calculated in the same manner as for the HFA, as described earlier.
The false-negative rate on the VSG system was calculated for each threshold at a particular location separately from other locations, and therefore only that individual threshold was rejected as unreliable in both normal and glaucomatous visual fields if it was 20% or higher. This method is less affected by overall visual field loss and yields more accurate estimates of the false-negative rate than that on HFA.23 The VSG system used a maximum-likelihood method that allows the slope of the psychometric function to vary across locations. When the slope of the psychometric function was less then 1.0, the datum for this particular location was considered unreliable. The slope reflected the overall reliability of the particular threshold and therefore thresholds with shallow slope were excluded as unreliable based on multiple criteria, except for one point of a patient with glaucoma. In the analysis of the age-related normative data, there were 16 points (6% of the data) with high false-negative rates, and in the data of the patients with glaucoma, there were 57 points (10% of the total).
The distribution parameter for each test location was the difference between the mean of the reversals (with steps of 0.15 log unit) and the threshold estimated by the maximum likelihood method. If stimulus contrast in most of the trials for a staircase is higher than the contrast at the 50% seen point of the psychometric function, threshold estimates become less accurate.23 The distribution parameter can be used to identify such staircases, and if it was higher than 0.15 log unit, then we considered that staircase unreliable because of failure to provide the appropriate distribution of stimulus contrasts. In the analysis of the age-related normative data, there were five points excluded based on this criterion.
Protocol
Data were collected at two visits, approximately 1 week apart, and were used to assess both learning effect and test-retest variability.24 The assessment was accomplished by modeling the learning effect as homogenous fluctuation and test-retest variability as heterogeneous fluctuation. Learning effects cause an overall increase in sensitivity from one test to the next (homogeneous fluctuation), as the person becomes more familiar with the task, and were modeled as the average change in sensitivity (in log units) across all locations tested. Test-retest variability is due to the psychophysical algorithm interacting with the subjects frequency-of-seeing curve and was modeled by allowing sensitivity at each location to vary independently (heterogeneous fluctuation).
During the first visit, contrast sensitivity and visual acuity were recorded using Pelli-Robson contrast sensitivity charts and ETDRS logMAR acuity charts.
During each visit, every young participant was tested on the VSG system with Gabor patches (eight and four reversals) and size III stimuli (four reversals), to estimate the test duration and test-retest variability for each set of reversals. Each older control subject and patient with glaucoma was tested on the VSG with Gabor stimuli (six reversals) and on the HFA with the SITA Standard 24-2 visual field test. At the end of the second visit, after the visual field test, patients with glaucoma had an optic nerve head assessment on the HRT II, as an index of structural abnormality.
Statistical Design
To assess variability with the new and standard perimetric stimuli, homogeneous and heterogeneous fluctuations were calculated. Heterogeneous fluctuation reflects variability of the individual threshold estimates and was calculated for every location as the absolute value of the difference in the log contrast sensitivities for the first and second visits. Homogeneous fluctuation reflects an overall change in sensitivity across two visits and was computed for each subject as the mean of differences in log contrast sensitivities between the first and second tests across all the locations. Homogeneous fluctuation sets an upper limit to the effect of change in the subjects criterion from one day to another, and hence sets a lower limit for heterogeneous fluctuation.
We asked several questions in the study:
Data were gathered from young control subjects to evaluate the effect of the number of reversals on test-retest variability. An F-test and two linear regressions were used, and statistical significance was set at P < 0.0125. The effect of the number of reversals on variability for Gabor patches with eight- and four-reversal staircases was assessed by using an F-test on heterogeneous fluctuation. To assess the dependence of variability on sensitivity, linear regression was performed on heterogeneous fluctuation versus sensitivity for the two staircases, and a z-score analysis was used to compare the slopes.
Homogeneous fluctuation was plotted against mean contrast sensitivity for each subject. The Bland-Altman analysis of agreement24 was used to compare learning effects for CSP and CAP.
To compare the variability of responses of patients with glaucoma to the Gabor and size III stimuli, we applied an F-test, linear regression, and slopes comparison; statistical significance was set at P < 0.0125. An F-test was performed on homogeneous fluctuation for the two devices. We assessed the relation between sensitivity and variability by performing linear regression on heterogeneous fluctuation, and the slopes were compared in z-score analyses.
Bland-Altman analysis of agreement was used to estimate whether the depth of defect obtained from testing patients on the VSG system is similar to that obtained from testing on the HFA. The two tests could not be compared on a point-by-point basis because different locations were tested on the two devices. The superior and inferior nasal quadrants were extensively tested on both of the tests and were used for this analysis. The sensitivity was converted into linear units and was averaged across all the points in the given quadrant for each patient for each test day, and this average was expressed in log units. Confidence limits were used to determine the range of fluctuations expected solely from test-retest variability. These were computed as 1.96 times the vector sum of standard deviations of homogeneous fluctuation for each test separately. To compensate for the difference in dynamic range on two devices, for CAP we assigned the values of depth of defect that were deeper than –1.34 log units the value of –1.34 log units. This value represented the deepest possible defect that could be measured with CSP on this particular group of subjects.
Bland-Altman analysis was used to determine the relations between contrast sensitivity and rim area. Contrast sensitivity from the two test dates was averaged across the nasal hemifields for each device and for every patient. Possible bias due to the different units used on different devices was minimized by expressing results as percentages of the mean normal.
The effect of eccentricity on contrast sensitivity in normal eyes was assessed using data from young and old control subjects tested with Gabor patches (Fig. 2) and size III stimuli with four-reversal staircases. Linear regression was performed on mean contrast sensitivity versus eccentricity, and z-scores were used to compare the slopes. Statistical significance was set at P < 0.05.
|
| Results |
|---|
|
|
|---|
|
|
|
|
|
|
| Discussion |
|---|
|
|
|---|
The present study extended the findings of Pan et al.6 by creating a clinically useful test. Pan et al. confirmed the findings of Pearson et al.5 with chromatic stimuli and extended their approach by using Gabor patches.6 However Pan et al. were restricted to using only eight locations, to allow the use of 12-reversal staircases without fatiguing the subjects. By retrospectively analyzing their data, they found that low variability could also have been achieved with Gabor stimuli with fewer reversals. In the present study, data from the young control subjects were used to assess the effect of the number of reversals on tests using a much larger number of locations than used by Pan et al. Although a large number of reversals can improve psychophysical threshold estimation, with the large number of locations tested in perimetry any increase in test duration can introduce fatigue and actually increase variability. Mean test duration (± SD) for the eight-reversal staircase was 14.0 ± 0.6 minutes, almost twice that for clinical perimetry, and did not yield significantly lower variability than the four-reversal staircase (5.0 ± 0.5 minutes). However, the dependence of variability on sensitivity was significantly lower in the eight-reversal staircase. To decrease test duration and to maintain low variability, older control subjects and patients with glaucoma were tested with a six-reversal staircase.
CAP with a size III stimulus has a dynamic range of 2.5 to 3.5 log units across locations and ages. With CSP, we found a dynamic range of 1.0 to 1.3 log units. There were two methods used to minimize potential statistical artifacts due to the difference in dynamic ranges. First, for the study of variability, with CAP all locations with sensitivities of 11 dB or smaller were scored as "not seen" and were excluded. Second, in the study of structure-function relations, all measures were converted to the percentage of mean normal.
A potential advantage of CSP over CAP is that the dynamic range appears to show less variation across subjects and locations. Figure 8 demonstrates the smaller eccentricity effects with the Gabor patches than with the size III stimuli, which is consistent with results in previous studies.25
With both CSP and CAP, homogeneous fluctuation showed a mean learning effect of less than 0.1 log unit (Fig. 5) , whereas heterogeneous fluctuation was often greater than 0.4 log unit (Fig. 3) . Because mean learning effects were small, we were readily able to detect the dependence of test-retest variability on sensitivity for the size III stimulus (Fig. 3) . Test-retest variability correlated with sensitivity for the size III stimulus but not for the Gabor patches (Fig. 3) . This finding cannot be attributed to learning effects, because mean test-retest variability with the Gabor stimuli was lower than with the size III stimuli.
The correlations between test-retest and sensitivity assessed in this study are consistent with those in a recent report by Artes et al.,7 who compared perimetric sensitivities for the 24-2 pattern of test locations with CAP and with the Humphrey FDT Matrix (Carl Zeiss Meditec, Inc.), which uses as a stimulus a 0.5-cyc/deg grating in a square region 5° across. They found that grating sensitivity had low variability without a significant dependence on sensitivity, just as we found with CSP (Fig. 5) . Artes et al. found that sensitivities to the two stimuli had a constant ratio over the first log unit of contrast sensitivity, and we found that the average effect of glaucomatous loss was similar on CSP and CAP (Fig. 6) . They found that when sensitivity was reduced by more than a log unit, sensitivity to a large sinusoidal stimulus was not reduced as much as sensitivity to the size III stimulus used in CAP, consistent with our results for patients with extensive CAP defects.
Garway-Heath et al.26 found a linear relation between area of the neuroretinal rim and perimetric sensitivity. The present study extended their finding by using the Bland-Altman analysis to quantify confidence limits and to assess systematic differences. The relation between rim area and perimetric sensitivity was similar for CAP and CSP for both upper and lower nasal quadrants. With both tests, the mean difference between rim and perimetric measured was close to zero when expressed as a percentage of mean normal. However, the slope of the regression line for Bland-Altman analysis was significantly shallower for CSP than for CAP. Our results indicate that with CAP but not with CSP the loss in sensitivity in the areas of severe damage was greater than the loss of rim area.
The overall test-retest variability was similar in CAP and CSP, with CSP having less dependence on sensitivity than CAP. With CAP, test-retest variability was higher in areas with low sensitivity, which becomes problematic for monitoring response to treatment. The reduced variability with CSP in damaged areas may reduce this problem. Further studies are needed to assess this finding.
There were several collateral findings in the study that cannot be adequately evaluated due to the small number of subjects; therefore, further study is needed. For the younger and older control subjects, evaluation of eccentricitys effects on CSP versus CAP was not the primary goal, but we did find strong evidence that the effects of eccentricity were smaller with the Gabor stimulus from CSP than with the size III stimulus from CAP. We also found that in severely damaged areas, defect depth could be much greater with CAP than with CSP. This study included patients with a wide range of severity of glaucomatous damage, and so a larger sample of patients with advanced or end-stage glaucoma would be needed to evaluate this finding.
| Conclusion |
|---|
|
|
|---|
Further research is necessary to evaluate findings that were not part of the primary focus of this study. Our primary finding is that CSP has strong potential to be used clinically to improve monitoring progression and response to treatment.
| Footnotes |
|---|
Submitted for publication September 17, 2007; revised February 28, 2008; accepted May 22, 2008.
Disclosure: A. Hot, None; M.W. Dul, None; W.H. Swanson, Carl Zeiss Meditec (C)
The publication costs of this article were defrayed in part by page charge payment. This article must therefore be marked "advertisement" in accordance with 18 U.S.C.
1734 solely to indicate this fact.
Corresponding author: Aliya Hot, SUNY State College of Optometry, 33 W 42nd Street, New York, NY 10036; alichkahot{at}gmail.com.
| References |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |