June 2007
Volume 48, Issue 6
Free
Glaucoma  |   June 2007
The Effect of Disc Size and Severity of Disease on the Diagnostic Accuracy of the Heidelberg Retina Tomograph Glaucoma Probability Score
Author Affiliations
  • Linda M. Zangwill
    From the Hamilton Glaucoma Center, Department of Ophthalmology, the
  • Sonia Jain
    Division of Biostatistics and Bioinformatics, Department of Family and Preventive Medicine, University of California, San Diego, La Jolla, California.
  • Lyne Racette
    From the Hamilton Glaucoma Center, Department of Ophthalmology, the
  • Karin B. Ernstrom
    Division of Biostatistics and Bioinformatics, Department of Family and Preventive Medicine, University of California, San Diego, La Jolla, California.
  • Christopher Bowd
    From the Hamilton Glaucoma Center, Department of Ophthalmology, the
  • Felipe A. Medeiros
    From the Hamilton Glaucoma Center, Department of Ophthalmology, the
  • Pamela A. Sample
    From the Hamilton Glaucoma Center, Department of Ophthalmology, the
  • Robert N. Weinreb
    From the Hamilton Glaucoma Center, Department of Ophthalmology, the
Investigative Ophthalmology & Visual Science June 2007, Vol.48, 2653-2660. doi:https://doi.org/10.1167/iovs.06-1314
  • Views
  • PDF
  • Share
  • Tools
    • Alerts
      ×
      This feature is available to authenticated users only.
      Sign In or Create an Account ×
    • Get Citation

      Linda M. Zangwill, Sonia Jain, Lyne Racette, Karin B. Ernstrom, Christopher Bowd, Felipe A. Medeiros, Pamela A. Sample, Robert N. Weinreb; The Effect of Disc Size and Severity of Disease on the Diagnostic Accuracy of the Heidelberg Retina Tomograph Glaucoma Probability Score. Invest. Ophthalmol. Vis. Sci. 2007;48(6):2653-2660. https://doi.org/10.1167/iovs.06-1314.

      Download citation file:


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

      ×
  • Supplements
Abstract

purpose. To compare the effect of disc size and disease severity on the Heidelberg Retina Tomograph (HRT) Glaucoma Probability Score (GPS) and the Moorfields Regression Analysis (MRA) for discriminating between glaucomatous and healthy eyes.

methods. Ninety-nine eyes with repeatable standard automated perimetry results showing glaucomatous damage and 62 normal eyes were included from the longitudinal Diagnostic Innovations in Glaucoma Study (DIGS). The severity of glaucomatous visual field defects ranged from early to severe (average [95% CI] pattern standard deviation [PSD] was 5.7 [5.0–6.5] dB). The GPS (HRTII ver. 3.0; Heidelberg Engineering, Heidelberg, Germany) utilizes two measures of peripapillary retinal nerve fiber layer shape (horizontal and vertical retinal nerve fiber layer curvature) and three measures of optic nerve head shape (cup depth, rim steepness, and cup size) as input into a relevance vector machine learning classifier that estimates a probability of having glaucoma. The MRA compares measured rim area with predicted rim area adjusted for disc size to categorize eyes as outside normal limits, borderline, or within normal limits. The effect of disc size and severity of disease on the diagnostic accuracy of both GPS and MRA was evaluated using the generalized estimating equation marginal logistic regression analysis.

results. Using the manufacturers’ suggested cutoffs for GPS global classification (>64% as outside normal limits), the sensitivity and specificity (95% CI) were 71.7% (62.2%–79.7%) and 82.3% (71.0%–89.8%), respectively. The sensitivity and specificity (95% CI) of the MRA result were 66.7% (58.0%–76.1%) and 88.7% (78.5%–94.34%), respectively. Likelihood ratios for regional GPS and MRA results outside normal limits ranged from 4.0 to 10.0, and 6.0 to infinity, respectively. Disc size and severity of disease were significantly associated with the sensitivity of both GPS and MRA.

conclusions. GPS tended to have higher sensitivities and somewhat lower specificities and lower likelihood ratios than MRA. These results suggest that in this population, GPS and MRA differentiate between glaucomatous and healthy eyes with good sensitivity and specificity. In addition, the likelihood ratios suggest that GPS may be most useful for confirming a normal disc, whereas MRA may be most helpful in confirming a suspicion of glaucoma. Larger disc size and more severe field loss were associated with improved diagnostic accuracy for both GPS and MRA.

