RESUMEN
Background: Post-acute sequelae of SARS-CoV-2, referred to as "long COVID", are a globally pervasive threat. While their many clinical determinants are commonly considered, their plausible social correlates are often overlooked. Methods: Here, we use data from a multinational prospective cohort study to compare social and clinical predictors of differences in quality of life with long COVID. We further measure the extent to which clinical intermediates may explain relationships between social variables and quality of life with long COVID. Findings: Beyond age, neuropsychological and rheumatological comorbidities, educational attainment, employment status, and female sex were important predictors of long COVID-associated quality of life days (long COVID QALDs). Furthermore, most of their associations could not be attributed to key long COVID-predicting comorbidities. In Norway, 90% (95% CI: 77%, 100%) of the adjusted association between belonging to the top two quintiles of educational attainment and long COVID QALDs was not explained by these clinical intermediates. The same was true for 86% (73%, 100%) and 93% (80%,100%) of the adjusted association between full-time employment and long COVID QALDs in the United Kingdom (UK) and Russia. Additionally, 77% (46%,100%) and 73% (52%, 94%) of the adjusted associations between female sex and long COVID QALDs in Norway and the UK were unexplained by the clinical mediators. Interpretation: Our findings highlight that socio-economic proxies and sex are key predictors of long COVID QALDs and that other (non-clinical) mechanisms drive their observed relationships. Importantly, we outline a multi-method, adaptable causal approach for evaluating the isolated contributions of social disparities to experiences with long COVID. Funding: UK Foreign, Commonwealth and Development Office; Wellcome Trust; Bill & Melinda Gates Foundation; Oxford COVID-19 Research Response Funding; UK National Institute for Health and Care Research; UK Medical Research Council; Public Health England; Liverpool Experimental Cancer Medicine Centre; Research Council of Norway; Vivaldi Invest A/S; South Eastern Norway Health Authority.
RESUMEN
BACKGROUND: The first licensed malaria vaccine, RTS,S/AS01E, confers moderate protection against symptomatic disease. Because many malaria infections are asymptomatic, we conducted a large-scale longitudinal parasite genotyping study of samples from a clinical trial exploring how vaccine dosing regimen affects vaccine efficacy. METHODS: Between Sept 28, 2017, and Sept 25, 2018, 1500 children aged 5-17 months were randomly assigned (1:1:1:1:1) to receive four different RTS,S/AS01E regimens or a rabies control vaccine in a phase 2b open-label clinical trial in Ghana and Kenya. Participants in the four RTS,S groups received two full doses at month 0 and month 1 and either full doses at month 2 and month 20 (group R012-20); full doses at month 2, month 14, month 26, and month 38 (group R012-14); fractional doses at month 2, month 14, month 26, and month 38 (group Fx012-14; early fourth dose); or fractional doses at month 7, month 20, and month 32 (group Fx017-20; delayed third dose). We evaluated the time to the first new genotypically detected infection and the total number of new infections during two follow-up periods (12 months and 20 months) in more than 36 000 dried blood spot specimens from 1500 participants. To study vaccine effects on time to the first new infection, we defined vaccine efficacy as one minus the hazard ratio (HR; RTS,S vs control) of the first new infection. We performed a post-hoc analysis of vaccine efficacy based on malaria infection status at first vaccination and force of infection by month 2. This trial (MAL-095) is registered with ClinicalTrials.gov, NCT03281291. FINDINGS: We observed significant and similar vaccine efficacy (25-43%; 95% CI union 9-53) against first new infection for all four RTS,S/AS01E regimens across both follow-up periods (12 months and 20 months). Each RTS,S/AS01E regimen significantly reduced the mean number of new infections in the 20-month follow-up period by 1·1-1·6 infections (95% CI union 0·6-2·1). Vaccine efficacy against first new infection was significantly higher in participants who were infected with malaria (68%; 95% CI 50-80) than in those who were uninfected (37%; 23-48) at the first vaccination (p=0·0053). INTERPRETATION: All tested dosing regimens blocked some infections to a similar degree. Improved vaccine efficacy in participants infected during vaccination could suggest new strategies for highly efficacious malaria vaccine development and implementation. FUNDING: GlaxoSmithKline Biologicals SA, PATH, Bill & Melinda Gates Foundation, and the German Federal Ministry of Education and Research.
Asunto(s)
Vacunas contra la Malaria , Malaria Falciparum , Humanos , Vacunas contra la Malaria/administración & dosificación , Vacunas contra la Malaria/inmunología , Ghana , Kenia/epidemiología , Lactante , Masculino , Femenino , Malaria Falciparum/prevención & control , Malaria Falciparum/epidemiología , Genotipo , Estudios Longitudinales , Eficacia de las Vacunas , Plasmodium falciparum/inmunología , Plasmodium falciparum/genética , Malaria/prevención & controlRESUMEN
A central goal of vaccine research is to characterize and validate immune correlates of protection (CoPs). In addition to helping elucidate immunological mechanisms, a CoP can serve as a valid surrogate endpoint for an infectious disease clinical outcome and thus qualifies as a primary endpoint for vaccine authorization or approval without requiring resource-intensive randomized, controlled phase 3 trials. Yet, it is challenging to persuasively validate a CoP, because a prognostic immune marker can fail as a reliable basis for predicting/inferring the level of vaccine efficacy against a clinical outcome, and because the statistical analysis of phase 3 trials only has limited capacity to disentangle association from cause. Moreover, the multitude of statistical methods garnered for CoP evaluation in phase 3 trials renders the comparison, interpretation, and synthesis of CoP results challenging. Toward promoting broader harmonization and standardization of CoP evaluation, this article summarizes four complementary statistical frameworks for evaluating CoPs in a phase 3 trial, focusing on the frameworks' distinct scientific objectives as measured and communicated by distinct causal vaccine efficacy parameters. Advantages and disadvantages of the frameworks are considered, dependent on phase 3 trial context, and perspectives are offered on how the frameworks can be applied and their results synthesized.
Asunto(s)
Eficacia de las Vacunas , Vacunas , Proyectos de Investigación , Biomarcadores/análisis , Causalidad , Ensayos Clínicos Controlados Aleatorios como AsuntoRESUMEN
Longitudinal modified treatment policies (LMTP) have been recently developed as a novel method to define and estimate causal parameters that depend on the natural value of treatment. LMTPs represent an important advancement in causal inference for longitudinal studies as they allow the non-parametric definition and estimation of the joint effect of multiple categorical, ordinal, or continuous treatments measured at several time points. We extend the LMTP methodology to problems in which the outcome is a time-to-event variable subject to a competing event that precludes observation of the event of interest. We present identification results and non-parametric locally efficient estimators that use flexible data-adaptive regression techniques to alleviate model misspecification bias, while retaining important asymptotic properties such as [Formula: see text]-consistency. We present an application to the estimation of the effect of the time-to-intubation on acute kidney injury amongst COVID-19 hospitalized patients, where death by other causes is taken to be the competing event.
Asunto(s)
Modelos Estadísticos , Análisis de Supervivencia , Humanos , Simulación por Computador , Estudios Longitudinales , Análisis de RegresiónRESUMEN
Background: The only licensed malaria vaccine, RTS,S/AS01 E , confers moderate protection against symptomatic disease. Because many malaria infections are asymptomatic, we conducted a large-scale longitudinal parasite genotyping study of samples from a clinical trial exploring how vaccine dosing regimen affects vaccine efficacy (VE). Methods: 1,500 children aged 5-17 months were randomized to receive four different RTS,S/AS01 E regimens or a rabies control vaccine in a phase 2b clinical trial in Ghana and Kenya. We evaluated the time to the first new genotypically detected infection and the total number of new infections during two follow-up periods in over 36K participant specimens. We performed a post hoc analysis of VE based on malaria infection status at first vaccination and force of infection. Results: We observed significant and comparable VE (25-43%, 95% CI union 9-53%) against first new infection for all four RTS,S/AS01 E regimens across both follow-up periods (12 and 20 months). Each RTS,S/AS01 E regimen significantly reduced the number of new infections in the 20-month follow-up period (control mean 4.1 vs. RTS,S/AS01 E mean 2.6-3.0). VE against first new infection was significantly higher in participants who were malaria-infected (68%; 95% CI, 50 to 80%) versus uninfected (37%; 95% CI, 23 to 48%) at the first vaccination (P=0.0053) and in participants experiencing greater force of infection between dose 1 and 3 (P=0.059). Conclusions: All tested dosing regimens blocked some infections to a similar degree. Improved VE in participants infected during vaccination could suggest new strategies for highly efficacious malaria vaccine development and implementation. ( ClinicalTrials.gov number, NCT03276962 ).
RESUMEN
The COVE trial randomized participants to receive two doses of mRNA-1273 vaccine or placebo on Days 1 and 29 (D1, D29). Anti-SARS-CoV-2 Spike IgG binding antibodies (bAbs), anti-receptor binding domain IgG bAbs, 50% inhibitory dilution neutralizing antibody (nAb) titers, and 80% inhibitory dilution nAb titers were measured at D29 and D57. We assessed these markers as correlates of protection (CoPs) against COVID-19 using stochastic interventional vaccine efficacy (SVE) analysis and principal surrogate (PS) analysis, frameworks not used in our previous COVE immune correlates analyses. By SVE analysis, hypothetical shifts of the D57 Spike IgG distribution from a geometric mean concentration (GMC) of 2737 binding antibody units (BAU)/mL (estimated vaccine efficacy (VE): 92.9% (95% CI: 91.7%, 93.9%)) to 274 BAU/mL or to 27,368 BAU/mL resulted in an overall estimated VE of 84.2% (79.0%, 88.1%) and 97.6% (97.4%, 97.7%), respectively. By binary marker PS analysis of Low and High subgroups (cut-point: 2094 BAU/mL), the ignorance interval (IGI) and estimated uncertainty interval (EUI) for VE were [85%, 90%] and (78%, 93%) for Low compared to [95%, 96%] and (92%, 97%) for High. By continuous marker PS analysis, the IGI and 95% EUI for VE at the 2.5th percentile (519.4 BAU/mL) vs. at the 97.5th percentile (9262.9 BAU/mL) of D57 Spike IgG concentration were [92.6%, 93.4%] and (89.2%, 95.7%) vs. [94.3%, 94.6%] and (89.7%, 97.0%). Results were similar for other D29 and D57 markers. Thus, the SVE and PS analyses additionally support all four markers at both time points as CoPs.
Asunto(s)
Vacuna nCoV-2019 mRNA-1273 , COVID-19 , Humanos , Anticuerpos Neutralizantes , Anticuerpos Antivirales , COVID-19/prevención & control , Inmunoglobulina G , Eficacia de las VacunasRESUMEN
BACKGROUND: Stochastic interventional vaccine efficacy (SVE) analysis is a new approach to correlate of protection (CoP) analysis of a phase III trial that estimates how vaccine efficacy (VE) would change under hypothetical shifts of an immune marker. METHODS: We applied nonparametric SVE methodology to the COVE trial of messenger RNA-1273 vs placebo to evaluate post-dose 2 pseudovirus neutralizing antibody (nAb) titer against the D614G strain as a CoP against COVID-19. Secondly, we evaluated the ability of these results to predict VE against variants based on shifts of geometric mean titers to variants vs D614G. Prediction accuracy was evaluated by 13 validation studies, including 12 test-negative designs. RESULTS: SVE analysis of COVE supported post-dose 2 D614G titer as a CoP: estimated VE ranged from 66.9% (95% confidence interval: 36.2, 82.8%) to 99.3% (99.1, 99.4%) at 10-fold decreased or increased titer shifts, respectively. The SVE estimates only weakly predicted variant-specific VE estimates (concordance correlation coefficient 0.062 for post 2-dose VE). CONCLUSION: SVE analysis of COVE supports nAb titer as a CoP for messenger RNA vaccines. Predicting variant-specific VE proved difficult due to many limitations. Greater anti-Omicron titers may be needed for high-level protection against Omicron vs anti-D614G titers needed for high-level protection against pre-Omicron COVID-19.
Asunto(s)
COVID-19 , Vacunas , Humanos , Anticuerpos Neutralizantes , Anticuerpos Antivirales , COVID-19/prevención & control , ARN Mensajero/genéticaRESUMEN
The covariance matrix plays a fundamental role in many modern exploratory and inferential statistical procedures, including dimensionality reduction, hypothesis testing, and regression. In low-dimensional regimes, where the number of observations far exceeds the number of variables, the optimality of the sample covariance matrix as an estimator of this parameter is well-established. High-dimensional regimes do not admit such a convenience. Thus, a variety of estimators have been derived to overcome the shortcomings of the canonical estimator in such settings. Yet, selecting an optimal estimator from among the plethora available remains an open challenge. Using the framework of cross-validated loss-based estimation, we develop the theoretical underpinnings of just such an estimator selection procedure. We propose a general class of loss functions for covariance matrix estimation and establish accompanying finite-sample risk bounds and conditions for the asymptotic optimality of the cross-validation selector. In numerical experiments, we demonstrate the optimality of our proposed selector in moderate sample sizes and across diverse data-generating processes. The practical benefits of our procedure are highlighted in a dimension reduction application to single-cell transcriptome sequencing data.
RESUMEN
The best assay or marker to define mRNA-1273 vaccine-induced antibodies as a correlate of protection (CoP) is unclear. In the COVE trial, participants received two doses of the mRNA-1273 COVID-19 vaccine or placebo. We previously assessed IgG binding antibodies to the spike protein (spike IgG) or receptor binding domain (RBD IgG) and pseudovirus neutralizing antibody 50 or 80% inhibitory dilution titer measured on day 29 or day 57, as correlates of risk (CoRs) and CoPs against symptomatic COVID-19 over 4 months after dose. Here, we assessed a new marker, live virus 50% microneutralization titer (LV-MN50), and compared and combined markers in multivariable analyses. LV-MN50 was an inverse CoR, with a hazard ratio of 0.39 (95% confidence interval, 0.19 to 0.83) at day 29 and 0.51 (95% confidence interval, 0.25 to 1.04) at day 57 per 10-fold increase. In multivariable analyses, pseudovirus neutralization titers and anti-spike binding antibodies performed best as CoRs; combining antibody markers did not improve correlates. Pseudovirus neutralization titer was the strongest independent correlate in a multivariable model. Overall, these results supported pseudovirus neutralizing and binding antibody assays as CoRs and CoPs, with the live virus assay as a weaker correlate in this sample set. Day 29 markers performed as well as day 57 markers as CoPs, which could accelerate immunogenicity and immunobridging studies.
Asunto(s)
Vacuna nCoV-2019 mRNA-1273 , COVID-19 , Humanos , Eficacia de las Vacunas , COVID-19/prevención & control , Anticuerpos Neutralizantes , Inmunoglobulina G , Anticuerpos AntiviralesRESUMEN
In the phase 3 trial of the AZD1222 (ChAdOx1 nCoV-19) vaccine conducted in the U.S., Chile, and Peru, anti-spike binding IgG concentration (spike IgG) and pseudovirus 50% neutralizing antibody titer (nAb ID50) measured four weeks after two doses were assessed as correlates of risk and protection against PCR-confirmed symptomatic SARS-CoV-2 infection (COVID-19). These analyses of SARS-CoV-2 negative participants were based on case-cohort sampling of vaccine recipients (33 COVID-19 cases by 4 months post dose two, 463 non-cases). The adjusted hazard ratio of COVID-19 was 0.32 (95% CI: 0.14, 0.76) per 10-fold increase in spike IgG concentration and 0.28 (0.10, 0.77) per 10-fold increase in nAb ID50 titer. At nAb ID50 below the limit of detection (< 2.612 IU50/ml), 10, 100, and 270 IU50/ml, vaccine efficacy was -5.8% (-651%, 75.6%), 64.9% (56.4%, 86.9%), 90.0% (55.8%, 97.6%) and 94.2% (69.4%, 99.1%). These findings provide further evidence towards defining an immune marker correlate of protection to help guide regulatory/approval decisions for COVID-19 vaccines.
RESUMEN
In the PREVENT-19 phase 3 trial of the NVX-CoV2373 vaccine (NCT04611802), anti-spike binding IgG concentration (spike IgG), anti-RBD binding IgG concentration (RBD IgG), and pseudovirus 50% neutralizing antibody titer (nAb ID50) measured two weeks post-dose two are assessed as correlates of risk and as correlates of protection against COVID-19. Analyses are conducted in the U.S. cohort of baseline SARS-CoV-2 negative per-protocol participants using a case-cohort design that measures the markers from all 12 vaccine recipient breakthrough COVID-19 cases starting 7 days post antibody measurement and from 639 vaccine recipient non-cases. All markers are inversely associated with COVID-19 risk and directly associated with vaccine efficacy. In vaccine recipients with nAb ID50 titers of 50, 100, and 7230 international units (IU50)/ml, vaccine efficacy estimates are 75.7% (49.8%, 93.2%), 81.7% (66.3%, 93.2%), and 96.8% (88.3%, 99.3%). The results support potential cross-vaccine platform applications of these markers for guiding decisions about vaccine approval and use.
Asunto(s)
COVID-19 , Humanos , Anticuerpos Neutralizantes , Anticuerpos Antivirales , COVID-19/prevención & control , Vacunas contra la COVID-19 , Inmunoglobulina G , SARS-CoV-2 , Eficacia de las Vacunas , Ensayos Clínicos Fase III como AsuntoRESUMEN
Inverse-probability-weighted estimators are the oldest and potentially most commonly used class of procedures for the estimation of causal effects. By adjusting for selection biases via a weighting mechanism, these procedures estimate an effect of interest by constructing a pseudopopulation in which selection biases are eliminated. Despite their ease of use, these estimators require the correct specification of a model for the weighting mechanism, are known to be inefficient, and suffer from the curse of dimensionality. We propose a class of nonparametric inverse-probability-weighted estimators in which the weighting mechanism is estimated via undersmoothing of the highly adaptive lasso, a nonparametric regression function proven to converge at nearly n - 1 / 3 $ n^{-1/3}$ -rate to the true weighting mechanism. We demonstrate that our estimators are asymptotically linear with variance converging to the nonparametric efficiency bound. Unlike doubly robust estimators, our procedures require neither derivation of the efficient influence function nor specification of the conditional outcome model. Our theoretical developments have broad implications for the construction of efficient inverse-probability-weighted estimators in large statistical models and a variety of problem settings. We assess the practical performance of our estimators in simulation studies and demonstrate use of our proposed methodology with data from a large-scale epidemiologic study.
Asunto(s)
Modelos Estadísticos , Probabilidad , Simulación por Computador , Sesgo de Selección , CausalidadRESUMEN
The widespread availability of high-dimensional biological data has made the simultaneous screening of many biological characteristics a central problem in computational and high-dimensional biology. As the dimensionality of datasets continues to grow, so too does the complexity of identifying biomarkers linked to exposure patterns. The statistical analysis of such data often relies upon parametric modeling assumptions motivated by convenience, inviting opportunities for model misspecification. While estimation frameworks incorporating flexible, data adaptive regression strategies can mitigate this, their standard variance estimators are often unstable in high-dimensional settings, resulting in inflated Type-I error even after standard multiple testing corrections. We adapt a shrinkage approach compatible with parametric modeling strategies to semiparametric variance estimators of a family of efficient, asymptotically linear estimators of causal effects, defined by counterfactual exposure contrasts. Augmenting the inferential stability of these estimators in high-dimensional settings yields a data adaptive approach for robustly uncovering stable causal associations, even when sample sizes are limited. Our generalized variance estimator is evaluated against appropriate alternatives in numerical experiments, and an open source R/Bioconductor package, biotmle, is introduced. The proposal is demonstrated in an analysis of high-dimensional DNA methylation data from an observational study on the epigenetic effects of tobacco smoking.
Asunto(s)
Biología , Proyectos de Investigación , Tamaño de la Muestra , CausalidadRESUMEN
Causal mediation analysis has historically been limited in two important ways: (i) a focus has traditionally been placed on binary exposures and static interventions and (ii) direct and indirect effect decompositions have been pursued that are only identifiable in the absence of intermediate confounders affected by exposure. We present a theoretical study of an (in)direct effect decomposition of the population intervention effect, defined by stochastic interventions jointly applied to the exposure and mediators. In contrast to existing proposals, our causal effects can be evaluated regardless of whether an exposure is categorical or continuous and remain well-defined even in the presence of intermediate confounders affected by exposure. Our (in)direct effects are identifiable without a restrictive assumption on cross-world counterfactual independencies, allowing for substantive conclusions drawn from them to be validated in randomized controlled trials. Beyond the novel effects introduced, we provide a careful study of nonparametric efficiency theory relevant for the construction of flexible, multiply robust estimators of our (in)direct effects, while avoiding undue restrictions induced by assuming parametric models of nuisance parameter functionals. To complement our nonparametric estimation strategy, we introduce inferential techniques for constructing confidence intervals and hypothesis tests, and discuss open-source software, the $\texttt{medshift}$$\texttt{R}$ package, implementing the proposed methodology. Application of our (in)direct effects and their nonparametric estimators is illustrated using data from a comparative effectiveness trial examining the direct and indirect effects of pharmacological therapeutics on relapse to opioid use disorder.
Asunto(s)
Análisis de Mediación , Modelos Estadísticos , Humanos , Modelos Teóricos , CausalidadRESUMEN
Measuring immune correlates of disease acquisition and protection in the context of a clinical trial is a prerequisite for improved vaccine design. We analysed binding and neutralizing antibody measurements 4 weeks post vaccination as correlates of risk of moderate to severe-critical COVID-19 through 83 d post vaccination in the phase 3, double-blind placebo-controlled phase of ENSEMBLE, an international randomized efficacy trial of a single dose of Ad26.COV2.S. We also evaluated correlates of protection in the trial cohort. Of the three antibody immune markers we measured, we found most support for 50% inhibitory dilution (ID50) neutralizing antibody titre as a correlate of risk and of protection. The outcome hazard ratio was 0.49 (95% confidence interval 0.29, 0.81; P = 0.006) per 10-fold increase in ID50; vaccine efficacy was 60% (43%, 72%) at non-quantifiable ID50 (<2.7 IU50 ml-1) and increased to 89% (78%, 96%) at ID50 = 96.3 IU50 ml-1. Comparison of the vaccine efficacy by ID50 titre curves for ENSEMBLE-US, the COVE trial of the mRNA-1273 vaccine and the COV002-UK trial of the AZD1222 vaccine supported the ID50 titre as a correlate of protection across trials and vaccine types.
Asunto(s)
Ad26COVS1 , COVID-19 , Humanos , COVID-19/prevención & control , ChAdOx1 nCoV-19 , Vacuna nCoV-2019 mRNA-1273 , Eficacia de las Vacunas , Anticuerpos NeutralizantesRESUMEN
Anti-spike IgG binding antibody, anti-receptor binding domain IgG antibody, and pseudovirus neutralizing antibody measurements four weeks post-vaccination were assessed as correlates of risk of moderate to severe-critical COVID-19 outcomes through 83 days post-vaccination and as correlates of protection following a single dose of Ad26.COV2.S COVID-19 vaccine in the placebo-controlled phase of ENSEMBLE, an international, randomized efficacy trial. Each marker had evidence as a correlate of risk and of protection, with strongest evidence for 50% inhibitory dilution (ID50) neutralizing antibody titer. The outcome hazard ratio was 0.49 (95% confidence interval 0.29, 0.81; p=0.006) per 10-fold increase in ID50; vaccine efficacy was 60% (43, 72%) at nonquantifiable ID50 (< 2.7 IU50/ml) and rose to 89% (78, 96%) at ID50 = 96.3 IU50/ml. Comparison of the vaccine efficacy by ID50 titer curves for ENSEMBLE-US, the COVE trial of the mRNA-1273 vaccine, and the COV002-UK trial of the AZD1222 vaccine supported consistency of the ID50 titer correlate of protection across trials and vaccine types.
RESUMEN
Several recently developed methods have the potential to harness machine learning in the pursuit of target quantities inspired by causal inference, including inverse weighting, doubly robust estimating equations and substitution estimators like targeted maximum likelihood estimation. There are even more recent augmentations of these procedures that can increase robustness, by adding a layer of cross-validation (cross-validated targeted maximum likelihood estimation and double machine learning, as applied to substitution and estimating equation approaches, respectively). While these methods have been evaluated individually on simulated and experimental data sets, a comprehensive analysis of their performance across real data based simulations have yet to be conducted. In this work, we benchmark multiple widely used methods for estimation of the average treatment effect using ten different nutrition intervention studies data. A nonparametric regression method, undersmoothed highly adaptive lasso, is used to generate the simulated distribution which preserves important features from the observed data and reproduces a set of true target parameters. For each simulated data, we apply the methods above to estimate the average treatment effects as well as their standard errors and resulting confidence intervals. Based on the analytic results, a general recommendation is put forth for use of the cross-validated variants of both substitution and estimating equation estimators. We conclude that the additional layer of cross-validation helps in avoiding unintentional over-fitting of nuisance parameter functionals and leads to more robust inferences.
Asunto(s)
Aprendizaje Automático , Proyectos de Investigación , Causalidad , Simulación por Computador , Humanos , Funciones de Verosimilitud , Modelos Estadísticos , Análisis de RegresiónRESUMEN
BACKGROUND: Hundreds of thousands of biodigesters have been constructed in Nepal. These household-level systems use human and animal waste to produce clean-burning biogas used for cooking, which can reduce household air pollution from woodburning cookstoves and prevent respiratory illnesses. The biodigesters, typically operated by female caregivers, require the handling of animal waste, which may increase domestic fecal contamination, exposure to diarrheal pathogens, and the risk of enteric infections, especially among young children. OBJECTIVE: We estimated the effect of daily reported biogas cookstove use on incident diarrhea among children <5y old in the Kavrepalanchok District of Nepal. Secondarily, we assessed effect measure modification and statistical interaction of individual- and household-level covariates (child sex, child age, birth order, exclusive breastfeeding, proof of vaccination, roof type, sanitation, drinking water treatment, food insecurity) as well as recent 14-d acute lower respiratory infection (ALRI) and season. METHODS: We analyzed 300,133 person-days for 539 children in an observational prospective cohort study to estimate the average effect of biogas stove use on incident diarrhea using cross-validated targeted maximum likelihood estimation (CV-TMLE). RESULTS: Households reported using biogas cookstoves in the past 3 d for 23% of observed person-days. The adjusted relative risk of diarrhea for children exposed to biogas cookstove use was 1.31 (95% confidence interval (CI): 1.00, 1.71) compared to unexposed children. The estimated effect of biogas stove use on diarrhea was stronger among breastfed children (2.09; 95% CI: 1.35, 3.25) than for nonbreastfed children and stronger during the dry season (2.03; 95% CI: 1.17, 3.53) than in the wet season. Among children exposed to biogas cookstove use, those with a recent ALRI had the highest mean risk of diarrhea, estimated at 4.53 events (95% CI: 1.03, 8.04) per 1,000 person-days. DISCUSSION: This analysis provides new evidence that child diarrhea may be an unintended health risk of biogas cookstove use. Additional studies are needed to identify exposure pathways of fecal pathogen contamination associated with biodigesters to improve the safety of these widely distributed public health interventions. https://doi.org/10.1289/EHP9468.