2025-11-12T07:49:15.324568

Poisson Network SIR Epidemic Model

Wairimu, Gothard, Rempala

We extend the classical Susceptible-Infected-Recovered (SIR) model to a network-based framework where the degree distribution of nodes follows a Poisson distribution. This extension incorporates an additional parameter representing the mean node degree, allowing for the inclusion of heterogeneity in contact patterns. Using this enhanced model, we analyze epidemic data from the 2018-20 Ebola outbreak in the Democratic Republic of the Congo, employing a survival approach combined with the Hamiltonian Monte Carlo method. Our results suggest that network-based models can more effectively capture the heterogeneity of epidemic dynamics compared to traditional compartmental models, without introducing unduly overcomplicated compartmental framework.

academic

Poisson Network SIR Epidemic Model

基本信息

论文ID: 2501.00187
标题: Poisson Network SIR Epidemic Model
作者: Josephine K. Wairimu (University of Nairobi), Andrew Gothard (The Ohio State University), Grzegorz A. Rempala (The Ohio State University)
分类: q-bio.PE (Quantitative Biology - Populations and Evolution), math.DS (Mathematics - Dynamical Systems), physics.soc-ph (Physics - Physics and Society), stat.AP (Statistics - Applications)
发表时间: 2024年12月30日 (arXiv提交)
论文链接: https://arxiv.org/abs/2501.00187

摘要

本文将经典的易感-感染-康复(SIR)模型扩展到基于网络的框架，其中节点的度分布遵循泊松分布。这种扩展引入了一个表示平均节点度的额外参数，允许纳入接触模式的异质性。使用这个增强模型，作者分析了2018-20年刚果民主共和国埃博拉疫情的流行病数据，采用生存分析方法结合哈密顿蒙特卡洛方法。结果表明，基于网络的模型能够比传统的舱室模型更有效地捕捉流行病动力学的异质性，且不会引入过度复杂的舱室框架。

研究背景与动机

问题定义

传统的SIR舱室模型虽然在流行病学建模中广泛应用，但存在几个重要局限性：

质量作用定律假设：假设人群中的社会接触模式是均匀的，忽略了个体行为的内在异质性
确定性模型的局限：忽略了疫情早期阶段的随机事件，此时感染人数较少，随机性对疾病传播起重要作用
复杂性增长问题：当涉及多种类型的感染者或易感者时，模型变得越来越复杂，难以解释和分析

研究重要性

人际传播是许多传染病的主要传播方式，准确建模接触网络对于：

预测疾病动力学
评估各种干预措施的影响
指导公共卫生策略制定至关重要。

研究动机

本文旨在通过引入基于网络的随机SIR模型来解决上述问题，特别关注泊松度分布网络，因为它具有精确的配对级闭包性质，使模型在保持数学可处理性的同时能够捕捉接触模式的异质性。

核心贡献

理论贡献：建立了经典SIR模型与泊松网络SIR模型之间的数学等价关系，证明了在特定条件下两者的近似性
方法创新：提出了基于动态生存分析(DSA)的参数估计方法，结合哈密顿蒙特卡洛(HMC)采样
实际应用：成功应用于2018-20年刚果民主共和国埃博拉疫情数据分析，展示了网络模型的实用价值
数据处理：开发了处理不完整流行病数据的似然函数方法，无需数据插补即可处理缺失的感染或康复时间

方法详解

任务定义

构建一个能够在保持数学可处理性的同时捕捉接触网络异质性的流行病传播模型，并开发相应的统计推断方法来估计模型参数和预测疫情发展。

模型架构

1. 网络SIR过程

在大小为n的网络上，SIR过程描述为：

初始时随机选择m个个体作为感染者
每个感染个体的感染期服从参数为γ的指数分布
在感染期内，个体以强度β的泊松过程接触其邻居
如果被接触的邻居是易感者，则立即被感染
感染期结束后，个体康复并获得免疫

2. 配对模型及其闭包

配对模型方程组为：

Ṡ = -β[SI]
İ = β[SI] - γ[I]  
Ṙ = γ[I]
[SI]˙ = -γ[SI] + β([SSI] - [ISI]) - β[SI]
[SS]˙ = -2β[SSI]

其中A、AB、ABC分别表示网络中单个、配对和三元组的计数。

3. 泊松网络的精确闭包

对于泊松度分布pk = exp(-μ)μk/k!，系统可以精确闭包：

[ASI] ≃ κ[AS][SI]/[S]

其中κ = 1（泊松分布的特殊性质）。

4. 简化的泊松网络SIR方程

最终得到的简化方程组：

Ṡ = -β̃S(1 + ρ - S + R̃₀⁻¹logS)
İ = -Ṡ - γI
R = 1 + ρ - S - I

其中R̃₀ = β̃/γ̃是网络基本再生数。

技术创新点

1. 精确配对闭包

利用泊松分布的特殊性质（κ = 1），实现了网络模型的精确配对级闭包，避免了高阶项的复杂性。

2. 动态生存分析框架

将网络SIR过程重新解释为生存分析问题：

S(t)作为不恰当的生存函数
感染时间TI的条件密度：fτT(t) = -Ṡ(t)/τT
康复时间的密度通过卷积获得

3. 缺失数据处理

开发了处理不完整数据的似然函数：

完整数据(ti, ri)：L(θ|ti, ri, T) = fτT(ti)γʷⁱe⁻ᵞ⁽ʳⁱ∧ᵀ⁻ᵗⁱ⁾
仅有感染时间：L(θ|ti, ◦, T) = fτT(ti)
仅有康复时间：L(θ|◦, ri, T) = gτT(ri)

实验设置

数据集

2018-20年刚果民主共和国埃博拉疫情数据：

时间范围：2019年5月27日至9月12日
样本数量：1,069名确诊EVD患者
数据内容：症状发作时间和住院时间
数据完整性：约30%的病例记录不完整，缺少感染或康复时间

评价指标

后验参数估计的均值和95%置信区间
基本再生数R̃₀
最终疫情规模K̂∞
模型拟合度（通过最大化似然函数评估）

对比方法

经典SIR模型（基于DSA方法的分析结果，来自文献11）

实现细节

采样方法：哈密顿蒙特卡洛(HMC)
软件工具：STAN统计软件，通过R的Rstan库调用
MCMC设置：2条链，每条3000次迭代，1000次预热
先验分布：
- β̃ ∈ (0.1, ∞)
- γ̃ ∈ (0, β̃)
- ρ ∈ (0, 0.01)
收敛诊断：使用Rubin's R统计量

实验结果

主要结果

参数估计对比（后验均值和95%置信区间）：

参数	经典SIR	网络SIR
β̃	0.235 (0.218, 0.253)	0.229 (0.209, 0.259)
γ̃	0.214 (0.199, 0.230)	0.215 (0.197, 0.242)
ρ	0.0067 (0.0055, 0.0081)	0.0055 (0.0046, 0.0073)
μ	-	39.48 (7.93, 93.00)
R̃₀	1.098 (1.061, 1.135)	1.071 (1.034, 1.109)
K̂∞	3481.41 (2877.416, 4155.878)	3773.37 (3373.245, 4226.315)