July 2015
Volume 56, Issue 8
Free
Glaucoma  |   July 2015
Revalidating the Usefulness of a “Sector-Wise Regression” Approach to Predict Glaucomatous Visual Function Progression
Author Affiliations & Notes
  • Kazunori Hirasawa
    Department of Ophthalmology The University of Tokyo, Graduate School of Medicine, Tokyo, Japan
    Orthopic and Visual Science, Department of Rehabilitation, School of Allied Health Sciences, Kitasato University, Kanagawa, Japan
  • Hiroshi Murata
    Department of Ophthalmology The University of Tokyo, Graduate School of Medicine, Tokyo, Japan
  • Ryo Asaoka
    Department of Ophthalmology The University of Tokyo, Graduate School of Medicine, Tokyo, Japan
  • Correspondence: Ryo Asaoka, Department of Ophthalmology, The University of Tokyo Graduate School of Medicine, 7-3-1 Hongo, Bunkyo-ku, Tokyo, 113-8655 Japan; rasaoka-tky@umin.ac.jp
Investigative Ophthalmology & Visual Science July 2015, Vol.56, 4332-4335. doi:https://doi.org/10.1167/iovs.15-16694
  • Views
  • PDF
  • Share
  • Tools
    • Alerts
      ×
      This feature is available to authenticated users only.
      Sign In or Create an Account ×
    • Get Citation

      Kazunori Hirasawa, Hiroshi Murata, Ryo Asaoka; Revalidating the Usefulness of a “Sector-Wise Regression” Approach to Predict Glaucomatous Visual Function Progression. Invest. Ophthalmol. Vis. Sci. 2015;56(8):4332-4335. https://doi.org/10.1167/iovs.15-16694.

      Download citation file:


      © ARVO (1962-2015); The Authors (2016-present)

      ×
  • Supplements
Abstract

Purpose: We previously developed a method to measure visual function progression using “Progression sectors,” which were derived from a clustering analysis of pointwise rates of glaucomatous visual field (VF) progression; however, only short series of VFs, where pointwise linear regression (PLR) is typically not reliable, were analyzed. The purpose of the current study is to further investigate the usefulness of a sector-wise approach in longer series of VFs where the accuracy of PLR is inherently greater.

Methods: We analyzed 16 VFs from each of 71 eyes in 45 patients. Sectorial averages of total deviation (TD) values were calculated at each “Progression sector,” consisting of 23 sectors as well as at “Nouri-Mahdavi” (NM) sectors, consisting of 10 sectors. Linear regression then was applied to predict sectorial averages in the 11th to 16th VFs and absolute prediction errors were compared to PLR.

Results: Absolute prediction errors associated with linear regression of the “Progression sectors” were significantly smaller than those from PLR (P < 0.05) when predicting values in the 14th to 16th VFs. Conversely, prediction errors with NM sectors were significantly larger than those from PLR and “Progression sectors” (P < 0.01) when predicting values in the 11th to 16th VFs.

Conclusions: It is useful to use “Progression sectors” to predict future VF progression in short (1–10 VFs) and long (>10 VFs) VF series, compared to PLR and NM sectors.

An accurate measurement of glaucomatous visual function progression is extremely important for clinical glaucoma management, not least because unnecessary medical and surgical IOP reduction interventions can be associated with various ocular and general complications.13 At the clinical setting, the rate of visual field (VF) progression generally is calculated by applying linear regression to VF measurements. For example, the rate of mean deviation (MD) change can be calculated using the Guided Progression Analysis (GPA) software on the Humphrey Visual Field Analyzer (HFA; Carl Zeiss Meditec AG, Dublin, CA, USA), while changes in pointwise sensitivity can be measured using the PROGRESSOR software (Medisoft, Inc., London, UK). Earlier detection of focal glaucoma progression can be achieved by using a sectorwise or pointwise linear regression (PLR) analysis415 compared to an event analysis,4,6,10 trend analysis with global indices,5,7,8 or pattern deviation.11 However, poor reproducibility of pointwise VF sensitivity measurements16 can hamper early detection17,18 and reduce the accuracy of predictions of future visual function damage. 
In our recent study, we reported that accurate predictions of visual function damage can be achieved by measuring progression in 23 small VF sectors (“Progression sectors,” see Fig. 1),19 because the effect of pointwise variability is reduced. The “Progression sectors” were developed based on clustering of progression rates in longitudinal series of VFs (n = 10) from 412 eyes in 412 open-angle glaucoma patients. This approach is in contrast to many other previous methods in which VF sectors were decided according to cross-sectional intertest point relationships.15,2023 However, Nouri-Mahdavi et al.24 also developed VF sectors consisting of 10 sectors based on pointwise progression speed. Our previous report suggested that the 10 VF sectors from the report of Nouri-Mahdavi et al.24 (“Nouri-Mahdavi” [NM] sectors) were not useful to accurately predict future progression compared to our proposed “Progression sectors.”19 One of the problems with the clustering technique is how to adequately choose the number of clusters.25 
Figure 1
 
