图片缩略图

fitVirusXX

版本2.0.3 (7.54 MB) 米兰巴蒂斯塔
基于logistic模型的冠状病毒流行评估

22下载

更新2020年7月13日

查看许可证

该函数拟合一组延迟logistic曲线与Covid-19疫情数据。数据的大小限制了曲线的数量。该模型是为日常监测流行病发展而设计的,不是为长期预测而设计的。生成的图表是该流行病的每日横截面图。对于任何预测,都需要监测图表的时间发展情况;最重要的是,疫情的最终规模。该模型是现象学的,在流行病持续蔓延的情况下优于单曲线模型(SI模型、SIR模型)。

该功能通过执行runMe.m函数调用。此功能默认情况下,呼吁世界中的数据。要获取数据斯洛文尼亚,作为一个例子,你需要在runMe.m功能变化@getDataWorld与@getDataSlovenia。

有关程序选项,请参见fitVirusXXdoc。mlx在doc文件夹中。

目前不同国家的数据可以从两个来源获得:
1.运行importTotalCases.m功能。此功能将读取数据<https://ourworldindata.org/coronavirus-source-data>并将其存储在数据文件夹。
2.运行importTotalCasesWM.m功能。此功能将读取数据<https://www.worldometers.info/coronavirus/>并将其存储在数据文件夹。

注1:该模型的预测取决于如何很好的模型模拟疫情的过程中,并在数据有多可靠。

注2:预测与新的或更改的数据变化。预测成为对泛滥的局面更加可靠,但在一开始不可靠的。

注3。在某些情况下,模型可能会失败。特别是,该模型可能不适合于模拟流行病的早期阶段或新一波疫情的出现。在这些情况下,我们可以得到不现实的指数增长,而这种增长通常会随着新的数据而减弱。如果最后一波呈指数级增长,那么我们需要减少波的数量,或者减少可能受到影响的波的数量(默认情况下,这个数字是7e9)。

注3。如果在最后一波中病例数迅速增加,波形较窄,给人一种波即将结束的感觉(logistic曲线相对于峰值是对称的)。在这种情况下,必须监测疫情规模的每日变化,然后对其持续时间作出结论。

注4.程序指示当数据足以计算初始近似值的开始日期。

注5。默认情况下,该函数尝试将数据拟合到四条逻辑逻辑曲线(参见“model”选项)。它还通过减去模型从数据中预测的极限大小来自动检测可能的下一波的开始。其他的逻辑曲线与这个额外的波动相匹配(参见‘w3’选项)。

注6。模型的参数是通过最小化目标函数得到的,目标函数是值的残差平方和。优化工具箱函数fminsearch用于计算未知模型参数的最优值。利用启发式方法得到了拟最优解。如果计算失败,则只绘制数据。

注7。4+波的模型可能会很费时间。对于一个180个数据的样本,四波模型所花费的时间约为180秒。

注8。根据该模型作出的每日预测载于
<https://www.researchgate.net/publication/339912313_Forecasting_of_final_COVID-19_epidemic_size_200616>

在疫情评估图上,将各地区划分为疫情的不同阶段(不规范,为方便起见,可任意选择):
红色——一个快速增长阶段(大约12%到88%的最终情况下)
黄色-平衡过渡阶段(98%)
绿化-最后阶段(高原阶段)

日例生长因子图两行1%(绿色)和5%
(红色)仅供方向之用。

免责声明:本软件和数据仅用于教育目的,不用于医疗或商业用途。它们按原样提供,所有明示或默示保证,包括适销性和适合某一特定用途的默示保证,均被拒绝。在某些情况下,模型可能会失败。特别是,模型可能不合适,或者模型可能在早期阶段失败。你可以自行决定如何使用它。

fitVirusXX释放不支持。万博1manbetx

引用作为

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

