全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 数据分析与数据挖掘
20307 4
2020-02-06


2月5日凌晨更新了SEIRS模型。

2月1日凌晨,根据国家卫建委官网公布的数据实时计算的R0值徘徊在了3.09(Tg=8.4)和3.61附近!R0值第一次出现了下降!详见下图表。

实时更新2019-nCoV基本再生数(R0)值:

基本再生数(R0)的定义:在发病初期,当所有人均为易感者时,一个病人在其平均患病期内所传染的人数。

基于SEIR模型计算的R0值公式如下:



,

v2-cce1a1eba86db4a6a5b19af5d570676e_hd.jpg

R0值Python实现代码如下:


复制代码



今天是2020年1月27日。根据国家卫生健康委的数据,截至2020年1月26日24时,中国30个省(区、市)累积报告确诊病例2744例,重症病例461例,累计死亡病例80例,累计治愈出院51例,疑似病例5794例,累计追踪密切接触者32799人,当日解除医学观察583人,现有30453人正在接受医学观察。

随着1月23日武汉市宣布全面封城后,各地确诊病例不断出现,2020年的新年因为新型冠状病毒(2019-nCoV)的传播而黯然失色。

本文尝试使用经典传染病模型对2019-nCoV进行建模,预测未来疾病走势。因传染病模型研究属于传染病动力学研究方向,不是本人的工作范围,因此,本人只是将模型中的微分方程,用Python进行实现,想起到抛砖引玉的作用。


(具体各个模型的理论细节,请移步其他文章)

模型一:SI-Model


复制代码

模型二:SIS-Model


复制代码

模型三:SIR-Model


复制代码

模型四:SIRS-Model


复制代码

v2-d416aa3e66d8afce1b1ed5499dcd0d51_hd.jpg



附件列表
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2020-2-6 09:40:13

模型五:SEIR-Model


复制代码

模型六:SEIRS-Model


复制代码

关于2019-nCoV的SEIR模型参数设置不具备参考性,所以做了删除处理。后续时间可能会更改重新放在这里,也可能更新其他模型(如考虑常量输出输入的SEIR模型、SQIR模型等。。。)

因本人以前从未研究过传染病动力学,只是一时感兴趣,所以不免存在问题,欢迎大家指正~

原文链接:https://zhuanlan.zhihu.com/p/104 ... KPQZb8ddyR1vmpyY%3D
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2020-2-6 17:45:33
优秀
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2020-2-10 16:49:19
thanks for sharing
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2020-3-1 10:07:22
非常好,点赞
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群