Oliveira Carlos R, Shapiro Eugene D, Weinberger Daniel M
Department of Pediatrics, Yale University School of Medicine, New Haven, CT, USA.
Department of Biostatistics, Yale University School of Public Health, New Haven, CT, USA.
Clin Epidemiol. 2022 Oct 18;14:1167-1175. doi: 10.2147/CLEP.S378039. eCollection 2022.
Vaccine effectiveness (VE) studies are often conducted after the introduction of new vaccines to ensure they provide protection in real-world settings. Control of confounding is often needed during the analyses, which is most efficiently done through multivariable modeling. When many confounders are being considered, it can be challenging to know which variables need to be included in the final model. We propose an intuitive Bayesian model averaging (BMA) framework for this task.
Data were used from a matched case-control study that aimed to assess the effectiveness of the Lyme vaccine post-licensure. Cases were residents of Connecticut, 15-70 years of age with confirmed Lyme disease. Up to 2 healthy controls were matched to each case subject by age. All participants were interviewed, and medical records were reviewed to ascertain immunization history and evaluate potential confounders. BMA was used to systematically search for potential models and calculate the weighted average VE estimate from the top subset of models. The performance of BMA was compared to three traditional single-best-model-selection methods: two-stage selection, stepwise elimination, and the leaps and bounds algorithm.
The analysis included 358 cases and 554 matched controls. VE ranged between 56% and 73% and 95% confidence intervals crossed zero in <5% of all candidate models. Averaging across the top 15 models, the BMA VE was 69% (95% CI: 18-88%). The two-stage, stepwise, and leaps and bounds algorithm yielded VE of 71% (95% CI: 21-90%), 73% (95% CI: 26-90%), and 74% (95% CI: 27-91%), respectively.
This paper highlights how the BMA framework can be used to generate transparent and robust estimates of VE. The BMA-derived VE and confidence intervals were similar to those estimated using traditional methods. However, by incorporating model uncertainty into the parameter estimation, BMA can lend additional rigor and credibility to a well-designed study.
新疫苗引入后通常会开展疫苗效力(VE)研究,以确保其在实际环境中提供保护。分析过程中常常需要控制混杂因素,而通过多变量建模能最有效地实现这一点。当考虑众多混杂因素时,确定最终模型中应纳入哪些变量可能具有挑战性。我们针对此任务提出了一个直观的贝叶斯模型平均(BMA)框架。
数据来自一项匹配病例对照研究,旨在评估莱姆病疫苗获批后的效力。病例为康涅狄格州15至70岁确诊莱姆病的居民。按年龄为每个病例对象匹配至多2名健康对照。对所有参与者进行访谈,并查阅病历以确定免疫史并评估潜在混杂因素。使用BMA系统搜索潜在模型,并从模型的顶级子集中计算加权平均VE估计值。将BMA的性能与三种传统的单最佳模型选择方法进行比较:两阶段选择、逐步排除和跳跃搜索算法。
分析纳入了358例病例和554名匹配对照。VE在56%至73%之间,所有候选模型中<5%的95%置信区间跨过零点。在前15个模型中进行平均,BMA的VE为69%(95%CI:18 - 88%)。两阶段、逐步和跳跃搜索算法得出的VE分别为71%(95%CI:21 - 90%)、73%(95%CI:26 - 90%)和74%(95%CI:27 - 91%)。
本文强调了BMA框架如何用于生成透明且稳健的VE估计值。BMA得出的VE和置信区间与使用传统方法估计的结果相似。然而,通过将模型不确定性纳入参数估计,BMA可为精心设计的研究增添额外的严谨性和可信度。