Python实现新型冠状病毒传播模型及预测代码实例


Posted in Python onFebruary 05, 2020

1.传染及发病过程

一个健康人感染病毒后进入潜伏期(时间长度为Q天),潜伏期之后进入发病期(时间长度为D天),发病期之后该患者有三个可能去向,分别是自愈、接收隔离、死亡。

2.模型假设

潜伏期Q=7天,根据报道潜伏期为2~14天,取中间值;发病期D=10天,根据文献报告,WHO认定SARS发病期为10天,假设武汉肺炎与此相同;潜伏期的患者不具有将病毒传染给他人的能力;发病期的患者具有将病毒传染给他人的能力;患者在发病期之后不再具有将病毒传染他人的能力;假设处于发病期的患者平均每天密切接触1人,致使该人患病的概率为γ最初只有一个人类感染者;病情自然发展,没有外部干扰。

3.模型公式

Python实现新型冠状病毒传播模型及预测代码实例

Python实现新型冠状病毒传播模型及预测代码实例:人类感染该病毒的天数

Python实现新型冠状病毒传播模型及预测代码实例: 第N天感染该病毒并且处于发病期的患者数量

Python实现新型冠状病毒传播模型及预测代码实例 : Q天前新被感染患者,当日进入发病期的数量

Python实现新型冠状病毒传播模型及预测代码实例: 当日发病期满,不再具有传染能力的患者数量

4.模型初始值

根据假设,最初只有一个人类感染者,所以:

Python实现新型冠状病毒传播模型及预测代码实例

5.实际疫情数据

人类感染病毒且发病的初始日期:根据财新网的报道,官方通报首例不明原因肺炎是在12月8日,考虑到确诊之前肯定已经尝试过各种治疗方案无效后认定为不明原因肺炎,所以有理由认为该名患者在12月8日已经处于发病期末端,根据假设发病期为10天,所以可以假设该名患者在11月29号发病,即N=1对应11月29日。

日新增发病数 Python实现新型冠状病毒传播模型及预测代码实例

近期疫情防控大事记:

Python实现新型冠状病毒传播模型及预测代码实例

从以上信息可以判断核酸检测试剂是在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日新增确诊病例:1287-830=457

考虑到:

1月20日之后发病即检测确诊检测用时2日国外新增病例在个位数,且不能保证发病即检测等因素,暂不考虑新增病例应该递增,所以1月22日新增数据异常,舍去

使用上溯每日新增数据,同时考虑到确诊需要2天, 可以得到:

Python实现新型冠状病毒传播模型及预测代码实例, Python实现新型冠状病毒传播模型及预测代码实例, Python实现新型冠状病毒传播模型及预测代码实例

Python实现新型冠状病毒传播模型及预测代码实例 对应11月29日 )

6.拟合确定 Python实现新型冠状病毒传播模型及预测代码实例

根据近期每日新增数据、模型初始值及模型公式,用最小二乘拟合得到 Python实现新型冠状病毒传播模型及预测代码实例

Python实现新型冠状病毒传播模型及预测代码实例

7.预测

患者数量:根据上文确定的模型及参数,从11月29日(N=1)至1月27日(N=60)人群中累计处于发病期的人数如下图所示:

Python实现新型冠状病毒传播模型及预测代码实例

根据模型,近期人群中患者数量计算如下:

Python实现新型冠状病毒传播模型及预测代码实例

人群中感染了病毒并处于发病期的患者数量,注意Pn一般需要延期2日才能确诊

每日新增患者数量:

根据模型,近期人群中每日新增患者数量计算如下:

Python实现新型冠状病毒传播模型及预测代码实例

注意:图中是人群中新增发病患者数量,可与晚2日的政府发布新增数量进行对比。即22日新增患者数量可与24日政府发布的新增病例进行对比。截至目前模型计算22日新增为369人,政府公布的24日新增病例457人

根据以上模型预计未来几日的情况如下:

1月25日将新确诊432例,人群中发病患者为4068人;

1月26日将新确诊505例,人群中发病患者为4759人;

1月27日将新确诊590例,人群中发病患者为5568人;

1月28日将新确诊691例,人群中发病患者为6514人;

1月29日将新确诊809例,人群中发病患者为7621人。

由于1月20日之后采取的各种措施,将导致发病期D下降,感染概率 Python实现新型冠状病毒传播模型及预测代码实例 下降,1月29日之后日均增长势头会减弱。

8.模型代码

import numpy as np
import matplotlib.pyplot as plt
 
gamma = 0.55
Q = 7
D = 10
P = np.zeros(300, dtype=np.float)
Psum = np.zeros(300, dtype=np.float)
for i in range(Q):
  P[i] = 1
