Using Parametric g-Computation to Estimate the Effect of Long-Term Exposure to Air Pollution on Mortality Risk and Simulate the Benefits of Hypothetical Policies: The Canadian Community Health Survey Cohort (2005 to 2015)
Authors: Chen Chen, Hong Chen, Aaron van Donkelaar, Richard T. Burnett, Randall V. Martin, Li Chen, Michael Tjepkema, Megan Kirby-McGregor, Yi Li, Jay S. Kaufman, Tarik Benmarhnia
Journal: Environmental Health Perspectives (2023)

Background: Numerous epidemiological studies have documented the adverse health impact of long-term exposure to fine particulate matter [particulate matter ≤2.5μm in aerodynamic diameter (PM2.5)] on mortality even at relatively low levels. However, methodological challenges remain to consider potential regulatory intervention's complexity and provide actionable evidence on the predicted benefits of interventions. We propose the parametric g-computation as an alternative analytical approach to such challenges. Method: We applied the parametric g-computation to estimate the cumulative risks of nonaccidental death under different hypothetical intervention strategies targeting long-term exposure to PM2.5 in the Canadian Community Health Survey cohort from 2005 to 2015. On both relative and absolute scales, we explored the benefits of hypothetical intervention strategies compared with the natural course that a) set the simulated exposure value at each follow-up year to a threshold value if exposure was above the threshold (8.8 μg/m3, 7.04 μg/m3, 5 μg/m3, and 4 μg/m3), and b) reduced the simulated exposure value by a percentage (5% and 10%) at each follow-up year. We used the 3-y average PM2.5 concentration with 1-y lag at the postal code of respondents' annual mailing addresses as their long-term exposure to PM2.5. We considered baseline and time-varying confounders, including demographics, behavior characteristics, income level, and neighborhood socioeconomic status. We also included the R syntax for reproducibility and replication. Results: All hypothetical intervention strategies explored led to lower 11-y cumulative mortality risks than the estimated value under the natural course without intervention, with the smallest reduction of 0.20 per 1,000 participants (95% CI: 0.06, 0.34) under the threshold of 8.8 μg/m3, and the largest reduction of 3.40 per 1,000 participants (95% CI: −0.23, 7.03) under the relative reduction of 10% per interval. The reductions in cumulative risk, or numbers of deaths that would have been prevented if the intervention was employed instead of maintaining the status quo, increased over time but flattened toward the end of the follow-up period. Estimates among those ≥65 years of age were greater with a similar pattern. Our estimates were robust to different model specifications. Discussion: We found evidence that any intervention further reducing the long-term exposure to PM2.5 would reduce the cumulative mortality risk, with greater benefits in the older population, even in a population already exposed to low levels of ambient PM2.5. The parametric g-computation used in this study provides flexibilities in simulating real-world interventions, accommodates time-varying exposure and confounders, and estimates adjusted survival curves with clearer interpretation and more information than a single hazard ratio, making it a valuable analytical alternative in air pollution epidemiological research. https://doi.org/10.1289/EHP11095