Pooled RRs with corresponding 95% CIs were calculated to evaluate the strength of relationship between body size and AMD incidence. Heterogeneity of effect sizes in the overall aggregations was assessed using Cochran
Q test and
I2 statistic, which denotes the proportion of between-study heterogeneity. The RRs were pooled using the random-effects model when the
P value for heterogeneity was lower than 0.1 or
I2 values were greater than 50%; otherwise, fixed-effects model was used to compute summary effect. Stratified analyses were performed to investigate the effect of AMD subtypes, BMI and AMD ascertainments, duration of follow-up, source country, AMD classification criteria, and smoking status on the relationship between body size and AMD. Dose–response meta-analysis was performed by the method described by Greenland and Longnecker
20 and Orsini et al.
21 to evaluate a potential trend between BMI levels and AMD risk. The estimated midpoint of each BMI category was utilized if mean or median was not reported in the study. The open-ended categories were assumed to have the same width as their adjacent categories. Because the references of exposure were different across studies, we used centered dose levels (each nonreference dose minus the reference dose within a study) for summarizing dose–response relation.
22 In this two-stage analysis, we first estimated restricted cubic spline model with three knots in settled percentiles (25%, 50%, and 75%) of the distribution, assuming the fixed-effects model. Then, the GLST command with the generalized least-squares regression, which required the mid values of BMI in each category, the number of patients and participants, and logarithms of RRs with 95% CIs, was used to carry out the dose–response meta-analysis.
P value for nonlinearity was calculated by testing the null hypothesis that the regression coefficients of the second-order spline transformations were equal to zero. Sensitivity analyses were performed by removing one study at a time and recalculating the pooled RR estimates for the remaining studies to evaluate whether the results were robust. Also, to ensure the robustness of the dose–response meta-analysis to different value assignment methods for the open-ended category in each study, we performed sensitivity analysis by assuming the width of open-ended categories to be 1.2, 1.5, and 2.0 times that of adjacent categories. This range of values is broad enough that the robustness of the results of this sensitivity analysis provides convincing evidence for the dose–response relationship. Publication bias was tested by the combination of a funnel plot–based method and use of Begg's test and Egger's test to estimate the number of missing studies and to calculate a corrected RR as if these studies had been present.
23,24 All analyses were performed with STATA version 12.0 (Stata Corp., College Station, TX, USA). A 2-tailed
P value of less than 0.05 was considered statistically significant.