Variables potentially associated with baseline CIGTS VF score were explored by using regression analysis. Variables tested included age, race, sex, education, baseline IOP, visual acuity, mean baseline reliability score, smoking (never smoked, former cigarette smoker, current cigarette smoker, current pipe or cigar smoker), family history of glaucoma, diabetes, hypertension and/or cardiovascular disease, evidence of cataract, right or left eye, and the alertness subscale from the SIP. In addition, CIGTS clinical center and technician effects were evaluated.
Several factors were found to affect baseline VF score significantly, as shown in
Table 2 . All variables together explained 23% of the variance in VF scores. A 1-point increase in the reliability score (worse VF reliability) was associated with a 2.7-point average increase in the VF score. VF scores of the men were higher than those of the women by 0.7 on average, and blacks were higher by 1.5 units than whites. A 10-letter decrease in VA score (equivalent to two Snellen lines) was associated with an average 1.5-unit increase in VF score. Cardiovascular disease was associated with a 1.1-unit increase in mean VF score, whereas diabetes was associated with a 1.8-unit decrease. Because patients with diabetes who had any evidence of diabetic retinopathy (≥10 microaneurysms in their retina) were excluded from the study, the eligible diabetic patients had higher visual function on average than nondiabetic patients. Thus, the diabetic effect is likely to be an artifact of this eligibility criterion and is not considered further.
IOP had a more complicated relationship with VF score, partly due to the study eligibility criteria that required qualifying VF loss for patients with IOPs less than 30 mm Hg (and later, 27 mm Hg), but did not require VF loss in patients with IOPs of 30 mm Hg or higher (later, ≥27 mm Hg) if glaucomatous optic disc damage was present. In the model, increasing IOP up to 30 mm Hg was associated with decreasing VF score, whereas an increase in IOP beyond 30 mm Hg was associated with increasing VF score. The IOP effects under 30 mm Hg are probably artifacts of the eligibility criteria, because patients with no qualifying VF defect entered only if they had elevated IOP. However, the magnitude of IOP effects at pressures higher than 30 mm Hg is not related to eligibility criteria. In that range, an increase of 10 mm Hg was associated with an increase of 1.6 units in VF score. An increase of 5 units on the alertness subscale (indicating more problems with attention) was associated with an increase of 1.2 units in VF score.
CIGTS clinical center effects were also significantly associated with baseline VF scores (P < 0.0001), although technicians within clinical centers were not significantly different (P = 0.48). The proportion of variation explained increased from 16% to 23% after including clinical center in the model. The distribution of center effects (after adjusting for all other effects) was fairly normal with a SD of 1.0 VF unit; the three centers that varied the most from the mean had deviations of 2.0, −1.4, and −1.1 VF units from the adjusted mean of all centers. Although the regression assumption of normally distributed residual errors was not met because of the floor effect in the VF measurements, no transformation of the data could adequately correct the problem.
Similar modeling was performed to find variables associated with MD. Starting with the same initial list of variables used in the CIGTS VF score models, we arrived at a final model with the same variables as in the CIGTS VF score models shown in
Table 2 . The MD model explained 18% of the variation without including CIGTS clinical center effects, and 24% including center. All coefficients were in the same direction (taking the reversed scaling into account) and had magnitudes similar to those in
Table 2 . The result is not surprising, given the high correlation between MD and the CIGTS VF score. One advantage of the MD model was that the residual errors were quite normally distributed.