数学中国

 找回密码
 注册
搜索
热搜: 活动 交友 discuz
查看: 4421|回复: 6

武汉肺炎感染人数应大致符合 Logistic 增长阻滞模型

[复制链接]
发表于 2020-1-29 08:51 | 显示全部楼层 |阅读模式
本帖最后由 Ysu2008 于 2021-5-5 20:59 编辑





从图上大致可看出目前正处于加速上升期……
谁来建个模?估算一下啥时候见顶。

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?注册

x

本帖被以下淘专辑推荐:

  • · 好貼|主题: 366, 订阅: 6
发表于 2020-1-29 14:33 | 显示全部楼层
自作聪明,模型早就有了。 SEIR 房室模型(compartmental model)
References
Nowak, M. and May, R. M. Virus dynamics: mathematical principles of immunology and virology. Oxford university press, 2000.
回复 支持 反对

使用道具 举报

 楼主| 发表于 2020-1-29 17:43 | 显示全部楼层
Nicolas2050 发表于 2020-1-29 14:33
自作聪明,模型早就有了。 SEIR 房室模型(compartmental model)
References
Nowak, M. and May, R. ...

,模型不唯一吧。
所有模型都是错的,只不过有些模型是有用的。
回复 支持 反对

使用道具 举报

发表于 2020-1-29 18:12 | 显示全部楼层
# -*- coding: utf-8 -*-
"""
Created on Tue Jan 21 22:24:00 2020
One is mever too old to learn.
@author: NicholasTU
一个感染者在传染期内在易感人群中引起的新感染者的期望人数,为传染病的传播速率,
有人称之为基本增生速率(Basic reproduction rate, R0)。
2020年1月25日
Dr. Eric Feigl-Ding, at Harvard University’s School of Public Health,
修正R0系数5.47;
2019 12 08是起始日期;
"""

import numpy as np
import matplotlib.pyplot as plt
plt.rcParams['font.sans-serif'] = ['SimHei'] #显示中文标签
plt.rcParams['axes.unicode_minus'] = False
import math

#武汉肺炎预测模型
gamma=1/14
#恢复率,44人
bb=0.05741
#感染率

n=13*10000*10000
#全国总人口
inf=1
rec=0
sus=n
t=1
inff=[]
#存放感染人数
recc=[]
#恢复健康的人数
suss=[]
#易感人群
time=[]
#从第一个病例发作开始计算天数
inff.append(inf)
recc.append(rec)
suss.append(sus)
time.append(t)
while True:       
     t=t+1
     if t<=48:
         k=5.47#哈佛大学推荐值
     if t>48:
           k=1.35#第48天,封城显现隔离效果,K减小;
     if t>200:
         break      
     beta=k*bb#beta感染率   
     a=inf   
     b=rec
     c=sus
     rec=b+gamma*a       
     inf=a+beta*a*c/n-gamma*a
     sus=n-rec-inf
     inff.append(round(inf))
     recc.append(rec)
     suss.append(round(sus))
     time.append(t)
for        i in range(200):
        print("第",time[i],"天全国累计的总感染人数=",n-suss[i])
   
#print(suss)
#print(inff)   
#print(recc)      
#plt.plot(time,suss,'b-',label='suspect')
plt.plot(time,inff,'r-',label='infect')
#plt.plot(time,recc,'g-',label='recover')
plt.xlabel('day')
plt.ylabel('population')
plt.xlim(0,201)
plt.title('武汉新型冠状病毒全国感染肺炎预测模型_2 by@NicholasTU')#
plt.legend(loc='upper right')
plt.show()
回复 支持 反对

使用道具 举报