意见及评分(20.

米兰巴蒂斯塔

你可以运行fitVirusXX物流功能。例如,如果你运行AA = fitVirusXX(@getDataWorld),则C =物流(10:100,AA,B)将预测的情况下为10倍的载体:100。

阿迪尔·马吉德

有没有一种方法,我们可以将生成的MATLAB图数据转换到Excel?如果是的话那么它如何可能呢?请帮忙

米兰巴蒂斯塔

它运行在我的电脑没有任何错误。

董gretch

这里的runMeXX的内容
关闭所有
目录(“数据”)
目录(“fitVirusXX”)
抽搐
aa = fitVirusXX (@getDataPhilippines,……
“回声”,“上”、“模型”,4,“ndat”,[],“拯救”,1,‘jpres’,‘r300’,‘w3’,3);
toc

米兰巴蒂斯塔

我没有得到任何错误。你能把你的runmex的内容显示出来吗?

董gretch

嗨米兰我仍然拥有最新版本的这个错误。你能帮我解决这个问题吗?
> > runMeXX
错误使用fitVirusXX / paramEst2
方法“paramEst2”没有为类“fitVirusXX”定义,或者从MATLAB的搜索路径中删除。

错误fitVirusXX(第215行)
paramEst2(OBJ, '回响',echof, '模型',型号, 'ndat',ndat, 'n最大',n最大);

runMeXX中的错误(第5行)
aa = fitVirusXX (@getDataPhilippines,……

米兰巴蒂斯塔

我没有得到任何错误当前版本。我用Matlab R2019b。

SangSup曹

谢谢你的关心。我用的是韩国的数据,即AA = fitVirusXX(@getDataSouthKorea,...
“回声”,“上”、“模型”,4,“ndat”,[],“拯救”,1,‘jpres’,‘r300’,‘w3’,1),我得到这个问题。
在我看来,数据处理功能好像出了问题。

米兰巴蒂斯塔

你可以指定你尝试运行,其作用是什么数据?

SangSup曹

谢谢你的好意。我得到这个错误。你能解释一下怎么修吗?
错误:fitVirusXX/plot(第206行)
出租车= datetime (datestr(趋势+ obj.time0));
datetime(第640行)

米兰巴蒂斯塔

如果“w3”选项不是0或1,v07版本将产生错误。所以你应该改变'w3',3到'w3',1:
aa = fitVirusXX (@getDataPhilippines,……
“回声”,“上”、“模型”,3,“ndat”,[],“拯救”,1,‘jpres’,‘r300’,‘w3 ', 1);

在当前版本v08中,这个问题得到了解决,所以你最多可以指定3个额外的wave搜索。

董gretch

嗨,米兰。我有这个错误。你能帮我纠正一下吗?
误差在runMeXX(第4行)
aa = fitVirusXX (@getDataPhilippines,……

米兰巴蒂斯塔

是的情况下,当疫情拖累。

ravib1996

谢谢你的努力,真正有助于可视化预测世界各地。这个模型比你之前的模型(fitVirusCV19)有更好的精度/RMSE吗?

不幸的是,这个链接没有用,因为它是用另一种语言写的。

米兰巴蒂斯塔

顺便说说。我使用了多个逻辑的模型,你释放你的medRxiv之前预印本https://www.rtvslo.si/zdravje/novi-koronavirus/stevilo-potrjenih-okuzb-preseglo-tisocico-epidemija-se-lahko-obrne-v-dve-smeri/519650

米兰巴蒂斯塔

我不知道一个数学公式怎么能获得版权,我没有。版权为程序所有,由Matlab自动插入。此外,logistic函数的求和与指数函数的求和(即傅里叶级数;版权吗?)。如果有类似的物流研究,而不是指数函数,你应该明确地提到它,以便我能引用它。不管怎样,这个主题已经有100年的历史了:
《逻辑曲线的总和》,《皇家统计学会学报》,第90卷,第4期(1927年),第729-746页

你是第一个提出研究员使用多(或您给出的实例为三个)物流(指数)函数的想法?还有谁提出并实施了使用多个(或三个在他的预印本)乙状结肠的分支的另一位研究员(指数)函数个月前发布,但您还没有提到这项工作;相反,你给你的“求和法”的版权?使用三个或多个指数函数的总和的想法已经被他人版权。即使更改求和公式有点像你一样,仍然需要你给的中心思想所需的全部学分的前业主。但是,你没有。

更新

2.0.3

改变形象

2.0.2

添加更多的选项迭代。更新初始猜测功能。

2.0.1

更新描述

2.0.0

主要的修改。波的数量现在只受数据大小的限制。添加从worldmeters.info导入数据的功能。

1.3.0

添加第四波。

1.2.0

更新“w3”选项。取消时间延迟的下限。

1.1.2

微小的修改

1.1.1

改变形象

1.1.0

自动检测多达4个波

1.0.1

微小的修改

MATLAB版本兼容性
创建R2020a
与任何版本兼容
平台兼容性
窗户 macOS Linux的

fitVirusXX_v13

fitVirusXX_v13 /数据

fitVirusXX_v13 / fitVirusXX

fitVirusXX_v13 / fitVirusXX / @ fitVirusXX

fitVirusXX_v13 /验证