Participant characteristics were presented by status of incident POAG. The distributions of the continuous variables were visually inspected for symmetry using density plots. All plots were adequately symmetric, and continuous variables were summarized with mean, standard deviation (SD), and range and compared using linear model ANOVA. Categorical variables were summarized by number and proportion and compared using the χ2 test.
Cox regression was used to assess the effect of SBP on the risk of incident POAG. Blood pressure was included as a time-updated variable to account for updated measurements collected through additional subsequent measurements taken during repeat visits. Participants were considered at risk from the initial assessment center visit until either diagnosis of POAG or censorship. Participants were censored by death, loss to follow-up, or end of available follow-up. The end of available follow-up was defined as the date of the latest available Showcase data update of hospital inpatient admission data (March 31, 2021).
Analyses with blood pressure represented as categorical predictors were conducted for SBP (<120, 120–130, 130–140, 140–150, ≥150 mmHg), DBP (<70, 70–80, 80–90, 90–100, ≥100 mmHg), MAP (<90, 90–100, 100–110, 110–120, ≥120 mmHg), and PP (<40, 40–50, 50–60, 60–70, ≥70 mmHg). Blood pressure was included in the model as a time-updated variable. This updates the value of the blood pressure measure for participants who had routine blood pressure recorded at subsequent study visits. Multivariable models included adjustment for age, gender, ethnicity (white, asian, black, other), BMI (underweight, healthy weight, overweight, obese), and education (secondary, tertiary, other, none of the above), Townsend deprivation index quintiles (index of material deprivation based on census variables describing car and house ownership, overcrowding, and unemployment), smoking status (never, previous, current), alcohol status (never, previous, current), diagnosis of diabetes at baseline, and prevalent cataract surgery at baseline. These covariates were included as baseline variables only and were not time-updated. Results are presented as hazard ratios (HRs) with 95% confidence intervals (CIs), including a P value for the linear trend of ordered categories for blood pressure parameters.
The second analysis included blood pressure parameters as continuous variables. SBP, DBP, MAP, and PP were included in the aforementioned multivariable model as time-updated continuous predictors modeled as restricted cubic splines with four knots placed at the 5th, 35th, 65th, and 95th percentiles, defined a priori.
21 The relative hazard ratios and 95% CIs for POAG were plotted on the
y-axis as a function of SBP, DBP, MAP, or PP on the
x-axis. The normal blood pressure values of SBP = 120 mmHg, DBP = 80 mmHg, MAP = 93.3 mmHg, and PP = 40 mmHg were used as reference values. Dashed vertical lines represent the knot locations for the restricted cubic splines.
The proportional hazards assumption was tested graphically using Schoenfeld residuals. No violations were present. A two-tailed
P value was set at 0.05 for statistical significance. Analyses were performed using R 4.1.1 (R Foundation for Statistical Computing, Vienna, Austria) with the packages
rms (v6.2-0)
22 and
survival (v3.2-11).
23
We conducted a further sensitivity analysis excluding participants who had taken antihypertensive medication anytime during the study. The antihypertensive medication classes included angiotensin-converting enzyme inhibitors, angiotensin receptor blockers, calcium channel blockers, and diuretics (thiazide and thiazide-like diuretics, loop diuretics, and potassium-sparing diuretics). The codes used to identify these medications are listed in the
Supplementary Material.