作为一名码农,我们用数据说话,下面为大家分析目前的新型冠状病毒真的严重吗?
声明:
目的:本文探讨的主要目的是,有没有一种算法可以根据已有的数据对疫情进行预测和评估,可以方便在疫情出现的时候,根据预测的数据来判断当前的疫情状况。
文章分为如下五个部分:
PS : 心急的朋友可以直接拖到末尾看结论。
由于当前新型冠状病毒的数据并不多(从 11 号-26 号),所以仅仅是新型冠状病毒的数据来做研究是不够的,必须还要一组参考验证组。(当然是越多越好,但是精力有限,所以只找了一组)。
首先明确一下判断疫情的严重情况的指标,目前疫情一般会公告这四个数据:
由于我们把地区限定在了中国(虽然各个地区不一样),所这个地区的指标就不用了,死亡人数在疫情的前中后期都不太一样,这个也不用。
疑似病例这个和各地的收诊情况有关系,比如香港的疑似病例就非常多,确诊很少,所以这个指标也不采用,剩下就以确诊病例作为唯一的参考指标。
猜想 1:假设病毒的传播和地区,气候以及病毒自身的特征有非常大的关系。
那么只要找到一个和新型冠状病毒这三个属性非常相似的病毒,只要算法能够对这种病毒的确诊病例进行预估和判断,那么算法也可以对新型冠状病毒进行预估和判断。
目前和新型冠状病毒最像的应该就是 SARS:
那么基于上述的条件,采用了 SARS 作为参考组。
猜想 2:可以给定一个 Δ。
如果可以拿到 SARS 在 1 月份的数据就好了,可惜没有一个比较完整的地方有记录,基本上卫生部发布的好像是从四月份开始。于是只能采用四月份的数据了。
评估的算法:采用了 LSTM(长短期记忆网络),搭建的话 Keras。
LSTM 已经是个老算法了,经常用于量化投资(简单点来说就是炒股之类,预测股票的增长),非常擅长于时间序列的预测。
最终调教后的算法模型是,根据过去三天的数据,可以预估最新一天的数据。
实验的环境:数据量比较少,直接租了一台 16 核 8G 的阿里云的主机进行测试。
猜想验证部分直接全部使用 WHO 的 SARS 的数据进行验证。采用的字段是 Cumulative number of case(s),累计病例数,我这里就当它是确诊人数了。
SARS 的数据分别有两部分,一部分是中国大陆的,一部分是中国香港的。在 5 月份后基本已经趋于平缓,所以只取了前两个月(约 60 天左右的数据)。
可以看到数据涵盖了病毒爆发期,病毒平稳增长期,以及病毒被消灭期,比较合适用来验证算法是否可以拟合整体的曲线。
猜想 1:这里采用前 50 天的数据做训练样本,训练的结果用来预测 60 天的情况。
预测和实际的对比如下:
从图中我们可以看到,在使用了 90% 的数据量下,数据是拟合的比较好的,肉眼看上去香港的数据基本可以用算法完全预估。
这里给一下数据误差的对比情况:
说人话就是,模型预测的数据和实际的数据最多差了 6%。效果是非常理想的。这里我们可以证明第一个猜想,即 LSTM 可以用来拟合 SARS 确诊人数的变化情况。
猜想 2:实际上虽然猜想一得到了证实,但是依旧没啥用,因为这里的模型是使用了 90% 的数据,基本已经到了传播的后期。
那么前期是否可以使用同一个模型进行验证呢? 这里只取中国大陆的数据进行试验。
根据三个阶段(这三个阶段是肉眼判断截取的,是否合理当时没有考虑)截取对应的数据进行建模,建模后的预测效果如下:
在前期的数据的预测中,到后面两个阶段实际的数据都比预测的数据比要大,说明在疫情开始了爆发。
在中期的数据的预测中,其他两个阶段预测的数据基本都比实际的数据要大, 说明在疫情在前期的感染速度较低,后期的感染速度也变低 。
根据上面两个数据,是可以推测是存在一个 Δ 去判断当前的疫情状况的。但是 Δ 究竟是多少呢?
这里倒是可以根据好几个因素去设置:
这里就偷懒了,可以明显看到各个时期都有一个或者多个转折点(消亡期不算)。
直接取最接近转折点的误差率来计算:
分别取前 6 天的数据(到 20 号),取前面 9 天的数据(到 25 号)的数据代入模型进行预估:
为什么选 20 号呢,20 号主要是央视有采访钟南山大佬,而且第一次说明了人传人的存在,以及 14 名医护被感染。
我觉得应该算是开始重视这个病毒的一个点,随后各地开始上报疫情,发布以及响应等。
再往前的时间点就数据就比较缺失了,模型精准度已经失去了使用的价值,所以采取了 20 号的这个时间点作为分隔来建造模型。
如果以到 20 号的模型来估算的话,第一次误差偏离 0.17 的时间是 23 号,误差为 0.19。23 号,广东,浙江启动了一级响应,确诊的病人增多。
但个人感觉应该是检测的方法变强导致的,侧面来看,模型的确检测除了异常值,不过并不能够说明病毒变强了。
如果以到 24 号的模型来看,预估的 25 号的数据是 1877,实际的数据 1972,误差率在 0.04,未到偏离值。
所以根据这个模型推断,当前疫情还在持续扩散中,但是并没有进一步恶化,可能国家的高度重视和很多省份都开始参与其中有关系。
依照当前的模型,明天的确诊人数的预估值是 2770,加上 0.17 的误差的话:
再次声明:以上的观点和结论都属于依据当前的确诊人数所作的推断,毫无参考的意义和价值。如果数据刚好对了,不能说明什么,如果数据错了,证明这个模型真的不行。
那么,如果我们带入 SARS 最严重的的时候,会怎么样?
还是按照误差 0.17 来算,在 23 号,误差达到了 0.18,也就是 3 天前,按照这个模型的数据推算,新型冠状的病毒传染性比 SARS 要强。 但很可能是和春运有关系,毕竟数据的时间不一致。
再再次声明:以上的观点和结论都属于依据当前的确诊人数所作的推断,毫无参考的意义和价值。如果数据刚好对了,不能说明什么,如果数据错了,证明这个模型真的不行。
写这个的目的纯粹是为了探讨和研究,实际上,真的要写一个这种疫情预测和评估的系统,需要非常多的参考的数值。
比如医护人员的数量,检测手法准确率,时间,病毒的潜伏期等等,需要非常多的人力物力才可以编写完成,像本文这种只靠单一维度的输入去判断和预估的模型和玩具没什么区别。
但是在编写的过程中也发现很多问题,例如政府之前的披露的数据真的非常少,疫情的相关数据还是多方寻找才找到一部分 。
希望在这次的疫情中保留比较详细的数据:
这些都可以作为系统的的输入,借助这些数据,打造一个强大完整的疫情监控和评估系统,甚至不是疫情,而是重大事件的监控和评估系统,待下次再出现疫情、灾情,可以快速的调度全国的资源进行处理和援助。
最后,祝大家春节快乐,多注意个人卫生,带好口罩,不要放松警惕。
建站 什么 虚拟主机 够用?这要看搭建的是什么类型的网站。比如个人博客类型的网...
信息化2.0时代提出开展智慧教育创新发展行动。2019年2月,中共中央、国务院印发...
摘要 元旦期间 订单业务线 告知 推送系统 无法正常收发消息,作为推送系统维护者...
在Python语言中有如下3种方法: 成员方法 类方法(classmethod) 静态方法(staticm...
从 10.0.0 版开始,异步迭代器就出现在 Node 中了,在本文中,我们将讨论异步迭...
Docker生成新镜像版本的两种方式 There are two ways Docker can generate new m...
前提条件 请您在购买前确保已完成注册和充值。详细操作请参见 如何注册公有云管...
【51CTO.com快译】 数据可视化工具不断发展,提供更强大的功能,同时改善可访问...
本文整理自直播《Hologres 数据导入/导出实践-王华峰(继儒)》 视频链接: https:/...
2021年3月24日,主题为《数据的世界,世界的数据》的星环科技2021春季新品发布会...