Department of Biostatistics, Boston University School of Public Health, Boston, Massachusetts, United States of America.
Department of Mathematics & Statistics, Boston University, Boston, Massachusetts, United States of America.
PLoS Comput Biol. 2022 Sep 1;18(9):e1010434. doi: 10.1371/journal.pcbi.1010434. eCollection 2022 Sep.
The reproductive number is an important metric that has been widely used to quantify the infectiousness of communicable diseases. The time-varying instantaneous reproductive number is useful for monitoring the real-time dynamics of a disease to inform policy making for disease control. Local estimation of this metric, for instance at a county or city level, allows for more targeted interventions to curb transmission. However, simultaneous estimation of local reproductive numbers must account for potential sources of heterogeneity in these time-varying quantities-a key element of which is human mobility. We develop a statistical method that incorporates human mobility between multiple regions for estimating region-specific instantaneous reproductive numbers. The model also can account for exogenous cases imported from outside of the regions of interest. We propose two approaches to estimate the reproductive numbers, with mobility data used to adjust incidence in the first approach and to inform a formal priori distribution in the second (Bayesian) approach. Through a simulation study, we show that region-specific reproductive numbers can be well estimated if human mobility is reasonably well approximated by available data. We use this approach to estimate the instantaneous reproductive numbers of COVID-19 for 14 counties in Massachusetts using CDC case report data and the human mobility data collected by SafeGraph. We found that, accounting for mobility, our method produces estimates of reproductive numbers that are distinct across counties. In contrast, independent estimation of county-level reproductive numbers tends to produce similar values, as trends in county case-counts for the state are fairly concordant. These approaches can also be used to estimate any heterogeneity in transmission, for instance, age-dependent instantaneous reproductive number estimates. As people are more mobile and interact frequently in ways that permit transmission, it is important to account for this in the estimation of the reproductive number.
繁殖数是一个重要的指标,已被广泛用于量化传染病的传染性。时变瞬时繁殖数可用于监测疾病的实时动态,为疾病控制提供决策依据。例如在县或市一级对该指标进行局部估计,可进行更有针对性的干预,以遏制传播。然而,对这些时变数量的局部繁殖数进行同时估计,必须考虑到这些数量潜在的异质性来源,其中一个关键因素是人类流动性。我们开发了一种统计方法,该方法将多个区域之间的人类流动纳入其中,用于估计特定区域的瞬时繁殖数。该模型还可以考虑从感兴趣区域以外输入的外生病例。我们提出了两种估计繁殖数的方法,第一种方法使用流动数据来调整发病率,第二种方法(贝叶斯方法)则利用流动数据来告知先验分布。通过模拟研究,我们表明,如果人类流动可以通过现有数据合理地近似,那么可以很好地估计特定区域的繁殖数。我们使用这种方法,利用疾病预防控制中心的病例报告数据和 SafeGraph 收集的人类流动数据,对马萨诸塞州的 14 个县的 COVID-19 瞬时繁殖数进行了估计。我们发现,考虑到流动性,我们的方法产生的繁殖数估计值在各个县之间存在差异。相比之下,独立估计县一级的繁殖数往往会产生相似的值,因为该州各县的病例计数趋势相当一致。这些方法还可用于估计任何传播的异质性,例如,年龄相关的瞬时繁殖数估计。由于人们的流动性更大,并且经常以允许传播的方式相互作用,因此在估计繁殖数时必须考虑到这一点。