Re: [討論] <最終確診數> 簡易推估

作者: perceval (摸魚中)   2020-03-16 23:24:48
※ 引述《hugoyo (阿佑)》之銘言:
: 標題: [討論] <最終確診數> 簡易推估
: 時間: Mon Mar 16 21:04:52 2020
:
: 前言:
:
: 首先,確診數≠感染人數,檢測病毒本來就受限於醫療、政治等等因素。 本方法只是純數
: 學+各種簡化,大概拿來預測一下最終大概會有多少確診。
: 補充:
: 確 診 數 不是感染數!!! 要看該國要不要驗,要不要公布,以下只有純數字分析。
: 討論:
: 一個足夠多確診的群體的確診人數,多會呈現所謂的S型曲線,一開始少少的,中間突然暴
: 增,最後趨於緩和。 數學上可以用error function來模擬。
:
: error function: https://i.imgur.com/xg1NGiC.png
不是公衛專長,不過這種族群成長模型比較常看到用logistic function描述
https://en.wikipedia.org/wiki/Logistic_function
一時找不到演算法fit參數,所以沒辦法比較
不過估算出的最大確診人數可能會差不多
但是目前歐洲國家其實還沒有辦法用這個模型去描述
應該都還沒有達到反曲點(inflection point 中國稱為拐點)
現在應該都還在指數型成長的階段
可以參考worldmeters網站上的資料,選定感興趣的國家(假設是義大利)
找到 Total Coronavirus Cases in Italy
https://www.worldometers.info/coronavirus/country/italy/
圖片上選取logarithmic,你會發現確診數取log後和時間是線性的關係
這表示目前還是指數形成長
有興趣的可以用wolfram alpha去估算
exp fit{}指令去fit參數
假設我取最近五筆資料 exp fit{12462,15113,17660,21157,24747}
得到的模型是 10662.4 e^(0.169288 x)
估計今天總確診數會是 29443
不過受限於每日檢驗量,應該沒辦法達到這個數字
但是其它歐洲國家(西班牙、德國、法國)應該都可以得到不錯的估計結果

Links booklink

Contact Us: admin [ a t ] ucptt.com