Journal of Glaucoma:
Combining Information From 3 Anatomic Regions in the Diagnosis of Glaucoma With Time-Domain Optical Coherence Tomography
Wang, Mingwu MD, PhD*; Lu, Ake Tzu-Hui PhD†; Varma, Rohit MD, MPH†; Schuman, Joel S. MD‡; Greenfield, David S. MD§; Huang, David MD, PhD∥; Advanced Imaging for Glaucoma Study Group
*Department of Ophthalmology and Vision Science, University of Arizona College of Medicine, Tucson, AZ
†Doheny Eye Institute and the Department of Ophthalmology, Keck School of Medicine, University of Southern California, Los Angeles, CA
‡University of Pittsburgh Medical Center, Pittsburgh, PA
§Bascom Palmer Eye Institute, University of Miami, Miami, FL
∥Casey Eye Institute, Oregon Health & Science University, Portland, OR
For the Advanced Imaging for Glaucoma Study Group, see http://http://www.AIGStudy.net.
Disclosure: Supported by NIH grants R01 EY013516 and P30 EY03040 and by a grant from Research to Prevent Blindness. D.H. and J.S.S. receive royalties from the Massachusetts Institute of Technology derived from an optical coherence tomography patent licensed to Carl Zeiss Meditec Inc. D.H. has a significant financial interest in Optovue Inc., a company that may have a commercial interest in the results of this research and technology. This potential individual conflict of interest has been reviewed and managed by Oregon Health & Science University. Other authors do not have financial interest in the subject of this article.
Reprints: David Huang, MD, PhD, Casey Eye Institute, Oregon Health & Science University, 3375S.W. Terwilliger Blvd., Portland, OR 97239-4197 (e-mail: email@example.com).
Received December 15, 2011
Accepted May 29, 2012
Purpose: To improve the diagnosis of glaucoma by combining time-domain optical coherence tomography (TD-OCT) measurements of the optic disc, circumpapillary retinal nerve fiber layer (RNFL), and macular retinal thickness.
Patients and Methods: Ninety-six age-matched normal and 96 perimetric glaucoma participants were included in this observational, cross-sectional study. Or-logic, support vector machine, relevance vector machine, and linear discrimination function were used to analyze the performances of combined TD-OCT diagnostic variables.
Results: The area under the receiver-operating curve (AROC) was used to evaluate the diagnostic accuracy and to compare the diagnostic performance of single and combined anatomic variables. The best RNFL thickness variables were the inferior (AROC=0.900), overall (AROC=0.892), and superior quadrants (AROC=0.850). The best optic disc variables were horizontal integrated rim width (AROC=0.909), vertical integrated rim area (AROC=0.908), and cup/disc vertical ratio (AROC=0.890). All macular retinal thickness variables had AROCs of 0.829 or less. Combining the top 3 RNFL and optic disc variables in optimizing glaucoma diagnosis, support vector machine had the highest AROC, 0.954, followed by or-logic (AROC=0.946), linear discrimination function (AROC=0.946), and relevance vector machine (AROC=0.943). All combination diagnostic variables had significantly larger AROCs than any single diagnostic variable. There are no significant differences among the combination diagnostic indices.
Conclusions: With TD-OCT, RNFL and optic disc variables had better diagnostic accuracy than macular retinal variables. Combining top RNFL and optic disc variables significantly improved diagnostic performance. Clinically, or-logic classification was the most practical analytical tool with sufficient accuracy to diagnose early glaucoma.
In glaucoma, between 30% and 50% of the ganglion cells may be lost before abnormalities appear in perimetric testing.1,2 Glaucomatous changes in the optic disc and retinal nerve fiber layer (RNFL) take place over time. Because less than optimal agreement has been reported in subjective assessment of optic disc photographs performed by different observers and even among glaucoma specialists,3 objective, quantitative measurement of ocular structures implicated in glaucoma is valuable.4 Optical coherence tomography (OCT) is a commonly used posterior segment imaging technology that has sufficient resolution to measure macular retinal thickness, RNFL thickness, and optic disc dimensions.5 Many studies have demonstrated the value of measuring these anatomic regions with OCT for the diagnosis and monitoring of glaucoma.6–9
Methodologies are still being refined to enhance the sensitivity and specificity for the detection of glaucomatous changes by OCT.10–15 With data from age-matched participants enrolled in the Advanced Imaging for Glaucoma Study (AIGS), we found that in distinguishing normal eyes from glaucomatous eyes the sensitivity and specificity were optimized with an “or-logic” combination of the 3 best RNFL variables.16 Using the area under the receiver-operating curve (AROC) and linear discriminant function (LDF), Medeiros et al17 demonstrated that a combination of optic disc and RNFL variables obtained by OCT improved the diagnostic accuracy for glaucoma detection.
In the present study, we used available AIGS data to explore the relative strengths of or-logic combination, support vector machine (SVM),18,19 relevance vector machine (RVM),20,21 and LDF approaches in analyzing OCT diagnostic variables. Measurements from 3 anatomic regions, the peripapillary RNFL, optic disc, and macular retina, were combined to increase the diagnostic accuracy of time-domain OCT (TD-OCT) in detecting early glaucomatous change. To the best of our knowledge, this is the first time in a single study that various statistical processing strategies have been tested in such an application.
MATERIALS AND METHODS
Study Population and Database
The AIGS is a multicenter bioengineering partnership and clinical study sponsored by the National Eye Institute. The designs and goals of AIGS were described previously.16,22 A more detailed description of the study protocols can be found in the AIGS Manual of Procedures, which can be downloaded from the Web site http://www.AIGStudy.net.
All study procedures adhered to the principles outlined in the Declaration of Helsinki for research involving human subjects. Written informed consent was obtained from all participants. Institutional review board and ethics committee approval was obtained in all participating institutions. Study participants are classified into normal (N) and perimetric glaucoma (PG) groups, according to criteria described previously.16,22 Only age-matched N and PG participants are used in this study. The ages of all participants were between 40 and 79 years at baseline. We performed age-matching selection, as previously described,16 to prevent bias from age differences between the N and PG groups.
One aim of the AIGS is to evaluate glaucoma diagnostic accuracy using quantitative imaging instruments. The diagnostic variables evaluated in this paper were fast peripapillary RNFL thickness, optic disc variables (described below), and macular retinal thickness obtained by TD-OCT (Stratus OCT system; software version 4.0; Carl Zeiss Meditec Inc., Dublin, CA). Fast Stratus OCT scans of the optic disc and macular retina were generated from six 6 mm linear scans in a spoke-like radial configuration, each of which were 30 degrees apart. Optic disc data were generated from the automated determination of the disc margin as defined by the Stratus OCT software, which included horizontal integrated rim width (HIRW), vertical integrated rim area (VIRA), vertical cup/disc ratio (VCDR), horizontal cup/disc ratio, rim area, cup/disc area ratio, and cup area. By convention of the Stratus OCT software, the central 6 mm diameter macular region (20.8 degrees) was divided into 4 quadrants: superior, inferior, nasal, and temporal. Each quadrant was further subdivided into inner and outer regions by a circumferential demarcation line 1.5 mm from the center of the macula. The mean macular retinal thickness was calculated as the weighted average of the sectoral macular thickness measurements given automatically by the Stratus. Performance of peripapillary RNFL thickness scans was as described previously.16
The fast scans of peripapillary RNFL thickness, optic disc, and macular retinal thickness were performed twice on the same day and the resulting measurements averaged. The photographers operating the Stratus OCT system were instructed to obtain all scans with signal strength scores of >7 if possible. Scans with signal strength <6 were excluded from analysis. Data from both eyes were used. Right and left eye clock-hour data were analyzed together based on the assumption of mirror-image symmetry.
Statistical Classification Analysis
We sought to combine diagnostic variables across anatomic regions to boost the diagnostic power by 2 types of classifiers: or-logic and machine learning classifiers. Both eyes of participants were used in the analysis whereas the intereye correlation was appropriately handled using available statistical approaches as we previously reported.16 In brief, the t tests were adjusted by the generalized estimating equation approach,23 and estimates of AROC, sensitivity, and variance were incorporated with robust variance-covariance estimates (Huber-White sandwich estimator).24,25
The data were analyzed to characterize the participants and assess the diagnostic power of OCT single and combined variables. Participant characteristics were compared between the N and PG groups using 2-tailed t tests for continuous variables and χ2 tests for categorical variables. The performance of OCT variables was assessed by AROC for all variables along with sensitivity at 95% specificity for combined variables. Estimates and comparisons of AROC were computed based on the Obuchowski method26 accounting for intereye correlation. This method is a generalization of the Delong method that has generally been used in other eye studies.27 For comparison of sensitivities, we used a generalized McNemar test to account for intereye correlation in estimation.28
A diagnosis of glaucoma was made if any of the component variables was abnormal based on or-logic classification. In our previous study we derived a composite score based on the minimum standardized deviate values to construct the AROC for or-logic classification.16 Because this method is only valid for Gaussian distributed variables, we had to generalize our approach so that it could be applied to any distribution. The generalization was done by transforming the value of each variable into the cumulative density function (CDF) based on the distribution of the N eyes. Thus, the composite score was formulated by the minimum CDF rather than the minimum standardized deviate. For example, overall, inferior, and superior RNFL thicknesses fit Gaussian distributions. If the composite score (based on previous derivation) of a particular eye had standardized deviate values of −1.75, −1.65, and −1.55 from the N group, then the CDF values would be 0.04, 0.05, and 0.06, respectively. The minimum 0.04 was used as the or-logic classifier to construct the AROC, and the eye would be diagnosed as a PG eye at the level of 95% specificity. The CDF calculation is valid for any parametric distribution and is more straightforward to classify an eye for any given cutoff value. Each OCT diagnostic variable was identified with well-known statistical distributions such as Gaussian, gamma, beta, and others to construct the CDF in or-logic classifier.
Machine Learning Classifiers
As they have been widely used in other eye studies,17,29–31 we selected 3 machine learning classifiers, LDF, SVM, and RVM, to classify N and PG eyes. To assess the performance of each approach, we used 16-fold cross-validation to train and test all classifiers. In k-fold cross-validation, each fold was tested by the model trained by the other (k−1) fold and the entire procedure was repeated k times. This way, the composite score of each eye was generated by the model constructed from independent observations to yield unbiased assessment. Such cross-validation was not necessary for or-logic classifier because it does not involve any optimized procedure.
In the LDF approach, several diagnostic variables were multiplied by different weight coefficients and then formulated through a linear combination. The weights were computed so that the summed variable was optimized for discrimination between N and PG groups. To allow nonlinear combination, SVM and RVM with a Gaussian kernel were used in this study. SVM sought to draw a boundary to maximize the “separation” of diagnostic variables between N and PG eyes and searched supporting vectors located on the boundary for model fitting. Even in relaxing the constraint of linear relationship, it is common to observe nonseparable data in high dimensional data. In this case, SVM assigns a parameter (generally denoted by C) to penalize the errors. With regard to the penalty for outliers and the relaxation of linear relationship, this approach yields more advantages over LDF. RVM used a model of identical functional form used under the SVM along with Bayesian framework to identify relevance vectors rather than the supporting vectors. The classifier generated a 0 to 1 composite score reflecting glaucoma probability, which would render to clinicians a more intuitive interpretation for glaucoma diagnosis. In addition, the Bayesian machine learning classifier has been used in the HRT3 scanning laser tomography machine to generate a glaucoma probability score.32,33
To account for the intereye correlation, the LDF was generalized with robust variance-covariance estimates.16 One eye was randomly selected from each participant in the training process of cross-validation under the SVM and RVM approaches as they required independent eyes to generate the support and relevance vectors, respectively. Then the vectors were applied to both eyes of participants to generate the composite scores in the test process. A set of predictors was selected for classification analysis based on the set finalized in the or-logic classification because our previous study successfully selected a sufficient set of variables in combining the RNFL variables.
Statistical significance was accepted at P<0.05. All analyses were done in SAS 9.1 and MATLAB 7.0. The MATLAB codes were freely available from http://asi.insa-rouen.fr/enseignants/~arakotom/toolbox/index.html for SVM and from http://www.miketipping.com/index.php?page=rvm for RVM.
Analyses were performed on 96 N (184 eyes) and 96 PG (139 eyes) age-matched participants that were selected from a database of 111 N (214 eyes) and 130 PG (188 eyes) participants. There were no significant differences between the N and PG groups in terms of age and sex; however, there were more whites in the N group (Table 1). In the PG group, 95 eyes (68.4%) had a mean deviation (MD) ≥−6.0 dB (early glaucoma), 31 eyes (22.3%) had a MD between −6.01 and −12.0 dB (moderate glaucoma), and 13 eyes (9.3%) had a MD≤−12 dB (advanced glaucoma).
All of the diagnostic variables were significantly different between N and PG groups (Table 2). All P values determined by the generalized estimating equation t tests were <0.0001 except 2 macular scans that had P values ≤0.01.
We analyzed AROCs for diagnostic variables from optic disc, peripapillary RNFL, and macular retina (Table 3). We included only RNFL variables from the superior quadrant, inferior quadrant, and overall average RNFL thickness because the nasal and temporal quadrants have low diagnostic accuracy.16 The 5 best diagnostic variables are optic disc HIRW, optic disc VIRA, inferior quadrant RNFL thickness, overall average RNFL thickness, and optic disc VCDR. The differences of AROCs between the top 5 variables were small and not significant. All macular retinal variables had relatively poor diagnostic accuracy. Among them, even the inferior outer macular thickness, which has the highest AROC value, performed significantly worse than the top 4 variables (P≤0.01) from optic disc and RNFL thickness.
To evaluate the diagnostic performance of or-logic classifier, we chose the top 3 RNFL variables (superior quadrant, inferior quadrant, and overall average) and the top 3 optic disc variables (HIRW, VIRA, and VCDR). The highest AROC values were obtained when applying or-logic combination to all of the top 6 variables (AROC=0.946, Table 4). That is, an eye would be diagnosed as having glaucoma if any of the 6 variables were abnormal. Such or-logic combination has a higher AROC than applying the same classifier to the top 3 RNFL variables (AROC=0.928, P=0.04) or the top 3 optic disc variables (AROC=0.916, P=0.07). The or-logic combinations that required any 2 variables to be abnormal for an abnormal classification did not perform as well as or-logic that required only 1 variable to be abnormal, although the difference was not statistically significant. Requiring more component variables to be abnormal further decreased AROC values (not listed in table).
Table 5 lists AROC values along with sensitivity at 95% specificity for the best combinations of optic disc and RNFL diagnostic variables. Machine learning classifiers and or-logic classification were used to compare with the single diagnostic variable HIRW that had the highest AROC. To compute the sensitivity of the or-logic, we fixed the specificity for each of the 6 variables at 99.3% to reach an overall 95% specificity. Compared with the best single variable HIRW, combining diagnostic variables achieved a significantly higher glaucoma diagnostic accuracy as measured by the AROC. Indeed, all of the 19 single variables (Table 3) had significantly (P<0.05) worse AROCs than any of the combination variables (Table 5). The diagnostic sensitivity of the best single variable, HIRW, was 0.648, significantly lower (P<0.02) than the combination variables (Table 5), which all achieved sensitivities >0.73. Although SVM achieved the highest diagnostic power, there were no statistically significant differences among the machine learning and or-logic classifiers.
Or-logic analysis was also performed on the subgroup of PG subject with early glaucoma (MD>−6 dB). We obtained an AROC (SE) of 0.927 (0.017), with sensitivity (SE) of 0.642 (0.053) at 95% specificity.
With inherently faster acquisition time and better image resolution and repeatability,34–37 Fourier-domain OCT (FD-OCT) is rapidly replacing TD-OCT in both clinical and research settings. However, TD-OCT machines are still widely used clinically. Although comparative studies have demonstrated good correlation between measurements obtained from TD-OCT and FD-OCT, systemic differences in the 2 generations of instruments do not allow their output to be used interchangeably.38–42 TD-OCT is capable of obtaining objective, quantitative, and reproducible measurements of certain anatomic regions of the eye, including RNFL thickness, optic disc topography, and macular retinal thickness.43–46 Many studies have shown that there are no statistically significant differences between FD-OCT and TD-OCT in discriminating normal from glaucomatous eyes.35,41,47–50 Thus it is still worthwhile to find the best way to clinically utilize the diagnostic information obtained from TD-OCT by using the large data set accumulated in the AIGS. The AIGS also utilizes FD-OCT systems, and that data will be presented separately.
The AIGS previously found that a simple or-logic combination of the 3 best RNFL thickness variables (overall, superior quadrant, inferior quadrant) worked significantly better than single variable, and provided a simple practical approach to improving glaucoma diagnosis.16 Using an LDF approach, Medeiros et al17 demonstrated that a combination of selected optic disc and RNFL variables, when applied to an independent sample group, resulted in an AROC of 0.97 for glaucoma detection using the Stratus OCT. Bowd et al29 obtained enhanced differentiation of healthy from glaucomatous eyes by using RVM and SVM classifiers trained on RNFL thickness measurements obtained by scanning laser polarimetry. In this article, we used a larger data set and considered optic disc variables in addition to the RNFL variables. In addition, we used more sophisticated machine-learning classifiers to evaluate the diagnostic power of combining various OCT-based diagnostic variables. The selection of variables was based on a post hoc approach in addition to information obtained from our previous study.16
We found that optic disc variables had the highest diagnostic power (AROCs) for discriminating glaucomatous from healthy eyes, followed by RNFL and macular retinal thickness variables. For our combined diagnostic variables, we chose as components the single variables with AROCs>0.85, with the exception of cup/disc area ratio. We excluded that particular ratio because it was highly collinear with the VCDR. Highly collinear variables work poorly in LDF and other discriminant function approaches.51 Because of low AROCs, all macular variables were excluded in the combination. For FD-OCT, such combination would have likely been different because macular ganglion cell complex mapping with FD-OCT achieves high diagnostic accuracy.52
Similar to previous articles, we found that combining multiple diagnostic variables significantly improved diagnostic accuracy over any single variable. The most surprising finding of the present study is that simple or-logic worked as well as the much more sophisticated machine learning and statistical classifiers, even in a relatively large combination of 6 variables from both the disc and RNFL regions. The fact that or-logic worked so well reflects the heterogeneous patterns of anatomic damage in glaucoma, some having more damage superiorly, some inferiorly, and some in a diffuse global pattern. This finding has great practical implication. In contrast to SVM, RVM, and LDF approaches, the or-logic approach does not require special software, does not require large training and evaluation data sets for validation, and does not require additional Food and Drug Administration approval for such software and data sets. Using the or-logic approach, the practicing physician can identify at risk or glaucomatous eyes by simply looking at the most important anatomic variables measured by OCT. In the present study, the or-logic combination increased the AROC value by 0.037 and diagnostic sensitivity by 9%, compared with the best single variable, HIRW (P=0.03). When applying or-logic to the group with early glaucoma, we obtained an excellent AROC value (0.927), suggesting that this approach is robust and applicable to the early diagnosis of glaucoma.
We have advocated for the use of or-logic combinations in a previous article that studied only RNFL variables.16 The addition of 3 disc variables to the or-logic combination significantly improved AROC, showing that combining more variables reduced specificity only slightly but greatly increased the sensitivity of glaucoma detection. This worked particularly well with the diagnostic threshold of the component variables set near the first percentile level. Thus we recommend a practical diagnostic approach where glaucoma is suspected when any one of the 3 disc variables (HIRW, VIRA, VCDR) or RNFL thickness variables (overall, inferior, superior quadrant averages) is abnormal at the first percentile level. Of course, clinicians should not solely rely on an OCT printout to make the diagnosis of glaucoma given the limitations of OCT and the wealth of other historical and clinical information that need to be considered in glaucoma evaluation. The clinician should also assess the reliability of OCT variables by looking at the signal strength, RNFL and disc boundary detection by automated software, and the correlation of anatomic loss patterns with visual fields.
In summary, our present study demonstrated that the use of TD-OCT diagnostic information can be improved by combining the top diagnostic variables from the optic disc and peripapillary RNFL regions. Simple or-logic classification worked equally well as more sophisticated machine learning and statistical classifier approaches, and deserves further study in other advanced imaging modalities.
1. Quigley HA, Addicks EM, Green WR.Optic nerve damage in human glaucoma. III. Quantitative correlation of nerve fiber loss and visual field defect in glaucoma, ischemic neuropathy, papilledema, and toxic neuropathy.Arch Ophthalmol.1982;100:135–146.
2. Mikelberg FS, Yidegiligne HM, Schulzer M.Optic nerve axon count and axon diameter in patients with ocular hypertension and normal visual fields.Ophthalmology.1995;102:342–348.
3. Altangerel U, Bayer A, Henderer JD, et al..Knowledge of chronology of optic disc stereophotographs influences the determination of glaucomatous change.Ophthalmology.2005;112:40–43.
4. Vessani RM, Moritz R, Batis L, et al..Comparison of quantitative imaging devices and subjective optic nerve head assessment by general ophthalmologists to differentiate normal from glaucomatous eyes.J Glaucoma.2009;18:253–261.
5. Huang D, Swanson EA, Lin CP, et al..Optical coherence tomography.Science.1991;254:1178–1181.
6. Medeiros FA, Zangwill LM, Alencar LM, et al..Detection of glaucoma progression with stratus OCT retinal nerve fiber layer, optic nerve head, and macular thickness measurements.Invest Ophthalmol Vis Sci.2009;50:5741–5748.
7. Paunescu LA, Schuman JS, Price LL, et al..Reproducibility of nerve fiber thickness, macular thickness, and optic nerve head measurements using Stratus OCT.Invest Ophthalmol Vis Sci.2004;45:1716–1724.
8. Budenz DL, Michael A, Chang RT, et al..Sensitivity and specificity of the StratusOCT for perimetric glaucoma.Ophthalmology.2005;112:3–9.
9. Naithani P, Sihota R, Sony P, et al..Evaluation of optical coherence tomography and Heidelberg retinal tomography parameters in detecting early and moderate glaucoma.Invest Ophthalmol Vis Sci.2007;48:3138–3145.
10. Schulze A, Lamparter J, Pfeiffer N, et al..Diagnostic ability of retinal ganglion cell complex, retinal nerve fiber layer, and optic nerve head measurements by Fourier-domain optical coherence tomography.Graefes Arch Clin Exp Ophthalmol.2011;249:1039–1045.
11. Bowd C, Hao J, Tavares IM, et al..Bayesian machine learning classifiers for combining structural and functional measurements to classify healthy and glaucomatous eyes.Invest Ophthalmol Vis Sci.2008;49:945–953.
12. Mwanza JC, Oakley JD, Budenz DL, et al..Ability of cirrus HD-OCT optic nerve head parameters to discriminate normal from glaucomatous eyes.Ophthalmology.2011;118:241–248.
13. Kaushik S, Singh Pandav S, Ichhpujani P, et al..Retinal nerve fiber layer measurement and diagnostic capability of spectral-domain versus time-domain optical coherence tomography.Eur J Ophthalmol.2011;21:566–572.
14. Bizios D, Heijl A, Hougaard JL, et al..Machine learning classifiers for glaucoma diagnosis based on classification of retinal nerve fibre layer thickness parameters measured by Stratus OCT.Acta Ophthalmol.2010;88:44–52.
15. Huang JY, Pekmezci M, Mesiwala N, et al..Diagnostic power of optic disc morphology, peripapillary retinal nerve fiber layer thickness, and macular inner retinal layer thickness in glaucoma diagnosis with Fourier-domain optical coherence tomography.J Glaucoma.2011;20:87–94.
16. Lu AT, Wang M, Varma R, et al..Combining nerve fiber layer parameters to optimize glaucoma diagnosis with optical coherence tomography.Ophthalmology.2008;115:1352–1357.
17. Medeiros FA, Zangwill LM, Bowd C, et al..Evaluation of retinal nerve fiber layer, optic nerve head, and macular thickness measurements for glaucoma detection using optical coherence tomography.Am J Ophthalmol.2005;139:44–55.
18. Boser BE, Guyon IM, Vapnik VN.A Training algorithm for optimal margin classifiers. In Proceedings of the 5th Annual ACM Workshop on Computational Learning Theory; ACM Press, Pittsburgh, PA:144–152.
19. Vapnik V.Statistical Learning Theory.1998.New York:Wiley.
20. Tipping M.Sparse Bayesian learning and the relevance vector machine.J Mach Learn Res.2001;1:211–244.
21. Tipping MSolla SA, Leen TK, Müller K-R.The relevance vector machine.Advances in Neural Information Processing Systems.2000.Cambridge, MA:MIT Press: Vol 12;652–658.
22. Tan O, Li G, Lu AT, et al..Mapping of macular substructures with optical coherence tomography for glaucoma diagnosis.Ophthalmology.2008;115:949–956.
23. Liang KY, Zeger SL.Longitudinal data analysis using generalized linear models.Biometrika.1986;73:13–22.
24. White H.A heteroskedasticity-consistent covariance matrix estimator and a direct test for heteroskedasticity.Econometrica.1980;48:817–830.
25. Williams RL.A note on robust variance estimation for cluster-correlated data.Biometrics.2000;56:645–646.
26. Obuchowski NA.Nonparametric analysis of clustered ROC curve data.Biometrics.1997;53:567–578.
27. DeLong ER, DeLong DM, Clarke-Pearson DL.Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach.Biometrics.1988;44:837–845.
28. Obuchowski NA.On the comparison of correlated proportions for clustered data.Stat Med.1998;17:1495–1507.
29. Bowd C, Medeiros FA, Zhang Z, et al..Relevance vector machine and support vector machine classifier analysis of scanning laser polarimetry retinal nerve fiber layer measurements.Invest Ophthalmol Vis Sci.2005;46:1322–1329.
30. Zangwill LM, Chan K, Bowd C, et al..Heidelberg retina tomograph measurements of the optic disc and parapapillary retina for detecting glaucoma analyzed by machine learning classifiers.Invest Ophthalmol Vis Sci.2004;45:3144–3151.
31. Townsend KA, Wollstein G, Danks D, et al..Heidelberg retina tomograph 3 machine learning classifiers for glaucoma detection.Br J Ophthalmol.2008;92:814–818.
32. Swindale NV, Stjepanovic G, Chin A, et al..Automated analysis of normal and glaucomatous optic nerve head topography images.Invest Ophthalmol Vis Sci.2000;41:1730–1742.
33. Coops A, Henson DB, Kwartz AJ, et al..Automated analysis of Heidelberg retina tomograph optic disc images by glaucoma probability score.Invest Ophthalmol Vis Sci.2006;47:5348–5355.
34. Savini G, Carbonelli M, Barboni P.Spectral-domain optical coherence tomography for the diagnosis and follow-up of glaucoma.Curr Opin Ophthalmol.2011;22:115–123.
35. Leung CK, Cheung CY, Weinreb RN, et al..Retinal nerve fiber layer imaging with spectral-domain optical coherence tomography: a variability and diagnostic performance study.Ophthalmology.2009;116:1257–12631263e1-e2.
36. Gabriele ML, Ishikawa H, Wollstein G, et al..Optical coherence tomography scan circle location and mean retinal nerve fiber layer measurement variability.Invest Ophthalmol Vis Sci.2008;49:2315–2321.
37. Mwanza JC, Chang RT, Budenz DL, et al..Reproducibility of peripapillary retinal nerve fiber layer thickness and optic nerve head parameters measured with cirrus HD-OCT in glaucomatous eyes.Invest Ophthalmol Vis Sci.2010;51:5724–5730.
38. Vizzeri G, Weinreb RN, Gonzalez-Garcia AO, et al..Agreement between spectral-domain and time-domain OCT for measuring RNFL thickness.Br J Ophthalmol.2009;93:775–781.
39. Gonzalez-Garcia AO, Vizzeri G, Bowd C, et al..Reproducibility of RTVue retinal nerve fiber layer thickness and optic disc measurements and agreement with Stratus optical coherence tomography measurements.Am J Ophthalmol.2009;147:1067–1074.
40. Sung KR, Kim DY, Park SB, et al..Comparison of retinal nerve fiber layer thickness measured by cirrus HD and Stratus optical coherence tomography.Ophthalmology.2009;116:1264–1270.
41. Cho JW, Sung KR, Hong JT, et al..Detection of glaucoma by spectral domain-scanning laser ophthalmoscopy/optical coherence tomography (SD-SLO/OCT) and time domain optical coherence tomography.J Glaucoma.2011;20:15–20.
42. Knight OJ, Chang RT, Feuer WJ, et al..Comparison of retinal nerve fiber layer measurements using time domain and spectral domain optical coherent tomography.Ophthalmology.2009;116:1271–1277.
43. Budenz DL, Fredette MJ, Feuer WJ, et al..Reproducibility of peripapillary retinal nerve fiber thickness measurements with stratus OCT in glaucomatous eyes.Ophthalmology.2008;115:661–666.
44. Leung CK, Cheung CY, Weinreb RN, et al..Evaluation of retinal nerve fiber layer progression in glaucoma: a comparison between the fast and the regular retinal nerve fiber layer scans.Ophthalmology.2011;118:763–767.
45. Bengtsson B, Andersson S, Heijl A.Performance of time-domain and spectral-domain optical coherence tomography for glaucoma screening.Acta Ophthalmol.2012;90:301–305.
46. Lee EJ, Kim TW, Weinreb RN, et al..Trend-based analysis of retinal nerve fiber layer thickness measured by optical coherence tomography in eyes with localized nerve fiber layer defects.Invest Ophthalmol Vis Sci.2011;52:1138–1144.
47. Kim NR, Kim JH, Kim CY, et al..Comparison of the optic nerve imaging by time-domain optical coherence tomography and Fourier-domain optical coherence tomography in distinguishing normal eyes from those with glaucoma.J Glaucoma.2011[Epub ahead of print].
48. Moreno-Montanes J, Olmo N, Alvarez A, et al..Cirrus high-definition optical coherence tomography compared with Stratus optical coherence tomography in glaucoma diagnosis.Invest Ophthalmol Vis Sci.2010;51:335–343.
49. Sehi M, Grewal DS, Sheets CW, et al..Diagnostic ability of Fourier-domain vs time-domain optical coherence tomography for glaucoma detection.Am J Ophthalmol.2009;148:597–605.
50. Jeoung JW, Park KH.Comparison of cirrus OCT and Stratus OCT on the ability to detect localized retinal nerve fiber layer defects in preperimetric glaucoma.Invest Ophthalmol Vis Sci.2010;51:938–945.
51. Hastie T, Buja A, Tibshirani R.Penalized discriminant analysis.Ann Stat.1995;23:73–102.
52. Tan O, Chopra V, Lu AT, et al..Detection of macular ganglion cell loss in glaucoma by Fourier-domain optical coherence tomography.Ophthalmology.2009;116:2305–2314.
optical coherence tomography; glaucoma; imaging; image processing
© 2014 by Lippincott Williams & Wilkins.
Highlight selected keywords in the article text.