- エボラの流行 2014年 アフリカ PART4
19 ：名無しさん＠お腹いっぱい。[]：2014/09/07(日) 11:07:07.36 ID:lCgbujxr - 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
|
|