Patterns of functional vision loss in glaucoma determined with archetypal analysis

Glaucoma is an optic neuropathy accompanied by vision loss which can be mapped by visual field (VF) testing revealing characteristic patterns related to the retinal nerve fibre layer anatomy. While detailed knowledge about these patterns is important to understand the anatomic and genetic aspects of glaucoma, current classification schemes are typically predominantly derived qualitatively. Here, we classify glaucomatous vision loss quantitatively by statistically learning prototypical patterns on the convex hull of the data space. In contrast to component-based approaches, this method emphasizes distinct aspects of the data and provides patterns that are easier to interpret for clinicians. Based on 13 231 reliable Humphrey VFs from a large clinical glaucoma practice, we identify an optimal solution with 17 glaucomatous vision loss prototypes which fit well with previously described qualitative patterns from a large clinical study. We illustrate relations of our patterns to retinal structure by a previously developed mathematical model. In contrast to the qualitative clinical approaches, our results can serve as a framework to quantify the various subtypes of glaucomatous visual field loss.

[1]  Chris A Johnson,et al.  Classification of visual field abnormalities in the ocular hypertension treatment study. , 2000, Archives of ophthalmology.

[2]  J. Paetzold,et al.  A mathematical description of nerve fiber bundle trajectories and their variability in the human retina , 2009, Vision Research.

[3]  N. Oden,et al.  Morphology of the normal visual field in a population-based random sample: principal components analysis. , 1992, Statistics in medicine.

[4]  Donald C. Hood,et al.  A framework for comparing structural and functional measures of glaucomatous damage , 2007, Progress in Retinal and Eye Research.

[5]  Robert N Weinreb,et al.  Using unsupervised learning with variational bayesian mixture of factor analysis to identify patterns of glaucomatous visual field defects. , 2004, Investigative ophthalmology & visual science.

[6]  D B Henson,et al.  Spatial classification of glaucomatous visual field loss. , 1996, The British journal of ophthalmology.

[7]  G. Lindgren,et al.  Normal variability of static perimetric threshold values across the central visual field. , 1987, Archives of ophthalmology.

[8]  Robert N Weinreb,et al.  Patterns of glaucomatous visual field loss in sita fields automatically identified using independent component analysis. , 2009, Transactions of the American Ophthalmological Society.

[9]  Douglas J. Rhee,et al.  Soluble Guanylate Cyclase α1–Deficient Mice: A Novel Murine Model for Primary Open Angle Glaucoma , 2013, PloS one.

[10]  Terrence J. Sejnowski,et al.  Variational Learning of Clusters of Undercomplete Nonsymmetric Independent Components , 2003, J. Mach. Learn. Res..

[11]  Nomdo M. Jansonius,et al.  A mathematical model for describing the retinal nerve fiber bundle trajectories in the human eye: average course, variability, and influence of refraction, optic disc size and optic disc position. , 2012, Experimental eye research.

[12]  G. Trick,et al.  Assessing the utility of reliability indices for automated visual fields. Testing ocular hypertensives. , 1989, Ophthalmology.

[13]  Yuko Ohno,et al.  Properties of perimetric threshold estimates from Full Threshold, SITA Standard, and SITA Fast strategies. , 2002, Investigative ophthalmology & visual science.

[14]  Robert N Weinreb,et al.  Using unsupervised learning with independent component analysis to identify patterns of glaucomatous visual field defects. , 2005, Investigative ophthalmology & visual science.

[15]  J. A. Hartigan,et al.  A k-means clustering algorithm , 1979 .

[16]  Manuel J. A. Eugster,et al.  From Spider-man to Hero - archetypal analysis in R , 2009 .

[17]  Michael V. Boland,et al.  Risk factors and open-angle glaucoma: classification and application. , 2007, Journal of glaucoma.

[18]  Kevin P. Murphy,et al.  Machine learning - a probabilistic perspective , 2012, Adaptive computation and machine learning series.

[19]  Jonathan Denniss,et al.  An anatomically customizable computational model relating the visual field to the optic nerve head in individual eyes. , 2012, Investigative ophthalmology & visual science.

[20]  Georg Lindgren,et al.  A package for the statistical analysis of visual fields , 1987 .

[21]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[22]  Bill Bynum,et al.  Lancet , 2015, The Lancet.

[23]  J. Tukey,et al.  Multiple-Factor Analysis , 1947 .

[24]  J. A. López del Val,et al.  Principal Components Analysis , 2018, Applied Univariate, Bivariate, and Multivariate Statistics Using Python.

[25]  Janey L. Wiggs,et al.  The p53 Codon 72 PRO/PRO Genotype May Be Associated with Initial Central Visual Field Defects in Caucasians with Primary Open Angle Glaucoma , 2012, PloS one.

[26]  J Katz,et al.  Comparison of analytic algorithms for detecting glaucomatous visual field loss. , 1991, Archives of ophthalmology.

[27]  E. Oja,et al.  Independent Component Analysis , 2013 .

[28]  Trevor Hastie,et al.  The Elements of Statistical Learning , 2001 .

[29]  Robert N Weinreb,et al.  Association of CAV1/CAV2 genomic variants with primary open-angle glaucoma overall and by gender and pattern of visual field loss. , 2014, Ophthalmology.

[30]  Donald C. Hood,et al.  Glaucomatous Damage of the Macula as Revealed by Averaging Optical Coherence Tomography Data , 2012 .

[31]  H. Rootzén,et al.  A new generation of algorithms for computerized threshold perimetry, SITA. , 2009, Acta ophthalmologica Scandinavica.

[32]  C. Ji An Archetypal Analysis on , 2005 .

[33]  Chris A Johnson,et al.  Staging functional damage in glaucoma: review of different classification methods. , 2007, Survey of ophthalmology.

[34]  A. Azzalini A class of distributions which includes the normal ones , 1985 .

[35]  Sterling C. Hilton,et al.  Classifying visual field data. , 1996, Statistics in medicine.