for j in range(300-Q):
  P[j+Q] = P[j+Q-1]+P[j]*gamma
  if j+Q-D-Q >= 0:
    P[j+Q] -= P[j+Q-D-Q]*gamma
  if j+Q == D:
    P[j+Q] -= 1
plt.xlabel("N")
plt.ylabel("PN")
plt.rcParams['font.sans-serif'] = ['SimHei']
plt.rcParams['axes.unicode_minus'] = False
plt.plot(range(1, 61), P[0:60])
plt.grid()
plt.show()

以上就是三水点靠木小编给大家整理的全部内容,如果大家有任何补充可以联系小编。

Python 相关文章推荐
python读取Android permission文件
Nov 01 Python
pymongo实现控制mongodb中数字字段做加法的方法
Mar 26 Python
设置python3为默认python的方法
Oct 31 Python
Python实现的旋转数组功能算法示例
Feb 23 Python
Python生成器的使用方法和示例代码
Mar 04 Python
python pandas模块基础学习详解
Jul 03 Python
python实现在函数图像上添加文字和标注的方法
Jul 08 Python
如何在keras中添加自己的优化器(如adam等)
Jun 19 Python
完美解决TensorFlow和Keras大数据量内存溢出的问题
Jul 03 Python
Python3爬虫关于代理池的维护详解
Jul 30 Python
Python Selenium操作Cookie的实例方法
Feb 28 Python
Python Pytorch查询图像的特征从集合或数据库中查找图像
Apr 09 Python
基于Tensorflow批量数据的输入实现方式
Feb 05 #Python
Python操作注册表详细步骤介绍
Feb 05 #Python
Python类继承和多态原理解析
Feb 05 #Python
Python模块 _winreg操作注册表
Feb 05 #Python
python3操作注册表的方法(Url protocol)
Feb 05 #Python
Python tkinter模版代码实例
Feb 05 #Python
Python Scrapy框架第一个入门程序示例
Feb 05 #Python
You might like
浅谈Laravel模板实体转义带来的坑
2019/10/22 PHP
读jQuery之八 包装事件对象
2011/06/21 Javascript
Javascript图像处理—为矩阵添加常用方法
2012/12/27 Javascript
纯JavaScript实现的兼容各浏览器的添加和移除事件封装
2015/03/28 Javascript
jquery获得当前html页面源码的方法
2015/07/14 Javascript
js实现的奥运倒计时时钟效果代码
2015/12/09 Javascript
js判断手机访问或者PC的几个例子(常用于手机跳转)
2015/12/15 Javascript
JS组件中bootstrap multiselect两大组件较量
2016/01/26 Javascript
Webwork 实现文件上传下载代码详解
2016/02/02 Javascript
JavaScript判断数组是否存在key的简单实例
2016/08/03 Javascript
xmlplus组件设计系列之下拉刷新(PullRefresh)(6)
2017/05/03 Javascript
强大的JavaScript响应式图表Chartist.js的使用
2017/09/13 Javascript
node实现的爬虫功能示例
2018/05/04 Javascript
原生JS实现的轮播图功能详解
2018/08/06 Javascript
vue项目中跳转到外部链接的实例讲解
2018/09/20 Javascript
在vue-cli3中使用axios获取本地json操作
2020/07/30 Javascript
[05:02]2014DOTA2 TI中国区预选赛精彩TOPPLAY第三弹
2014/06/25 DOTA
Python中getattr函数和hasattr函数作用详解
2016/06/14 Python
python机器学习库常用汇总
2017/11/15 Python
python 哈希表实现简单python字典代码实例
2019/09/27 Python
Python-openCV读RGB通道图实例
2020/01/17 Python
Python使用uuid库生成唯一标识ID
2020/02/12 Python
Html5调用手机摄像头并实现人脸识别的实现
2018/12/21 HTML / CSS
Kenneth Cole官网:纽约时尚优雅品牌
2016/11/14 全球购物
国际贸易专业推荐信
2013/11/15 职场文书
模具设计与制造专业推荐信
2014/02/16 职场文书
宿舍违规用电检讨书
2014/02/16 职场文书
终止合同协议书
2014/04/17 职场文书
幼儿园教师演讲稿
2014/05/06 职场文书
关于爱国的演讲稿
2014/05/07 职场文书
国际贸易专业求职信
2014/06/04 职场文书
小学竞选班干部演讲稿
2014/08/20 职场文书
使用qt quick-ListView仿微信好友列表和聊天列表的示例代码
2021/06/13 Python
MySQL系列之十一 日志记录
2021/07/02 MySQL
Nginx反向代理学习实例教程
2021/10/24 Servers
如何使用注解方式实现 Redis 分布式锁
2022/07/23 Redis