fortran66のブログ

fortran について書きます。

【やっつけ】武漢新型肺炎死者数グラフ

中共政府公式データ

最初の内は武漢地方政府の発表で、途中から中共全国の結果に代わっています。 犠牲者数が少なかったころは、症例を細かに書いていましたが、今は数字だけに。 休みのせいか入れ替えの関係か、令和二年一月十九日分のデータがありません。

http://www.nhc.gov.cn/xcs/yqtb/list_gzbd.shtml

グラフの前半は、データ隠ぺいの反動で増え方が急なのだろうと仮定して、後半の傾きが緩くなったところを目分量で延長してみました。新型肺炎患者数と死者数が描かれています。データ欠損部分が醜くなっていますw 横軸も糞ですがw三月の頭頃まで伸ばしてみました。患者数二千万人、死者数五十万くらいになります。

後半の傾きが緩くなっているのが、病院等が飽和して十分に数を数えらていない可能性もありますが。

f:id:fortran66:20200201224332p:plain

bukan.xy

15 41 2
16 45 2
17 62 2
18 198 3
19 0 0
20 291 6
21 440 9
22 571 17
23 830 25
24 1287 41
25 1975 56
26 2744 80
27 4515 106
28 5974 132
29 7711 170
30 9692 213
31 11791 259
open(10, file='bukan.xy')      
print *, 'start'
do i = 1, 10000    
  read(10, *, end = 99) id, ipatient, ideath
  write(9, *) id, ipatient, ideath, 4.0*exp(0.25*i), 200.0 * exp(0.25*i)
end do
99 continue
do j = 1, 30
  write(9, *) id + j, 0, 0, 4.0*exp(0.25*(i + j - 1)), 200.0 * exp(0.25*(i + j - 1))
end do
close(10)
close(9)
stop 'normal end'
end
 start
STOP normal end
%fig:
x, y1, y2, y3, y4 = np.loadtxt('fort.9', unpack = True) 
plt.grid(which='both')
plt.yscale('log')
plt.plot(x, y1)
plt.plot(x, y2)
plt.plot(x, y3)
plt.plot(x, y4)
plt.title('Coronavirus death toll', fontsize = 25)
plt.xlabel('date', fontsize = 10)
plt.ylabel('number of person', fontsize = 10)