发表于 2020-1-29 18:12 | 显示全部楼层
第 1 天全国累计的总感染人数= 0
第 2 天全国累计的总感染人数= 1
第 3 天全国累计的总感染人数= 2
第 4 天全国累计的总感染人数= 2
第 5 天全国累计的总感染人数= 3
第 6 天全国累计的总感染人数= 4
第 7 天全国累计的总感染人数= 4
第 8 天全国累计的总感染人数= 6
第 9 天全国累计的总感染人数= 7
第 10 天全国累计的总感染人数= 9
第 11 天全国累计的总感染人数= 11
第 12 天全国累计的总感染人数= 14
第 13 天全国累计的总感染人数= 17
第 14 天全国累计的总感染人数= 22
第 15 天全国累计的总感染人数= 27
第 16 天全国累计的总感染人数= 33
第 17 天全国累计的总感染人数= 42
第 18 天全国累计的总感染人数= 52
第 19 天全国累计的总感染人数= 64
第 20 天全国累计的总感染人数= 80
第 21 天全国累计的总感染人数= 99
第 22 天全国累计的总感染人数= 124
第 23 天全国累计的总感染人数= 154
第 24 天全国累计的总感染人数= 191
第 25 天全国累计的总感染人数= 237
第 26 天全国累计的总感染人数= 295
第 27 天全国累计的总感染人数= 367
第 28 天全国累计的总感染人数= 456
第 29 天全国累计的总感染人数= 566
第 30 天全国累计的总感染人数= 704
第 31 天全国累计的总感染人数= 875
第 32 天全国累计的总感染人数= 1087
第 33 天全国累计的总感染人数= 1351
第 34 天全国累计的总感染人数= 1679
第 35 天全国累计的总感染人数= 2086
第 36 天全国累计的总感染人数= 2592
第 37 天全国累计的总感染人数= 3221
第 38 天全国累计的总感染人数= 4003
第 39 天全国累计的总感染人数= 4974
第 40 天全国累计的总感染人数= 6181
第 41 天全国累计的总感染人数= 7680
第 42 天全国累计的总感染人数= 9543
第 43 天全国累计的总感染人数= 11859
第 44 天全国累计的总感染人数= 14736
第 45 天全国累计的总感染人数= 18311
第 46 天全国累计的总感染人数= 22753
第 47 天全国累计的总感染人数= 28273
第 48 天全国累计的总感染人数= 35132
第 49 天全国累计的总感染人数= 37236
第 50 天全国累计的总感染人数= 39352
第 51 天全国累计的总感染人数= 41481
第 52 天全国累计的总感染人数= 43623
第 53 天全国累计的总感染人数= 45778
第 54 天全国累计的总感染人数= 47946
第 55 天全国累计的总感染人数= 50127
第 56 天全国累计的总感染人数= 52322
第 57 天全国累计的总感染人数= 54530
第 58 天全国累计的总感染人数= 56751
第 59 天全国累计的总感染人数= 58986
第 60 天全国累计的总感染人数= 61234
第 61 天全国累计的总感染人数= 63496
第 62 天全国累计的总感染人数= 65772
第 63 天全国累计的总感染人数= 68061
第 64 天全国累计的总感染人数= 70364
第 65 天全国累计的总感染人数= 72682
第 66 天全国累计的总感染人数= 75013
第 67 天全国累计的总感染人数= 77359
第 68 天全国累计的总感染人数= 79718
第 69 天全国累计的总感染人数= 82093
第 70 天全国累计的总感染人数= 84481
第 71 天全国累计的总感染人数= 86884
第 72 天全国累计的总感染人数= 89302
第 73 天全国累计的总感染人数= 91734
第 74 天全国累计的总感染人数= 94181
第 75 天全国累计的总感染人数= 96643
第 76 天全国累计的总感染人数= 99120
第 77 天全国累计的总感染人数= 101612
第 78 天全国累计的总感染人数= 104118
第 79 天全国累计的总感染人数= 106641
第 80 天全国累计的总感染人数= 109178
第 81 天全国累计的总感染人数= 111731
第 82 天全国累计的总感染人数= 114299
第 83 天全国累计的总感染人数= 116883
第 84 天全国累计的总感染人数= 119483
第 85 天全国累计的总感染人数= 122098
第 86 天全国累计的总感染人数= 124729
第 87 天全国累计的总感染人数= 127377
第 88 天全国累计的总感染人数= 130040
第 89 天全国累计的总感染人数= 132719
第 90 天全国累计的总感染人数= 135415
第 91 天全国累计的总感染人数= 138127
第 92 天全国累计的总感染人数= 140855
第 93 天全国累计的总感染人数= 143600
第 94 天全国累计的总感染人数= 146362
第 95 天全国累计的总感染人数= 149140
第 96 天全国累计的总感染人数= 151936
第 97 天全国累计的总感染人数= 154748
第 98 天全国累计的总感染人数= 157577
第 99 天全国累计的总感染人数= 160423
第 100 天全国累计的总感染人数= 163287
第 101 天全国累计的总感染人数= 166168
第 102 天全国累计的总感染人数= 169067
第 103 天全国累计的总感染人数= 171983
第 104 天全国累计的总感染人数= 174916
第 105 天全国累计的总感染人数= 177868
第 106 天全国累计的总感染人数= 180837
第 107 天全国累计的总感染人数= 183825
第 108 天全国累计的总感染人数= 186830
第 109 天全国累计的总感染人数= 189854
第 110 天全国累计的总感染人数= 192896
第 111 天全国累计的总感染人数= 195957
第 112 天全国累计的总感染人数= 199036
第 113 天全国累计的总感染人数= 202133
第 114 天全国累计的总感染人数= 205250
第 115 天全国累计的总感染人数= 208385
第 116 天全国累计的总感染人数= 211540
第 117 天全国累计的总感染人数= 214713
第 118 天全国累计的总感染人数= 217906
第 119 天全国累计的总感染人数= 221118
第 120 天全国累计的总感染人数= 224350
第 121 天全国累计的总感染人数= 227601
第 122 天全国累计的总感染人数= 230872
第 123 天全国累计的总感染人数= 234162
第 124 天全国累计的总感染人数= 237473
第 125 天全国累计的总感染人数= 240803
第 126 天全国累计的总感染人数= 244154
第 127 天全国累计的总感染人数= 247525
第 128 天全国累计的总感染人数= 250917
第 129 天全国累计的总感染人数= 254329
第 130 天全国累计的总感染人数= 257762
第 131 天全国累计的总感染人数= 261215
第 132 天全国累计的总感染人数= 264690
第 133 天全国累计的总感染人数= 268185
第 134 天全国累计的总感染人数= 271702
第 135 天全国累计的总感染人数= 275240
第 136 天全国累计的总感染人数= 278799
第 137 天全国累计的总感染人数= 282380
第 138 天全国累计的总感染人数= 285983
第 139 天全国累计的总感染人数= 289607
第 140 天全国累计的总感染人数= 293254
第 141 天全国累计的总感染人数= 296922
第 142 天全国累计的总感染人数= 300613
第 143 天全国累计的总感染人数= 304326
第 144 天全国累计的总感染人数= 308061
第 145 天全国累计的总感染人数= 311819
第 146 天全国累计的总感染人数= 315600
第 147 天全国累计的总感染人数= 319404
第 148 天全国累计的总感染人数= 323231
第 149 天全国累计的总感染人数= 327081
第 150 天全国累计的总感染人数= 330954
第 151 天全国累计的总感染人数= 334851
第 152 天全国累计的总感染人数= 338771
第 153 天全国累计的总感染人数= 342715
第 154 天全国累计的总感染人数= 346683
第 155 天全国累计的总感染人数= 350675
第 156 天全国累计的总感染人数= 354691
第 157 天全国累计的总感染人数= 358731
第 158 天全国累计的总感染人数= 362796
第 159 天全国累计的总感染人数= 366885
第 160 天全国累计的总感染人数= 370999
第 161 天全国累计的总感染人数= 375138
第 162 天全国累计的总感染人数= 379302
第 163 天全国累计的总感染人数= 383492
第 164 天全国累计的总感染人数= 387706
第 165 天全国累计的总感染人数= 391946
第 166 天全国累计的总感染人数= 396212
第 167 天全国累计的总感染人数= 400504
第 168 天全国累计的总感染人数= 404821
第 169 天全国累计的总感染人数= 409165
第 170 天全国累计的总感染人数= 413534
第 171 天全国累计的总感染人数= 417931
第 172 天全国累计的总感染人数= 422354
第 173 天全国累计的总感染人数= 426803
第 174 天全国累计的总感染人数= 431280
第 175 天全国累计的总感染人数= 435783
第 176 天全国累计的总感染人数= 440314
第 177 天全国累计的总感染人数= 444872
第 178 天全国累计的总感染人数= 449458
第 179 天全国累计的总感染人数= 454071
第 180 天全国累计的总感染人数= 458713
第 181 天全国累计的总感染人数= 463382
第 182 天全国累计的总感染人数= 468080
第 183 天全国累计的总感染人数= 472806
第 184 天全国累计的总感染人数= 477560
第 185 天全国累计的总感染人数= 482344
第 186 天全国累计的总感染人数= 487156
第 187 天全国累计的总感染人数= 491997
第 188 天全国累计的总感染人数= 496868
第 189 天全国累计的总感染人数= 501768
第 190 天全国累计的总感染人数= 506697
第 191 天全国累计的总感染人数= 511657
第 192 天全国累计的总感染人数= 516646
第 193 天全国累计的总感染人数= 521665
第 194 天全国累计的总感染人数= 526715
第 195 天全国累计的总感染人数= 531795
第 196 天全国累计的总感染人数= 536906
第 197 天全国累计的总感染人数= 542048
第 198 天全国累计的总感染人数= 547221
第 199 天全国累计的总感染人数= 552425
第 200 天全国累计的总感染人数= 557661
回复 支持 反对

