Deterministic compartmental models are predominantly used in the modeling of infectious diseases, though stochastic models are considered more realistic, yet are complicated to estimate due to missing data. In this paper we present a novel algorithm for estimating the stochastic SIR/SEIR epidemic model within a Bayesian framework, which can be readily extended to more complex stochastic compartmental models. Specifically, based on the infinitesimal conditional independence properties of the model, we are able to find a proposal distribution for a Metropolis algorithm which is very close to the correct posterior distribution. As a consequence, rather than perform a Metropolis step updating one missing data point at a time, as in the current benchmark Markov chain Monte Carlo (MCMC) algorithm, we are able to extend our proposal to the entire set of missing observations. This improves the MCMC methods dramatically and makes the stochastic models now a viable modeling option. A number of real data illustrations and the necessary mathematical theory supporting our results are presented.
翻译:传染性疾病模型中主要使用决定性的区划模型,尽管人们认为随机模型比较现实,但由于缺少数据而难以估计。在本文中,我们提出了一个新的算法,用于在巴伊西亚框架内估计随机SIR/SEIR流行病模型,这个算法可以很容易地推广到更复杂的随机区划模型。具体地说,根据模型的无限的有条件独立特性,我们能够找到一个大都会算法的分布建议,该算法非常接近正确的后院分布。因此,我们不象目前基准的Markov链Monte Carlo(MCMC)算法那样,一次更新一个缺失的数据点,而是进行大都会步骤,我们能够将我们的建议扩大到全部缺失的观察。这极大地改进了MCMC方法,使这些随机模型现在成为一个可行的模型。我们提出了一些真实的数据说明和必要的数学理论,以支持我们的结果。