The “Progression sectors” identified by clustering based on the rate of progression. Progression sectors19 consist of 23 sectors, and are illustrated in the 24-2 VF test pattern for a right eye (left). “NM sectors”24 consist of 10 sectors, which are based on hierarchical cluster analysis on pointwise rates of progression based on linear regression analysis of the threshold sensitivity at individual test locations over time They also are illustrated for comparison (right). The blind spot is filled in black.
Figure 1
 
The “Progression sectors” identified by clustering based on the rate of progression. Progression sectors19 consist of 23 sectors, and are illustrated in the 24-2 VF test pattern for a right eye (left). “NM sectors”24 consist of 10 sectors, which are based on hierarchical cluster analysis on pointwise rates of progression based on linear regression analysis of the threshold sensitivity at individual test locations over time They also are illustrated for comparison (right). The blind spot is filled in black.
Previous research suggests that at least five,26 eight, or even more2729 VFs are needed to obtain a reliable forecast of future VF progression. Thus, our original study merely validated the usefulness of “Progression sectors” to predict future damage (a patient's 10th VF) using short series of VFs (3–9 VFs measured approximately every 6 months). The purpose of the current study is to further investigate the usefulness of the “Progression sectors” approach, in comparison with PLR, to predict progression further into the future (11th–16th VFs) using a longer series of VFs (10 VFs) where PLR is inherently more reliable. 
Methods
Study Parameters
The study was approved by the Research Ethics Committee of the Graduate School of Medicine and Faculty of Medicine at the University of Tokyo. Informed consent was obtained from all subjects. This study was performed according to the tenets of the Declaration of Helsinki. 
Visual field data were obtained from 71 eyes in 45 patients with glaucoma who were followed up in the glaucoma clinic at the University of Tokyo Hospital. These subjects are a subset of those in our previous report.19 Visual field data were obtained retrospectively if a subject had more than 10 VFs when our previous study was conducted, otherwise VF data were obtained prospectively in the period of 2013 to 2015. Table 1 shows the demographic and ocular characteristics of the patients. Measurements of VF were performed using the HFA with the 30-2 or 24-2 program and the Swedish Interactive Threshold Algorithm Standard. When VFs were obtained with the 30-2 test pattern, only the 52 test locations overlapping with the 24-2 test pattern were used in the analysis. Patients' first VFs were excluded from the analysis. Other inclusion criteria in this study were best corrected visual acuity better than 6/12, refraction within ±6 diopter ametropia, no previous ocular surgery except for cataract extraction and intraocular lens implantation, and no other posterior segment eye disease. Reliability criteria for VFs were applied: fixation losses less than 20% and false-positive responses less than 15%; the false-negative rate was not used following a previous report.30 A glaucomatous VF was defined if a VF met any of the following criteria31: (1) a cluster of ≥3 points in the pattern deviation plot in a single hemifield (superior/inferior) with P < 0.05, one of which must have been P < 0.01, excluding the outermost test point of Humphrey Field Analyzer 30-2 program; (2) glaucoma hemifield test (GHT) result outside of normal limits; or (3) abnormal PSD with P < 0.05. Patients who underwent intraocular surgical treatments during the observed period were excluded. The VF of a right eye was mirror-imaged to that of a left eye. 
Table 1
 
Demographic and Ocular Characteristics of the Patients
Table 1
 
Demographic and Ocular Characteristics of the Patients
Prediction accuracy was investigated in the 11th to 16th VFs of a patient using their initial 10 VFs. This is in contrast to our previous study in which the 10th VF was predicted using a patient's previous VFs. The calculation of prediction accuracy was performed as follows: first, in each cluster, the average of total deviation (TD) values from all test locations within respective cluster was calculated at each sector (23 “Progression sectors,” see Figs. 1 and 10 “NM sectors”24) and a linear regression was fit for each cluster location in the VF using the average TD values at the respective cluster location in the first 10 VF follow-ups. The visual function of 11th to 16th VFs were predicted by extrapolating the linear regression. Then, the absolute value of the difference between the predicted sectorial TD values and the actual pointwise TD values as average absolute prediction error was calculated. Finally, prediction errors from the two sector-wise methods were compared to the prediction error associated with PLR. 
Statistical Analysis
All statistical analyses were done using the statistical programming language R (ver. 2.15.0; The R Foundation for Statistical Computing, Vienna, Austria). Absolute prediction errors were compared using the Wilcoxon test where P values were adjusted using Benjamini-Hochberg procedure for multiple testing.32 
Results
The demographics of the patients in this study are summarized in Table 1. The mean of all TD values in the initial VF was −7.5 ± 5.4 dB (mean ± SD) and −11.3 ± 6.3 dB in the 16th VF. Patients' initial age was 50.9 ± 12.5 years. The global progression rate of the mean of TD values was −0.41 ± 0.88 dB/y. The interval between the first and 10th VF was 4.9 ± 0.8 years and the interval between the 10th and 16th VF was 3.8 ± 1.2 years. 
Figure 2 and Table 2 summarize the absolute prediction errors associated with PLR and the sector-wise methods. The prediction errors associated with the “Progression sectors” were significantly smaller than those with PLR (P < 0.05, Wilcoxon test adjusted with Benjamini-Hochberg procedure for multiple testing32) when predicting the 14th to 16th VFs, although the prediction error of the method was significantly larger than PLR when predicting the 11th VF (P < 0.05, Wilcoxon test adjusted with Benjamini-Hochberg procedure for multiple testing32). Absolute prediction errors with NM sectors were significantly larger than those from PLR and “Progression sectors” (P < 0.01, Wilcoxon test adjusted with Benjamini-Hochberg procedure for multiple testing32) when predicting all future (11th–16th) VFs. 
Figure 2
 
Mean absolute prediction errors when predicting the 11th to 16th VF using the initial 10 VFs. Prediction errors for regression based on PLR, “NM sectors,”24 and “Progression sectors” proposed in the current study. Statistically significant difference at **P < 0.01 and *P < 0.05, respectively, are indicated.
Figure 2
 
Mean absolute prediction errors when predicting the 11th to 16th VF using the initial 10 VFs. Prediction errors for regression based on PLR, “NM sectors,”24 and “Progression sectors” proposed in the current study. Statistically significant difference at **P < 0.01 and *P < 0.05, respectively, are indicated.
Table 2
 
Mean Absolute Prediction Errors Predicting the 11th to 16th VF Using Initial 10 VFs
Table 2
 
Mean Absolute Prediction Errors Predicting the 11th to 16th VF Using Initial 10 VFs
Discussion
In our previous study we reported that the prediction of future visual function progression could be improved by dividing the VF into small sectors.19 In the current study, the usefulness of these “Progression sectors” was further investigated using this longer series of VF data. As a result, accurate predictions were observed with the PLR and “Progression sectors” approach than the “NM sectors” approach in predicting a patient's 11th to 13th VFs, and more accurate predictions were observed with the “Progression sectors” approach than with PLR and also “NM sectors” in predicting a patient's 14th to 16th VFs. 
It is widely acknowledged that PLR is relatively unreliable until a considerable number of VFs are accumulated.2628 Indeed, our previous studies strongly supported this notion.29 In the current study a larger number of VFs (n = 10) was used to predict future visual function damage, and prediction accuracy was once again observed to fall as the interval to the predicted VF lengthened. This could be caused by the unavoidable short- and long-term fluctuation of VF measurements,33,34 including poor fixation/inattention or fatigue. Indeed, significant eye movements can occur during an otherwise reliable VF test.35,36 In the current study, 10 VFs were acquired over approximately 5 years, yet PLR remained comparatively unreliable when predicting distant future progression. This is clinically very important because treatments are intensified according to estimates of VF progression rates. The current data were obtained from a standard hospital glaucoma clinic where VF tests were conducted in 6-monthly intervals; however, clinicians often are required to make treatment decisions with less frequently measured VFs.37 As shown in the current study, future visual function progression was more accurately predicted in the 14th to 16th VFs by applying the “Progression sectors” approach compared to PLR. Conversely, prediction accuracy of the “Progression sectors” method was worse for the immediate (11th) VF compared to PLR; however, the average increase in prediction error for the sectors method for the range of VFs predicted (between 11th and 16th VF) was merely 1.0 dB, which was half of that of PLR. 
Nouri-Mahdavi et al.24 also have reported a sector-based method for measuring VF progression. Glaucomatous VF damage follows the structural pattern of retinal nerve fibers,3840 and our “Progression sectors” and “NM sectors” mirror the structure-function relationship; however, in the current study, “NM sectors” led to a significantly worse prediction of future visual function damage compared to PLR. Glaucomatous progression usually starts focally and, hence, the accuracy of progression detection methods is a balance between reducing the variability of pointwise VF sensitivity by averaging test points and not masking focal damage by averaging too many VF points in sectors. Thus, the poor performance of “NM sectors” for prediction, observed in the current study, could be attributed to the difference in the size of the clusters, since the VF is divided into much larger sectors using the approach of Nouri-Mahdavi et al.24 
As shown in Figure 1, our “Progression sectors” are smaller in the central area, such as sectors 12 and 15; this result is in agreement with previous reports that suggest these areas tend to progress independently from other test locations.41 This is probably because of the rich retinal nerve fiber layer (RNFL) distribution in this area. Indeed, our group previously suggested that an even finer sector can be identified in this area,42 which is smaller than the 6° interval of the 24-2 HFA VF. A future study should be performed using series of VFs from patients with grids of 10-2 and 24-2, to investigate whether it is useful to measure progression in smaller sectors in the central area. 
In conclusion, predicting future visual function damage is difficult due to the inherent variability of these measurements; however, the accuracy for future visual function damage using our sector-based approach is equal or improving in short and long VF series compared to PLR and NM sectors. 
Acknowledgments
Supported in part by JST-CREST and JSPS KAKENHI Grant 26462679. 
Disclosure: K. Hirasawa, None; H. Murata, None; R. Asaoka, None 
References
Shields M. B. Textbook of Glaucoma. Baltimore, MD: William & Wilkins; 1997.
Aihara M, Shirato S, Sakata R. Incidence of deepening of the upper eyelid sulcus after switching from latanoprost to bimatoprost. Jpn J Ophthalmol. 2011; 55: 600–604.
Yamada Y, Takayanagi R, Tsuchiya K, et al. Assessment of systemic adverse reactions induced by ophthalmic beta-adrenergic receptor antagonists. J Ocul Pharmacol Ther. 2001; 17: 235–248.
McNaught AI, Crabb DP, Fitzke FW, Hitchings RA. Visual field progression: comparison of Humphrey Statpac2 and pointwise linear regression analysis. Graefes Arch Clin Exp Ophthalmol. 1996; 234: 411–418.
Smith SD, Katz J, Quigley HA. Analysis of progressive change in automated visual fields in glaucoma. Invest Ophthalmol Vis Sci. 1996; 37: 1419–1428.
Viswanathan AC, Fitzke FW, Hitchings RA. Early detection of visual field progression in glaucoma: a comparison of PROGRESSOR and STATPAC 2. Br J Ophthalmol. 1997; 81: 1037–1042.
Wild JM, Hutchings N, Hussey MK, Flanagan JG, Trope GE. Pointwise univariate linear regression of perimetric sensitivity against follow-up time in glaucoma. Ophthalmology. 1997; 104: 808–815.
Nouri-Mahdavi K, Brigatti L, Weitzman M, Caprioli J. Comparison of methods to detect visual field progression in glaucoma. Ophthalmology. 1997; 104: 1228–1236.
Mayama C, Araie M, Suzuki Y, et al. Statistical evaluation of the diagnostic accuracy of methods used to determine the progression of visual field defects in glaucoma. Ophthalmology. 2004; 111: 2117–2125.
Nouri-Mahdavi K, Caprioli J, Coleman AL, Hoffman D, Gaasterland D. Pointwise linear regression for evaluation of visual field outcomes and comparison with the advanced glaucoma intervention study methods. Arch Ophthalmol. 2005; 123: 193–199.
Manassakorn A, Nouri-Mahdavi K, Koucheki B, Law SK, Caprioli J. Pointwise linear regression analysis for detection of visual field progression with absolute versus corrected threshold sensitivities. Invest Ophthalmol Vis Sci. 2006; 47: 2896–2903.
Nouri-Mahdavi K, Hoffman D, Ralli M, Caprioli J. Comparison of methods to predict visual field progression in glaucoma. Arch Ophthalmol. 2007; 125: 1176–1181.
Nassiri N, Moghimi S, Coleman AL, Law SK, Caprioli J, Nouri-Mahdavi K. Global and pointwise rates of decay in glaucoma eyes deteriorating according to pointwise event analysis. Invest Ophthalmol Vis Sci. 2013; 54: 1208–1213.
Katz J, Sommer A, Gaasterland DE, Anderson DR. Comparison of analytic algorithms for detecting glaucomatous visual field loss. Arch Ophthalmol. 1991; 109: 1684–1689.
Mandava S, Zulauf M, Zeyen T, Caprioli J. An evaluation of clusters in the glaucomatous visual field. Am J Ophthalmol. 1993; 116: 684–691.
Artes PH, Iwase A, Ohno Y, Kitazawa Y, Chauhan BC. Properties of perimetric threshold estimates from Full Threshold SITA Standard, and SITA Fast strategies. Invest Ophthalmol Vis Sci. 2002; 43: 2654–2659.
De Moraes CG, Liebmann CA, Susanna R,Jr, Ritch R, Liebmann JM. Examination of the performance of different pointwise linear regression progression criteria to detect glaucomatous visual field change. Clin Experiment Ophthalmol. 2012; 40: e190–e196.
Jansonius NM. On the accuracy of measuring rates of visual field change in glaucoma. Br J Ophthalmol. 2010; 94: 1404–1405.
Hirasawa K, Murata H, Hirasawa H, Mayama C, Asaoka R. Clustering visual field test points based on rates of progression to improve the prediction of future damage. Invest Ophthalmol Vis Sci. 2014; 55: 7681–7685.
Suzuki Y, Kitazawa Y, Araie M, et al. Mathematical and optimal clustering of test points of the central 30-degree visual field of glaucoma. J Glaucoma. 2001; 10: 121–128.
Asman P, Heijl A. Arcuate cluster analysis in glaucoma perimetry. J Glaucoma. 1993; 2: 13–20.
Chauhan BC, Drance SM, Lai C. A cluster analysis for threshold perimetry. Graefes Arch Clin Exp Ophthalmol. 1989; 227: 216–220.
Suzuki Y, Araie M, Ohashi Y. Sectorization of the central 30 degrees visual field in glaucoma. Ophthalmology. 1993; 100: 69–75.
Nouri-Mahdavi K, Mock D, Hosseini H, et al. Pointwise rates of visual field progression cluster according to retinal nerve fiber layer bundles. Invest Ophthalmol Vis Sci. 2012; 53: 2390–2394.
Rokach L. Clustering methods. In: Maimon O, Rokach L, eds. Data Mining and Knowledge Discovery Handbook. Berlin: Springer; 2010: 321–352.
Krakau CE. A statistical trap in the evaluation of visual field decay. Acta Ophthalmol Suppl. 1985; 173: 19–21.
Holmin C, Krakau CE. Regression analysis of the central visual field in chronic glaucoma cases. A follow-up study using automatic perimetry. Acta Ophthalmol (Copenh). 1982; 60: 267–274.
Spry PG, Bates AB, Johnson CA, Chauhan BC. Simulation of longitudinal threshold visual field data. Invest Ophthalmol Vis Sci. 2000; 41: 2192–2200.
Taketani Y, Murata H, Fujino Y, Mayama C, Asaoka R. How many visual fields are required to precisely predict future test results in glaucoma patients when using different trend analyses? Invest Ophthalmol Vis Sci. In press.
Bengtsson B, Heijl A. False-negative responses in glaucoma perimetry: indicators of patient performance or test reliability? Am J Ophthalmol. 2000; 130: 689.
Anderson DR, Patella VM. Automated Static Perimetry, 2nd ed. St. Louis, MO: Mosby; 1999.
Benjamini YYH. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Stat Soc Series B. 1995; 57: 289–300.
Flammer J, Drance SM, Fankhauser F, Augustiny L. Differential light threshold in automated static perimetry. Factors influencing short-term fluctuation. Arch Ophthalmol. 1984; 102: 876–879.
Flammer J, Drance SM, Zulauf M. Differential light threshold. Short- and long-term fluctuation in patients with glaucoma normal controls, and patients with suspected glaucoma. Arch Ophthalmol. 1984; 102: 704–706.
Demirel S, Vingrys AJ. Eye movements during perimetry and the effect that fixational instability has on perimetric outcomes. J Glaucoma. 1994; 3: 28–35.
Demirel S, Vingrys AJ. Fixational instability during perimetry and the blindspot monitor. In: Mills RP, ed. Perimetry Update, 1992/1993. Amsterdam: Kugler Publications; 1992: 515–520.
Malik R, Baker H, Russell RA, Crabb DP. A survey of attitudes of glaucoma subspecialists in England and Wales to visual field test intervals in relation to NICE guidelines. BMJ Open. 2013; 3.
Garway-Heath DF, Poinoosawmy D, Fitzke FW, Hitchings RA. Mapping the visual field to the optic disc in normal tension glaucoma eyes. Ophthalmology. 2000; 107: 1809–1815.
Weber J, Dannheim F, Dannheim D. The topographical relationship between optic disc and visual field in glaucoma. Acta Ophthalmol (Copenh). 1990; 68: 568–574.
Asman P, Heijl A. Glaucoma Hemifield Test. Automated visual field evaluation. Arch Ophthalmol. 1992; 110: 812–819.
Zalta AH. Use of a central 10 degrees field and size V stimulus to evaluate and monitor small central islands of vision in end stage glaucoma. Br J Ophthalmol. 1991; 75: 151–154.
Asaoka R. Mapping glaucoma patients' 30-2 and 10-2 visual fields reveals clusters of test points damaged in the 10-2 grid that are not sampled in the sparse 30-2 grid. PLoS One. 2014; 9: e98525.
Figure 1
 
The “Progression sectors” identified by clustering based on the rate of progression. Progression sectors19 consist of 23 sectors, and are illustrated in the 24-2 VF test pattern for a right eye (left). “NM sectors”24 consist of 10 sectors, which are based on hierarchical cluster analysis on pointwise rates of progression based on linear regression analysis of the threshold sensitivity at individual test locations over time They also are illustrated for comparison (right). The blind spot is filled in black.
Figure 1
 
The “Progression sectors” identified by clustering based on the rate of progression. Progression sectors19 consist of 23 sectors, and are illustrated in the 24-2 VF test pattern for a right eye (left). “NM sectors”24 consist of 10 sectors, which are based on hierarchical cluster analysis on pointwise rates of progression based on linear regression analysis of the threshold sensitivity at individual test locations over time They also are illustrated for comparison (right). The blind spot is filled in black.
Figure 2
 
Mean absolute prediction errors when predicting the 11th to 16th VF using the initial 10 VFs. Prediction errors for regression based on PLR, “NM sectors,”24 and “Progression sectors” proposed in the current study. Statistically significant difference at **P < 0.01 and *P < 0.05, respectively, are indicated.
Figure 2
 
Mean absolute prediction errors when predicting the 11th to 16th VF using the initial 10 VFs. Prediction errors for regression based on PLR, “NM sectors,”24 and “Progression sectors” proposed in the current study. Statistically significant difference at **P < 0.01 and *P < 0.05, respectively, are indicated.
Table 1
 
Demographic and Ocular Characteristics of the Patients
Table 1
 
Demographic and Ocular Characteristics of the Patients
Table 2
 
Mean Absolute Prediction Errors Predicting the 11th to 16th VF Using Initial 10 VFs
Table 2
 
Mean Absolute Prediction Errors Predicting the 11th to 16th VF Using Initial 10 VFs
×
×

This PDF is available to Subscribers Only

Sign in or purchase a subscription to access this content. ×

You must be signed into an individual account to use this feature.

×