图像缩略图

fitVirus

2.1.2版(3.47 MB​​)通过 米兰巴蒂斯塔
由Logistic模型冠状COVID-19疫情量估计

203下载

更新2020年4月14日

查看许可证

编者按:这个文件被选为MATLAB中心本周的挑选

该功能fitVirus03实现了从日常预测epidemy最终尺寸的估计的逻辑模型。该模型是数据驱动的,所以它的预测是一样好数据。此外,假定模型是一个阶段流行的一个合理的说明。但是,如果疫情发展到第二阶段的模型变得无用。该模型也没用初始流行阶段。

的贡献包含冠状奥地利,比利时,中国,克罗地亚,丹麦,德国,匈牙利,法国,伊朗,意大利,伦巴第,挪威,荷兰,纽约州,葡萄牙,斯洛文尼亚,韩国,西班牙,瑞士,英国,美国数据和数据在中国以外(最多24.Mar.2020)

回归收敛失败可能纯的初始猜测或小的数据集。因此,该方法并不适用于流行的早期阶段。此外,结果是无用的,如果回归统计不符合最低标准,说R ^ 2> 0.8,p值<0.05。

在流行病学评价图上区域颜色分开流行病学阶段(这些不是标准,但为了方便随意选择):
红 - 快速增长阶段
黄色 - 过渡到稳态阶段
绿色 - 收尾阶段(稳定期)

产生的第二数字是每日epidemy尺寸的评价。如果这些值不收敛于一个常数,然后疫情可能是尚未稳定。

更详细的描述可以在
https://www.researchgate.net/publication/339240777_Estimation_of_the_final_size_of_coronavirus_epidemic_by_the_logistic_model
例子中可以找到
https://www.researchgate.net/publication/339912313_Forecasting_of_final_COVID-19_epidemic_size_200318

基于SIR模型的新版本可在
//www.tianjin-qmedu.com/matlabcentral/fileexchange/74658-fitviruscovid19

对于其他国家的数据都可以从
https://ourworldindata.org/coronavirus-source-data

免责声明。软件和数据是用于教育,而不是用于医疗或商业用途。

引用作为

