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

看板nCoV2019作者 (摸魚中)時間4年前 (2020/03/16 23:24), 編輯推噓1(101)
留言2則, 1人參與, 4年前最新討論串2/3 (看更多)
※ 引述《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 不過受限於每日檢驗量,應該沒辦法達到這個數字 但是其它歐洲國家(西班牙、德國、法國)應該都可以得到不錯的估計結果 -- The road to hell is paved with good intentions - -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.116.21.126 (臺灣) ※ 文章網址: https://www.ptt.cc/bbs/nCoV2019/M.1584372292.A.B4E.html

03/19 07:11, 4年前 , 1F
還是覺得分布函數積分後用erf誤差比較小
03/19 07:11, 1F

03/19 07:14, 4年前 , 2F
只是Logistic在計算上比較好看
03/19 07:14, 2F
文章代碼(AID): #1URvf4jE (nCoV2019)
文章代碼(AID): #1URvf4jE (nCoV2019)