Applying "Lasso" Regression to Predict Future Visual Field Progression in Glaucoma Patients.

PURPOSE We evaluated the usefulness of various regression models, including least absolute shrinkage and selection operator (Lasso) regression, to predict future visual field (VF) progression in glaucoma patients. METHODS Series of 10 VFs (Humphrey Field Analyzer 24-2 SITA-standard) from each of 513 eyes in 324 open-angle glaucoma patients, obtained in 4.9 ± 1.3 years (mean ± SD), were investigated. For each patient, the mean of all total deviation values (mTD) in the 10th VF was predicted using varying numbers of prior VFs (ranging from the first three VFs to all previous VFs) by applying ordinary least squares linear regression (OLSLR), M-estimator robust regression (M-robust), MM-estimator robust regression (MM-robust), skipped regression (Skipped), deepest regression (Deepest), and Lasso regression. Absolute prediction errors then were compared. RESULTS With OLSLR, prediction error varied between 5.7 ± 6.1 (using the first three VFs) and 1.2 ± 1.1 dB (using all nine previous VFs). Prediction accuracy was not significantly improved with M-robust, MM-robust, Skipped, or Deepest regression in almost all VF series; however, a significantly smaller prediction error was obtained with Lasso regression even with a small number of VFs (using first 3 VFs, 2.0 ± 2.2; using all nine previous VFs, 1.2 ± 1.1 dB). CONCLUSIONS Prediction errors using OLSLR are large when only a small number of VFs are included in the regression. Lasso regression offers much more accurate predictions, especially in short VF series.

[1]  M. Hubert,et al.  The Deepest Regression Method , 2002 .

[2]  V. Yohai.  HIGH BREAKDOWN-POINT AND HIGH EFFICIENCY ROBUST ESTIMATES FOR REGRESSION , 1987 .

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

[4]  M. Araie,et al.  A placebo-controlled 3-year study of a calcium blocker on visual field and ocular circulation in glaucoma with low-normal pressure. , 2008, Ophthalmology.

[5]  M. Aihara,et al.  Incidence of deepening of the upper eyelid sulcus on treatment with a tafluprost ophthalmic solution , 2014, Japanese Journal of Ophthalmology.

[6]  P. Rousseeuw,et al.  Robustness of Deepest Regression , 2000 .

[7]  H. Jampel,et al.  Natural history of normal-tension glaucoma. , 2001, Ophthalmology.

[8]  F. Fankhauser,et al.  Differential light threshold in automated static perimetry. Factors influencing short-term fluctuation. , 1984, Archives of ophthalmology.

[9]  Y. Sawada,et al.  Assessment of systemic adverse reactions induced by ophthalmic beta-adrenergic receptor antagonists. , 2001, Journal of ocular pharmacology and therapeutics : the official journal of the Association for Ocular Pharmacology and Therapeutics.

[10]  D P Crabb,et al.  Interobserver agreement on visual field progression in glaucoma: a comparison of methods , 2003, The British journal of ophthalmology.

[11]  C. Krakau.  A statistical trap in the evaluation of visual field decay , 1985, Acta ophthalmologica. Supplement.

[12]  R. Wilcox.  Introduction to Robust Estimation and Hypothesis Testing , 1997 .

[13]  C. Johnson,et al.  Simulation of longitudinal threshold visual field data. , 2000, Investigative ophthalmology & visual science.

[14]  M. Bruce Shields,et al.  Textbook of glaucoma , 1987 .

[15]  C. Mayama,et al.  The 5-year incidence of bleb-related infection and its risk factors after filtering surgeries with adjunctive mitomycin C: collaborative bleb-related infection incidence and treatment study 2. , 2014, Ophthalmology.

[16]  M. Aihara,et al.  Incidence of deepening of the upper eyelid sulcus after switching from latanoprost to bimatoprost , 2011, Japanese Journal of Ophthalmology.

[17]  Tzyy-Ping Jung,et al.  Progression of patterns (POP): a machine classifier algorithm to identify glaucoma progression in visual fields. , 2012, Investigative ophthalmology & visual science.

[18]  N M Jansonius,et al.  On the accuracy of measuring rates of visual field change in glaucoma , 2010, British Journal of Ophthalmology.

[19]  C. Krakau,et al.  REGRESSION ANALYSIS OF THE CENTRAL VISUAL FIELD IN CHRONIC GLAUCOMA CASES , 1982, Acta ophthalmologica.

[20]  Douglas R. Anderson.  Automated Static Perimetry , 1992 .

[21]  D B Henson,et al.  Influence of fixation accuracy on threshold variability in patients with open angle glaucoma. , 1996, Investigative ophthalmology & visual science.

[22]  Richard A. Russell,et al.  On alternative methods for measuring visual field decay: Tobit linear regression. , 2011, Investigative ophthalmology & visual science.

[23]  J Flammer,et al.  Differential light threshold. Short- and long-term fluctuation in patients with glaucoma, normal controls, and patients with suspected glaucoma. , 1984, Archives of ophthalmology.

[24]  P. J. Huber.  Robust Regression: Asymptotics, Conjectures and Monte Carlo , 1973 .

[25]  H. Theil.  A Rank-Invariant Method of Linear and Polynomial Regression Analysis , 1992 .

[26]  H. Jampel.  Target pressure in glaucoma therapy. , 1997, Journal of glaucoma.

[27]  F. Fitzke,et al.  Analysis of visual field progression in glaucoma. , 1996, The British journal of ophthalmology.

[28]  Kouros Nouri-Mahdavi,et al.  A method to measure and predict rates of regional visual field decay in glaucoma. , 2011, Investigative ophthalmology & visual science.

[29]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .

[30]  Trevor Hastie,et al.  Regularization Paths for Generalized Linear Models via Coordinate Descent. , 2010, Journal of statistical software.

[31]  B. Bengtsson,et al.  False-negative responses in glaucoma perimetry: indicators of patient performance or test reliability? , 2000, American journal of ophthalmology.

[32]  Anders Heijl,et al.  Prediction of glaucomatous visual field loss by extrapolation of linear trends. , 2008, Archives of ophthalmology.

[33]  Richard A. Russell,et al.  Detecting Changes in Retinal Function: Analysis with Non-Stationary Weibull Error Regression and Spatial Enhancement (ANSWERS) , 2014, PloS one.

[34]  Ted Maddess,et al.  Locally countable properties and the perceptual salience of textures. , 2013, Journal of the Optical Society of America. A, Optics, image science, and vision.

[35]  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.

[36]  Huseyin Seker,et al.  A hybrid dynamic Bayesian network approach for modelling temporal associations of gene expressions for hypertension diagnosis , 2014, 2014 36th Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[37]  Z. Bai,et al.  Asymptotic distributions of the maximal depth estimators for regression and multivariate location , 1999 .

[38]  Ryo Asaoka,et al.  Measuring Visual Field Progression in the Central 10 Degrees Using Additional Information from Central 24 Degrees Visual Fields and ‘Lasso Regression’ , 2013, PloS one.

[39]  Ryo Asaoka,et al.  A new approach to measure visual field progression in glaucoma patients using variational bayes linear regression. , 2014, Investigative ophthalmology & visual science.

[40]  R. Tibshirani.  Regression Shrinkage and Selection via the Lasso , 1996 .