米兰巴蒂斯塔(2020)。fitVirus(//www.tianjin-qmedu.com/matlabcentral/fileexchange/74411-fitvirus),MATLAB中央文件交换。检索

意见及评级(43

您可以添加科索沃数据?

非常感谢您更新自己的程序和固定的失败。

罗伯托,你已经安装了统计工具箱?

您好,我有一个问题想运行fitVirus03功能和Matlab的礼物我这个错误:

> > fitVirus03 (@getDataGermany);
****估计德国疫情规模
初始猜测K = 126123 R = 0.287012 A = 18645
使用错误optimoptions(线105)
“SpecifyObjectiveGradient”不是LSQCURVEFIT的选项。
可以在lsqefit文档页面上找到选项列表。

fitVirus03中的错误(第50行)
选择= optimoptions (“lsqcurvefit”,“显示”,“关闭”,…

他们可以帮助我知道我做错了。
谢谢你!

谢谢!

请更新此代码与图形从fitVirusCOVID19。

谢谢!

卢马克

贡纳尔

很好,但是如果modell也能在GNU Octave中运行的话,它会有更大的影响(也就是说,optimoptions和嵌套函数需要兼容的版本)。

有没有人试图将Matlab代码移植到其他解决方案,比如GNU Octave?

我将巴西的数据发给了do Mr. Batista,以便他可以添加到报告中。

对所有。SIR模型改进了收敛性和初始猜测计算。我认为它比logistic模型更有效,但是更健壮。
//www.tianjin-qmedu.com/matlabcentral/fileexchange/74658-fitviruscovid19

亲爱的克劳迪奥,感谢您的建议。请记住,物流模式是非常简单的。每日预测可以是非常不错的。我对斯洛文尼亚预测是由3月19日百分之几。但在这一天,我们有一个地方爆发(跳跃)。这样的事件后,预测是没用几天,因为每天的预测值低于实际的。这种情况在几天之内改变(如长安2月12日)。SIR模型也有类似的问题。

亲爱的米兰,我一直在使用你的函数在智利和过去三天,预测是相当不错的。我加了95%的置信区间为“第二天”的预测。既然你已经在使用SML工具箱,它可能是更多的用户非常有用。这里是我的代码行:

[betaNL,RNL,JNL] = nlinfit(samplaTime(1:N),sampleC(1:N),@乐趣,COEF);
[Ypred,δ= nlpredci(@fun,[samplaTime(端)+1]”,betaNL,RNL, '雅可比行列式',JNL);
T =表(samplaTime(结束)+ 2,轮(Ypred),圆形(三角), 'VariableNames',{ '日', '预测', 'CI'})

再次感谢您的分享。

很酷!

伟大的工作米兰。对于葡萄牙,我建议要检查这个源(葡萄牙卫生部门):
https://covid19.min-saude.pt/ponto-de-situacao-atual-em-portugal/

“fitnlm”需要统计和机器学习工具箱,您应该更新需求。

对于一个更新的版本与冷凝的代码(在一个单一的.m文件)和自动数据检索用于COVID-19和每个国家:
//www.tianjin-qmedu.com/matlabcentral/fileexchange/74603-covid19modeling-fitvirus-adaptation

米兰伟大的工作!

嗨,你是什么意思你的最后一句话是什么意思?

sjcosta

你好,

首先感谢大家的MATLAB模型。这似乎很好地工作。我与葡萄牙的情况下更新,这似乎是完美也预测。如何更新葡萄牙号码?

Rolf Boelens教授为荷兰和美国提供数据和脚本。

优秀的典范。使用它的每一天。感谢您的辛勤工作。任何想法时,我们可以期待一个美国模式?

谢谢。该计划的预期目标是帮助人们评估时的疫情就会结束了,如果措施得力估计。现在,我在发布上述网址每日报告。

谢谢你的更新!伟大的工作!
我们可以期待新的图表,隔日?
在这里或其他地方?

彼得Graat

不错,但需要优化工具箱

ilPlus30

我尝试了模型更新的意大利数据。做得好。感谢分享。

连续回归使用MATLAB功能lsqcurvefit它有没有统计输出。于是又发作了与MATLAB功能fitnlm制造。结果可能会有所不同(小数据集) - 我不知道为什么 - 因此警告只是提醒一句要小心结果的解释。

绝对不可思议的工作和模式,现在只是在几个国家进行尝试。在德国,虽然回归模型看起来很合适,但我得到了如下信息:“***警告:lsqefit和fitnlm的结果存在显著差异。”
Knlm / Klsq = 358.476
rnlm / rlsq = 0.998951
Anlm / Alsq = 357.319"

我无法理解根本原因,因为最初的猜测是成功的……

我有没有Github上的经验,但我会努力做到你建议什么。

@milan我想把代码放在Github/Gitlab上,这样其他人就可以发送pull request等等。

好工作米兰!感谢分享。

他们可以做出自己的贡献MATLAB和随意添加fitVirus脚本。你有什么建议?

这是一个极好的主意,并为学生学习了一下有关这两个MATLAB和统计数据的大好机会。

什么是贡献数据文件对其他国家的最佳方式?

维车队

更新

2.1.2

改变形象

2.1.1

添加RMSE的图。更新数据。

魅惑

更改图形布局。更新数据。

2.0.25

禁止无效的预测时,案件的实际数量比实际更大

2.0.24

正确的描述

2.0.23

添加数据源

2.0.22

添加数据,丹麦,匈牙利,挪威,纽约州

2.0.21

提高初始猜测。比利时,克罗地亚,英国添加数据。升级数据荷兰(感谢罗尔夫Boelens)

2.0.20

新增约fitVirusCOVID19程序说明

2.0.19

正确的链接到新版本

2.0.18

正确的猜测(感谢Nikolas Wernecke)。将实际的每日用例添加到图中。病例总数和感染率的R2被添加到概要中。数据包括葡萄牙。

2.0.17

添加关于程序新版本的说明(感谢Joshua McGee)

2.0.16

更新的要求

2.0.15

数据与荷兰和美国的脚本添加(感谢罗尔夫Boelens)

2.0.14

更新说明

2.0.13

添加示例链接

2.0.12

为奥地利添加数据

)

正确的数据

2.0.10

更新数据。添加法国、瑞士的数据

2.0.9

更改摘要

2.0.8

为西班牙添加数据

2.0.7

更新说明

2.0.6

添加数据为德国

2.0.5

更新数据。加入总结报告现场脚本。

2.0.4

去掉上限。

2.0.3

更新数据。小的变化。R2现在包含在表中。

2.0.2

微小的变化

2.0.1

更改imaga

2.0.0

主要的修改。剔除威布尔回归,剔除高峰时间曲线图,合并曲线图进行流行评估及其发生率。增加流行持续时间和结束日期报告。用C0替换回归变量A。添加斯洛文尼亚的数据。

1.0.14

微小的修改

1.0.13

添加伊朗和中国以外的数据。威布尔回归现在是可选的。

1.0.12

将数据添加到2020年3月7日

1.0.11

韩国正确的数据

1.0.10

添加数据5.Mar.2020

1.0.9

每日预测的相对误差的正确计算

1.0.8

为意大利添加数据

1.0.7

正确的描述

1.0.6

添加数据4.Mar.2020

1.0.5

正确的描述

1.0.4

3添加数据2020年3月

1.0.3

更新说明

1.0.2

例如更新

1.0.1

添加图片

MATLAB版本兼容性
创建与R2019b
兼容任何版本
平台兼容性
视窗 苹果系统 Linux的