Re: [讨论] <最终确诊数> 简易推估

楼主: perceval (摸鱼中)   2020-03-16 23:24:48
※ 引述《hugoyo (阿佑)》之铭言:
: 标题: [讨论] <最终确诊数> 简易推估
: 时间: Mon Mar 16 21:04:52 2020
:
: 前言:
:
: 首先,确诊数≠感染人数,检测病毒本来就受限于医疗、政治等等因素。 本方法只是纯数
: 学+各种简化,大概拿来预测一下最终大概会有多少确诊。
: 补充:
: 确 诊 数 不是感染数!!! 要看该国要不要验,要不要公布,以下只有纯数字分析。
: 讨论:
: 一个足够多确诊的群体的确诊人数,多会呈现所谓的S型曲线,一开始少少的,中间突然暴
: 增,最后趋于缓和。 数学上可以用error function来模拟。
:
: error function: https://i.imgur.com/xg1NGiC.png
不是公卫专长,不过这种族群成长模型比较常看到用logistic function描述
https://en.wikipedia.org/wiki/Logistic_function
一时找不到算法fit参数,所以没办法比较
不过估算出的最大确诊人数可能会差不多
但是目前欧洲国家其实还没有办法用这个模型去描述
应该都还没有达到反曲点(inflection point 中国称为拐点)
现在应该都还在指数型成长的阶段
可以参考worldmeters网站上的资料,选定感兴趣的国家(假设是意大利)
找到 Total Coronavirus Cases in Italy
https://www.worldometers.info/coronavirus/country/italy/
图片上选取logarithmic,你会发现确诊数取log后和时间是线性的关系
这表示目前还是指数形成长
有兴趣的可以用wolfram alpha去估算
exp fit{}指令去fit参数
假设我取最近五笔资料 exp fit{12462,15113,17660,21157,24747}
得到的模型是 10662.4 e^(0.169288 x)
估计今天总确诊数会是 29443
不过受限于每日检验量,应该没办法达到这个数字
但是其它欧洲国家(西班牙、德国、法国)应该都可以得到不错的估计结果

Links booklink

Contact Us: admin [ a t ] ucptt.com