中共政府公式データ
最初の内は武漢地方政府の発表で、途中から中共全国の結果に代わっています。 犠牲者数が少なかったころは、症例を細かに書いていましたが、今は数字だけに。 休みのせいか入れ替えの関係か、令和二年一月十九日分のデータがありません。
http://www.nhc.gov.cn/xcs/yqtb/list_gzbd.shtml
グラフの前半は、データ隠ぺいの反動で増え方が急なのだろうと仮定して、後半の傾きが緩くなったところを目分量で延長してみました。新型肺炎患者数と死者数が描かれています。データ欠損部分が醜くなっていますw 横軸も糞ですがw三月の頭頃まで伸ばしてみました。患者数二千万人、死者数五十万くらいになります。
後半の傾きが緩くなっているのが、病院等が飽和して十分に数を数えらていない可能性もありますが。
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)