A database of the visual field test results (SITA 24-2 algorithm, Humphrey Visual Field Analyzer; Carl Zeiss Meditec Inc., Dublin, CA) from 6696 eyes of 3586 patients with suspicious/diagnosed glaucoma (age 66.0 ± 13.0 years old, collected at Manchester Royal Eye Hospital) were classified into seven perimetric stages with the Brusini staging method.
19 Eyes with stage 0 were defined as normal group (
n = 2344) and eyes with stage 2 and 3 were selected as a defective group (
n = 2222) for further analysis. Eyes graded as borderline and stage 1 were excluded from the analysis, as there is considerable overlap in these grades between normal and early glaucomatous eyes. Eyes with advanced visual field loss (stages 4 and 5) were also excluded from the analysis on the basis that this level of loss does not offer a suitable diagnostic challenge. We included all eyes that met these criteria and did not exclude eyes on the basis of reliability indices.
Positive predictive value (PPV) is a statistical parameter used to estimate the performance of a diagnostic test. It is defined as the ratio between the true positives (TP) and total positive calls (false positives [FP] plus TP).
To reduce any bias that might result from a single sample of visual field data, we used a bootstrap method (Matlab, version 2008; Mathworks Inc., Natick, MA) to generate 200 datasets from the original sample. For each dataset, 4566 visual field tests were randomly selected, with replacement, from the original dataset. For each of the 200 datasets the PPV of the 54 test locations was calculated using a cutoff criterion of pattern deviation probability less than 0.01 to define a location as being normal or defective. For each of the 200 datasets, the location with the highest PPV was then identified and the most frequent location with the highest PPV was selected for the optimized test pattern. Visual field results of eyes in which this abnormal location was missing were then removed in each of the 200 datasets and the PPV of the residual sample calculated. This process was repeated until all defective eyes in the set had been detected (see
Fig. 1).
The diagnostic performance of the optimized test patterns (sensitivity and specificity) were then calculated at each step (i.e., for 1, 2, 3… n test locations), from the original dataset (n = 2344 normal and 2222 defective eyes). The sensitivity was defined as the proportion of the abnormal eyes detected with defects at the n test locations and the specificity was the proportion of normal eyes detected with no defects at the n test locations.
In addition, five random test patterns were generated to establish the benefits of using optimized distributions. The five random series of test locations were generated with a Matlab program (version 2008; Mathworks). Receiver operating characteristics (ROC) curves were plotted for both the PPV optimized pattern and randomized test patterns based on the sensitivity obtained at a specificity level of approximately 95%, 90%, 85%, 80%, 75%, 70%, and 60%.
To analyze the characteristics of visual field defects detected with increasing number of test locations, the average and SD of mean deviation (MD), pattern standard deviation (PSD), and total number of defective locations, were calculated for 1–10, 11–20, 21–30, and 31–40 optimized test locations from the original dataset. One-way analysis of variance (SPSS 16.0 for Windows; SPSS Inc., Chicago, IL) was used to compare means of those parameters between the defined location groups.