使用道具 举报

发表于 2020-2-10 18:42 | 显示全部楼层
各地疫情缓和日期预测

laiyindata  莱茵数据在线

疫情在全国范围内已经肆虐了大半个月,看着每天上升的确诊人数,14亿人的感受想必完全一致:担忧!紧接着,我们又会不禁发问:疫情何时结束?往日的繁华何时可以恢复?这些日子以来,不同的平台都给出各种预测,答案也是五花八门。2月2日,诺贝尔奖得主Michael Levitt就预测疫情很快会结束;可是一些专家也认为疫情的肆虐可能还需要持续很长一段时间。在这样的背景下,大家可能也很关心自己的家乡和工作地何时迎来平静的生活。针对不同地区的精细分析还十分缺乏。我们团队利用Logistic回归和各地已经公布的数据做了统计回归,预测了各个省级行政区的疫情增长情况,给出了各地疫情明显缓和的日期。


图1.基于全国数据的Logistics拟合

Logistic 回归,译作对数几率回归或者罗吉斯回归,是一种对数几率模型,是离散选择法模型之一,是生物统计学、临床、数量心理学、计量经济学、市场营销等统计实证分析中最常用的方法之一。获得结果尽管不能尽善尽美,但参考价值很高。

我们的模型显示,全国总体疫情将于2月22日缓解(确诊病例不再明显增长)。各地疫情会在此前相继落幕。最受大家关心的湖北,我们预言会在2月22日出现平台期,即单日增加确诊人数小于50人。而浙江、广东和河南三个疫情肆虐的省份会在2月11日、2月14日和2月16日进入缓和阶段。具体地的疫情缓和日期预测,我们详细地列入了表1。供大家参考,快来看看吧。特别说明:部分地区的数据不足以用来建立模型,我们就没有进行预测。

大家一定要有信心,困难一定会过去。在党和政府的领导下,我们一定可以迎来曙光!有那么医护人员和科研人员在为大家奋斗,我们一定可以迎来曙光!

表1:各地疫情缓和(确诊病例不再明显增长)时间预测

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?注册

x
回复 支持 反对

使用道具 举报

 楼主| 发表于 2020-2-10 20:28 | 显示全部楼层
英雄所见略同。
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|手机版|小黑屋|数学中国 ( 京ICP备05040119号 )

GMT+8, 2024-3-29 03:47 , Processed in 0.081055 second(s), 18 queries .

Powered by Discuz! X3.4

Copyright © 2001-2020, Tencent Cloud.

快速回复 返回顶部 返回列表