【エボラ出血熱】エボラ熱死者1900人超える [9/4] ★5

このエントリーをはてなブックマークに追加
756名無しさん@0新周年@転載は禁止
y'(t)=k*y(t)の解の平行移動であるy(t)=a*10^b+cに8/31までのWHO発表感染者数データを最小二乗法でフィッティングして
246.391 + 199.68*10^(0.014121 t)
を得た
感染者が死亡率Pでd日後に死ぬと仮定すると、
(246.391 + 199.68 10^(0.014121 (-d + t))) P
このモデルを死者数データに最小二乗法でフィッティングして
0.830483 (246.391 + 199.68*10^(0.014121 (-15.8181 + t)))
を得た

最新版8/31付

累積感染者数・死者数および予測グラフ
通常プロット http://imgur.com/1kFvnHY.png
片対数プロットhttp://imgur.com/FAV7Q1U.png

これらのフィッティング式を微分して、1日あたりの感染者数増加・死者数増加を得ると

通常プロット http://i.imgur.com/iak2XSo.png
片対数プロットhttp://imgur.com/GFYt6be.png

死亡率は83.0%、感染してから死亡するまでの生存期間は15.8日であり
http://imgur.com/oja4eQa.png
残差のlogを2次元プロットし妥当性を確認した

このPとdを生データに適用すると、全期間に渡り矛盾なく死亡率を説明することが出来る
http://imgur.com/FfdvHaG.png

ちなみにグラフはこんな風に作ってる
http://i.imgur.com/mikYjDf.png