Python实现新型冠状病毒传播模型及预测代码实例实现新型冠状病毒传播模型及预测代码实例
在本篇文章里小编给大家整理的是关于Python实现新型冠状病毒传播模型及预测代码内容,有兴趣的朋友们可
以学习下。
1.传染及发病过程传染及发病过程
一个健康人感染病毒后进入潜伏期(时间长度为Q天),潜伏期之后进入发病期(时间长度为D天),发病期之后该患者有三
个可能去向,分别是自愈、接收隔离、死亡。
2.模型假设模型假设
潜伏期Q=7天,根据报道潜伏期为2~14天,取中间值;发病期D=10天,根据文献报告,WHO认定SARS发病期为10天,假设
武汉肺炎与此相同;潜伏期的患者不具有将病毒传染给他人的能力;发病期的患者具有将病毒传染给他人的能力;患者在发病
期之后不再具有将病毒传染他人的能力;假设处于发病期的患者平均每天密切接触1人,致使该人患病的概率为γ最初只有一
个人类感染者;病情自然发展,没有外部干扰。
3.模型公式模型公式
:人类感染该病毒的天数
: 第N天感染该病毒并且处于发病期的患者数量
: Q天前新被感染患者,当日进入发病期的数量
: 当日发病期满,不再具有传染能力的患者数量
4.模型初始值模型初始值
根据假设,最初只有一个人类感染者,所以:
5.实际疫情数据实际疫情数据
人类感染病毒且发病的初始日期人类感染病毒且发病的初始日期:根据财新网的报道,官方通报首例不明原因肺炎是在12月8日,考虑到确诊之前肯定已经尝
试过各种治疗方案无效后认定为不明原因肺炎,所以有理由认为该名患者在12月8日已经处于发病期末端,根据假设发病期为
10天,所以可以假设该名患者在11月29号发病,即N=1对应11月29日。
日新增发病数日新增发病数
近期疫情防控大事记:
从以上信息可以判断核酸检测试剂是在1月16日、17日大幅使用的,18日、19日确诊大量病例,因此1月20日之前的确诊病例
数对模型参考意义不大。1月20日之后,可以认为新发病例,发病即检测。
截至1月20日24时,国家卫健委公告累计确诊病例291
截至1月21日24时,国家卫健委公告累计确诊病例440
截至1月22日24时,国家卫健委公告累计确诊病例571
截至1月23日24时,国家卫健委公告累计确诊病例830
截至1月24日24时,国家卫健委公告累计确诊病例1287
所以:
1月21日新增确诊病例:440-291=1491月22日新增确诊病例:571-440=1311月23日新增确诊病例:830-571=2591月24日新
评论0