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 相关文章推荐
使用grappelli为django admin后台添加模板
Nov 18 Python
Python中用于检查英文字母大写的isupper()方法
May 19 Python
Python中Iterator迭代器的使用杂谈
Jun 20 Python
Python探索之ModelForm代码详解
Oct 26 Python
python机器学习实战之最近邻kNN分类器
Dec 20 Python
TensorFlow 实战之实现卷积神经网络的实例讲解
Feb 26 Python
基于python OpenCV实现动态人脸检测
May 25 Python
python中多个装饰器的调用顺序详解
Jul 16 Python
python多进程重复加载的解决方式
Dec 13 Python
python解析多层json操作示例
Dec 30 Python
Python基于execjs运行js过程解析
Nov 27 Python
python 指定源路径来解决import问题的操作
Mar 04 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
php学习 函数 课件
2008/06/15 PHP
PHP文件上传之多文件上传的实现思路
2016/01/27 PHP
Laravel中七个非常有用但很少人知道的Carbon方法
2017/09/21 PHP
Laravel 登录后清空COOKIE的操作方法
2019/10/14 PHP
关于火狐(firefox)及ie下event获取的两种方法
2012/12/27 Javascript
nodejs开发环境配置与使用
2014/11/17 NodeJs
基于javascript实现九九乘法表
2016/03/27 Javascript
浅谈bootstrap源码分析之tab(选项卡)
2016/06/06 Javascript
vue router路由嵌套不显示问题的解决方法
2017/06/17 Javascript
Iphone手机、安卓手机浏览器控制默认缩放大小的方法总结(附代码)
2017/08/18 Javascript
JS 实现百度搜索功能
2018/02/01 Javascript
Taro集成Redux快速上手的方法示例
2018/06/21 Javascript
微信小程序去除左上角返回键的实现方法
2020/03/06 Javascript
Python批量按比例缩小图片脚本分享
2015/05/21 Python
Django 忘记管理员或忘记管理员密码 重设登录密码的方法
2018/05/30 Python
使用批处理脚本自动生成并上传NuGet包(操作方法)
2019/11/19 Python
python实现批量处理将图片粘贴到另一张图片上并保存
2019/12/12 Python
tensorflow如何继续训练之前保存的模型实例
2020/01/21 Python
Python importlib动态导入模块实现代码
2020/04/16 Python
tensorflow实现将ckpt转pb文件的方法
2020/04/22 Python
vscode配置anaconda3的方法步骤
2020/08/08 Python
Python requests上传文件实现步骤
2020/09/15 Python
open_basedir restriction in effect. 原因与解决方法
2021/03/14 PHP
YOOX台湾:意大利奢侈品电商
2018/10/13 全球购物
改变生活的男士内衣:SAXX Underwear
2019/08/28 全球购物
护士专业推荐信
2013/11/02 职场文书
学院书画协会部门职责
2013/11/28 职场文书
弘扬雷锋精神演讲稿
2014/05/10 职场文书
奥巴马胜选演讲稿
2014/05/15 职场文书
大学生自我评价范文
2015/03/03 职场文书
Python中Cookies导出某站用户数据的方法
2021/05/17 Python
Java内存模型之happens-before概念详解
2021/06/13 Java/Android
使用canvas仿Echarts实现金字塔图的实例代码
2021/11/11 HTML / CSS
win10双系统怎么删除一个系统?win10电脑有两个系统删除一个的操作方法
2022/07/15 数码科技
Nginx使用ngx_http_upstream_module实现负载均衡功能示例
2022/08/05 Servers
Windows server 2016服务器基本设置
2022/08/14 Servers