1
|
Spiliopoulos L, Sørensen AIV, Bager P, Nielsen NM, Hansen JV, Koch A, Meder IK, Videbech P, Ethelberg S, Hviid A. Postacute symptoms 4 months after SARS-CoV-2 infection during the Omicron period: a nationwide Danish questionnaire study. Am J Epidemiol 2024; 193:1106-1114. [PMID: 37981717 PMCID: PMC11299024 DOI: 10.1093/aje/kwad225] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2022] [Revised: 10/27/2023] [Accepted: 11/14/2023] [Indexed: 11/21/2023] Open
Abstract
Postacute symptoms are not uncommon after severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection with pre-Omicron variants. How the Omicron variant and coronavirus disease 2019 (COVID-19) booster vaccination influence the risk of postacute symptoms is less clear. We analyzed data from a nationwide Danish questionnaire study, EFTER-COVID, comprising 36 109 individuals aged ≥15 years who were tested between July 2021 and January 2022, to evaluate the associations of the Omicron variant and COVID-19 booster vaccination with postacute symptoms and new-onset general health problems 4 months after infection with SARS-CoV-2. Risk differences (RDs) were estimated by comparing Omicron cases with controls, comparing Omicron cases with Delta cases, and comparing Omicron cases vaccinated with 3 doses with those vaccinated with 2 doses, adjusting for age, sex, body mass index, self-reported chronic diseases, Charlson comorbidity index, health-care occupation, and vaccination status. Four months after testing for SARS-CoV-2 during the Omicron period, cases experienced substantial postacute symptoms and new-onset health problems in comparison with controls; the largest RD was observed for memory issues (RD = 7.4%; 95% CI, 6.4-8.3). However, risks were generally lower than those in the Delta period, particularly for dysosmia (RD = -15.0%; 95% CI, -17.0 to -13.2) and dysgeusia (RD = -11.2%; 95% CI, -13.2 to -9.5). Booster vaccination was associated with fewer postacute symptoms and new-onset health problems 4 months after Omicron infection as compared with 2 doses of COVID-19 vaccine.
Collapse
Affiliation(s)
- Lampros Spiliopoulos
- Department of Epidemiology Research, Statens Serum Institut, 2300 Copenhagen, Denmark
| | - Anna Irene Vedel Sørensen
- Department of Infectious Disease Epidemiology and Prevention, Statens Serum Institut, 2300 Copenhagen, Denmark
| | - Peter Bager
- Department of Epidemiology Research, Statens Serum Institut, 2300 Copenhagen, Denmark
- Department of Infectious Disease Epidemiology and Prevention, Statens Serum Institut, 2300 Copenhagen, Denmark
| | - Nete Munk Nielsen
- Department of Epidemiology Research, Statens Serum Institut, 2300 Copenhagen, Denmark
- Focused Research Unit in Neurology, Department of Neurology, Hospital of Southern Jutland, University of Southern Denmark, 6200 Aabenraa, Denmark
| | - Jørgen Vinsløv Hansen
- Department of Epidemiology Research, Statens Serum Institut, 2300 Copenhagen, Denmark
| | - Anders Koch
- Department of Infectious Disease Epidemiology and Prevention, Statens Serum Institut, 2300 Copenhagen, Denmark
- Department of Infectious Diseases, Rigshospitalet University Hospital, 2100 Copenhagen, Denmark
- Global Health Section, Department of Public Health, Faculty of Health and Medical Sciences, University of Copenhagen, 1014 Copenhagen, Denmark
| | - Inger Kristine Meder
- Department of Epidemiology Research, Statens Serum Institut, 2300 Copenhagen, Denmark
| | - Poul Videbech
- Centre for Neuropsychiatric Depression Research, Mental Health Centre Glostrup, 2600 Glostrup, Denmark
- Institute for Clinical Medicine, Faculty of Health and Medical Sciences, University of Copenhagen, 2200 Copenhagen, Denmark
| | - Steen Ethelberg
- Department of Infectious Disease Epidemiology and Prevention, Statens Serum Institut, 2300 Copenhagen, Denmark
- Global Health Section, Department of Public Health, Faculty of Health and Medical Sciences, University of Copenhagen, 1014 Copenhagen, Denmark
| | - Anders Hviid
- Department of Epidemiology Research, Statens Serum Institut, 2300 Copenhagen, Denmark
- Pharmacovigilance Research Centre, Department of Drug Design and Pharmacology, University of Copenhagen, 2100 Copenhagen, Denmark
| |
Collapse
|
2
|
Loh WW, Ren D, West SG. Parametric g-formula for Testing Time-Varying Causal Effects: What It Is, Why It Matters, and How to Implement It in Lavaan. MULTIVARIATE BEHAVIORAL RESEARCH 2024:1-24. [PMID: 38963381 DOI: 10.1080/00273171.2024.2354228] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/05/2024]
Abstract
Psychologists leverage longitudinal designs to examine the causal effects of a focal predictor (i.e., treatment or exposure) over time. But causal inference of naturally observed time-varying treatments is complicated by treatment-dependent confounding in which earlier treatments affect confounders of later treatments. In this tutorial article, we introduce psychologists to an established solution to this problem from the causal inference literature: the parametric g-computation formula. We explain why the g-formula is effective at handling treatment-dependent confounding. We demonstrate that the parametric g-formula is conceptually intuitive, easy to implement, and well-suited for psychological research. We first clarify that the parametric g-formula essentially utilizes a series of statistical models to estimate the joint distribution of all post-treatment variables. These statistical models can be readily specified as standard multiple linear regression functions. We leverage this insight to implement the parametric g-formula using lavaan, a widely adopted R package for structural equation modeling. Moreover, we describe how the parametric g-formula may be used to estimate a marginal structural model whose causal parameters parsimoniously encode time-varying treatment effects. We hope this accessible introduction to the parametric g-formula will equip psychologists with an analytic tool to address their causal inquiries using longitudinal data.
Collapse
Affiliation(s)
- Wen Wei Loh
- Department of Methodology and Statistics, Faculty of Health, Medicine and Life Sciences (FHML), Maastricht University, Maastricht, The Netherlands
- Department of Quantitative Theory and Methods, Emory University, Atlanta, GA, USA
| | - Dongning Ren
- Department of Work and Social Psychology, Maastricht University, Maastricht, The Netherlands
| | - Stephen G West
- Department of Psychology, Arizona State University, Tempe, AZ, USA
| |
Collapse
|
3
|
Chiu YH, Wen L, McGrath S, Logan R, Dahabreh IJ, Hernán MA. Evaluating Model Specification When Using the Parametric G-Formula in the Presence of Censoring. Am J Epidemiol 2023; 192:1887-1895. [PMID: 37338985 PMCID: PMC11043789 DOI: 10.1093/aje/kwad143] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2022] [Revised: 03/13/2023] [Accepted: 06/15/2023] [Indexed: 06/22/2023] Open
Abstract
The noniterative conditional expectation (NICE) parametric g-formula can be used to estimate the causal effect of sustained treatment strategies. In addition to identifiability conditions, the validity of the NICE parametric g-formula generally requires the correct specification of models for time-varying outcomes, treatments, and confounders at each follow-up time point. An informal approach for evaluating model specification is to compare the observed distributions of the outcome, treatments, and confounders with their parametric g-formula estimates under the "natural course." In the presence of loss to follow-up, however, the observed and natural-course risks can differ even if the identifiability conditions of the parametric g-formula hold and there is no model misspecification. Here, we describe 2 approaches for evaluating model specification when using the parametric g-formula in the presence of censoring: 1) comparing factual risks estimated by the g-formula with nonparametric Kaplan-Meier estimates and 2) comparing natural-course risks estimated by inverse probability weighting with those estimated by the g-formula. We also describe how to correctly compute natural-course estimates of time-varying covariate means when using a computationally efficient g-formula algorithm. We evaluate the proposed methods via simulation and implement them to estimate the effects of dietary interventions in 2 cohort studies.
Collapse
Affiliation(s)
- Yu-Han Chiu
- Correspondence to Dr. Yu-Han Chiu, Department of Public Health Sciences, College of Medicine, Pennsylvania State University, Hershey, PA 17033 (e-mail: )
| | | | | | | | | | | |
Collapse
|
4
|
Malekifar P, Nedjat S, Abdollahpour I, Nazemipour M, Malekifar S, Mansournia MA. Impact of Alcohol Consumption on Multiple Sclerosis Using Model-based Standardization and Misclassification Adjustment Via Probabilistic Bias Analysis. ARCHIVES OF IRANIAN MEDICINE 2023; 26:567-574. [PMID: 38310413 PMCID: PMC10862089 DOI: 10.34172/aim.2023.83] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/07/2023] [Accepted: 08/06/2023] [Indexed: 02/05/2024]
Abstract
BACKGROUND The etiology of multiple sclerosis (MS) is still not well-demonstrated, and assessment of some risk factors like alcohol consumption has problems like confounding and measurement bias. To determine the causal effect of alcohol consumption on MS after adjusting for alcohol consumption misclassification bias and confounders. METHODS In a population-based incident case-control study, 547 patients with MS and 1057 healthy people were recruited. A minimally sufficient adjustment set of confounders was derived using the causal directed acyclic graph. The probabilistic bias analysis method (PBAM) using beta, logit-logistic, and triangular probability distributions for sensitivity/specificity to adjust for misclassification bias in self-reporting alcohol consumption and model-based standardization (MBS) to estimate the causal effect of alcohol consumption were used. Population attributable fraction (PAF) estimates with 95% Monte Carlo sensitivity analysis (MCSA) intervals were calculated using PBAM and MBS analysis. Bootstrap was used to deal with random errors. RESULTS The adjusted risk ratio (95% MCSA interval) from the probabilistic bias analysis and MBS between alcohol consumption and MS using the three distribution was in the range of 1.93 (1.07 to 4.07) to 2.02 (1.15 to 4.69). The risk difference (RD) in all three scenarios was 0.0001 (0.0000 to 0.0005) and PAF was in the range of 0.15 (0.010 to 0.50) to 0.17 (0.001 to 0.47). CONCLUSION After adjusting for measurement bias, confounding, and random error alcohol consumption had a positive causal effect on the incidence of MS.
Collapse
Affiliation(s)
- Pooneh Malekifar
- Department of Epidemiology and Biostatistics, School of Public Health, Tehran University of Medical Sciences, Tehran, Iran
| | - Saharnaz Nedjat
- Department of Epidemiology and Biostatistics, School of Public Health, Tehran University of Medical Sciences, Tehran, Iran
| | - Ibrahim Abdollahpour
- Child Growth and Development Research Center, Research Institute for Primordial Prevention of Non-Communicable Disease, Isfahan University of Medical Science, Isfahan, Iran
| | - Maryam Nazemipour
- Department of Epidemiology and Biostatistics, School of Public Health, Tehran University of Medical Sciences, Tehran, Iran
| | - Saeed Malekifar
- Department of Computer Engineering, Iran University of Science and Technology, Tehran, Iran
| | - Mohammad Ali Mansournia
- Department of Epidemiology and Biostatistics, School of Public Health, Tehran University of Medical Sciences, Tehran, Iran
| |
Collapse
|
5
|
Jun I, Cohen SA, Ser SE, Marini S, Lucero RJ, Bian J, Prosperi M. Optimizing Dynamic Antibiotic Treatment Strategies against Invasive Methicillin-Resistant Staphylococcus Aureus Infections using Causal Survival Forests and G-Formula on Statewide Electronic Health Record Data. PROCEEDINGS OF MACHINE LEARNING RESEARCH 2023; 218:98-115. [PMID: 37854935 PMCID: PMC10584043] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Figures] [Subscribe] [Scholar Register] [Indexed: 10/20/2023]
Abstract
Developing models for individualized, time-varying treatment optimization from observational data with large variable spaces, e.g., electronic health records (EHR), is problematic because of inherent, complex bias that can change over time. Traditional methods such as the g-formula are robust, but must identify critical subsets of variables due to combinatorial issues. Machine learning approaches such as causal survival forests have fewer constraints and can provide fine-tuned, individualized counterfactual predictions. In this study, we aimed to optimize time-varying antibiotic treatment -identifying treatment heterogeneity and conditional treatment effects- against invasive methicillin-resistant Staphylococcus Aureus (MRSA) infections, using statewide EHR data collected in Florida, USA. While many previous studies focused on measuring the effects of the first empiric treatment (i.e., usually vancomycin), our study focuses on dynamic sequential treatment changes, comparing possible vancomycin switches with other antibiotics at clinically relevant time points, e.g., after obtaining a bacterial culture and susceptibility testing. Our study population included adult individuals admitted to the hospital with invasive MRSA. We collected demographic, clinical, medication, and laboratory information from the EHR for these patients. Then, we followed three sequential antibiotic choices (i.e., their empiric treatment, subsequent directed treatment, and final sustaining treatment), evaluating 30-day mortality as the outcome. We applied both causal survival forests and g-formula using different clinical intervention policies. We found that switching from vancomycin to another antibiotic improved survival probability, yet there was a benefit from initiating vancomycin compared to not using it at any time point. These findings show consistency with the empiric choice of vancomycin before confirmation of MRSA and shed light on how to manage switches on course. In conclusion, this application of causal machine learning on EHR demonstrates utility in modeling dynamic, heterogeneous treatment effects that cannot be evaluated precisely using randomized clinical trials.
Collapse
Affiliation(s)
- Inyoung Jun
- Department of Epidemiology, College of Public Health and Health Professions, University of Florida, Gainesville, FL 32610, USA
| | - Scott A Cohen
- Department of Epidemiology, College of Public Health and Health Professions, University of Florida, Gainesville, FL 32610, USA
| | - Sarah E Ser
- Department of Epidemiology, College of Public Health and Health Professions, University of Florida, Gainesville, FL 32610, USA
| | - Simone Marini
- Department of Epidemiology, College of Public Health and Health Professions, University of Florida, Gainesville, FL 32610, USA
| | - Robert J Lucero
- School of Nursing, University of California, Los Angeles, Los Angeles, CA 90095, USA
| | - Jiang Bian
- Department of Health Outcomes and Biomedical Informatics, College of Medicine, University of Florida, Gainesville, FL 32610, USA
| | - Mattia Prosperi
- Department of Epidemiology, College of Public Health and Health Professions, University of Florida, Gainesville, FL 32610, USA
| |
Collapse
|
6
|
Takeuchi M, Shinkawa K, Yanagita M, Kawakami K. Hypothetical intervention of targeted systolic blood pressure control of < 120 mmHg on renal prognosis for persons with stage 3-4 chronic kidney disease: an application of parametric g-formula using health checkup data in Japan. Clin Exp Nephrol 2023; 27:542-547. [PMID: 36947305 DOI: 10.1007/s10157-023-02341-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2021] [Accepted: 03/10/2023] [Indexed: 03/23/2023]
Abstract
BACKGROUND The revised KDIGO guidelines recommend maintaining systolic blood pressure (sBP) < 120 mmHg in patients with chronic kidney disease (CKD), based on cardiovascular and survival benefits. However, the renal benefit of this strategy remains less clear. METHODS We used data of routine health checkups in Japan. Persons whose estimated glomerular filtration rate (eGFR) was < 60 mL/min/1.73 m2 in 2015 without end-stage disease were followed until 2020. We estimated the 5-year benefit of hypothetical targeted sBP control using parametric g-formula modeling, accounting for both time-fixed and time-varying confounding variables. Four sensitivity analyses, including analysis using a marginal structural model (MSM) and positive control outcome analysis, were also done. RESULTS We enrolled 28,972 persons with CKD (median age: 54 years, male: 69%, baseline eGFR [median]: 56 mL/min/1.73m2). As compared with the natural course without a targeted intervention, there was no renoprotective effect of targeted sBP control, with a 5-year difference in eGFR of 0.65 mL/min/1.73 m2 (95% confidence interval - 0.42 to 1.65 mL/min/1.73 m2). MSM analysis found a similar result. In contrast, the positive control analysis using the cardiovascular outcome showed that targeted sBP control would reduce the cardiovascular disease incidence by 6.0% over 5 years. CONCLUSIONS A targeted sBP control strategy maintaining < 120 mmHg may not yield a renoprotective effect for patients with stage 3-4 CKD, although it was expected to offer a cardiovascular benefit. Future research may be warranted in higher-risk populations, such as elderly people or those with more advanced kidney disease.
Collapse
Affiliation(s)
- Masato Takeuchi
- Department of Pharmacoepidemiology, Graduate School of Medicine and Public Health, Kyoto University, Yoshidakonoe-cho, Sakyo-ku, Kyoto, 606-8501, Japan.
| | - Kanna Shinkawa
- Department of Pharmacoepidemiology, Graduate School of Medicine and Public Health, Kyoto University, Yoshidakonoe-cho, Sakyo-ku, Kyoto, 606-8501, Japan
- Department of Nephrology, Graduate School of Medicine, Kyoto University, Kyoto, Japan
| | - Motoko Yanagita
- Department of Nephrology, Graduate School of Medicine, Kyoto University, Kyoto, Japan
- Institute for the Advanced Study of Human Biology (WPI-ASHBi), Kyoto University, Kyoto, Japan
| | - Koji Kawakami
- Department of Pharmacoepidemiology, Graduate School of Medicine and Public Health, Kyoto University, Yoshidakonoe-cho, Sakyo-ku, Kyoto, 606-8501, Japan
| |
Collapse
|
7
|
Smith TJS, Keil AP, Buckley JP. Estimating Causal Effects of Interventions on Early-life Environmental Exposures Using Observational Data. Curr Environ Health Rep 2023; 10:12-21. [PMID: 36418665 DOI: 10.1007/s40572-022-00388-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 10/03/2022] [Indexed: 11/27/2022]
Abstract
PURPOSE OF REVIEW We discuss how epidemiologic studies have used observational data to estimate the effects of potential interventions on early-life environmental exposures. We summarize the value of posing questions about interventions, how a group of techniques known as "g-methods" can provide advantages for estimating intervention effects, and how investigators have grappled with the strong assumptions required for causal inference. RECENT FINDINGS We identified nine studies that estimated health effects of hypothetical interventions on early-life environmental exposures. Of these, six examined air pollution. Interventions evaluated by these studies included setting exposure levels at a specific value, shifting exposure distributions, and limiting exposure levels to less than a threshold value. Only one study linked exposure contrasts to a specific intervention on an exposure source, however. There is growing interest in estimating intervention effects of early-life environmental exposures, in part because intervention effects are directly related to possible public health actions. Future studies can build on existing work by linking research questions to specific hypothetical interventions that could reduce exposure levels. We discuss how framing questions around interventions can help overcome some of the barriers to causal inference and how advances related to machine learning may strengthen studies by sidestepping the overly restrictive assumptions of parametric regression models. By leveraging advancements in causal inference and exposure science, an intervention framework for environmental epidemiology can guide actionable solutions to improve children's environmental health.
Collapse
Affiliation(s)
- Tyler J S Smith
- Department of Environmental Health & Engineering, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD, USA
| | - Alexander P Keil
- Department of Epidemiology, University of North Carolina Gillings School of Global Public Health, Chapel Hill, NC, USA
| | - Jessie P Buckley
- Departments of Environmental Health & Engineering and Epidemiology, Johns Hopkins Bloomberg School of Public Health, 615 North Wolfe Street, W7515, Baltimore, MD, 21205, USA.
| |
Collapse
|
8
|
Ross RK, Cole SR, Westreich D, Edwards JK, Musonda P, Vwalika B, Kasaro MP, Price JT, Stringer JSA. Different effects for different questions: An illustration using short cervix and the risk of preterm birth. Int J Gynaecol Obstet 2023; 160:842-849. [PMID: 35899762 PMCID: PMC11155393 DOI: 10.1002/ijgo.14372] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2021] [Revised: 07/05/2022] [Accepted: 07/20/2022] [Indexed: 11/07/2022]
Abstract
OBJECTIVE To illustrate the difference between exposure effects and population attributable effects. METHODS We examined the effect of mid-pregnancy short cervical length (<25 mm) on preterm birth using data from a prospective cohort of pregnant women in Lusaka, Zambia. Preterm birth was live birth or stillbirth before 37 weeks of pregnancy. For estimation, we used multivariable regression and parametric g-computation. RESULTS Among 1409 women included in the analysis, short cervix was rare (2.4%); 13.6% of births were preterm. Exposure effect estimates were large (marginal risk ratio 2.86, 95% confidence interval [CI] 1.80-4.54), indicating that the preterm birth risk was substantially higher among women with a short cervix compared with women without a short cervix. However, the population attributable effect estimates were close to the null (risk ratio 1.06, 95% CI 1.02-1.10), indicating that an intervention to counteract the impact of short cervix on preterm birth would have minimal effect on the population risk of preterm birth. CONCLUSION Although authors often refer to "the" effect, there are actually different types of effects, as we have illustrated here. In planning research, it is important to consider which effect to estimate to ensure that the estimate aligns with the research objective.
Collapse
Affiliation(s)
- Rachael K Ross
- Department of Epidemiology, Gillings School of Global Public Health, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, USA
| | - Stephen R Cole
- Department of Epidemiology, Gillings School of Global Public Health, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, USA
| | - Daniel Westreich
- Department of Epidemiology, Gillings School of Global Public Health, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, USA
| | - Jessie K Edwards
- Department of Epidemiology, Gillings School of Global Public Health, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, USA
| | - Patrick Musonda
- Department of Epidemiology and Biostatistics, School of Public Health, University of Zambia, Lusaka, Zambia
| | - Bellington Vwalika
- Department of Obstetrics and Gynecology, School of Medicine, University of Zambia, Lusaka, Zambia
| | | | - Joan T Price
- Department of Obstetrics and Gynecology, School of Medicine, University of Zambia, Lusaka, Zambia
- University of North Carolina Global Projects Zambia, Lusaka, Zambia
- Department of Obstetrics and Gynecology, School of Medicine, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, USA
| | - Jeffrey S A Stringer
- Department of Epidemiology, Gillings School of Global Public Health, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, USA
- Department of Obstetrics and Gynecology, School of Medicine, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, USA
| |
Collapse
|
9
|
Chen C, Chen H, van Donkelaar A, Burnett RT, Martin RV, Chen L, Tjepkema M, Kirby-McGregor M, Li Y, Kaufman JS, Benmarhnia T. 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). ENVIRONMENTAL HEALTH PERSPECTIVES 2023; 131:37010. [PMID: 36920446 PMCID: PMC10016347 DOI: 10.1289/ehp11095] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
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.
Collapse
Affiliation(s)
- Chen Chen
- Scripps Institution of Oceanography, University of California, San Diego, La Jolla, California, USA
| | - Hong Chen
- Environmental Health Science and Research Bureau, Health Canada, Ottawa, Ontario, Canada
- Public Health Ontario, Toronto, Ontario, Canada
- ICES, Toronto, Ontario, Canada
- Dalla Lana School of Public Health, University of Toronto, Toronto, Ontario, Canada
| | - Aaron van Donkelaar
- Department of Energy, Environmental, and Chemical Engineering, Washington University in St. Louis, St. Louis, Missouri, USA
| | - Richard T. Burnett
- Environmental Health Science and Research Bureau, Health Canada, Ottawa, Ontario, Canada
| | - Randall V. Martin
- Department of Energy, Environmental, and Chemical Engineering, Washington University in St. Louis, St. Louis, Missouri, USA
| | - Li Chen
- Environmental Health Science and Research Bureau, Health Canada, Ottawa, Ontario, Canada
| | - Michael Tjepkema
- Health Analysis Division, Statistics Canada, Ottawa, Ontario, Canada
| | - Megan Kirby-McGregor
- Department of Epidemiology, Biostatistics and Occupational Health, McGill University, Montreal, Quebec, Canada
| | - Yi Li
- Department of Epidemiology, Biostatistics and Occupational Health, McGill University, Montreal, Quebec, Canada
| | - Jay S. Kaufman
- Department of Epidemiology, Biostatistics and Occupational Health, McGill University, Montreal, Quebec, Canada
| | - Tarik Benmarhnia
- Scripps Institution of Oceanography, University of California, San Diego, La Jolla, California, USA
| |
Collapse
|
10
|
Saarela O, Stephens DA, Moodie EEM. The Role of Exchangeability in Causal Inference. Stat Sci 2023. [DOI: 10.1214/22-sts879] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/24/2023]
Affiliation(s)
- Olli Saarela
- Olli Saarela is Associate Professor, Dalla Lana School of Public Health, University of Toronto, 155 College Street, Toronto, Ontario M5T 3M7, Canada
| | - David A. Stephens
- David A. Stephens is Professor, Department of Mathematics and Statistics, McGill University, Burnside Hall, 805 Sherbrooke Street West, Montreal, Quebec H3A 0B9, Canada
| | - Erica E. M. Moodie
- Erica E. M. Moodie is Professor, Department of Epidemiology and Biostatistics, McGill University, 2001 McGill College Ave, Montreal, Quebec H3A 1G1, Canada
| |
Collapse
|
11
|
Randomization, design and analysis for interdependency in aging research: no person or mouse is an island. NATURE AGING 2022; 2:1101-1111. [PMID: 37063472 PMCID: PMC10099485 DOI: 10.1038/s43587-022-00333-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Investigators traditionally use randomized designs and corresponding analysis procedures to make causal inferences about the effects of interventions, assuming independence between an individual's outcome and treatment assignment and the outcomes of other individuals in the study. Often, such independence may not hold. We provide examples of interdependency in model organism studies and human trials and group effects in aging research and then discuss methodologic issues and solutions. We group methodologic issues as they pertain to (1) single-stage individually randomized trials; (2) cluster-randomized controlled trials; (3) pseudo-cluster-randomized trials; (4) individually randomized group treatment; and (5) two-stage randomized designs. Although we present possible strategies for design and analysis to improve the rigor, accuracy and reproducibility of the science, we also acknowledge real-world constraints. Consequences of nonadherence, differential attrition or missing data, unintended exposure to multiple treatments and other practical realities can be reduced with careful planning, proper study designs and best practices.
Collapse
|
12
|
Kühne F, Schomaker M, Stojkov I, Jahn B, Conrads-Frank A, Siebert S, Sroczynski G, Puntscher S, Schmid D, Schnell-Inderst P, Siebert U. Causal evidence in health decision making: methodological approaches of causal inference and health decision science. GERMAN MEDICAL SCIENCE : GMS E-JOURNAL 2022; 20:Doc12. [PMID: 36742460 PMCID: PMC9869404 DOI: 10.3205/000314] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Figures] [Subscribe] [Scholar Register] [Received: 12/10/2021] [Indexed: 02/07/2023]
Abstract
Objectives Public health decision making is a complex process based on thorough and comprehensive health technology assessments involving the comparison of different strategies, values and tradeoffs under uncertainty. This process must be based on best available evidence and plausible assumptions. Causal inference and health decision science are two methodological approaches providing information to help guide decision making in health care. Both approaches are quantitative methods that use statistical and modeling techniques and simplifying assumptions to mimic the complexity of the real world. We intend to review and lay out both disciplines with their aims, strengths and limitations based on a combination of textbook knowledge and expert experience. Methods To help understanding and differentiating the methodological approaches of causal inference and health decision science, we reviewed both methods with the focus on aims, research questions, methods, assumptions, limitations and challenges, and software. For each methodological approach, we established a group of four experts from our own working group to carefully review and summarize each method, followed by structured discussion rounds and written reviews, in which the experts from all disciplines including HTA and medicine were involved. The entire expert group discussed objectives, strengths and limitations of both methodological areas, and potential synergies. Finally, we derived recommendations for further research and provide a brief outlook on future trends. Results Causal inference methods aim for drawing causal conclusions from empirical data on the relationship of pre-specified interventions on a specific target outcome and apply a counterfactual framework and statistical techniques to derive causal effects of exposures or interventions from these data. Causal inference is based on a causal diagram, more specifically, a directed acyclic graph (DAG), which encodes the assumptions regarding the causal relations between variables. Depending on the type of confounding and selection bias, traditional statistical methods or more complex g-methods are needed to derive valid causal effects. Besides the correct specification of the DAG and the statistical model, assumptions such as consistency, positivity, and exchangeability must be checked when aiming at causal inference. Health decision science aims for guiding policy decision making regarding health interventions considering and balancing multiple competing objectives of a decision based on data from multiple sources and studies, for example prevalence studies, clinical trials and long-term observational routine effectiveness studies, and studies on preferences and costs. It involves decision analysis, a systematic, explicit and quantitative framework to guide decisions under uncertainty. Decision analyses are based on decision-analytic models to mimic the course of disease as well as aspects and consequences of the intervention in order to quantitatively optimize the decision. Depending on the type of decision problem, decision trees, state-transition models, discrete event simulation models, dynamic transmission models, or other model types are applied. Models must be validated against observed data, and comprehensive sensitivity analyses must be performed to assess uncertainty. Besides the appropriate choice of the model type and the valid specification of the model structure, it must be checked if input parameters of effects can be interpreted as causal parameters in the model. Otherwise results will be biased. Conclusions Both causal inference and health decision science aim for providing best causal evidence for informed health decision making. The strengths and limitations of both methods differ and a good understanding of both methods is essential for correct application but also for correct interpretation of findings from the described methods. Importantly, decision-analytic modeling should be combined with causal inference when developing guidance and recommendations regarding decisions on health care interventions.
Collapse
Affiliation(s)
- Felicitas Kühne
- Institute of Public Health, Medical Decision Making and Health Technology Assessment, Department of Public Health, Health Services Research and Health Technology Assessment, UMIT TIROL – University for Health Sciences, Medical Informatics and Technology, Hall i.T., Austria
| | - Michael Schomaker
- Institute of Public Health, Medical Decision Making and Health Technology Assessment, Department of Public Health, Health Services Research and Health Technology Assessment, UMIT TIROL – University for Health Sciences, Medical Informatics and Technology, Hall i.T., Austria
- Centre for Infectious Disease Epidemiology & Research, University of Cape Town, South Africa
| | - Igor Stojkov
- Institute of Public Health, Medical Decision Making and Health Technology Assessment, Department of Public Health, Health Services Research and Health Technology Assessment, UMIT TIROL – University for Health Sciences, Medical Informatics and Technology, Hall i.T., Austria
| | - Beate Jahn
- Institute of Public Health, Medical Decision Making and Health Technology Assessment, Department of Public Health, Health Services Research and Health Technology Assessment, UMIT TIROL – University for Health Sciences, Medical Informatics and Technology, Hall i.T., Austria
- Division of Health Technology Assessment, ONCOTYROL – Center for Personalized Cancer Medicine, Innsbruck, Austria
| | - Annette Conrads-Frank
- Institute of Public Health, Medical Decision Making and Health Technology Assessment, Department of Public Health, Health Services Research and Health Technology Assessment, UMIT TIROL – University for Health Sciences, Medical Informatics and Technology, Hall i.T., Austria
| | - Silke Siebert
- Institute of Public Health, Medical Decision Making and Health Technology Assessment, Department of Public Health, Health Services Research and Health Technology Assessment, UMIT TIROL – University for Health Sciences, Medical Informatics and Technology, Hall i.T., Austria
| | - Gaby Sroczynski
- Institute of Public Health, Medical Decision Making and Health Technology Assessment, Department of Public Health, Health Services Research and Health Technology Assessment, UMIT TIROL – University for Health Sciences, Medical Informatics and Technology, Hall i.T., Austria
| | - Sibylle Puntscher
- Institute of Public Health, Medical Decision Making and Health Technology Assessment, Department of Public Health, Health Services Research and Health Technology Assessment, UMIT TIROL – University for Health Sciences, Medical Informatics and Technology, Hall i.T., Austria
| | - Daniela Schmid
- Institute of Public Health, Medical Decision Making and Health Technology Assessment, Department of Public Health, Health Services Research and Health Technology Assessment, UMIT TIROL – University for Health Sciences, Medical Informatics and Technology, Hall i.T., Austria
| | - Petra Schnell-Inderst
- Institute of Public Health, Medical Decision Making and Health Technology Assessment, Department of Public Health, Health Services Research and Health Technology Assessment, UMIT TIROL – University for Health Sciences, Medical Informatics and Technology, Hall i.T., Austria
| | - Uwe Siebert
- Institute of Public Health, Medical Decision Making and Health Technology Assessment, Department of Public Health, Health Services Research and Health Technology Assessment, UMIT TIROL – University for Health Sciences, Medical Informatics and Technology, Hall i.T., Austria
- Division of Health Technology Assessment, ONCOTYROL – Center for Personalized Cancer Medicine, Innsbruck, Austria
- Center for Health Decision Science, Departments of Epidemiology and Health Policy & Management, Harvard T.H. Chan School of Public Health, Boston, MA, USA
- Program on Cardiovascular Research, Institute for Technology Assessment and Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA
| |
Collapse
|
13
|
Impact of lowering fine particulate matter from major emission sources on mortality in Canada: A nationwide causal analysis. Proc Natl Acad Sci U S A 2022; 119:e2209490119. [PMID: 36442082 PMCID: PMC9894124 DOI: 10.1073/pnas.2209490119] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022] Open
Abstract
Emissions of fine particulate matter (PM2.5) from human activities have been linked to substantial disease burdens, but evidence regarding how reducing PM2.5 at its sources would improve public health is sparse. We followed a population-based cohort of 2.7 million adults across Canada from 2007 through 2016. For each participant, we estimated annual mean concentrations of PM2.5 and the fractional contributions to PM2.5 from the five leading anthropogenic sources at their residential address using satellite observations in combination with a global atmospheric chemistry transport model. For each source, we estimated the causal effects of six hypothetical interventions on 10-y nonaccidental mortality risk using the parametric g-formula, a structural causal model. We conducted stratified analyses by age, sex, and income. This cohort would have experienced tangible health gains had contributions to PM2.5 from any of the five sources been reduced. Compared with no intervention, a 10% annual reduction in PM2.5 contributions from transportation and power generation, Canada's largest and fifth-largest anthropogenic sources, would have prevented approximately 175 (95%CI: 123-226) and 90 (95%CI: 63-117) deaths per million by 2016, respectively. A more intensive 50% reduction per year in PM2.5 contributions from the two sources would have averted 360 and 185 deaths per million, respectively, by 2016. The potential health benefits were greater among men, older adults, and low-income earners. In Canada, where PM2.5 levels are among the lowest worldwide, reducing PM2.5 contributions from anthropogenic sources by as little as 10% annually would yield meaningful health gains.
Collapse
|
14
|
Kuehne F, Arvandi M, Hess LM, Faries DE, Matteucci Gothe R, Gothe H, Beyrer J, Zeimet AG, Stojkov I, Mühlberger N, Oberaigner W, Marth C, Siebert U. Causal analyses with target trial emulation for real-world evidence removed large self-inflicted biases: systematic bias assessment of ovarian cancer treatment effectiveness. J Clin Epidemiol 2022; 152:269-280. [PMID: 36252741 DOI: 10.1016/j.jclinepi.2022.10.005] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2022] [Revised: 08/17/2022] [Accepted: 10/03/2022] [Indexed: 11/05/2022]
Abstract
BACKGROUND AND OBJECTIVES Drawing causal conclusions from real-world data (RWD) poses methodological challenges and risk of bias. We aimed to systematically assess the type and impact of potential biases that may occur when analyzing RWD using the case of progressive ovarian cancer. METHODS We retrospectively compared overall survival with and without second-line chemotherapy (LOT2) using electronic medical records. Potential biases were determined using directed acyclic graphs. We followed a stepwise analytic approach ranging from crude analysis and multivariable-adjusted Cox model up to a full causal analysis using a marginal structural Cox model with replicates emulating a reference randomized controlled trial (RCT). To assess biases, we compared effect estimates (hazard ratios [HRs]) of each approach to the HR of the reference trial. RESULTS The reference trial showed an HR for second line vs. delayed therapy of 1.01 (95% confidence interval [95% CI]: 0.82-1.25). The corresponding HRs from the RWD analysis ranged from 0.51 for simple baseline adjustments to 1.41 (95% CI: 1.22-1.64) accounting for immortal time bias with time-varying covariates. Causal trial emulation yielded an HR of 1.12 (95% CI: 0.96-1.28). CONCLUSION Our study, using ovarian cancer as an example, shows the importance of a thorough causal design and analysis if one is expecting RWD to emulate clinical trial results.
Collapse
Affiliation(s)
- Felicitas Kuehne
- Department of Public Health, Health Services Research and Health Technology Assessment, Institute of Public Health, Medical Decision Making and Health Technology Assessment, UMIT TIROL - University for Health Sciences, Medical Informatics and Technology, Hall i.T., Austria
| | - Marjan Arvandi
- Department of Public Health, Health Services Research and Health Technology Assessment, Institute of Public Health, Medical Decision Making and Health Technology Assessment, UMIT TIROL - University for Health Sciences, Medical Informatics and Technology, Hall i.T., Austria
| | - Lisa M Hess
- Eli Lilly and Company, Indianapolis, IN, USA
| | | | - Raffaella Matteucci Gothe
- Department of Public Health, Health Services Research and Health Technology Assessment, Institute of Public Health, Medical Decision Making and Health Technology Assessment, UMIT TIROL - University for Health Sciences, Medical Informatics and Technology, Hall i.T., Austria
| | - Holger Gothe
- Department of Public Health, Health Services Research and Health Technology Assessment, Institute of Public Health, Medical Decision Making and Health Technology Assessment, UMIT TIROL - University for Health Sciences, Medical Informatics and Technology, Hall i.T., Austria; Chair of Health Sciences/Public Health, Medical Faculty "Carl Gustav Carus", Technical University Dresden, Dresden, Germany
| | | | - Alain Gustave Zeimet
- Department of Obstetrics and Gynecology, Innsbruck Medical University, Innsbruck, Austria
| | - Igor Stojkov
- Department of Public Health, Health Services Research and Health Technology Assessment, Institute of Public Health, Medical Decision Making and Health Technology Assessment, UMIT TIROL - University for Health Sciences, Medical Informatics and Technology, Hall i.T., Austria
| | - Nikolai Mühlberger
- Department of Public Health, Health Services Research and Health Technology Assessment, Institute of Public Health, Medical Decision Making and Health Technology Assessment, UMIT TIROL - University for Health Sciences, Medical Informatics and Technology, Hall i.T., Austria
| | - Willi Oberaigner
- Department of Public Health, Health Services Research and Health Technology Assessment, Institute of Public Health, Medical Decision Making and Health Technology Assessment, UMIT TIROL - University for Health Sciences, Medical Informatics and Technology, Hall i.T., Austria; Institute for Clinical Epidemiology, Cancer Registry Tyrol, Tirol Kliniken, Innsbruck, Austria
| | - Christian Marth
- Department of Obstetrics and Gynecology, Innsbruck Medical University, Innsbruck, Austria
| | - Uwe Siebert
- Department of Public Health, Health Services Research and Health Technology Assessment, Institute of Public Health, Medical Decision Making and Health Technology Assessment, UMIT TIROL - University for Health Sciences, Medical Informatics and Technology, Hall i.T., Austria; Center for Health Decision Science and Departments of Epidemiology and Health Policy & Management, Harvard T.H. Chan School of Public Health, Boston, MA, USA; Institute for Technology Assessment and Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA.
| |
Collapse
|
15
|
Grembi JA, Rogawski McQuade ET. Introducing riskCommunicator: An R package to obtain interpretable effect estimates for public health. PLoS One 2022; 17:e0265368. [PMID: 35849588 PMCID: PMC9292119 DOI: 10.1371/journal.pone.0265368] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2022] [Accepted: 06/24/2022] [Indexed: 12/02/2022] Open
Abstract
Common statistical modeling methods do not necessarily produce the most relevant or interpretable effect estimates to communicate risk. Overreliance on the odds ratio and relative effect measures limit the potential impact of epidemiologic and public health research. We created a straightforward R package, called riskCommunicator, to facilitate the presentation of a variety of effect measures, including risk differences and ratios, number needed to treat, incidence rate differences and ratios, and mean differences. The riskCommunicator package uses g-computation with parametric regression models and bootstrapping for confidence intervals to estimate effect measures in time-fixed data. We demonstrate the utility of the package using data from the Framingham Heart Study to estimate the effect of prevalent diabetes on the 24-year risk of cardiovascular disease or death. The package promotes the communication of public-health relevant effects and is accessible to a broad range of epidemiologists and health researchers with little to no expertise in causal inference methods or advanced coding.
Collapse
Affiliation(s)
- Jessica A. Grembi
- Division of Infectious Diseases and Global Health, Department of Medicine, Stanford University School of Medicine, Stanford, CA, United States of America
| | | |
Collapse
|
16
|
Li Y, Wang Q, Liang S, Feng C, Yang H, Yu H, Yuan D, Yang S. Effect of Switching Antiretroviral Treatment Regimen in Patients With Drug-Resistant HIV-1 Infection: Retrospective Observational Cohort Study. JMIR Public Health Surveill 2022; 8:e33429. [PMID: 35749212 PMCID: PMC9270715 DOI: 10.2196/33429] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2021] [Revised: 04/16/2022] [Accepted: 05/10/2022] [Indexed: 11/13/2022] Open
Abstract
BACKGROUND Evidence on the efficacy of antiretroviral therapy (ART) regimen switches on the mortality of patients with HIV drug resistance (HIVDR) is limited. OBJECTIVE We aim to provide policy guidance for ART regimen selection and evaluate the effectiveness of ART regime switches for people living with HIV and HIV-1 drug resistance. METHODS This retrospective observational cohort study included 179 people living with HIV and HIV-1 drug resistance from 2011 to 2020. The time that participants switched treatment regimens either to protease inhibitor (PI)-based ART regimens (PIs) or nonnucleoside reverse transcriptase inhibitor (NNRTI)-based ART regimens (NNRTIs) was taken as an observation starting point and followed up every 12 months. The parametric g-formula was used to estimate the 5-year risk of mortality under the situations of (1) natural course, (2) immediate switch to NNRTIs, (3) immediate switch to PIs, and (4) if CD4(+) T cells<200 switched to PIs. RESULTS The follow-up time of the 179 patients ranged from 30 to 119 months. The median follow-up time was 90 months. During a follow-up of 15,606 person-months, 27 individuals died in the cohort. The estimated 5-year risk of mortality under natural course, immediate switch to NNRTIs, immediate switch to PIs, and if CD4(+), and switch to PIs if T cells<200 were 11.62% (95% CI 7.82-17.11), 31.88% (95% CI 20.79-44.94), 2.87% (95% CI 0.32-7.07), and 5.30% (95% CI 2.07-10.21), respectively. The risk ratios (RRs) of immediate switch to NNRTIs, immediate switch to PIs, and switch to PIs if CD4(+) T cells<200, compared with natural course mortality rate, were 2.74 (95% CI 2.01-3.47), 0.25 (95% CI: 0.04-0.54), and 0.46 (95% CI 0.22-0.71), respectively. The risk differences were 20.26% (95% CI 10.96-28.61), -8.76% (95% CI -13.34 to -5.09) and -6.32% (95% CI -9.75 to -3.11), respectively. CONCLUSIONS Our study found that a PI-based ART regimen was beneficial for reducing mortality in people living with HIV and HIV-1 drug resistance. More effort should be given to find HIV-1 drug resistance earlier to ensure a timely adjustment to PI-based ART, thereby maximizing the benefit of early switch treatment for people living with HIV and HIV-1 drug resistance.
Collapse
Affiliation(s)
- Yiping Li
- Center for AIDS/STD Control and Prevention, Sichuan Center for Disease Control and Prevention, Chengdu, China
| | - Qinjian Wang
- West China School of Public Health and West China Fourth Hospital, Sichuan University, Chengdu, China
| | - Shu Liang
- Center for AIDS/STD Control and Prevention, Sichuan Center for Disease Control and Prevention, Chengdu, China
| | - Chuanteng Feng
- Institute for Disaster Management and Reconstruction, Sichuan University, Chengdu, China
| | - Hong Yang
- Center for AIDS/STD Control and Prevention, Sichuan Center for Disease Control and Prevention, Chengdu, China
| | - Hang Yu
- Center for AIDS/STD Control and Prevention, Sichuan Center for Disease Control and Prevention, Chengdu, China
| | - Dan Yuan
- Center for AIDS/STD Control and Prevention, Sichuan Center for Disease Control and Prevention, Chengdu, China
| | - Shujuan Yang
- West China School of Public Health and West China Fourth Hospital, Sichuan University, Chengdu, China.,International Institute of Spatial Lifecourse Health, Wuhan University, Wuhan, China
| |
Collapse
|
17
|
Tai AS, Huang YT, Yang HI, Lan LV, Lin SH. G-Computation to Causal Mediation Analysis With Sequential Multiple Mediators-Investigating the Vulnerable Time Window of HBV Activity for the Mechanism of HCV Induced Hepatocellular Carcinoma. Front Public Health 2022; 9:757942. [PMID: 35071157 PMCID: PMC8779208 DOI: 10.3389/fpubh.2021.757942] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2021] [Accepted: 11/30/2021] [Indexed: 12/09/2022] Open
Abstract
Regression-based approaches are widely used in causal mediation analysis. The presence of multiple mediators, however, increases the complexity and difficulty of mediation analysis. In such cases, regression-based approaches cannot efficiently address estimation issues. Hence, a flexible approach to mediation analysis is needed. Therefore, we developed a method for using g-computation algorithm to conduct causal mediation analysis in the presence of multiple ordered mediators. Compared to regression-based approaches, the proposed simulation-based approach increases flexibility in the choice of models and increases the range of the outcome scale. The Taiwanese Cohort Study dataset was used to evaluate the efficacy of the proposed approach for investigating the mediating role of early and late HBV viral load in the effect of HCV infection on hepatocellular carcinoma (HCC) in HBV seropositive patients (n = 2,878; HCV carrier n = 123). Our results indicated that early HBV viral load had a negative mediating role in HCV-induced HCC. Additionally, early exposure to a low HBV viral load affected HCC through a lag effect on HCC incidence [OR = 0.873, 95% CI = (0.853, 0.893)], and the effect of early exposure to a low HBV viral load on HCC incidence was slightly larger than that of a persistently low viral load on HCC incidence [OR = 0.918, 95% CI = (0.896, 0.941)].
Collapse
Affiliation(s)
- An-Shun Tai
- Institute of Statistics, National Yang Ming Chiao Tung University, Hsinchu, Taiwan
| | - Yen-Tsung Huang
- Institute of Statistical Science, Academia Sinica, Taipei, Taiwan
| | - Hwai-I Yang
- Genomics Research Center, Academia Sinica, Taipei, Taiwan
| | - Lauren V Lan
- Institute of Statistics, National Yang Ming Chiao Tung University, Hsinchu, Taiwan.,Department of Biostatistics, Johns Hopkins University, Baltimore, MD, United States
| | - Sheng-Hsuan Lin
- Institute of Statistics, National Yang Ming Chiao Tung University, Hsinchu, Taiwan
| |
Collapse
|
18
|
Rudolph JE, Cartus A, Bodnar LM, Schisterman EF, Naimi AI. The Role of the Natural Course in Causal Analysis. Am J Epidemiol 2022; 191:341-348. [PMID: 34643230 DOI: 10.1093/aje/kwab248] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2021] [Accepted: 10/05/2021] [Indexed: 11/13/2022] Open
Abstract
The average causal effect compares counterfactual outcomes if everyone had been exposed versus if everyone had been unexposed, which can be an unrealistic contrast. Alternatively, we can target effects that compare counterfactual outcomes against the factual outcomes observed in the sample (i.e., we can compare against the natural course). Here, we demonstrate how the natural course can be estimated and used in causal analyses for model validation and effect estimation. Our example is an analysis assessing the impact of taking aspirin on pregnancy, 26 weeks after randomization, in the Effects of Aspirin in Gestation and Reproduction trial (United States, 2006-2012). To validate our models, we estimated the natural course using g-computation and then compared that against the observed incidence of pregnancy. We observed good agreement between the observed and model-based natural courses. We then estimated an effect that compared the natural course against the scenario in which participants assigned to aspirin always complied. If participants had always complied, there would have been 5.0 (95% confidence interval: 2.2, 7.8) more pregnancies per 100 women than was observed. It is good practice to estimate the natural course for model validation when using parametric models, but whether one should estimate a natural course contrast depends on the underlying research questions.
Collapse
|
19
|
Sudharsanan N, Bijlsma MJ. Educational note: causal decomposition of population health differences using Monte Carlo integration and the g-formula. Int J Epidemiol 2022; 50:2098-2107. [PMID: 34999885 PMCID: PMC8743135 DOI: 10.1093/ije/dyab090] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 04/19/2021] [Indexed: 11/14/2022] Open
Abstract
One key objective of the population health sciences is to understand why one social group has different levels of health and well-being compared with another. Whereas several methods have been developed in economics, sociology, demography, and epidemiology to answer these types of questions, a recent method introduced by Jackson and VanderWeele (2018) provided an update to decompositions by anchoring them within causal inference theory. In this paper, we demonstrate how to implement the causal decomposition using Monte Carlo integration and the parametric g-formula. Causal decomposition can help to identify the sources of differences across populations and provide researchers with a way to move beyond estimating inequalities to explaining them and determining what can be done to reduce health disparities. Our implementation approach can easily and flexibly be applied for different types of outcome and explanatory variables without having to derive decomposition equations. We describe the concepts of the approach and the practical steps and considerations needed to implement it. We then walk through a worked example in which we investigate the contribution of smoking to sex differences in mortality in South Korea. For this example, we provide both pseudocode and R code using our package, cfdecomp. Ultimately, we outline how to implement a very general decomposition algorithm that is grounded in counterfactual theory but still easy to apply to a wide range of situations.
Collapse
Affiliation(s)
| | - Maarten J Bijlsma
- Laboratory of Population Health, Max Planck Institute for Demographic Research, Germany
- Groningen Research Institute of Pharmacy, Unit Pharmacotherapy, -Epidemiology & -Economics (PTEE), University of Groningen, the Netherlands
| |
Collapse
|
20
|
Mooney SJ, Shev AB, Keyes KM, Tracy M, Cerdá M. G-Computation and Agent-Based Modeling for Social Epidemiology: Can Population Interventions Prevent Posttraumatic Stress Disorder? Am J Epidemiol 2022; 191:188-197. [PMID: 34409437 PMCID: PMC8897987 DOI: 10.1093/aje/kwab219] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2020] [Revised: 08/02/2021] [Accepted: 08/11/2021] [Indexed: 11/13/2022] Open
Abstract
Agent-based modeling and g-computation can both be used to estimate impacts of intervening on complex systems. We explored each modeling approach within an applied example: interventions to reduce posttraumatic stress disorder (PTSD). We used data from a cohort of 2,282 adults representative of the adult population of the New York City metropolitan area from 2002-2006, of whom 16.3% developed PTSD over their lifetimes. We built 4 models: g-computation, an agent-based model (ABM) with no between-agent interactions, an ABM with violent-interaction dynamics, and an ABM with neighborhood dynamics. Three interventions were tested: 1) reducing violent victimization by 37.2% (real-world reduction); 2) reducing violent victimization by100%; and 3) supplementing the income of 20% of lower-income participants. The g-computation model estimated population-level PTSD risk reductions of 0.12% (95% confidence interval (CI): -0.16, 0.29), 0.28% (95% CI: -0.30, 0.70), and 1.55% (95% CI: 0.40, 2.12), respectively. The ABM with no interactions replicated the findings from g-computation. Introduction of interaction dynamics modestly decreased estimated intervention effects (income-supplement risk reduction dropped to 1.47%), whereas introduction of neighborhood dynamics modestly increased effectiveness (income-supplement risk reduction increased to 1.58%). Compared with g-computation, agent-based modeling permitted deeper exploration of complex systems dynamics at the cost of further assumptions.
Collapse
Affiliation(s)
- Stephen J Mooney
- Correspondence to Dr. Stephen Mooney, 1959 NE Pacific Street, Health Sciences Building F-262, Box 357236, Seattle, WA 98195 (e-mail: )
| | | | | | | | | |
Collapse
|
21
|
Cartus AR, Jarlenski MP, Himes KP, James AE, Naimi AI, Bodnar LM. Adverse Cardiovascular Events Following Severe Maternal Morbidity. Am J Epidemiol 2022; 191:126-136. [PMID: 34343230 DOI: 10.1093/aje/kwab208] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2021] [Revised: 07/16/2021] [Accepted: 07/16/2021] [Indexed: 11/13/2022] Open
Abstract
Severe maternal morbidity (SMM) affects 50,000 women annually in the United States, but its consequences are not well understood. We aimed to estimate the association between SMM and risk of adverse cardiovascular events during the 2 years postpartum. We analyzed 137,140 deliveries covered by the Pennsylvania Medicaid program (2016-2018), weighted with inverse probability of censoring weights to account for nonrandom loss to follow-up. SMM was defined as any diagnosis on the Centers for Disease Control and Prevention list of SMM diagnoses and procedures and/or intensive care unit admission occurring at any point from conception through 42 days postdelivery. Outcomes included heart failure, ischemic heart disease, and stroke/transient ischemic attack up to 2 years postpartum. We used marginal standardization to estimate average treatment effects. We found that SMM was associated with increased risk of each adverse cardiovascular event across the follow-up period. Per 1,000 deliveries, relative to no SMM, SMM was associated with 12.1 (95% confidence interval (CI): 6.2, 18.0) excess cases of heart failure, 6.4 (95% CI: 1.7, 11.2) excess cases of ischemic heart disease, and 8.2 (95% CI: 3.2, 13.1) excess cases of stroke/transient ischemic attack at 26 months of follow-up. These results suggest that SMM identifies a group of women who are at high risk of adverse cardiovascular events after delivery. Women who survive SMM may benefit from more comprehensive postpartum care linked to well-woman care.
Collapse
|
22
|
Sheth AN, Adimora AA, Golub ET, Kassaye SG, Rana A, Westreich D, Cyriaque JW, Parish C, Konkle-Parker D, Jones DL, Kempf MC, Ofotokun I, Kanthula RM, Donohue J, Raccamarich P, Tisdale T, Ramirez C, Warren-Jeanpiere L, Tien PC, Alcaide ML. Study of Treatment and Reproductive Outcomes Among Reproductive-Age Women With HIV Infection in the Southern United States: Protocol for a Longitudinal Cohort Study. JMIR Res Protoc 2021; 10:e30398. [PMID: 34932006 PMCID: PMC8726043 DOI: 10.2196/30398] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2021] [Revised: 06/30/2021] [Accepted: 07/16/2021] [Indexed: 01/26/2023] Open
Abstract
BACKGROUND Nearly a quarter of the 1.1 million individuals with HIV in the United States are women. Racial and ethnic minority women in the Southern United States are disproportionately impacted. Reproductive-age women with HIV are prone to poor HIV outcomes but remain underrepresented in HIV research. We will answer contemporary questions related to the health outcomes in this population by enrolling a prospective cohort of reproductive-age women with and without HIV in the Southern United States. OBJECTIVE The Study of Treatment and Reproductive Outcomes (STAR) will enroll and retain 2000 reproductive-age women with and without HIV. The STAR will leverage the infrastructure of the US-based Multicenter AIDS Cohort Study (MACS)/Women's Interagency HIV Study (WIHS) Combined Cohort Study, comprising the WIHS (a cohort of women with and at risk for HIV, which began in 1993), and the MACS (a cohort of gay and bisexual men with and at risk for HIV, which began in 1984). Although the advancing age of the participants enrolled in the MACS/WIHS Combined Cohort Study provides an opportunity to address the questions related to HIV and aging, the research questions pertinent to the reproductive years must also be addressed. The STAR will conduct high-priority scientific research in key areas with the overall aim of addressing the unique needs of reproductive-age women with HIV. METHODS The STAR is a prospective, observational cohort study that will be conducted at 6 sites in the United States-Atlanta, Georgia; Birmingham, Alabama; Jackson, Mississippi; Chapel Hill, North Carolina; Miami, Florida; and Washington, District of Columbia. Visits will occur semiannually for 2 years, with additional visits for up to 5 years. At each visit, the participating women will complete a structured interview for collecting key demographic, psychosocial, and clinical variables, and undergo biospecimen collection for laboratory testing and repositing (blood, urine, hair, vaginal, anal, and oral specimens). Pregnant women and infants will undergo additional study assessments. The initial scientific focus of the STAR is to understand the roles of key social determinants of health, depression, reproductive health, and oral health on HIV and pregnancy outcomes across the reproductive life span. RESULTS Enrollment in the STAR commenced in February 2021 and is ongoing. CONCLUSIONS Through in-depth, longitudinal data and biospecimen collection, the newly initiated STAR cohort will create a platform to answer scientific questions regarding reproductive-age women with and without HIV. STAR will be uniquely positioned to enable investigators to conduct high-impact research relevant to this population. Building on the legacy of the MACS and WIHS cohorts, the STAR is designed to foster multidisciplinary collaborations to galvanize scientific discoveries to improve the health of reproductive-age women with HIV and ameliorate the effects of the HIV epidemic in this population in the United States.
Collapse
Affiliation(s)
- Anandi N Sheth
- Division of Infectious Diseases, Department of Medicine, Emory University School of Medicine, Atlanta, GA, United States
- Infectious Diseases Program, Grady Health System, Atlanta, GA, United States
| | - Adaora A Adimora
- Division of Infectious Diseases, Department of Medicine, University of North Carolina-Chapel Hill, Chapel Hill, NC, United States
- Department of Epidemiology, University of North Carolina-Chapel Hill, Chapel Hill, NC, United States
| | - Elizabeth Topper Golub
- Department of Epidemiology, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD, United States
| | - Seble G Kassaye
- Department of Medicine, Georgetown University, Washington, DC, United States
| | - Aadia Rana
- Division of Infectious Diseases, Department of Medicine, University of Alabama-Birmingham Heersink School of Medicine, Birmingham, AL, United States
| | - Daniel Westreich
- Department of Epidemiology, University of North Carolina-Chapel Hill, Chapel Hill, NC, United States
| | - Jennifer Webster Cyriaque
- Division of Oral and Craniofacial Sciences, Department of Microbiology and Immunology, University of North Carolina - Chapel Hill, Chapel Hill, NC, United States
| | - Carrigan Parish
- Department of Sociomedical Sciences, Mailman School of Public Health, Columbia University, NY, United States
| | - Deborah Konkle-Parker
- Schools of Nursing, Medicine, and Population Health Sciences, University of Mississippi Medical Center, Jackson, MS, United States
| | - Deborah L Jones
- Department of Psychiatry and Behavioral Sciences, University of Miami Miller School of Medicine, Miami, FL, United States
| | - Mirjam-Colette Kempf
- Division of Infectious Diseases, Department of Medicine, University of Alabama-Birmingham Heersink School of Medicine, Birmingham, AL, United States
- Departments of Epidemiology and Health Behavior, University of Alabama-Birmingham Ryals School of Public Health, Birmingham, AL, United States
- Department of Nursing Family, Community & Health Systems, University of Alabama-Birmingham School of Nursing, Birmingham, AL, United States
| | - Igho Ofotokun
- Division of Infectious Diseases, Department of Medicine, Emory University School of Medicine, Atlanta, GA, United States
- Infectious Diseases Program, Grady Health System, Atlanta, GA, United States
| | - Ruth M Kanthula
- Department of Pediatrics, Georgetown University, Washington, DC, United States
| | - Jessica Donohue
- Department of Epidemiology, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD, United States
| | - Patricia Raccamarich
- Division of Infectious Diseases, Department of Medicine, University of Miami Miller School of Medicine, Miami, FL, United States
| | - Tina Tisdale
- Division of Infectious Diseases, Department of Medicine, Emory University School of Medicine, Atlanta, GA, United States
| | - Catalina Ramirez
- Division of Infectious Diseases, Department of Medicine, University of North Carolina-Chapel Hill, Chapel Hill, NC, United States
| | | | - Phyllis C Tien
- Department of Medicine, University of California San Francisco, San Francisco, CA, United States
- Department of Veteran Affairs Medical Center, San Francisco, CA, United States
| | - Maria L Alcaide
- Division of Infectious Diseases, Department of Medicine, University of Miami Miller School of Medicine, Miami, FL, United States
- Department of Obstetrics and Gynecology, University of Miami Miller School of Medicine, Miami, FL, United States
| |
Collapse
|
23
|
Cerdá M, Jalali MS, Hamilton AD, DiGennaro C, Hyder A, Santaella-Tenorio J, Kaur N, Wang C, Keyes KM. A Systematic Review of Simulation Models to Track and Address the Opioid Crisis. Epidemiol Rev 2021; 43:147-165. [PMID: 34791110 DOI: 10.1093/epirev/mxab013] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2020] [Revised: 10/20/2021] [Accepted: 11/04/2021] [Indexed: 01/04/2023] Open
Abstract
The opioid overdose crisis is driven by an intersecting set of social, structural, and economic forces. Simulation models offer a tool to help us understand and address this complex, dynamic, and nonlinear social phenomenon. We conducted a systematic review of the literature on simulation models of opioid use and overdose up to September 2019. We extracted modeling types, target populations, interventions, and findings. Further, we created a database of model parameters used for model calibration, and evaluated study transparency and reproducibility. Of the 1,398 articles screened, we identified 88 eligible articles. The most frequent types of models were compartmental (36%), Markov (20%), system dynamics (16%), and Agent-Based models (16%). Over a third evaluated intervention cost-effectiveness (40%), and another third (39%) focused on treatment and harm reduction services for people with opioid use disorder (OUD). More than half (61%) discussed calibrating their models to empirical data, and 31% discussed validation approaches used in their modeling process. From the 63 studies that provided model parameters, we extracted the data sources on opioid use, OUD, OUD treatment, cessation/relapse, emergency medical services, and mortality parameters. This database offers a tool that future modelers can use to identify potential model inputs and evaluate comparability of their models to prior work. Future applications of simulation models to this field should actively tackle key methodological challenges, including the potential for bias in the choice of parameter inputs, investment in model calibration and validation, and transparency in the assumptions and mechanics of simulation models to facilitate reproducibility.
Collapse
Affiliation(s)
- Magdalena Cerdá
- Center for Opioid Epidemiology and Policy, Division of Epidemiology, Department of Population Health, New York University School of Medicine, New York, New York
| | | | - Ava D Hamilton
- Department of Epidemiology, Mailman School of Public Health, Columbia University, New York, New York
| | | | - Ayaz Hyder
- Division of Environmental Health Sciences, College of Public Health, and Translational Data Analytics Institute, The Ohio State University, Columbus, Ohio
| | - Julian Santaella-Tenorio
- Center for Opioid Epidemiology and Policy, Division of Epidemiology, Department of Population Health, New York University School of Medicine, New York, New York
| | - Navdep Kaur
- Department of Epidemiology, Mailman School of Public Health, Columbia University, New York, New York
| | - Christina Wang
- Center for Opioid Epidemiology and Policy, Division of Epidemiology, Department of Population Health, New York University School of Medicine, New York, New York
| | - Katherine M Keyes
- Department of Epidemiology, Mailman School of Public Health, Columbia University, New York, New York
| |
Collapse
|
24
|
Stoner MCD, Edwards JK, Westreich D, Kilburn K, Ahern J, Lippman SA, Gómez-Olivé FX, Kahn K, Pettifor A. Modeling Cash Plus Other Psychosocial and Structural Interventions to Prevent HIV Among Adolescent Girls and Young Women in South Africa (HPTN 068). AIDS Behav 2021; 25:133-143. [PMID: 33475905 PMCID: PMC8292433 DOI: 10.1007/s10461-021-03158-3] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 01/08/2021] [Indexed: 11/25/2022]
Abstract
Poverty alleviation programs can reduce HIV incidence but may have greater impacts when combined with other psychosocial interventions. We modeled the change in HIV incidence among South African adolescent girls and young women (AGYW) associated with combining a cash transfer (the South African Child Support Grant (CSG)) with other structural and behavioral interventions. We modeled observational data from the HPTN 068 study where 2328 HIV negative AGYW (13–20 years) were followed for 4 years. In a Monte Carlo simulation based on this cohort (N = 10,000), CSG receipt was not independently associated with HIV incidence. Providing the CSG combined with increasing caregiver care and reducing adolescent depression had the largest reduction in HIV incidence with the fewest number of combined interventions (RD − 3.0%; (95% CI − 5.1%, − 0.9%). Combining a monthly grant with interventions to increase caregiver care and reduce adolescent depression could substantially reduce HIV incidence above the provision of cash alone.
Collapse
Affiliation(s)
- Marie C D Stoner
- Women's Global Health Imperative, RTI International, Berkeley, CA, USA.
- Carolina Population Center, University of North Carolina, Chapel Hill, NC, USA.
| | - Jessie K Edwards
- Department of Epidemiology, University of North Carolina, Chapel Hill, NY, USA
| | - Daniel Westreich
- Department of Epidemiology, University of North Carolina, Chapel Hill, NY, USA
| | - Kelly Kilburn
- Duke Global Health Institute, Duke University, Durham, NC, USA
| | - Jennifer Ahern
- Division of Epidemiology & Biostatistics, School of Public Health, University of California, Berkeley, USA
| | - Sheri A Lippman
- MRC/Wits Rural Public Health and Health Transitions Research Unit (Agincourt), School of Public Health, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa
- Department of Medicine, University of California, San Francisco, USA
| | - F Xavier Gómez-Olivé
- MRC/Wits Rural Public Health and Health Transitions Research Unit (Agincourt), School of Public Health, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa
- INDEPTH Network, Accra, Ghana
| | - Kathleen Kahn
- MRC/Wits Rural Public Health and Health Transitions Research Unit (Agincourt), School of Public Health, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa
- INDEPTH Network, Accra, Ghana
- Epidemiology and Global Health Unit, Department of Public Health and Clinical Medicine, Umeå University, Umeå, Sweden
| | - Audrey Pettifor
- Carolina Population Center, University of North Carolina, Chapel Hill, NC, USA
- Department of Epidemiology, University of North Carolina, Chapel Hill, NY, USA
- MRC/Wits Rural Public Health and Health Transitions Research Unit (Agincourt), School of Public Health, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa
| |
Collapse
|
25
|
Stoner MCD, Westreich D, Ahern J, Edwards J, Gómez-Olivé FX, Tollman SM, Lippman S, Kahn K, Pettifor A. Modeling Combination Interventions to Prevent Human Immunodeficiency Virus in Adolescent Girls and Young Women in South Africa (HIV Prevention Trials Network 068). Clin Infect Dis 2021; 73:e1911-e1918. [PMID: 33325509 DOI: 10.1093/cid/ciaa1598] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2020] [Indexed: 11/14/2022] Open
Abstract
BACKGROUND Combination interventions may be an effective way to prevent human immunodeficiency virus (HIV) in adolescent girls and young women. However, current studies are not designed to understand which specific interventions and combinations will be most effective. We estimate the possible impacts of interventions on a combination of factors associated with HIV. METHODS We used the g-formula to model interventions on combinations of HIV risk factors to identify those that would prevent the most incident HIV infections, including low school attendance, intimate partner violence, depression, transactional sex, and age-disparate partnerships. We used data from the HIV Prevention Trials Network (HPTN) 068 study in rural South Africa from 2011 to 2017. We estimated HIV incidence under a potential intervention that reduced each risk factor and compared this to HIV incidence under the current distribution of these risk factors. RESULTS Although many factors had strong associations with HIV, potential intervention estimates did not always suggest large reductions in HIV incidence because the prevalence of risk factors was low. When modeling combination effects, an intervention to increase schooling, decrease depression, and decease transactional sex showed the largest reduction in incident infection (risk difference, -1.4%; 95% confidence interval [CI], -2.7% to -.2%), but an intervention on only transactional sex and depression still reduced HIV incidence by -1.3% (95% CI, -2.6% to -.2%). CONCLUSIONS To achieve the largest reductions in HIV, both prevalence of the risk factor and strength of association with HIV must be considered. Additionally, intervening on more risk factors may not necessarily result in larger reductions in HIV incidence.
Collapse
Affiliation(s)
- Marie C D Stoner
- Carolina Population Center, University of North Carolina, Chapel Hill, North Carolina, USA.,Women's Global Health Imperative, RTI International, San Francisco, California, USA
| | - Daniel Westreich
- Department of Epidemiology, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, USA
| | - Jennifer Ahern
- School of Public Health, Division of Epidemiology and Biostatistics, University of California, Berkeley, Berkeley, California, USA
| | - Jessie Edwards
- Department of Epidemiology, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, USA
| | - F Xavier Gómez-Olivé
- Medical Research Council/Wits Rural Public Health and Health Transitions Research Unit (Agincourt), School of Public Health, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa
| | - Stephen M Tollman
- Medical Research Council/Wits Rural Public Health and Health Transitions Research Unit (Agincourt), School of Public Health, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa.,INDEPTH Network, Accra, Ghana
| | - Sheri Lippman
- Medical Research Council/Wits Rural Public Health and Health Transitions Research Unit (Agincourt), School of Public Health, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa.,Department of Medicine, University of California, San Francisco, San Francisco, California, USA
| | - Kathleen Kahn
- Medical Research Council/Wits Rural Public Health and Health Transitions Research Unit (Agincourt), School of Public Health, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa.,INDEPTH Network, Accra, Ghana.,Epidemiology and Global Health Unit, Department of Public Health and Clinical Medicine, Umeå University, Umeå, Sweden
| | - Audrey Pettifor
- Carolina Population Center, University of North Carolina, Chapel Hill, North Carolina, USA.,Department of Epidemiology, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, USA.,Medical Research Council/Wits Rural Public Health and Health Transitions Research Unit (Agincourt), School of Public Health, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa
| |
Collapse
|
26
|
Filshtein TJ, Li X, Zimmerman SC, Ackley SF, Glymour MM, Power MC. Proof of Concept Example for Use of Simulation to Allow Data Pooling Despite Privacy Restrictions. Epidemiology 2021; 32:638-647. [PMID: 34183527 PMCID: PMC8338788 DOI: 10.1097/ede.0000000000001373] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Abstract
BACKGROUND Integrating results from multiple samples is often desirable, but privacy restrictions may preclude full data pooling, and most datasets do not include fully harmonized variable sets. We propose a simulation-based method leveraging partial information across datasets to guide creation of synthetic data based on explicit assumptions about the underlying causal structure that permits pooled analyses that adjust for all desired confounders in the context of privacy restrictions. METHODS This proof-of-concept project uses data from the Health and Retirement Study (HRS) and Atherosclerosis Risk in Communities (ARIC) study. We specified an estimand of interest and a directed acyclic graph (DAG) summarizing the presumed causal structure for the effect of glycated hemoglobin (HbA1c) on cognitive change. We derived publicly reportable statistics to describe the joint distribution of each variable in our DAG. These summary estimates were used as data-generating rules to create synthetic datasets. After pooling, we imputed missing covariates in the synthetic datasets and used the synthetic data to estimate the pooled effect of HbA1c on cognitive change, adjusting for all desired covariates. RESULTS Distributions of covariates and model coefficients and associated standard errors for our model estimating the effect of HbA1c on cognitive change were similar across cohort-specific original and preimputation synthetic data. The estimate from the pooled synthetic incorporates control for confounders measured in either original dataset. DISCUSSION Our approach has advantages over meta-analysis or individual-level pooling/data harmonization when privacy concerns preclude data sharing and key confounders are not uniformly measured across datasets.
Collapse
Affiliation(s)
- Teresa J. Filshtein
- Department of Epidemiology and Biostatistics, University of California San Francisco School of Public Health
| | - Xiang Li
- Department of Epidemiology, Milken Institute School of Public Health, George Washington University
| | - Scott C. Zimmerman
- Department of Epidemiology and Biostatistics, University of California San Francisco School of Public Health
| | - Sarah F. Ackley
- Department of Epidemiology and Biostatistics, University of California San Francisco School of Public Health
| | - M. Maria Glymour
- Department of Epidemiology and Biostatistics, University of California San Francisco School of Public Health
| | - Melinda C. Power
- Department of Epidemiology, Milken Institute School of Public Health, George Washington University
| |
Collapse
|
27
|
McManus KA, Christensen B, Nagraj VP, Furl R, Yerkes L, Swindells S, Weissman S, Rhodes A, Targonski P, Rogawski McQuade E, Dillingham R. Evidence From a Multistate Cohort: Enrollment in Affordable Care Act Qualified Health Plans' Association With Viral Suppression. Clin Infect Dis 2021; 71:2572-2580. [PMID: 31734691 PMCID: PMC7744983 DOI: 10.1093/cid/ciz1123] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2019] [Accepted: 11/12/2019] [Indexed: 11/23/2022] Open
Abstract
Background Healthcare delivery changes associated with viral suppression (VS) could contribute to the United States’ “Ending the HIV Epidemic” (EtHE) initiative. This study aims to determine whether Qualified Health Plans (QHPs) purchased by AIDS Drug Assistance Programs (ADAPs) are associated with VS for low-income people living with HIV (PLWH) across 3 states. Methods A multistate cohort of ADAP clients eligible for ADAP-funded QHPs were studied (2014–2015). A log-binomial model was used to estimate the association of demographics and healthcare delivery factors with QHP enrollment prevalence and 1-year risk of VS. A number needed to treat/enroll (NNT) for 1 additional person to achieve viral suppression was calculated. Results Of the cohort (n = 7776), 52% enrolled in QHPs. QHP enrollment in 2015 was associated with QHP coverage in 2014 (adjusted PR [aPR], 3.28; 95% confidence intervals [CIs], 3.06–3.53) and engagement in care in 2014 (aPR, 1.16; 1.04–1.28). PLWH who were engaged in care (n = 4597) and had QHPs had a higher VS rate than those who received medications from Direct ADAP (86.0% vs 80.2%). QHPs’ NNT for an additional person to achieve VS is 20 (14.1–34.5). Starting undetectable (adjusted risk ratio [aRR], 1.39; 1.28–1.52) and enrolling in QHPs in 2015 (aRR, 1.06; 0.99–1.14) was associated with VS. Conclusions Once enrolled in ADAP-funded QHPs, ADAP clients stay enrolled. Enrollment is associated with VS across states/demographic groups. ADAPs, especially in the South and in Medicaid nonexpansion states, should consider investing in QHPs because increased enrollment could improve VS rates. This evidence-based intervention could be part of EtHE.
Collapse
Affiliation(s)
- Kathleen A McManus
- Division of Infectious Diseases and International Health, University of Virginia, Charlottesville, Virginia, USA.,Center for Health Policy, University of Virginia, Charlottesville, Virginia, USA
| | | | - V Peter Nagraj
- School of Medicine Research Computing, University of Virginia, Charlottesville, Virginia, USA
| | - Renae Furl
- Division of Infectious Diseases, University of Nebraska Medical Center, Omaha, Nebraska, USA
| | - Lauren Yerkes
- Virginia Department of Health, Richmond, Virginia, USA
| | - Susan Swindells
- Division of Infectious Diseases, University of Nebraska Medical Center, Omaha, Nebraska, USA
| | - Sharon Weissman
- Division of Infectious Diseases, University of South Carolina, Columbia, South Carolina, USA
| | - Anne Rhodes
- Virginia Department of Health, Richmond, Virginia, USA
| | - Paul Targonski
- Public Health Sciences, University of Virginia, Charlottesville, Virginia, USA.,Department of Medicine, University of Virginia, Charlottesville, Virginia, USA
| | - Elizabeth Rogawski McQuade
- Division of Infectious Diseases and International Health, University of Virginia, Charlottesville, Virginia, USA.,Public Health Sciences, University of Virginia, Charlottesville, Virginia, USA
| | - Rebecca Dillingham
- Division of Infectious Diseases and International Health, University of Virginia, Charlottesville, Virginia, USA
| |
Collapse
|
28
|
Karaboyas A, Morgenstern H, Fleischer NL, Schaubel DE, Robinson BM. Replicating Randomized Trial Results with Observational Data Using the Parametric g-Formula: An Application to Intravenous Iron Treatment in Hemodialysis Patients. Clin Epidemiol 2020; 12:1249-1260. [PMID: 33204166 PMCID: PMC7667704 DOI: 10.2147/clep.s283321] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2020] [Accepted: 10/27/2020] [Indexed: 11/23/2022] Open
Abstract
BACKGROUND Reproducibility of clinical and epidemiologic research is important to generalize findings and has increasingly been scrutinized. A recently published randomized trial, PIVOTAL, evaluated high vs low intravenous iron dosing strategies to manage anemia in hemodialysis patients in the UK. Our objective was to assess the reproducibility of the PIVOTAL trial findings using data from a well-established cohort study, the Dialysis Outcomes and Practice Patterns Study (DOPPS). METHODS To overcome the absence of randomization in the DOPPS, we applied the parametric g-formula, an extension of standardization to longitudinal data. We estimated the effect of a proactive high-dose vs reactive low-dose iron supplementation strategy on all-cause mortality (primary outcome), hemoglobin, two measures of iron concentration (ferritin and TSAT), and erythropoiesis-stimulating agent dose over 12 months of follow-up in 6325 DOPPS patients. RESULTS Comparing high- vs low-iron dose strategies, the 1-year mortality risk difference was 0.020 (95% CI: 0.008, 0.031) and risk ratio was 1.20 (95% CI: 1.07, 1.33), compared with null 1-year findings in the PIVOTAL trial. Differences in secondary outcomes were directionally consistent but of lesser magnitude than in the PIVOTAL trial. CONCLUSION Our findings are somewhat consistent with the recent PIVOTAL trial, with discrepancies potentially attributable to model misspecification and differences between the two study populations. In addition to the importance of our results to nephrologists and hence hemodialysis patients, our analysis illustrates the utility of the parametric g-formula for generalizing results and comparing complex and dynamic treatment strategies using observational data.
Collapse
Affiliation(s)
- Angelo Karaboyas
- Arbor Research Collaborative for Health, Ann Arbor, MI, USA
- University of Michigan, Department of Epidemiology, Ann Arbor, MI, USA
| | - Hal Morgenstern
- University of Michigan, Departments of Epidemiology and Environmental Health Sciences, School of Public Health, and Department of Urology, Medical School, Ann Arbor, MI, USA
| | - Nancy L Fleischer
- University of Michigan, Department of Epidemiology, Ann Arbor, MI, USA
| | - Douglas E Schaubel
- University of Michigan, Department of Biostatistics, Ann Arbor, MI, USA
- University of Pennsylvania, Department of Biostatistics, Epidemiology and Informatics, Philadelphia, PA, USA
| | - Bruce M Robinson
- Arbor Research Collaborative for Health, Ann Arbor, MI, USA
- University of Michigan, Department of Internal Medicine, Ann Arbor, MI, USA
| |
Collapse
|
29
|
Tchetgen Tchetgen EJ, Fulcher IR, Shpitser I. Auto-G-Computation of Causal Effects on a Network. J Am Stat Assoc 2020; 116:833-844. [PMID: 34366505 PMCID: PMC8345318 DOI: 10.1080/01621459.2020.1811098] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2018] [Revised: 02/05/2020] [Accepted: 07/31/2020] [Indexed: 10/23/2022]
Abstract
Methods for inferring average causal effects have traditionally relied on two key assumptions: (i) the intervention received by one unit cannot causally influence the outcome of another; and (ii) units can be organized into nonoverlapping groups such that outcomes of units in separate groups are independent. In this article, we develop new statistical methods for causal inference based on a single realization of a network of connected units for which neither assumption (i) nor (ii) holds. The proposed approach allows both for arbitrary forms of interference, whereby the outcome of a unit may depend on interventions received by other units with whom a network path through connected units exists; and long range dependence, whereby outcomes for any two units likewise connected by a path in the network may be dependent. Under network versions of consistency and no unobserved confounding, inference is made tractable by an assumption that the networks outcome, treatment and covariate vectors are a single realization of a certain chain graph model. This assumption allows inferences about various network causal effects via the auto-g-computation algorithm, a network generalization of Robins' well-known g-computation algorithm previously described for causal inference under assumptions (i) and (ii). Supplementary materials for this article are available online.
Collapse
Affiliation(s)
| | - Isabel R. Fulcher
- Department of Global Health and Social Medicine, Harvard Medical School, Boston, MA
| | - Ilya Shpitser
- Department of Computer Science, Johns Hopkins Whiting School of Engineering, Baltimore, MD
| |
Collapse
|
30
|
Breger TL, Edwards JK, Cole SR, Westreich D, Pence BW, Adimora AA. Two-stage g-computation: Evaluating Treatment and Intervention Impacts in Observational Cohorts When Exposure Information Is Partly Missing. Epidemiology 2020; 31:695-703. [PMID: 32657953 PMCID: PMC8725064 DOI: 10.1097/ede.0000000000001233] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
Abstract
Illustrations of the g-computation algorithm to evaluate population average treatment and intervention effects have been predominantly implemented in settings with complete exposure information. Thus, worked examples of approaches to handle missing data in this causal framework are needed to facilitate wider use of these estimators. We illustrate two-stage g-computation estimators that leverage partially observed information on the full study sample and complete exposure information on a subset to estimate causal effects. In a hypothetical cohort of 1,623 human immunodeficiency virus (HIV)-positive women with 30% complete opioid prescription information, we illustrate a two-stage extrapolation g-computation estimator for the average treatment effect of shorter or longer duration opioid prescriptions; we further illustrate two-stage inverse probability weighting and imputation g-computation estimators for the average intervention effect of shortening the duration of prescriptions relative to the status quo. Two-stage g-computation estimators approximated the true risk differences for the population average treatment and intervention effects while g-computation fit to the subset of complete cases was biased. In 10,000 Monte Carlo simulations, two-stage approaches considerably reduced bias and mean squared error and improved the coverage of 95% confidence limits. Although missing data threaten validity and precision, two-stage g-computation designs offer principled approaches to handling missing information.
Collapse
Affiliation(s)
- Tiffany L. Breger
- Department of Epidemiology, University of North Carolina at Chapel Hill, Chapel Hill, NC
| | - Jessie K. Edwards
- Department of Epidemiology, University of North Carolina at Chapel Hill, Chapel Hill, NC
| | - Stephen R. Cole
- Department of Epidemiology, University of North Carolina at Chapel Hill, Chapel Hill, NC
| | - Daniel Westreich
- Department of Epidemiology, University of North Carolina at Chapel Hill, Chapel Hill, NC
| | - Brian W. Pence
- Department of Epidemiology, University of North Carolina at Chapel Hill, Chapel Hill, NC
| | - Adaora A. Adimora
- Department of Epidemiology, University of North Carolina at Chapel Hill, Chapel Hill, NC
- Department of Medicine, University of North Carolina at Chapel Hill, Chapel Hill, NC
| |
Collapse
|
31
|
Spieker AJ, Ko EM, Roy JA, Mitra N. Nested g-computation: A causal approach to analysis of censored medical costs in the presence of time-varying treatment. J R Stat Soc Ser C Appl Stat 2020; 69:1189-1208. [PMID: 34108743 DOI: 10.1111/rssc.12441] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Rising medical costs are an emerging challenge in policy decisions and resource allocation planning. When cumulative medical cost is the outcome, right-censoring induces informative missingness due to heterogeneity in cost accumulation rates across subjects. Inverse-weighting approaches have been developed to address the challenge of informative cost trajectories in mean cost estimation, though these approaches generally ignore post-baseline treatment changes. In post-hysterectomy endometrial cancer patients, data from a linked database of Medicare records and the Surveillance, Epidemiology, and End Results program of the National Cancer Institute reveal substantial within-subject variation in treatment over time. In such a setting, the utility of existing intent-to-treat approaches is generally limited. Estimates of population mean cost under a hypothetical time-varying treatment regime can better assist with resource allocation when planning for a treatment policy change; such estimates must inherently take time-dependent treatment and confounding into account. In this paper, we develop a nested g-computation approach to cost analysis to address this challenge, while accounting for censoring. We develop a procedure to evaluate sensitivity to departures from baseline treatment ignorability. We further conduct a variety of simulations and apply our nested g-computation procedure to two-year costs from endometrial cancer patients.
Collapse
Affiliation(s)
- Andrew J Spieker
- Department of Biostatistics, Vanderbilt University Medical Center, Nashville, TN 37203
| | - Emily M Ko
- Leonard Davis Institute of Health Economics, University of Pennsylvania, Philadelphia, PA 19104.,Department of Obstetrics and Gynecology, University of Pennsylvania Perelman School of Medicine, Philadelphia, PA 19104
| | - Jason A Roy
- Department of Biostatistics and Epidemiology, Rutgers University, Piscataway, NJ 08854
| | - Nandita Mitra
- Department of Biostatistics, Epidemiology, and Informatics, University of Pennsylvania, Philadelphia, PA 19104.,Leonard Davis Institute of Health Economics, University of Pennsylvania, Philadelphia, PA 19104
| |
Collapse
|
32
|
Kawahara T, Shinozaki T, Matsuyama Y. Doubly robust estimator of risk in the presence of censoring dependent on time-varying covariates: application to a primary prevention trial for coronary events with pravastatin. BMC Med Res Methodol 2020; 20:204. [PMID: 32736528 PMCID: PMC7395418 DOI: 10.1186/s12874-020-01087-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2020] [Accepted: 07/23/2020] [Indexed: 11/11/2022] Open
Abstract
Background In the presence of dependent censoring even after stratification of baseline covariates, the Kaplan–Meier estimator provides an inconsistent estimate of risk. To account for dependent censoring, time-varying covariates can be used along with two statistical methods: the inverse probability of censoring weighted (IPCW) Kaplan–Meier estimator and the parametric g-formula estimator. The consistency of the IPCW Kaplan–Meier estimator depends on the correctness of the model specification of censoring hazard, whereas that of the parametric g-formula estimator depends on the correctness of the models for event hazard and time-varying covariates. Methods We combined the IPCW Kaplan–Meier estimator and the parametric g-formula estimator into a doubly robust estimator that can adjust for dependent censoring. The estimator is theoretically more robust to model misspecification than the IPCW Kaplan–Meier estimator and the parametric g-formula estimator. We conducted simulation studies with a time-varying covariate that affected both time-to-event and censoring under correct and incorrect models for censoring, event, and time-varying covariates. We applied our proposed estimator to a large clinical trial data with censoring before the end of follow-up. Results Simulation studies demonstrated that our proposed estimator is doubly robust, namely it is consistent if either the model for the IPCW Kaplan–Meier estimator or the models for the parametric g-formula estimator, but not necessarily both, is correctly specified. Simulation studies and data application demonstrated that our estimator can be more efficient than the IPCW Kaplan–Meier estimator. Conclusions The proposed estimator is useful for estimation of risk if censoring is affected by time-varying risk factors.
Collapse
Affiliation(s)
- Takuya Kawahara
- Clinical Research Promotion Center, The University of Tokyo Hospital, 7-3-1, Hongo, Bunkyo-ku, Tokyo, 113-8655, Japan.
| | - Tomohiro Shinozaki
- Department of Information and Computer Technology, Graduate School of Engineering, Tokyo University of Science, Tokyo, Japan
| | - Yutaka Matsuyama
- Department of Biostatistics, School of Public Health, The University of Tokyo, Tokyo, Japan
| |
Collapse
|
33
|
Keil AP, Buckley JP, O'Brien KM, Ferguson KK, Zhao S, White AJ. A Quantile-Based g-Computation Approach to Addressing the Effects of Exposure Mixtures. ENVIRONMENTAL HEALTH PERSPECTIVES 2020; 128:47004. [PMID: 32255670 PMCID: PMC7228100 DOI: 10.1289/ehp5838] [Citation(s) in RCA: 641] [Impact Index Per Article: 160.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/17/2023]
Abstract
BACKGROUND Exposure mixtures frequently occur in data across many domains, particularly in the fields of environmental and nutritional epidemiology. Various strategies have arisen to answer questions about exposure mixtures, including methods such as weighted quantile sum (WQS) regression that estimate a joint effect of the mixture components. OBJECTIVES We demonstrate a new approach to estimating the joint effects of a mixture: quantile g-computation. This approach combines the inferential simplicity of WQS regression with the flexibility of g-computation, a method of causal effect estimation. We use simulations to examine whether quantile g-computation and WQS regression can accurately and precisely estimate the effects of mixtures in a variety of common scenarios. METHODS We examine the bias, confidence interval (CI) coverage, and bias-variance tradeoff of quantile g-computation and WQS regression and how these quantities are impacted by the presence of noncausal exposures, exposure correlation, unmeasured confounding, and nonlinearity of exposure effects. RESULTS Quantile g-computation, unlike WQS regression, allows inference on mixture effects that is unbiased with appropriate CI coverage at sample sizes typically encountered in epidemiologic studies and when the assumptions of WQS regression are not met. Further, WQS regression can magnify bias from unmeasured confounding that might occur if important components of the mixture are omitted from the analysis. DISCUSSION Unlike inferential approaches that examine the effects of individual exposures while holding other exposures constant, methods like quantile g-computation that can estimate the effect of a mixture are essential for understanding the effects of potential public health actions that act on exposure sources. Our approach may serve to help bridge gaps between epidemiologic analysis and interventions such as regulations on industrial emissions or mining processes, dietary changes, or consumer behavioral changes that act on multiple exposures simultaneously. https://doi.org/10.1289/EHP5838.
Collapse
Affiliation(s)
- Alexander P Keil
- Department of Epidemiology, University of North Carolina, Chapel Hill, North Carolina, USA
- Epidemiology Branch, National Institute of Environmental Health Sciences, National Institutes of Health, Department of Health and Human Services, Research Triangle Park, North Carolina, USA
| | - Jessie P Buckley
- Department of Environmental Health and Engineering, Johns Hopkins University, Baltimore, Maryland, USA
- Department of Epidemiology, Johns Hopkins University, Baltimore, Maryland, USA
| | - Katie M O'Brien
- Epidemiology Branch, National Institute of Environmental Health Sciences, National Institutes of Health, Department of Health and Human Services, Research Triangle Park, North Carolina, USA
| | - Kelly K Ferguson
- Epidemiology Branch, National Institute of Environmental Health Sciences, National Institutes of Health, Department of Health and Human Services, Research Triangle Park, North Carolina, USA
| | - Shanshan Zhao
- Biostatistics Branch, National Institute of Environmental Health Sciences, National Institutes of Health, Department of Health and Human Services, Research Triangle Park, North Carolina, USA
| | - Alexandra J White
- Epidemiology Branch, National Institute of Environmental Health Sciences, National Institutes of Health, Department of Health and Human Services, Research Triangle Park, North Carolina, USA
| |
Collapse
|
34
|
Waddy SP, Solomon AJ, Becerra AZ, Ward JB, Chan KE, Fwu CW, Norton JM, Eggers PW, Abbott KC, Kimmel PL. Racial/Ethnic Disparities in Atrial Fibrillation Treatment and Outcomes among Dialysis Patients in the United States. J Am Soc Nephrol 2020; 31:637-649. [PMID: 32079604 PMCID: PMC7062215 DOI: 10.1681/asn.2019050543] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2019] [Accepted: 12/10/2019] [Indexed: 12/28/2022] Open
Abstract
BACKGROUND Because stroke prevention is a major goal in the management of ESKD hemodialysis patients with atrial fibrillation, investigating racial/ethnic disparities in stroke among such patients is important to those who could benefit from strategies to maximize preventive measures. METHODS We used the United States Renal Data System to identify ESKD patients who initiated hemodialysis from 2006 to 2013 and then identified those with a subsequent atrial fibrillation diagnosis and Medicare Part A/B/D. Patients were followed for 1 year for all-cause stroke, mortality, prescription medications, and cardiovascular disease procedures. The survival mediational g-formula quantified the percentage of excess strokes attributable to lower use of atrial fibrillation treatments by race/ethnicity. RESULTS The study included 56,587 ESKD hemodialysis patients with atrial fibrillation. Black, white, Hispanic, and Asian patients accounted for 19%, 69%, 8%, and 3% of the population, respectively. Compared with white patients, black, Hispanic, or Asian patients were more likely to experience stroke (13%, 15%, and 16%, respectively) but less likely to fill a warfarin prescription (10%, 17%, and 28%, respectively). Warfarin prescription was associated with decreased stroke rates. Analyses suggested that equalizing the warfarin distribution to that in the white population would prevent 7%, 10%, and 12% of excess strokes among black, Hispanic, and Asian patients, respectively. We found no racial/ethnic disparities in all-cause mortality or use of cardiovascular disease procedures. CONCLUSIONS Racial/ethnic disparities in all-cause stroke among hemodialysis patients with atrial fibrillation are partially mediated by lower use of anticoagulants among black, Hispanic, and Asian patients. The reasons for these disparities are unknown, but strategies to maximize stroke prevention in minority hemodialysis populations should be further investigated.
Collapse
Affiliation(s)
- Salina P Waddy
- Department of Neurology, Atlanta Veterans Affairs Medical Center, Decatur, Georgia
| | - Allen J Solomon
- Division of Cardiology, Department of Medicine, George Washington University, Washington, DC
| | - Adan Z Becerra
- Department of Public Health Sciences, Social and Scientific Systems, Silver Spring, Maryland; and
| | - Julia B Ward
- Department of Public Health Sciences, Social and Scientific Systems, Silver Spring, Maryland; and
| | - Kevin E Chan
- Division of Kidney, Urologic, and Hematologic Diseases, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland
| | - Chyng-Wen Fwu
- Department of Public Health Sciences, Social and Scientific Systems, Silver Spring, Maryland; and
| | - Jenna M Norton
- Division of Kidney, Urologic, and Hematologic Diseases, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland
| | - Paul W Eggers
- Division of Kidney, Urologic, and Hematologic Diseases, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland
| | - Kevin C Abbott
- Division of Kidney, Urologic, and Hematologic Diseases, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland
| | - Paul L Kimmel
- Division of Kidney, Urologic, and Hematologic Diseases, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland
| |
Collapse
|
35
|
Breskin A, Westreich D, Hurt CB, Cole SR, Hudgens MG, Seaberg EC, Thio CL, Tien PC, Adimora AA. The Effects of Hepatitis C Treatment Eligibility Criteria on All-cause Mortality Among People With Human Immunodeficiency Virus. Clin Infect Dis 2019; 69:1613-1620. [PMID: 30615096 PMCID: PMC6792128 DOI: 10.1093/cid/ciz008] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2018] [Accepted: 01/04/2019] [Indexed: 01/26/2023] Open
Abstract
BACKGROUND The cost of direct-acting antivirals (DAAs) for hepatitis C virus (HCV) prompted many payers to restrict treatment to patients who met non-evidence-based criteria. These restrictions have implications for survival of people with HCV, especially for people with human immunodeficiency virus (HIV)/HCV coinfection who are at high risk for liver disease progression. The goal of this work was to estimate the effects of DAA access policies on 10-year all-cause mortality among people with HIV. METHODS The study population included 3056 adults with HIV in the Women's Interagency HIV Study and Multicenter AIDS Cohort Study from 1 October 1994 through 30 September 2015. We used the parametric g-formula to estimate 10-year all-cause mortality under DAA access policies that included treating (i) all people with HCV; (ii) only people with suppressed HIV; (iii) only people with severe fibrosis; and (iv) only people with HIV suppression and severe fibrosis. RESULTS The 10-year risk difference of treating all coinfected persons with DAAs compared with no treatment was -3.7% (95% confidence interval [CI], -9.1% to .6%). Treating only those with suppressed HIV and severe fibrosis yielded a risk difference of -1.1% (95% CI, -2.8% to .6%), with 51% (95% CI, 38%-59%) of coinfected persons receiving DAAs. Treating a random selection of 51% of coinfected persons at baseline decreased the risk by 1.9% (95% CI, -4.7% to .3%). CONCLUSIONS Restrictive DAA access policies may decrease survival compared to treating similar proportions of people with HIV/HCV coinfection with DAAs at random. These findings suggest that lives could be saved by thoughtfully revising access policies.
Collapse
Affiliation(s)
- Alexander Breskin
- Department of EpidemiologyInstitute for Global Health and Infectious Diseases, Baltimore, Maryland
| | - Daniel Westreich
- Department of EpidemiologyInstitute for Global Health and Infectious Diseases, Baltimore, Maryland
| | - Christopher B Hurt
- Institute for Global Health and Infectious Diseases, Baltimore, Maryland
| | - Stephen R Cole
- Department of EpidemiologyInstitute for Global Health and Infectious Diseases, Baltimore, Maryland
| | - Michael G Hudgens
- Department of Biostatistics, University of North Carolina at Chapel Hill, Baltimore, Maryland
| | - Eric C Seaberg
- Department Epidemiology, Johns Hopkins University, Baltimore, Maryland
| | - Chloe L Thio
- Department Medicine, Johns Hopkins University, Baltimore, Maryland
| | - Phyllis C Tien
- Department of Medicine, University of California, San Francisco
| | - Adaora A Adimora
- Department of EpidemiologyInstitute for Global Health and Infectious Diseases, Baltimore, Maryland
- Institute for Global Health and Infectious Diseases, Baltimore, Maryland
| |
Collapse
|
36
|
Kuehne F, Jahn B, Conrads-Frank A, Bundo M, Arvandi M, Endel F, Popper N, Endel G, Urach C, Gyimesi M, Murray EJ, Danaei G, Gaziano TA, Pandya A, Siebert U. Guidance for a causal comparative effectiveness analysis emulating a target trial based on big real world evidence: when to start statin treatment. J Comp Eff Res 2019; 8:1013-1025. [PMID: 31512926 DOI: 10.2217/cer-2018-0103] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Abstract
Aim: The aim of this project is to describe a causal (counterfactual) approach for analyzing when to start statin treatment to prevent cardiovascular disease using real-world evidence. Methods: We use directed acyclic graphs to operationalize and visualize the causal research question considering selection bias, potential time-independent and time-dependent confounding. We provide a study protocol following the 'target trial' approach and describe the data structure needed for the causal assessment. Conclusion: The study protocol can be applied to real-world data, in general. However, the structure and quality of the database play an essential role for the validity of the results, and database-specific potential for bias needs to be explicitly considered.
Collapse
Affiliation(s)
- Felicitas Kuehne
- Department of Public Health, Health Services Research & Health Technology Assessment, Institute of Public Health, Medical Decision Making & Health Technology Assessment, UMIT - University for Health Sciences, Medical Informatics & Technology, Hall iT, Austria
| | - Beate Jahn
- Department of Public Health, Health Services Research & Health Technology Assessment, Institute of Public Health, Medical Decision Making & Health Technology Assessment, UMIT - University for Health Sciences, Medical Informatics & Technology, Hall iT, Austria
| | - Annette Conrads-Frank
- Department of Public Health, Health Services Research & Health Technology Assessment, Institute of Public Health, Medical Decision Making & Health Technology Assessment, UMIT - University for Health Sciences, Medical Informatics & Technology, Hall iT, Austria
| | - Marvin Bundo
- Department of Public Health, Health Services Research & Health Technology Assessment, Institute of Public Health, Medical Decision Making & Health Technology Assessment, UMIT - University for Health Sciences, Medical Informatics & Technology, Hall iT, Austria
| | - Marjan Arvandi
- Department of Public Health, Health Services Research & Health Technology Assessment, Institute of Public Health, Medical Decision Making & Health Technology Assessment, UMIT - University for Health Sciences, Medical Informatics & Technology, Hall iT, Austria
| | | | - Niki Popper
- DEXHELPP, Vienna, Austria.,TU Wien, Research Unit of Information and Software Engineering, Austria.,dwh GmbH, Simulation Services & Technical Solutions, Austria
| | - Gottfried Endel
- Department EWG, Main Association of Austrian Social Security Institutions, Vienna, Austria
| | - Christoph Urach
- dwh GmbH, Simulation Services & Technical Solutions, Austria
| | | | - Eleanor J Murray
- Department of Epidemiology, Boston University School of Public Health, Boston, MA 02118, USA.,Department of Epidemiology, Harvard TH Chan School of Public Health, Boston, MA 02115, USA
| | - Goodarz Danaei
- Department of Epidemiology, Harvard TH Chan School of Public Health, Boston, MA 02115, USA.,Department of Global Health & Population, Harvard TH Chan School of Public Health, Boston, MA 02115, USA
| | - Thomas A Gaziano
- Department of Health Policy & Management, Center for Health Decision Science, Harvard TH Chan School of Public Health, Boston, MA 02115, USA.,Division of Cardiovascular Medicine, Brigham & Women's Hospital, Boston, MA 02115, USA
| | - Ankur Pandya
- Department of Health Policy & Management, Center for Health Decision Science, Harvard TH Chan School of Public Health, Boston, MA 02115, USA
| | - Uwe Siebert
- Department of Public Health, Health Services Research & Health Technology Assessment, Institute of Public Health, Medical Decision Making & Health Technology Assessment, UMIT - University for Health Sciences, Medical Informatics & Technology, Hall iT, Austria.,Department of Health Policy & Management, Center for Health Decision Science, Harvard TH Chan School of Public Health, Boston, MA 02115, USA.,Cardiovascular Research Program, Institute for Technology Assessment & Department of Radiology, Massachusetts General Hospital, Harvard Medical School, Boston, MA 02115, USA.,Division of Health Technology Assessment & Bioinformatics, ONCOTYROL - Center for Personalized Cancer Medicine, Innsbruck, Austria
| |
Collapse
|
37
|
Schomaker M, Luque-Fernandez MA, Leroy V, Davies MA. Using longitudinal targeted maximum likelihood estimation in complex settings with dynamic interventions. Stat Med 2019; 38:4888-4911. [PMID: 31436859 DOI: 10.1002/sim.8340] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2018] [Revised: 07/10/2019] [Accepted: 07/19/2019] [Indexed: 11/12/2022]
Abstract
Longitudinal targeted maximum likelihood estimation (LTMLE) has very rarely been used to estimate dynamic treatment effects in the context of time-dependent confounding affected by prior treatment when faced with long follow-up times, multiple time-varying confounders, and complex associational relationships simultaneously. Reasons for this include the potential computational burden, technical challenges, restricted modeling options for long follow-up times, and limited practical guidance in the literature. However, LTMLE has desirable asymptotic properties, ie, it is doubly robust, and can yield valid inference when used in conjunction with machine learning. It also has the advantage of easy-to-calculate analytic standard errors in contrast to the g-formula, which requires bootstrapping. We use a topical and sophisticated question from HIV treatment research to show that LTMLE can be used successfully in complex realistic settings, and we compare results to competing estimators. Our example illustrates the following practical challenges common to many epidemiological studies: (1) long follow-up time (30 months); (2) gradually declining sample size; (3) limited support for some intervention rules of interest; (4) a high-dimensional set of potential adjustment variables, increasing both the need and the challenge of integrating appropriate machine learning methods; and (5) consideration of collider bias. Our analyses, as well as simulations, shed new light on the application of LTMLE in complex and realistic settings: We show that (1) LTMLE can yield stable and good estimates, even when confronted with small samples and limited modeling options; (2) machine learning utilized with a small set of simple learners (if more complex ones cannot be fitted) can outperform a single, complex model, which is tailored to incorporate prior clinical knowledge; and (3) performance can vary considerably depending on interventions and their support in the data, and therefore critical quality checks should accompany every LTMLE analysis. We provide guidance for the practical application of LTMLE.
Collapse
Affiliation(s)
- M Schomaker
- Centre for Infectious Disease Epidemiology & Research, University of Cape Town, Cape Town, South Africa.,Institute of Public Health, Medical Decision Making and Health Technology Assessment, Department of Public Health, Health Services Research and Health Technology Assessment, UMIT - University for Health Sciences, Medical Informatics and Technology, Hall in Tirol, Austria
| | - M A Luque-Fernandez
- Biomedical Research Institute of Granada - Noncommunicable and Cancer Epidemiology Group, Andalusian School of Public Health, University of Granada, Granada, Spain.,Department of Epidemiology and Population Health, London School of Hygiene and Tropical Medicine, London, UK.,Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, Massachusetts
| | | | - M A Davies
- Centre for Infectious Disease Epidemiology & Research, University of Cape Town, Cape Town, South Africa
| |
Collapse
|
38
|
Albert JM, Cho JI, Liu Y, Nelson S. Generalized causal mediation and path analysis: Extensions and practical considerations. Stat Methods Med Res 2019; 28:1793-1807. [PMID: 29869589 PMCID: PMC6428612 DOI: 10.1177/0962280218776483] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
Causal mediation analysis seeks to decompose the effect of a treatment or exposure among multiple possible paths and provide casually interpretable path-specific effect estimates. Recent advances have extended causal mediation analysis to situations with a sequence of mediators or multiple contemporaneous mediators. However, available methods still have limitations, and computational and other challenges remain. The present paper provides an extended causal mediation and path analysis methodology. The new method, implemented in the new R package, gmediation (described in a companion paper), accommodates both a sequence (two stages) of mediators and multiple mediators at each stage, and allows for multiple types of outcomes following generalized linear models. The methodology can also handle unsaturated models and clustered data. Addressing other practical issues, we provide new guidelines for the choice of a decomposition, and for the choice of a reference group multiplier for the reduction of Monte Carlo error in mediation formula computations. The new method is applied to data from a cohort study to illuminate the contribution of alternative biological and behavioral paths in the effect of socioeconomic status on dental caries in adolescence.
Collapse
Affiliation(s)
- Jeffrey M. Albert
- Department of Population and Quantitative Health Sciences, Case Western Reserve University, Cleveland, OH, USA
| | - Jang Ik Cho
- Department of Population and Quantitative Health Sciences, Case Western Reserve University, Cleveland, OH, USA
| | - Yiying Liu
- Department of Population and Quantitative Health Sciences, Case Western Reserve University, Cleveland, OH, USA
| | - Suchitra Nelson
- Department of Community Dentistry, Case School of Dental Medicine, Cleveland, OH, USA
| |
Collapse
|
39
|
Abstract
The United States Environmental Protection Agency considers nutrient pollution in stream ecosystems one of the U.S.' most pressing environmental challenges. But limited independent replicates, lack of experimental randomization, and space- and time-varying confounding handicap causal inference on effects of nutrient pollution. In this paper the causal g-methods are extended to allow for exposures to vary in time and space in order to assess the effects of nutrient pollution on chlorophyll a - a proxy for algal production. Publicly available data from North Carolina's Cape Fear River and a simulation study are used to show how causal effects of upstream nutrient concentrations on downstream chlorophyll a levels may be estimated from typical water quality monitoring data. Estimates obtained from the parametric g-formula, a marginal structural model, and a structural nested model indicate that chlorophyll a concentrations at Lock and Dam 1 were influenced by nitrate concentrations measured 86 to 109 km upstream, an area where four major industrial and municipal point sources discharge wastewater.
Collapse
Affiliation(s)
- Bradley C Saul
- Department of Biostatistics, University of North Carolina Chapel Hill
| | - Michael G Hudgens
- Department of Biostatistics, University of North Carolina Chapel Hill
| | - Michael A Mallin
- Center for Marine Science, University of North Carolina Wilmington
| |
Collapse
|
40
|
Abstract
Epidemiologists often wish to estimate quantities that are easy to communicate and correspond to the results of realistic public health interventions. Methods from causal inference can answer these questions. We adopt the language of potential outcomes under Rubin's original Bayesian framework and show that the parametric g-formula is easily amenable to a Bayesian approach. We show that the frequentist properties of the Bayesian g-formula suggest it improves the accuracy of estimates of causal effects in small samples or when data are sparse. We demonstrate an approach to estimate the effect of environmental tobacco smoke on body mass index among children aged 4-9 years who were enrolled in a longitudinal birth cohort in New York, USA. We provide an algorithm and supply SAS and Stan code that can be adopted to implement this computational approach more generally.
Collapse
Affiliation(s)
- Alexander P. Keil
- Department of Epidemiology, University of North Carolina, Chapel Hill, USA
| | - Eric J. Daza
- Stanford Prevention Research Center, Stanford University School of Medicine, Palo Alto, USA
| | - Stephanie M. Engel
- Department of Epidemiology, University of North Carolina, Chapel Hill, USA
| | - Jessie P. Buckley
- Department of Epidemiology, University of North Carolina, Chapel Hill, USA
| | - Jessie K. Edwards
- Department of Epidemiology, University of North Carolina, Chapel Hill, USA
| |
Collapse
|
41
|
Rogawski ET, Liu J, Platts-Mills JA, Kabir F, Lertsethtakarn P, Siguas M, Khan SS, Praharaj I, Murei A, Nshama R, Mujaga B, Havt A, Maciel IA, Operario DJ, Taniuchi M, Gratz J, Stroup SE, Roberts JH, Kalam A, Aziz F, Qureshi S, Islam MO, Sakpaisal P, Silapong S, Yori PP, Rajendiran R, Benny B, McGrath M, Seidman JC, Lang D, Gottlieb M, Guerrant RL, Lima AAM, Leite JP, Samie A, Bessong PO, Page N, Bodhidatta L, Mason C, Shrestha S, Kiwelu I, Mduma ER, Iqbal NT, Bhutta ZA, Ahmed T, Haque R, Kang G, Kosek MN, Houpt ER. Use of quantitative molecular diagnostic methods to investigate the effect of enteropathogen infections on linear growth in children in low-resource settings: longitudinal analysis of results from the MAL-ED cohort study. LANCET GLOBAL HEALTH 2018; 6:e1319-e1328. [PMID: 30287125 PMCID: PMC6227248 DOI: 10.1016/s2214-109x(18)30351-6] [Citation(s) in RCA: 232] [Impact Index Per Article: 38.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/18/2018] [Revised: 06/28/2018] [Accepted: 07/11/2018] [Indexed: 12/15/2022]
Abstract
BACKGROUND Enteropathogen infections in early childhood not only cause diarrhoea but contribute to poor growth. We used molecular diagnostics to assess whether particular enteropathogens were associated with linear growth across seven low-resource settings. METHODS We used quantitative PCR to detect 29 enteropathogens in diarrhoeal and non-diarrhoeal stools collected from children in the first 2 years of life obtained during the Etiology, Risk Factors, and Interactions of Enteric Infections and Malnutrition and the Consequences for Child Health and Development (MAL-ED) multisite cohort study. Length was measured monthly. We estimated associations between aetiology-specific diarrhoea and subclinical enteropathogen infection and quantity and attained length in 3 month intervals, at age 2 and 5 years, and used a longitudinal model to account for temporality and time-dependent confounding. FINDINGS Among 1469 children who completed 2 year follow-up, 35 622 stool samples were tested and yielded valid results. Diarrhoeal episodes attributed to bacteria and parasites, but not viruses, were associated with small decreases in length after 3 months and at age 2 years. Substantial decrements in length at 2 years were associated with subclinical, non-diarrhoeal, infection with Shigella (length-for-age Z score [LAZ] reduction -0·14, 95% CI -0·27 to -0·01), enteroaggregative Escherichia coli (-0·21, -0·37 to -0·05), Campylobacter (-0·17, -0·32 to -0·01), and Giardia (-0·17, -0·30 to -0·05). Norovirus, Cryptosporidium, typical enteropathogenic E coli, and Enterocytozoon bieneusi were also associated with small decrements in LAZ. Shigella and E bieneusi were associated with the largest decreases in LAZ per log increase in quantity per g of stool (-0·13 LAZ, 95% CI -0·22 to -0·03 for Shigella; -0·14, -0·26 to -0·02 for E bieneusi). Based on these models, interventions that successfully decrease exposure to Shigella, enteroaggregative E coli, Campylobacter, and Giardia could increase mean length of children by 0·12-0·37 LAZ (0·4-1·2 cm) at the MAL-ED sites. INTERPRETATION Subclinical infection and quantity of pathogens, particularly Shigella, enteroaggregative E coli, Campylobacter, and Giardia, had a substantial negative association with linear growth, which was sustained during the first 2 years of life, and in some cases, to 5 years. Successfully reducing exposure to certain pathogens might reduce global stunting. FUNDING Bill & Melinda Gates Foundation.
Collapse
Affiliation(s)
- Elizabeth T Rogawski
- Division of Infectious Diseases and International Health, University of Virginia, Charlottesville, VA, USA; Department of Public Health Sciences, University of Virginia, Charlottesville, VA, USA
| | - Jie Liu
- Division of Infectious Diseases and International Health, University of Virginia, Charlottesville, VA, USA
| | - James A Platts-Mills
- Division of Infectious Diseases and International Health, University of Virginia, Charlottesville, VA, USA
| | | | | | | | - Shaila S Khan
- International Centre for Diarrhoeal Disease Research, Dhaka, Bangladesh
| | | | | | | | - Buliga Mujaga
- Kilimanjaro Clinical Research Institute, Moshi, Tanzania
| | | | | | - Darwin J Operario
- Division of Infectious Diseases and International Health, University of Virginia, Charlottesville, VA, USA
| | - Mami Taniuchi
- Division of Infectious Diseases and International Health, University of Virginia, Charlottesville, VA, USA
| | - Jean Gratz
- Division of Infectious Diseases and International Health, University of Virginia, Charlottesville, VA, USA
| | - Suzanne E Stroup
- Division of Infectious Diseases and International Health, University of Virginia, Charlottesville, VA, USA
| | - James H Roberts
- Department of Public Health Sciences, University of Virginia, Charlottesville, VA, USA
| | | | | | | | - M Ohedul Islam
- International Centre for Diarrhoeal Disease Research, Dhaka, Bangladesh
| | - Pimmada Sakpaisal
- Armed Forces Research Institute of Medical Sciences (AFRIMS), Bangkok, Thailand
| | - Sasikorn Silapong
- Armed Forces Research Institute of Medical Sciences (AFRIMS), Bangkok, Thailand
| | - Pablo P Yori
- Asociación Benéfica PRISMA, Iquitos, Peru; Bloomberg School of Public Health, Johns Hopkins University, Baltimore, MD, USA
| | | | | | - Monica McGrath
- Bloomberg School of Public Health, Johns Hopkins University, Baltimore, MD, USA; Fogarty International Center, National Institutes of Health, Bethesda, MD, USA
| | - Jessica C Seidman
- Fogarty International Center, National Institutes of Health, Bethesda, MD, USA
| | - Dennis Lang
- Foundation for the National Institutes of Health, Bethesda, MD, USA
| | - Michael Gottlieb
- Foundation for the National Institutes of Health, Bethesda, MD, USA
| | - Richard L Guerrant
- Division of Infectious Diseases and International Health, University of Virginia, Charlottesville, VA, USA
| | | | | | | | | | - Nicola Page
- National Institute for Communicable Diseases, Johannesburg, South Africa
| | - Ladaporn Bodhidatta
- Fogarty International Center, National Institutes of Health, Bethesda, MD, USA
| | - Carl Mason
- Armed Forces Research Institute of Medical Sciences (AFRIMS), Bangkok, Thailand
| | - Sanjaya Shrestha
- Walter Reed/AFRIMS Research Unit, Nepal, Kathmandu, Nepal; University of Bergen, Bergen, Norway
| | - Ireen Kiwelu
- Kilimanjaro Clinical Research Institute, Moshi, Tanzania
| | | | | | | | - Tahmeed Ahmed
- International Centre for Diarrhoeal Disease Research, Dhaka, Bangladesh
| | - Rashidul Haque
- International Centre for Diarrhoeal Disease Research, Dhaka, Bangladesh
| | | | - Margaret N Kosek
- Asociación Benéfica PRISMA, Iquitos, Peru; Bloomberg School of Public Health, Johns Hopkins University, Baltimore, MD, USA
| | - Eric R Houpt
- Division of Infectious Diseases and International Health, University of Virginia, Charlottesville, VA, USA.
| | | |
Collapse
|
42
|
Stoner MCD, Edwards JK, Miller WC, Aiello AE, Halpern CT, Julien A, Rucinski KB, Selin A, Twine R, Hughes JP, Wang J, Agyei Y, Gómez-Olivé FX, Wagner RG, Laeyendecker O, Macphail C, Kahn K, Pettifor A. Does Partner Selection Mediate the Relationship Between School Attendance and HIV/Herpes Simplex Virus-2 Among Adolescent Girls and Young Women in South Africa: An Analysis of HIV Prevention Trials Network 068 Data. J Acquir Immune Defic Syndr 2018; 79:20-27. [PMID: 29847479 PMCID: PMC6092209 DOI: 10.1097/qai.0000000000001766] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
Abstract
OBJECTIVE School attendance prevents HIV and herpes simplex virus-2 (HSV-2) in adolescent girls and young women, but the mechanisms to explain this relationship remain unclear. Our study assesses the extent to which characteristics of sex partners, partner age, and number mediate the relationship between attendance and risk of infection in adolescent girls and young women in South Africa. DESIGN We use longitudinal data from the HIV Prevention Trials Network 068 randomized controlled trial in rural South Africa, where girls were enrolled in early adolescence and followed in the main trial for more than 3 years. We examined older partners and the number of partners as possible mediators. METHODS We used the parametric g-formula to estimate 4-year risk differences for the effect of school attendance on the cumulative incidence of HIV/HSV-2 overall and the controlled direct effect (CDE) for mediation. We examined mediation separately and jointly for the mediators of interest. RESULTS We found that young women with high attendance in school had a lower cumulative incidence of HIV compared with those with low attendance (risk difference = -1.6%). Partner age difference (CDE = -1.2%) and the number of partners (CDE = -0.4%) mediated a large portion of this effect. In fact, when we accounted for the mediators jointly, the effect of schooling on HIV was almost removed, showing full mediation (CDE = -0.3%). The same patterns were observed for the relationship between school attendance and cumulative incidence of HSV-2 infection. CONCLUSION Increasing school attendance reduces the risk of acquiring HIV and HSV-2. Our results indicate the importance of school attendance in reducing partner number and partner age difference in this relationship.
Collapse
Affiliation(s)
- Marie C D Stoner
- Department of Epidemiology, University of North Carolina, Chapel Hill, NC
| | - Jessie K Edwards
- Department of Epidemiology, University of North Carolina, Chapel Hill, NC
| | - William C Miller
- Division of Epidemiology, The Ohio State University, Columbus, OH
| | - Allison E Aiello
- Department of Epidemiology, University of North Carolina, Chapel Hill, NC
| | - Carolyn T Halpern
- Department of Maternal and Child Health, University of North Carolina, Chapel Hill, NC
| | - Aimée Julien
- Department of Epidemiology, University of North Carolina, Chapel Hill, NC
- MRC/Wits Rural Public Health and Health Transitions Research Unit (Agincourt), School of Public Health, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa
| | | | - Amanda Selin
- Department of Epidemiology, University of North Carolina, Chapel Hill, NC
| | - Rhian Twine
- MRC/Wits Rural Public Health and Health Transitions Research Unit (Agincourt), School of Public Health, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa
| | - James P Hughes
- Department of Biostatistics, University of Washington, Seattle, WA
- Fred Hutchinson Cancer Research Center, Seattle, WA
| | - Jing Wang
- Fred Hutchinson Cancer Research Center, Seattle, WA
| | - Yaw Agyei
- School of Medicine, Johns Hopkins University, Baltimore, MD
| | - Francesc Xavier Gómez-Olivé
- MRC/Wits Rural Public Health and Health Transitions Research Unit (Agincourt), School of Public Health, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa
- INDEPTH Network, Accra, Ghana
| | - Ryan G Wagner
- MRC/Wits Rural Public Health and Health Transitions Research Unit (Agincourt), School of Public Health, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa
- Department of Public Health and Clinical Medicine, Epidemiology and Global Health Unit, Umeå University, Umeå, Sweden
| | - Oliver Laeyendecker
- Division of Intramural Research, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Baltimore, MD
| | - Catherine Macphail
- MRC/Wits Rural Public Health and Health Transitions Research Unit (Agincourt), School of Public Health, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa
- School of Health and Society, University of Wollongong, Wollongong, New South Wales, Australia
- Wits Reproductive Health and HIV Research Institute, University of the Witwatersrand, Johannesburg, South Africa
| | - Kathleen Kahn
- MRC/Wits Rural Public Health and Health Transitions Research Unit (Agincourt), School of Public Health, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa
- INDEPTH Network, Accra, Ghana
- Department of Public Health and Clinical Medicine, Epidemiology and Global Health Unit, Umeå University, Umeå, Sweden
| | - Audrey Pettifor
- Department of Epidemiology, University of North Carolina, Chapel Hill, NC
- MRC/Wits Rural Public Health and Health Transitions Research Unit (Agincourt), School of Public Health, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa
| |
Collapse
|
43
|
Keil AP, Edwards JK. A review of time scale fundamentals in the g-formula and insidious selection bias. CURR EPIDEMIOL REP 2018; 5:205-213. [PMID: 30555772 PMCID: PMC6289285] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
PURPOSE OF REVIEW We review recent examples of data analysis with the g-formula, a powerful tool for analyzing longitudinal data and survival analysis. Specifically, we focus on the common choices of time scale and review inferential issues that may arise. RECENT FINDINGS Researchers are increasingly engaged with questions that require time scales subject to left-truncation and right-censoring. The assumptions necessary for allowing right-censoring are well defined in the literature, whereas similar assumptions for left-truncation are not well defined. Policy and biologic considerations sometimes dictate that observational data must be analyzed on time scales that are subject to left-truncation, such as age. SUMMARY Further consideration of left-truncation is needed, especially when biologic or policy considerations dictate that age is the relevant time scale of interest. Methodologic development is needed to reduce potential for bias when left-truncation may occur.
Collapse
|
44
|
Keil AP, Richardson DB, Westreich D, Steenland K. Estimating the Impact of Changes to Occupational Standards for Silica Exposure on Lung Cancer Mortality. Epidemiology 2018; 29:658-665. [PMID: 29870429 PMCID: PMC6066423 DOI: 10.1097/ede.0000000000000867] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
Abstract
BACKGROUND Respiratory exposure to silica is associated with the risk of death owing to malignant and nonmalignant disease. 2.3 million US workers are exposed to silica. Occupational exposure limits for silica are derived from a number of lines of evidence, including observational studies. Observational studies may be subject to healthy worker survivor bias, which could result in underestimates of silica's impact on worker mortality and, in turn, bias risk estimates for occupational exposure limits. METHODS Using data on 65,999 workers pooled across multiple industries, we estimate the impacts of several hypothetical occupational exposure limits on silica exposure on lung cancer and all-cause mortality. We use the parametric g-formula, which can account for healthy worker survivor bias. RESULTS Assuming we could eliminate occupational exposure, we estimate that there would be 20.7 fewer deaths per 1,000 workers in our pooled study by age 80 (95% confidence interval = 14.5, 26.8), including 3.91 fewer deaths owing to lung cancer (95% CI = 1.53, 6.30). Less restrictive interventions demonstrated smaller but still substantial risk reductions. CONCLUSIONS Our results suggest that occupational exposure limits for silica can be further strengthened to reduce silica-associated mortality and illustrate how current risk analysis for occupational limits can be improved.
Collapse
Affiliation(s)
- Alexander P Keil
- From the Department of Epidemiology, University of North Carolina, Chapel Hill, NC
| | - David B Richardson
- From the Department of Epidemiology, University of North Carolina, Chapel Hill, NC
| | - Daniel Westreich
- From the Department of Epidemiology, University of North Carolina, Chapel Hill, NC
| | - Kyle Steenland
- Department of Environmental Health, Rollins School of Public Health, Emory University, Atlanta, GA
| |
Collapse
|
45
|
Edwards JK, Cole SR, Moore RD, Mathews WC, Kitahata M, Eron JJ. Sensitivity Analyses for Misclassification of Cause of Death in the Parametric G-Formula. Am J Epidemiol 2018; 187:1808-1816. [PMID: 29420696 DOI: 10.1093/aje/kwy028] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2017] [Accepted: 02/02/2018] [Indexed: 01/03/2023] Open
Abstract
Cause-specific mortality is an important outcome in studies of interventions to improve survival, yet causes of death can be misclassified. Here, we present an approach to performing sensitivity analyses for misclassification of cause of death in the parametric g-formula. The g-formula is a useful method to estimate effects of interventions in epidemiologic research because it appropriately accounts for time-varying confounding affected by prior treatment and can estimate risk under dynamic treatment plans. We illustrate our approach using an example comparing acquired immune deficiency syndrome (AIDS)-related mortality under immediate and delayed treatment strategies in a cohort of therapy-naive adults entering care for human immunodeficiency virus infection in the United States. In the standard g-formula approach, 10-year risk of AIDS-related mortality under delayed treatment was 1.73 (95% CI: 1.17, 2.54) times the risk under immediate treatment. In a sensitivity analysis assuming that AIDS-related death was measured with sensitivity of 95% and specificity of 90%, the 10-year risk ratio comparing AIDS-related mortality between treatment plans was 1.89 (95% CI: 1.13, 3.14). When sensitivity and specificity are unknown, this approach can be used to estimate the effects of dynamic treatment plans under a range of plausible values of sensitivity and specificity of the recorded event type.
Collapse
Affiliation(s)
- Jessie K Edwards
- Department of Epidemiology, Gillings School of Global Public Health, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina
| | - Stephen R Cole
- Department of Epidemiology, Gillings School of Global Public Health, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina
| | - Richard D Moore
- School of Medicine, Johns Hopkins University, Baltimore, Maryland
| | | | - Mari Kitahata
- Department of Medicine, University of Washington, Seattle, Washington
| | - Joseph J Eron
- School of Medicine, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina
| |
Collapse
|
46
|
Commentary: The Limits of Risk Factors Revisited: Is It Time for a Causal Architecture Approach? Epidemiology 2018; 28:1-5. [PMID: 27755276 DOI: 10.1097/ede.0000000000000578] [Citation(s) in RCA: 36] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
|
47
|
Keil AP, Edwards JK. A Review of Time Scale Fundamentals in the g-Formula and Insidious Selection Bias. CURR EPIDEMIOL REP 2018. [DOI: 10.1007/s40471-018-0153-0] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
48
|
Naimi AI, Cole SR, Kennedy EH. An introduction to g methods. Int J Epidemiol 2018; 46:756-762. [PMID: 28039382 DOI: 10.1093/ije/dyw323] [Citation(s) in RCA: 81] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 10/17/2016] [Indexed: 11/13/2022] Open
Abstract
Robins' generalized methods (g methods) provide consistent estimates of contrasts (e.g. differences, ratios) of potential outcomes under a less restrictive set of identification conditions than do standard regression methods (e.g. linear, logistic, Cox regression). Uptake of g methods by epidemiologists has been hampered by limitations in understanding both conceptual and technical details. We present a simple worked example that illustrates basic concepts, while minimizing technical complications.
Collapse
Affiliation(s)
| | - Stephen R Cole
- Department of Epidemiology, University of North Carolina at Chapel Hill and
| | | |
Collapse
|
49
|
Schomaker M, Leroy V, Wolfs T, Technau KG, Renner L, Judd A, Sawry S, Amorissani-Folquet M, Noguera-Julian A, Tanser F, Eboua F, Navarro ML, Chimbetete C, Amani-Bosse C, Warszawski J, Phiri S, N'Gbeche S, Cox V, Koueta F, Giddy J, Sygnaté-Sy H, Raben D, Chêne G, Davies MA. Optimal timing of antiretroviral treatment initiation in HIV-positive children and adolescents: a multiregional analysis from Southern Africa, West Africa and Europe. Int J Epidemiol 2018; 46:453-465. [PMID: 27342220 DOI: 10.1093/ije/dyw097] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/31/2016] [Indexed: 01/08/2023] Open
Abstract
Background There is limited knowledge about the optimal timing of antiretroviral treatment initiation in older children and adolescents. Methods A total of 20 576 antiretroviral treatment (ART)-naïve patients, aged 1-16 years at enrolment, from 19 cohorts in Europe, Southern Africa and West Africa, were included. We compared mortality and growth outcomes for different ART initiation criteria, aligned with previous and recent World Health Organization criteria, for 5 years of follow-up, adjusting for all measured baseline and time-dependent confounders using the g-formula. Results Median (1st;3rd percentile) CD4 count at baseline was 676 cells/mm 3 (394; 1037) (children aged ≥ 1 and < 5 years), 373 (172; 630) (≥ 5 and < 10 years) and 238 (88; 425) (≥ 10 and < 16 years). There was a general trend towards lower mortality and better growth with earlier treatment initiation. In children < 10 years old at enrolment, by 5 years of follow-up there was lower mortality and a higher mean height-for-age z-score with immediate ART initiation versus delaying until CD4 count < 350 cells/mm 3 (or CD4% < 15% or weight-for-age z-score < -2) with absolute differences in mortality and height-for-age z-score of 0.3% (95% confidence interval: 0.1%; 0.6%) and -0.08 (-0.09; -0.06) (≥ 1 and < 5 years), and 0.3% (0.04%; 0.5%) and -0.07 (-0.08; -0.05) (≥ 5 and < 10 years). In those aged > 10 years at enrolment we did not find any difference in mortality or growth with immediate ART initiation, with estimated differences of -0.1% (-0.2%; 0.6%) and -0.03 (-0.05; 0.00), respectively. Growth differences in children aged < 10 years persisted for treatment thresholds using higher CD4 values. Regular follow-up led to better height and mortality outcomes. Conclusions Immediate ART is associated with lower mortality and better growth for up to 5 years in children < 10 years old. Our results on adolescents were inconclusive.
Collapse
Affiliation(s)
- Michael Schomaker
- Centre for Infectious Disease Epidemiology and Research, University of Cape Town, Cape Town, South Africa
| | - Valeriane Leroy
- Inserm, U1027, Université Paul Sabatier Toulouse 3 Toulouse, France
| | - Tom Wolfs
- Children's Hospital/UMCU, Department of Infectious Diseases, Utrecht, The Netherlands
| | - Karl-Günter Technau
- Department of Paediatrics & Child Health, Rahima Moosa Mother and Child Hospital and University of the Witwatersrand, Johannesburg, South Africa.,Empilweni Services and Research Unit, Department of Paediatrics & Child Health, Rahima Moosa Mother and Child Hospital and University of the Witwatersrand, Johannesburg, South Africa
| | - Lorna Renner
- University of Ghana Medical School, Accra, Ghana
| | - Ali Judd
- MRC Clinical Trials Unit, University College London, London, UK
| | - Shobna Sawry
- University of the Witwatersrand, Wits Reproductive Health and HIV Institute, Chris Hani Baragwanath Academic Hospital, Soweto, South Africa
| | | | - Antoni Noguera-Julian
- Pediatrics Department, Hospital Sant Joan de Déu, Universitat de Barcelona, Barcelona, Spain
| | - Frank Tanser
- Africa Centre for Health and Population Studies, University of KwaZulu-Natal, Somkhele, South Africa.,School of Nursing and Public Health, University of KwaZulu-Natal, Durban, South Africa.,Centre for the AIDS Programme of Research in South Africa - CAPRISA, University of KwaZulu-Natal, Congella, South Africa
| | | | | | | | | | - Josiane Warszawski
- Centre de recherche en épidémiologie et santé des populations, 1018 Inserm, France
| | - Sam Phiri
- Lighthouse Trust Clinic, Kamuzu Central Hospital, Lilongwe, Malawi
| | - Sylvie N'Gbeche
- Centre de Prise en Charge de Recherche et de Formation Enfants, Abidjan, Côte d'Ivoire
| | - Vivian Cox
- Médecins Sans Frontiéres South Africa, Cape Town, South Africa
| | - Fla Koueta
- Charles de Gaulle University Hospital, Ouagadougou, Burkina Faso
| | - Janet Giddy
- Sinikithemba Clinic, McCord Hospital, Durban, South Africa
| | | | - Dorthe Raben
- Department of Infectious Diseases, Rigshospitalet, University of Copenhagen, Copenhagen, Denmark
| | - Geneviève Chêne
- University of Bordeaux Bordeaux, ISPED, Centre INSERM U1219-Bordeaux Population Health, F-33000 Bordeaux, France.,INSERM, ISPED, Centre INSERM U1219-Bordeaux Population Health, F-33000 Bordeaux, France.,CHU de Bordeaux, Pôle de santé publique, Service d information médicale, F-33000 Bordeaux, France
| | - Mary-Ann Davies
- Centre for Infectious Disease Epidemiology and Research, University of Cape Town, Cape Town, South Africa
| | | |
Collapse
|
50
|
Affiliation(s)
- Sandro Galea
- Boston University School of Public Health, 715 Albany Street, Talbot 301, Boston, MA, 02118, USA.
| |
Collapse
|