Glaucoma is a progressive optic neuropathy diagnosed by identifying characteristic optic nerve, retinal nerve fiber layer (RNFL), and visual field damage. For more than 10 years, confocal scanning laser ophthalmoscopy (CSLO), scanning laser polarimetry (SLP), and optical coherence tomography (OCT) have provided objective, reproducible measurements that, when used with other clinical information, can assist the clinician in differentiating between normal and glaucomatous eyes. Imaging instruments recently have incorporated normative databases into their statistical analyses so that an automatic assessment of whether the eye is outside normal limits (ONL) is provided to the clinician. The diagnostic accuracy of optic disc and RNFL measurements and the automatic classifications of these instruments have been evaluated extensively. 1 2 3 4 Although the overall diagnostic accuracy of the best parameters of each of these imaging instruments may be similar, the performance of diagnostic tests can vary among subgroups of patients with glaucoma, according to clinical and nonclinical characteristics. 2 3 4 5 6  
It is well established that the accuracy of tests increases with increasing severity of disease. It is therefore important to provide an estimate of diagnostic precision at various stages of disease. In addition, optic disc size has been shown to influence the diagnostic accuracy of imaging instruments, particularly the confocal scanning laser ophthalmoscope. 6 7 8 9 10 11 12 Recently Medeiros et al. 8 have used marginal logistic regression methods 13 14 15 for simultaneous evaluation of the effect of severity of disease and disc size on the diagnostic accuracy of imaging instruments. 
Each imaging instrument has specific advantages and limitations. 1 A limitation of a commercially available CSLO, the Heidelberg Retina Tomograph (HRT, Heidelberg Engineering, Heidelberg, Germany) has been its reliance on an operator to outline the disc margin before topographic optic disc parameters can be calculated. The outlining of the disc margin adds processing time, and differences in how it is completed can lead to interobserver variability in stereometric variables. 16 17 In addition, many of these topographic optic disc parameters are calculated with a reference plane. 
The recently released HRT software version 3.0 includes the Glaucoma Probability Score (GPS). As the GPS calculation is based on the overall shape of the optic nerve head and posterior pole and does not rely on the outlining of the disc margin for its calculation, 18 it may be less influenced by optic disc size than are conventional CSLO topographic optic disc parameters and the Moorfields Regression Analysis (MRA). The objective of this study was to compare the GPS with the MRA for discriminating between glaucomatous and healthy eyes and to evaluate the influence of disease severity and optic disc size on the diagnostic accuracy of these two classification systems. 
Methods
Subjects
One randomly selected eye from each of 99 patients with glaucoma and 62 normal subjects participating in the longitudinal Diagnostic Innovations in Glaucoma Study (DIGS) was included in the study. 
All participants underwent a complete ophthalmic examination including slit lamp biomicroscopy, intraocular pressure measurement, dilated stereoscopic fundus examination, and standard automated perimetry (SAP) using the Swedish Interactive Threshold Algorithm (SITA) and the 24-2 program (Humphrey Field Analyzer; Carl Zeiss Meditec, Inc., Dublin, CA). Visual fields were reliable (fixation losses and false positive and false negative responses <33%). To be included in DIGS, at study entry, all participants had open angles, a best corrected acuity of 20/40 or better, a spherical refraction within ±5.0 D, and cylinder correction within ±3.0 D. A family history of glaucoma was allowed. 
Participants were excluded from DIGS if they had a history of intraocular surgery except for uncomplicated cataract or glaucoma surgery. Participants were also excluded if there was evidence of secondary causes of elevated IOP (e.g., iridocyclitis, trauma), other intraocular eye disease, other diseases affecting the visual field (e.g., pituitary lesions, demyelinating diseases, HIV+ or AIDS, or diabetic retinopathy), or medications known to affect visual field sensitivity. 
For purposes of this analysis, patients were classified as having glaucoma if they had at least two consecutive standard automated perimetry examinations with either a pattern standard deviation (PSD) outside the 95% normal limits or a glaucoma hemifield test (GHT) result outside the 99% normal limits. At least one of the abnormal fields was obtained within 6 months of CSLO imaging. The appearance of the optic disc was not used as criteria for designation as glaucomatous. 
Normal control eyes had IOP <22 mm Hg with no history of elevated IOP and with normal visual field results, defined as a PSD within 95% confidence limits and a GHT result within normal limits (WNL). Normal control eyes also had no evidence of glaucomatous optic disc damage (no diffuse or focal rim thinning, or RNFL defects) as evaluated by clinical examination. 
The severity of visual field damage was assessed on a scale of 0 (no field loss) to 20 (end-stage glaucoma), according to the Advanced Glaucoma Intervention Study (AGIS) severity score. The AGIS score is based on the extent of damage measured by the total deviation plot at different visual field locations and has been described in detail elsewhere. 19  
The research adhered to the tenets of the Declaration of Helsinki. Informed consent was obtained from all participants and the University of California, San Diego, Human Subjects Committee approved all methodology. 
Confocal Scanning Laser Ophthalmoscopy
The HRTII provides topographical measures of the optic disc and peripapillary retina and has been discussed in detail elsewhere. 1 Three scans centered on the optic disc were automatically obtained for each test eye, and a mean topography was created. Magnification errors were corrected by using patients’ corneal curvature measurements. The optic disc margin was outlined on the mean topography image by trained technicians while they viewed simultaneous stereoscopic photographs of the optic disc. All images included in the analysis were reviewed for adequate centration, focus, and illumination; all mean topography images had a standard deviation of <50 μm. The scans were obtained with HRT software version 1.5.9.0 or earlier, but were analyzed with the recently released software version 3.0. 
HRT software version 3.0 includes improved alignment algorithms, a larger normative database, and the calculation of the GPS. 18 Two measures of peripapillary retinal nerve fiber layer shape (horizontal and vertical retinal nerve fiber layer curvature) and three measures of optic nerve head shape (cup depth, rim steepness, and cup size) are used as input into a relevance vector machine learning classifier to estimate the probability of having glaucoma as between 0% and 100%. Two mathematical functions are used to model the topography of the optic nerve head: (1) A Gaussian cumulative distribution function is used to model the optic disc, and (2) a quadratic (parabolic) surface is used to model the peripapillary retina. To parameterize the cup, a parabolic surface is fitted to the peripapillary region of each topograph. As outlined in Swindale et al., 18 the parabolic surface, which serves as a reference plane for estimating the cup parameters, is then subtracted from the topography. The average location of the deepest points in the difference topograph is used to identify the cup center. In contrast to Swindale et al., 18 the GPS constructs a cumulative Gaussian distribution of topograph heights to estimate the cup radius (r) such that p(radius ≤ r) = 0.5. The cup radius r serves as a cup margin. Thus, the cup area is computed as the area of circle of radius r and the mean cup depth is computed as the average height of the measurements inside the cup in the difference topograph. The rim steepness estimates are derived from the radial topograph–height gradients. No contour line or reference plane is used in the GPS calculation. GPS output is then automatically classified into three categories: outside normal limits (ONL; GPS > 64%), borderline (BL; GPS between 24% and 64%) and within normal limits (WNL; GPS < 24%). 
The MRA compares measured rim area to predicted rim area adjusted for disc size, to categorize eyes as ONL, BL, or WNL. 9 It relies on a contour line and the standard reference plane (50 μm below the mean height of the contour in the temporal sector between 350° and 356°) for its measurements. By using the HRT 3.0 software, both the GPS and MRA classify eyes as within normal limits (WNL), borderline (BL), or outside normal limits (ONL), according to the same normative database of 700 eyes of whites and 200 eyes of African-Americans. For this analysis, the white normative database was used because most of the DIGS participants are of European descent. The comparison to the normative database is provided in six regions (superior temporal, inferior temporal, temporal, superior nasal, inferior nasal, and nasal), and as an overall global classification (if any of the six regions are ONL, then the eye is classified as ONL). For analysis using the MRA and GPS as categorical variables (ONL versus WNL), BL values were considered WNL for estimates of the sensitivity, specificity, and likelihood ratio. 
In addition to estimating diagnostic accuracy by using the MRA and GPS as categorical variables, we evaluated the sensitivity of each at fixed specificities of 80% and 90%. We converted the MRA into a continuous variable by subtracting the predicted MRA from the actual MRA. This difference is used for determining whether the MRA is ONL. The difference between the predicted and actual MRA for each region was used to estimate sensitivity. For GPS the continuous variable used was the RVM output between 0 and 100. 
The area under the receiver operator characteristic curve (AUROC) was calculated for both the three-level categorical variables (WNL, BL, and ONL) and for the MRA and GPS continuous variables MRA (predicted minus actual) and GPS (relevance vector machine output between 0%–100%). 
Statistical Analysis
The sensitivity, specificity, and likelihood ratios of the MRA and GPS were compared for both global and regional results. The 95% CI for sensitivity and specificity were calculated by the Wilson Score Method without continuity correction for proportions. 20 21 Likelihood ratio confidence intervals (CI) were computed with the Simel et al. method, 22 and AUROC CIs were calculated using the method of Delong et al. 23  
We compared the influence of disc size and severity of glaucoma on the diagnostic accuracy of the MRA and GPS as categorical variables (ONL versus not ONL) by using a generalized estimating equation (GEE) marginal logistic regression modeling approach. This method directly compares the effect of covariates on several tests performed on the same group of subjects and adjusts for subject-specific correlation. 13 14 15 Sensitivity is the dependent variable in the GEE logistic regression model. The covariates, disc area, AGIS score, and test type (GPS or MRA) were entered into the model as independent variables, and first-order interaction terms were included. An exchangeable correlation structure among observations was assumed. This method has recently been applied to compare the influence of disc size and severity of glaucoma on the sensitivity of the GDxVCC (Carl Zeiss Meditec, Inc.), the HRT (Heidelberg Engineering), and the Stratus OCT (Carl Zeiss Meditec, Inc.) for glaucoma detection. 8  
In addition, as a confirmatory analysis of the categorical diagnostic findings, we evaluated the sensitivity for detection of glaucoma at two levels of specificity: 80% and 90%. For this analysis, the 80% and 90% cutoffs for specificity were determined by using the healthy control eyes and were applied to the glaucomatous eyes to determine the corresponding sensitivity cutoffs. This allowed sensitivity to be studied at the same level of specificity. GEE marginal logistic regression was then used to evaluate the influence of disc area and severity of glaucoma on diagnostic performance. The dependent variable, sensitivity, was dichotomized to 1 or 0 if the result of the diagnostic test was above or below the sensitivity cutoff. 
Statistical analyses were performed with commercial software (JMP ver. 6.0; SAS Institute, Cary, NC) and R (http://www.r-project.org/), version 2.1.1. 
Results
The demographic and ocular characteristics of the 99 subjects with glaucoma and the 62 healthy control subjects are presented in Table 1 . The average age (95% CI) of patients with glaucoma was significantly higher than the age of those with the healthy control eyes: 67.9 (65.8–70.0) and 57.7 (55.2–60.3) years, respectively. Mean deviation (95% CI) of the standard automated perimetry (SAP) test closest to the CSLO imaging date was −5.2 (−6.2 to −4.2) dB in glaucomatous eyes and −0.05 (−0.7 to −0.2) dB in control eyes. The severity of glaucomatous visual field defects ranged from early to severe (average [95% CI]) PSD was 5.7 (5.0–6.5) dB; and AGIS score was 3.9 (3.1–4.7). 
The AUROC curves are presented for the GPS and MRA as three-level categorical variables in Table 2and as continuous variables in Table 3 . In general, the AUROC curves estimated from the continuous variables tended to be higher for GPS than for MRA. However, these differences reached statistical significance only in the temporal superior region. The AUROC curve differences between the GPS and MRA were less consistent for estimates based on categorical variables. 
The analysis of the GPS and MRA as categorical variables (outside normal limits versus not outside normal limits) found that both GPS and MRA have relatively high specificity, at least 82% for all regions. In general, the MRA had a higher specificity (range, from 88.7% for global to 100% in the temporal superior region) than GPS (range from 82.3% for global to 93.6% in the temporal superior region). However, the GPS tended to have higher sensitivity (range, from 64.7% temporal superior to 72.7% in the temporal inferior region), than MRA (25.3% in the temporal superior region to 67.7% globally). These differences in diagnostic accuracy between MRA and GPS were found to be statistically significant in the logistic marginal regression models for each of the four regional analyses (variable test type, P < 0.001; Table 4 , Fig. 1 ). 
The MRA tended to have larger likelihood ratios of an ONL result (range, 5.99 to infinity) than the GPS (range, 4.04–10.02). Likelihood ratios can be interpreted by their effect on the posttest probability of disease, 24 with likelihood ratios between 0.5 and 2 considered as having an insignificant effect, ratios between 0.2 and 0.5 or between 2 and 5 as having a small effect, ratios between 0.1 and 0.2 or between 5 and 10 as having a moderate effect, and ratios <0.1 or >10 as having a large effect. The MRA likelihood ratios of an ONL result were all ≥6, indicating a moderate to large effect on the posttest probability of glaucoma. The GPS ONL value had a moderate effect on posttest probabilities of glaucoma, with values ranging from 4.04 to 10.02. The likelihood ratios for a WNL result were smaller for GPS than for MRA, with most GPS likelihood ratios having a moderate to large effect on posttest probabilities (range, 0.014–0.18), whereas MRA likelihood ratios had an insignificant effect (range, 0.20–0.59). BL values had an insignificant or weak effect on posttest probabilities, with likelihood ratios for MRA ranging from 1.06 to 3.54 and for GPS ranging from 0.66 to 1.19. 
Moreover, as shown in Table 4and Figure 1 , the GEE logistic marginal regression models indicate that for each region, the independent variables AGIS score and disc size influenced the diagnostic accuracy (normal versus ONL) of both GPS and MRA (P < 0.05). For the global parameter, however, only severity of glaucoma (AGIS) was positively associated with increased sensitivity (P = 0.007); disc size did not reach statistical significance (P = 0.081). 
The GEE marginal logistic model used was  
\[\mathrm{logit}(\mathrm{sensitivity})\ {=}\ {\beta}_{0}\ {+}\ {\beta}_{1}(\mathrm{AGIS})\ {+}\ {\beta}_{2}(\mathrm{disc\ area})\ {+}\ {\beta}_{3}(\mathrm{test-type}){+}\ {\beta}_{4}(\mathrm{AGIS}\ {\times}\ \mathrm{disc\ area})\ {+}\ {\beta}_{5}(\mathrm{AGIS}\ {\times}\ \mathrm{test-type})\ {+}\ {\beta}_{6}(\mathrm{disc\ area}\ {\times}\ \mathrm{test-type}),\]
where AGIS represents the AGIS scores for severity of visual field defect, disc area represents the optic disc area, test type is an indicator variable coding for type of test, MRA versus GPS (GPS is the reference test), and sensitivity is a dichotomous dependent variable. For the continuous case, sensitivity was dichotomized based on cutoffs determined at fixed specificities. Three first-order (pair-wise) interaction terms were also included. 
The parameter estimate of the slope suggests that as disc size increases, the sensitivity increases. Similarly, sensitivity increases with the severity of visual field damage. The interaction terms (disc size versus test type and AGIS score versus test type) were not statistically significant, suggesting that disc size and severity of glaucomatous visual field damage does not affect MRA differently than it affects GPS. 
In addition to analyzing the GPS and MRA as categorical variables, we modeled the sensitivity of the GPS and MRA for detecting glaucoma at two fixed specificities, 80% and 90%, to ensure that the analysis on GPS and MRA as categorical variables was not driven by the inherent tradeoff between sensitivity and specificity. At a fixed specificity of 90%, both severity of visual field damage, and disc size were significantly associated with GPS and MRA sensitivity for detecting glaucoma (Table 5 , Fig. 2 ). At both levels of fixed specificity, no significant differences were found between the sensitivity of the GPS and MRA. However, at a specificity of 80%, the influence of disc size on sensitivity reached statistical significance in the temporal superior region (P = 0.029), but not in the global (P = 0.093), temporal inferior (P = 0.098), nasal superior (P = 0.141), and nasal inferior (P = 0.158) regions. 
Discussion
We found that, when compared to the internal normative database the GPS results ONL tended to have higher sensitivity, but lower specificity than MRA results. The differences in sensitivity and specificity were larger in the regional output compared with the global MRA and GPS output. In addition, this study demonstrates that disc size and severity of visual field damage influence the diagnostic accuracy of both the GPS and MRA in a similar manner. Sensitivity improves with increasing disc size and severity of visual field damage. This finding suggests that even though outlining of the disc margin is not necessary for the GPS, the size of the disc still influences the diagnostic accuracy of the measurement. 
The finding that GPS sensitivity tends to be higher than that of the MRA are consistent with the sensitivity for detection of early glaucoma reported by Harizman et al., 25 72.3% and 59.6%, respectively. This study also confirms previous reports that the diagnostic accuracy of HRT parameters, particularly linear discriminant functions and the MRA improve with increasing disc size, 6 7 8 9 10 11 12 26 probably due to the difficulty in detecting neuroretinal rim loss in a small disc compared with a large disc. 8 Most of these studies evaluated the association of diagnostic accuracy with disc sizes using univariate stratified analysis or regression techniques that do not take into consideration the severity of disease. It is likely that disc size and other covariates have more of an effect on diagnostic accuracy in eyes with early compared with severe glaucomatous visual field damage. The advantage of the multivariate technique used in the present study is that it can evaluate and control for the severity of visual field damage and disc size in the same analysis. 
To compare more directly the results of the two classification parameters, we also chose to examine the sensitivity at two levels of specificity. At fixed specificities, the sensitivity of GPS and MRA were not significantly different. At 90% specificity, both disc size and severity of visual field damage were associated with the sensitivity of the measurements. However, at 80% specificity, the sensitivity was significantly associated with severity of visual field damage, but was not consistently associated with disc size (P ranged from 0.029 for the temporal superior region to 0.158 for the nasal superior region). It is unclear why disc size was not consistently associated with the sensitivity at 80% specificity. 
To complete the analysis at fixed specificity, we incorporated the MRA in the analysis as a continuous variable (actual minus predicted) and compared it directly to the GPS. This comparison of continuous variables utilizes the normative database to calculate the predicted values, and can be considered comparable to the categorical values provided to the clinician on an HRT printout. 
There are several ways to describe and summarize the ability of a diagnostic test to detect disease. The most common measures of diagnostic accuracy include sensitivity, specificity, and AUROC curves. The advantages and limitations of these measures have been described recently. 8 27 In brief, sensitivity, specificity, and AUROC curve provide important information about the overall diagnostic accuracy of a test and were therefore used to compare the diagnostic accuracy of the MRA and GPS parameters, and to evaluate the effects of glaucoma severity and disc size. However, sensitivity, specificity, and AUROC do not necessarily provide information in a form that is useful to the clinician or patient in clinical decision-making. For example, sensitivity and specificity are related (as sensitivity increases, specificity decreases and vice versa) and depend on the specific cutoffs used to define the disease. For this reason, we included a comparison of the sensitivity of the GPS and MRA at fixed specificities. By definition, at a high fixed specificity of 90%, 10% of normal eyes will be classified as glaucomatous. Unfortunately, it is difficult to apply this information to a specific patient. Similarly, the AUROC curve is important for comparing the diagnostic accuracy of two different diagnostic tests, but has little clinical meaning for making decisions regarding a particular patient. In contrast, the likelihood ratio provides this type of information; it expresses the magnitude by which the probability of a diagnosis in a given patient is modified by the results of the test. In another words, the likelihood ratio indicates how much a given diagnostic test result will raise or lower the pretest probability of the disease in question. We therefore reported the likelihood ratios for the three categorical outputs of the MRA and GPS: ONL, BL, and WNL. We found that an MRA output of ONL had a moderate to large effect on the posttest probability of glaucoma both globally and for each region. A GPS output of ONL had a moderate effect on posttest probability. A GPS output of WNL was much more strongly associated with the probability that the test result was normal than was the same output from the MRA. The results in this study population suggest that GPS provides better information for confirming a normal disc, whereas MRA is most helpful in confirming a suspicion of glaucoma. It should be noted that even small changes in posttest probability may be relevant, depending on other relevant clinical information and the pretest probability of disease. 
There are several possible limitations to the present study. First, the subjects with glaucoma were older than the normal subjects, which may lead to an overestimation of the diagnostic accuracy of the methods. However, the main objective of this study was to use the same population to compare the influence of disc size and severity of disease between the MRA and GPS, and since there was no relationship between age and disc size (R 2 = 0.006, P = 0.153) and a very weak relationship between age and AGIS score (R 2 = 0.017, P = 0.018), it is unlikely that the age difference influenced the comparison. Second, the sample size of this study was modest, leading to relatively wide confidence limits for estimates of AUROC, sensitivity, specificity, and likelihood ratios. Larger studies are needed to provide more precise estimates of the diagnostic accuracy. Finally, in our limited population, disc area was larger in the glaucomatous eyes than in the normal eyes, which may affect how disc size influenced the diagnostic accuracy of the MRA and GPS. However, the analysis of sensitivity at a fixed specificity using glaucomatous eyes only confirmed the relationship of better sensitivity with increasing disc size for both MRA and GPS. Other investigators also have found larger discs in subjects with glaucoma than in normal control subjects. 26 A likely explanation for this difference in disc size is sampling bias related to glaucoma’s being more difficult to detect in small optic discs. 8 If this is the case, then small discs are likely to be underrepresented and large discs overrepresented among patients with diagnosed glaucoma in glaucoma specialty clinics. 26  
In conclusion, GPS can differentiate between glaucomatous and healthy eyes with relatively good sensitivity and specificity. Using the manufacturer’s suggested cutoff values for classification as ONL, the GPS tended to have higher sensitivities and lower specificities and likelihood ratios than did the MRA. Disc size influences the diagnostic accuracy of both the GPS and MRA. 
 
Table 1.
 
Demographic and Ocular Characteristics of the Glaucoma and Healthy Control Study Groups
Table 1.
 
Demographic and Ocular Characteristics of the Glaucoma and Healthy Control Study Groups
Glaucoma Eyes (n = 99) Normal Eyes (n = 62) P
Age (y)* 67.9 (65.8, 70.0) 57.7 (55.2, 60.3) <0.0001
Gender (% male) 52% 26% 0.001
Visual field pattern standard deviation (dB)* 5.7 (5.0, 6.5) 1.5 (1.4, 1.6) <0.0001
Visual field mean deviation (dB)* −5.2 (−6.2, −4.2) −0.5 (−0.7,−0.2) <0.0001
AGIS Score* 3.9 (3.1, 4.7) .02 (−0.02, 0.04) <0.0001
Optic disc size* (mm2) 1.97 (1.86, 2.07) 1.76 (1.67, 1.85) 0.006
Table 2.
 
AUROC, Sensitivity, Specificity, and Likelihood Ratios of the MRA and GPS as Categorical Variables
Table 2.
 
AUROC, Sensitivity, Specificity, and Likelihood Ratios of the MRA and GPS as Categorical Variables
AUROC (95% CI) Sensitivity (95% CI) Specificity (95% CI) Likelihood Ratio
Outside Normal Limits (95% CI) Borderline (95% CI) Within Normal Limits (95% CI)
Global
 MRA 0.74 (0.66–0.83) 67.7 (58.0–76.1) 88.7 (78.5–94.4) 5.99 (2.94–12.20) 1.06 (0.52–2.17) 0.20 (0.12–0.34)
 GPS 0.70 (0.61–0.80) 71.7 (62.2–79.7) 82.3 (71.0–89.8) 4.04 (2.33–7.00) 0.91 (0.48–1.71) 0.14 (0.07–0.28)
Temporal superior
 MRA 0.67 (0.60–0.73) 25.3 (17.7–34.6) 100 (94.2–100.0) Inf NA 1.96 (0.89–4.33) 0.59 (0.48–0.72)
 GPS 0.68 (0.58–0.77) 64.7 (54.8–73.4) 93.6 (84.6–97.5) 10.02 (3.84–26.14) 0.80 (0.47–1.35) 0.18 (0.1–0.32)
Temporal inferior
 MRA 0.78 (0.71–0.85) 48.5 (38.9–58.2) 98.4 (91.4–99.7) 30.11 (4.19–218.80) 2.63 (1.04–1.61) 0.33 (0.24–0.45)
 GPS 0.69 (0.59–0.78) 72.7 (63.2–80.5) 88.7 (78.5–94.4) 6.44 (3.17–13.08) 0.66 (0.35–1.25) 0.18 (0.09–0.34)
Nasal superior
 MRA 0.72 (0.66–0.78) 35.4 (26.6–45.2) 98.4 (91.4–99.7) 21.96 (3.08–155.96) 3.54 (1.08–11.61) 0.50 (0.4–0.63)
 GPS 0.74 (0.65–0.83) 67.7 (58.0–76.1) 85.5 (74.7–92.2) 4.66 (2.51–8.66) 1.19 (0.61–2.3) 0.16 (0.09–0.29)
Nasal inferior
 MRA 0.74 (0.68–0.81) 43.4 (34.1–53.3) 96.8 (89.0–99.1) 13.45 (3.38–53.62) 2.66 (0.93–7.54) 0.43 (0.33–0.56)
 GPS 0.70 (0.60–0.79) 68.7 (59.0–77.0) 88.7 (78.5–94.4) 6.08 (2.99–12.38) 0.83 (0.45–1.5) 0.17 (0.09–0.3)
Table 3.
 
Sensitivity at 80% and 90% Specificity for the MRA and GPS
Table 3.
 
Sensitivity at 80% and 90% Specificity for the MRA and GPS
AUROC (95% CI) Sensitivity (95% CI)
At 80% Specificity At 90% Specificity
Global
 MRA 0.82 (0.75–0.88) 67% (59%–74%) 60% (52%–68%)
 GPS 0.86 (0.80–0.92) 70% (62%–77%) 67% (59%–74%)
Temporal superior
 MRA 0.79 (0.73–0.86) 60% (52%–68%) 45% (37%–53%)
 GPS 0.87 (0.81–0.93) 77% (70%–83%) 70% (62%–77%)
Temporal inferior
 MRA 0.86 (0.81–0.92) 76% (68%–82%) 68% (60%–75%)
 GPS 0.87 (0.81–0.93) 76% (68%–82%) 71% (63%–78%)
Nasal superior
 MRA 0.80 (0.73–0.87) 72% (64%–79%) 58% (50%–65%)
 GPS 0.86 (0.80–0.92) 75% (67%–81%) 64% (56%–71%)
Nasal superior
 MRA 0.81 (0.74–0.87) 74% (66%–81%) 68% (60%–75%)
 GPS 0.86 (0.80–0.92) 76% (68%–82%) 67% (59%–74%)
Table 4.
 
GEE Logistic Regression Modeling Sensitivity of MRA and GPS
Table 4.
 
GEE Logistic Regression Modeling Sensitivity of MRA and GPS
Variable Coef. Global Temporal Superior Temporal Inferior Nasal Superior Nasal Inferior
Est. SE P Est. SE P Est. SE P Est. SE P Est. SE P
Constant 1.11 0.27 <0.001 0.87 0.29 0.003 1.24 0.30 <0.001 0.91 0.26 <0.001 1.02 0.29 <0.001
AGIS β1 0.17 0.06 .007 0.17 0.07 0.015 0.21 0.08 .008 0.16 0.06 0.007 0.18 0.07 .007
Disc area β2 1.34 0.77 0.081 2.46 0.86 0.004 1.49 0.76 .051 1.65 0.78 0.035 1.85 0.84 .027
Test-type β3 −0.21 0.31 0.494 −2.11 0.34 <0.001 −1.27 0.30 <0.001 −1.60 0.32 <0.001 −1.29 0.29 <0.001
AGIS × disc area β4 −0.07 0.11 0.536 −0.03 0.09 0.738 0.05 0.10 0.572 −0.08 0.10 0.418 −0.02 0.11 .835
AGIS × test-type β5 −0.04 0.08 0.577 0.00 0.08 0.978 −0.09 0.08 .0264 −0.04 0.07 0.601 −0.05 0.08 .558
Disc area × test-type β6 0.35 0.79 0.655 −0.97 1.03 0.349 −0.53 0.65 0.409 −0.47 0.82 0.566 −0.49 0.74 .509
Figure 1.
 
After adjustment for disc size and AGIS score, the sensitivity of the GPS was significantly higher than the MRA when the internal normative database was used to classify eyes as ONL. With the temporal inferior results as an example, the illustration shows that MRA (A) and GPS (B) sensitivity both increased with increasing disc area and more severe visual field damage (higher AGIS scores) and that GPS sensitivity tended to be higher.
Figure 1.
 
After adjustment for disc size and AGIS score, the sensitivity of the GPS was significantly higher than the MRA when the internal normative database was used to classify eyes as ONL. With the temporal inferior results as an example, the illustration shows that MRA (A) and GPS (B) sensitivity both increased with increasing disc area and more severe visual field damage (higher AGIS scores) and that GPS sensitivity tended to be higher.
Table 5.
 
GEE Logistic Regression Modeling Sensitivity of MRA and GPS at 90% Specificity
Table 5.
 
GEE Logistic Regression Modeling Sensitivity of MRA and GPS at 90% Specificity
Variable Coef. Global Temporal Superior Temporal Inferior Nasal Superior Nasal Inferior
Est. SE P Est. SE P Est. SE P Est. SE P Est. SE P
Constant 0.87 0.25 0.001 1.10 0.30 <0.001 1.08 0.28 0.000 0.64 0.23 0.006 0.66 0.24 0.007
AGIS β1 0.18 0.06 0.004 0.19 0.09 0.026 0.18 0.07 0.005 0.12 0.05 0.019 0.13 0.05 0.015
Disc area β2 1.58 0.75 0.035 2.11 0.77 0.006 1.77 0.78 0.024 1.25 0.63 0.047 1.57 0.72 0.030
Test-type β3 −0.37 0.30 0.227 −1.23 0.29 <0.001 −0.21 0.29 0.472 −0.23 0.30 0.445 −0.01 0.30 0.986
AGIS × disc area β4 −0.05 0.11 0.655 0.11 0.13 0.406 0.02 0.12 0.853 −0.10 0.09 0.265 −0.11 0.10 0.279
AGIS × test type β5 −0.09 0.08 0.270 −0.03 0.08 0.685 −0.03 0.08 0.666 −0.07 0.07 0.271 −0.11 0.07 0.101
Disc area × test type β6 0.90 0.87 0.304 0.01 0.67 0.990 −0.79 0.75 0.290 0.90 0.77 0.241 0.92 0.80 0.247
Figure 2.
 
After adjustment for disc size and AGIS score, the sensitivity (at a fixed specificity of 90%) for glaucoma detection of the GPS was similar to that of the MRA. With the temporal inferior results as an example, the illustration shows that MRA (A) and GPS (B) sensitivity both increased with increasing disc area and more severe visual field damage (higher AGIS score values) and that the sensitivity was similar for both classification systems.
Figure 2.
 
After adjustment for disc size and AGIS score, the sensitivity (at a fixed specificity of 90%) for glaucoma detection of the GPS was similar to that of the MRA. With the temporal inferior results as an example, the illustration shows that MRA (A) and GPS (B) sensitivity both increased with increasing disc area and more severe visual field damage (higher AGIS score values) and that the sensitivity was similar for both classification systems.
ZangwillLM, MedeirosFA, BowdC, WeinrebRN. Optic nerve imaging devices: recent advances.GrehnF StampherR eds. Essentials in Ophthalmology: Glaucoma. 2004;63–91.Springer-Verlag Heidelberg, Germany.
GreaneyMJ, HoffmanDC, Garway-HeathDF, et al. Comparison of optic nerve imaging methods to distinguish normal eyes from those with glaucoma. Invest Ophthalmol Vis Sci. 2002;43:140–145. [PubMed]
MedeirosFA, ZangwillLM, BowdC, WeinrebRN. Comparison of the GDx VCC Scanning Laser Polarimeter, HRT II Confocal Scanning Laser Ophthalmoscope, and Stratus OCT Optical Coherence Tomograph for the detection of glaucoma. Arch Ophthalmol. 2004;122:827–837. [CrossRef] [PubMed]
ZangwillLM, BowdC, BerryCC, et al. Discriminating between normal and glaucomatous eyes using the Heidelberg Retina Tomograph, GDx Nerve Fiber Analyzer, and Optical Coherence Tomograph. Arch Ophthalmol. 2001;119:985–993. [CrossRef] [PubMed]
IesterM, JonasJB, MardinCY, BuddeWM. Discriminant analysis models for early detection of glaucomatous optic disc changes. Br J Ophthalmol. 2000;84:464–468. [CrossRef] [PubMed]
IesterM, MikelbergFS, DranceSM. The effect of optic disc size on diagnostic precision with the Heidelberg retina tomograph. Ophthalmology. 1997;104:545–548. [CrossRef] [PubMed]
BowdC, ZangwillLM, BlumenthalEZ, et al. Imaging of the optic disc and retinal nerve fiber layer: the effects of age, optic disc area, refractive error, and gender. J Opt Soc Am A Opt Image Sci Vis. 2002;19:197–207. [CrossRef] [PubMed]
MedeirosFA, ZangwillLM, BowdC, et al. Influence of disease severity and optic disc size on the diagnostic performance of imaging instruments in glaucoma. Invest Ophthalmol Vis Sci. 2006;47:1008–1015. [CrossRef] [PubMed]
WollsteinG, Garway-HeathDF, HitchingsRA. Identification of early glaucoma cases with the scanning laser ophthalmoscope. Ophthalmology. 1998;105:1557–1563. [CrossRef] [PubMed]
BathijaR, ZangwillL, BerryCC, et al. Detection of early glaucomatous structural damage with confocal scanning laser tomography. J Glaucoma. 1998;7:121–127. [PubMed]
FordBA, ArtesPH, McCormickTA, et al. Comparison of data analysis tools for detection of glaucoma with the Heidelberg Retina Tomograph. Ophthalmology. 2003;110:1145–1150. [CrossRef] [PubMed]
MardinCY, HornFK. Influence of optic disc size on the sensitivity of the Heidelberg Retina Tomograph. Graefes Arch Clin Exp Ophthalmol. 1998;236:641–645. [CrossRef] [PubMed]
LeisenringW, PepeMS, LongtonG. A marginal regression modelling framework for evaluating medical diagnostic tests. Stat Med. 1997;16:1263–1281. [CrossRef] [PubMed]
LiangKY, ZegerSL. Longitudinal data-analysis using generalized linear-models. Biometrika. 1986;73:13–22. [CrossRef]
MartusP, StrouxA, JunemannAM, et al. GEE approaches to marginal regression models for medical diagnostic tests. Stat Med. 2004;23:1377–1398. [CrossRef] [PubMed]
IesterM, MikelbergFS, CourtrightP, et al. Interobserver variability of optic disk variables measured by confocal scanning laser tomography. Am J Ophthalmol. 2001;132:57–62. [CrossRef] [PubMed]
MigliorS, AlbeE, GuareschiM, et al. Intraobserver and interobserver reproducibility in the evaluation of optic disc stereometric parameters by Heidelberg Retina Tomograph. Ophthalmology. 2002;109:1072–1077. [CrossRef] [PubMed]
SwindaleNV, StjepanovicG, ChinA, MikelbergFS. Automated analysis of normal and glaucomatous optic nerve head topography images. Invest Ophthalmol Vis Sci. 2000;41:1730–1742. [PubMed]
Advanced Glaucoma Intervention Study. 2. Visual field test scoring and reliability. Ophthalmology. 1994;101:1445–1455. [CrossRef] [PubMed]
NewcombeRG. Two-sided confidence intervals for the single proportion: comparison of seven methods. Stat Med. 1998;17:857–872. [CrossRef] [PubMed]
WilsonE. Probable inference, the law of succession, and statistical inference. J Am Stat Assoc. 1927;22:209–217. [CrossRef]
SimelDL, SamsaGP, MatcharDB. Likelihood ratios with confidence: sample size estimation for diagnostic test studies. J Clin Epidemiol. 1991;44:763–770. [CrossRef] [PubMed]
DeLongER, DeLongDM, Clarke-PearsonDL. Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach. Biometrics. 1988;44:837–845. [CrossRef] [PubMed]
JaeschkeR, GuyattGH, SackettDL. Users’ guides to the medical literature. III. How to use an article about a diagnostic test B—what are the results and will they help me in caring for my patients? The Evidence-Based Medicine Working Group. JAMA. 1994;271:703–707. [CrossRef] [PubMed]
HarizmanN, ZelefskyJR, IlitchevE, et al. Detection of glaucoma using operator-dependent versus operator-independent classification in the Heidelberg retinal tomograph-III. Br J Ophthalmol. 2006;90:1390–1392. [CrossRef] [PubMed]
CoopsA, HensonDB, KwartzAJ, ArtesPH. Automated analysis of heidelberg retina tomograph optic disc images by glaucoma probability score. Invest Ophthalmol Vis Sci. 2006;47:5348–5355. [CrossRef] [PubMed]
LanglotzCP. Fundamental measures of diagnostic examination performance: usefulness for clinical decision making and research. Radiology. 2003;228:3–9. [CrossRef] [PubMed]
Figure 1.
 
After adjustment for disc size and AGIS score, the sensitivity of the GPS was significantly higher than the MRA when the internal normative database was used to classify eyes as ONL. With the temporal inferior results as an example, the illustration shows that MRA (A) and GPS (B) sensitivity both increased with increasing disc area and more severe visual field damage (higher AGIS scores) and that GPS sensitivity tended to be higher.
Figure 1.
 
After adjustment for disc size and AGIS score, the sensitivity of the GPS was significantly higher than the MRA when the internal normative database was used to classify eyes as ONL. With the temporal inferior results as an example, the illustration shows that MRA (A) and GPS (B) sensitivity both increased with increasing disc area and more severe visual field damage (higher AGIS scores) and that GPS sensitivity tended to be higher.
Figure 2.
 
After adjustment for disc size and AGIS score, the sensitivity (at a fixed specificity of 90%) for glaucoma detection of the GPS was similar to that of the MRA. With the temporal inferior results as an example, the illustration shows that MRA (A) and GPS (B) sensitivity both increased with increasing disc area and more severe visual field damage (higher AGIS score values) and that the sensitivity was similar for both classification systems.
Figure 2.
 
After adjustment for disc size and AGIS score, the sensitivity (at a fixed specificity of 90%) for glaucoma detection of the GPS was similar to that of the MRA. With the temporal inferior results as an example, the illustration shows that MRA (A) and GPS (B) sensitivity both increased with increasing disc area and more severe visual field damage (higher AGIS score values) and that the sensitivity was similar for both classification systems.
Table 1.
 
Demographic and Ocular Characteristics of the Glaucoma and Healthy Control Study Groups
Table 1.
 
Demographic and Ocular Characteristics of the Glaucoma and Healthy Control Study Groups
Glaucoma Eyes (n = 99) Normal Eyes (n = 62) P
Age (y)* 67.9 (65.8, 70.0) 57.7 (55.2, 60.3) <0.0001
Gender (% male) 52% 26% 0.001
Visual field pattern standard deviation (dB)* 5.7 (5.0, 6.5) 1.5 (1.4, 1.6) <0.0001
Visual field mean deviation (dB)* −5.2 (−6.2, −4.2) −0.5 (−0.7,−0.2) <0.0001
AGIS Score* 3.9 (3.1, 4.7) .02 (−0.02, 0.04) <0.0001
Optic disc size* (mm2) 1.97 (1.86, 2.07) 1.76 (1.67, 1.85) 0.006
Table 2.
 
AUROC, Sensitivity, Specificity, and Likelihood Ratios of the MRA and GPS as Categorical Variables
Table 2.
 
AUROC, Sensitivity, Specificity, and Likelihood Ratios of the MRA and GPS as Categorical Variables
AUROC (95% CI) Sensitivity (95% CI) Specificity (95% CI) Likelihood Ratio
Outside Normal Limits (95% CI) Borderline (95% CI) Within Normal Limits (95% CI)
Global
 MRA 0.74 (0.66–0.83) 67.7 (58.0–76.1) 88.7 (78.5–94.4) 5.99 (2.94–12.20) 1.06 (0.52–2.17) 0.20 (0.12–0.34)
 GPS 0.70 (0.61–0.80) 71.7 (62.2–79.7) 82.3 (71.0–89.8) 4.04 (2.33–7.00) 0.91 (0.48–1.71) 0.14 (0.07–0.28)
Temporal superior
 MRA 0.67 (0.60–0.73) 25.3 (17.7–34.6) 100 (94.2–100.0) Inf NA 1.96 (0.89–4.33) 0.59 (0.48–0.72)
 GPS 0.68 (0.58–0.77) 64.7 (54.8–73.4) 93.6 (84.6–97.5) 10.02 (3.84–26.14) 0.80 (0.47–1.35) 0.18 (0.1–0.32)
Temporal inferior
 MRA 0.78 (0.71–0.85) 48.5 (38.9–58.2) 98.4 (91.4–99.7) 30.11 (4.19–218.80) 2.63 (1.04–1.61) 0.33 (0.24–0.45)
 GPS 0.69 (0.59–0.78) 72.7 (63.2–80.5) 88.7 (78.5–94.4) 6.44 (3.17–13.08) 0.66 (0.35–1.25) 0.18 (0.09–0.34)
Nasal superior
 MRA 0.72 (0.66–0.78) 35.4 (26.6–45.2) 98.4 (91.4–99.7) 21.96 (3.08–155.96) 3.54 (1.08–11.61) 0.50 (0.4–0.63)
 GPS 0.74 (0.65–0.83) 67.7 (58.0–76.1) 85.5 (74.7–92.2) 4.66 (2.51–8.66) 1.19 (0.61–2.3) 0.16 (0.09–0.29)
Nasal inferior
 MRA 0.74 (0.68–0.81) 43.4 (34.1–53.3) 96.8 (89.0–99.1) 13.45 (3.38–53.62) 2.66 (0.93–7.54) 0.43 (0.33–0.56)
 GPS 0.70 (0.60–0.79) 68.7 (59.0–77.0) 88.7 (78.5–94.4) 6.08 (2.99–12.38) 0.83 (0.45–1.5) 0.17 (0.09–0.3)
Table 3.
 
Sensitivity at 80% and 90% Specificity for the MRA and GPS
Table 3.
 
Sensitivity at 80% and 90% Specificity for the MRA and GPS
AUROC (95% CI) Sensitivity (95% CI)
At 80% Specificity At 90% Specificity
Global
 MRA 0.82 (0.75–0.88) 67% (59%–74%) 60% (52%–68%)
 GPS 0.86 (0.80–0.92) 70% (62%–77%) 67% (59%–74%)
Temporal superior
 MRA 0.79 (0.73–0.86) 60% (52%–68%) 45% (37%–53%)
 GPS 0.87 (0.81–0.93) 77% (70%–83%) 70% (62%–77%)
Temporal inferior
 MRA 0.86 (0.81–0.92) 76% (68%–82%) 68% (60%–75%)
 GPS 0.87 (0.81–0.93) 76% (68%–82%) 71% (63%–78%)
Nasal superior
 MRA 0.80 (0.73–0.87) 72% (64%–79%) 58% (50%–65%)
 GPS 0.86 (0.80–0.92) 75% (67%–81%) 64% (56%–71%)
Nasal superior
 MRA 0.81 (0.74–0.87) 74% (66%–81%) 68% (60%–75%)
 GPS 0.86 (0.80–0.92) 76% (68%–82%) 67% (59%–74%)
Table 4.
 
GEE Logistic Regression Modeling Sensitivity of MRA and GPS
Table 4.
 
GEE Logistic Regression Modeling Sensitivity of MRA and GPS
Variable Coef. Global Temporal Superior Temporal Inferior Nasal Superior Nasal Inferior
Est. SE P Est. SE P Est. SE P Est. SE P Est. SE P
Constant 1.11 0.27 <0.001 0.87 0.29 0.003 1.24 0.30 <0.001 0.91 0.26 <0.001 1.02 0.29 <0.001
AGIS β1 0.17 0.06 .007 0.17 0.07 0.015 0.21 0.08 .008 0.16 0.06 0.007 0.18 0.07 .007
Disc area β2 1.34 0.77 0.081 2.46 0.86 0.004 1.49 0.76 .051 1.65 0.78 0.035 1.85 0.84 .027
Test-type β3 −0.21 0.31 0.494 −2.11 0.34 <0.001 −1.27 0.30 <0.001 −1.60 0.32 <0.001 −1.29 0.29 <0.001
AGIS × disc area β4 −0.07 0.11 0.536 −0.03 0.09 0.738 0.05 0.10 0.572 −0.08 0.10 0.418 −0.02 0.11 .835
AGIS × test-type β5 −0.04 0.08 0.577 0.00 0.08 0.978 −0.09 0.08 .0264 −0.04 0.07 0.601 −0.05 0.08 .558
Disc area × test-type β6 0.35 0.79 0.655 −0.97 1.03 0.349 −0.53 0.65 0.409 −0.47 0.82 0.566 −0.49 0.74 .509
Table 5.
 
GEE Logistic Regression Modeling Sensitivity of MRA and GPS at 90% Specificity
Table 5.
 
GEE Logistic Regression Modeling Sensitivity of MRA and GPS at 90% Specificity
Variable Coef. Global Temporal Superior Temporal Inferior Nasal Superior Nasal Inferior
Est. SE P Est. SE P Est. SE P Est. SE P Est. SE P
Constant 0.87 0.25 0.001 1.10 0.30 <0.001 1.08 0.28 0.000 0.64 0.23 0.006 0.66 0.24 0.007
AGIS β1 0.18 0.06 0.004 0.19 0.09 0.026 0.18 0.07 0.005 0.12 0.05 0.019 0.13 0.05 0.015
Disc area β2 1.58 0.75 0.035 2.11 0.77 0.006 1.77 0.78 0.024 1.25 0.63 0.047 1.57 0.72 0.030
Test-type β3 −0.37 0.30 0.227 −1.23 0.29 <0.001 −0.21 0.29 0.472 −0.23 0.30 0.445 −0.01 0.30 0.986
AGIS × disc area β4 −0.05 0.11 0.655 0.11 0.13 0.406 0.02 0.12 0.853 −0.10 0.09 0.265 −0.11 0.10 0.279
AGIS × test type β5 −0.09 0.08 0.270 −0.03 0.08 0.685 −0.03 0.08 0.666 −0.07 0.07 0.271 −0.11 0.07 0.101
Disc area × test type β6 0.90 0.87 0.304 0.01 0.67 0.990 −0.79 0.75 0.290 0.90 0.77 0.241 0.92 0.80 0.247
×
×

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.

×