一個簡單的模型估算。貼子有點長,僅供討論。

來源: 簡體 2020-02-16 16:23:16 [] [博客] [舊帖] [給我悄悄話] 本文已被閱讀: 次 (10875 bytes)

這是個很簡單的模型。

 

初始設定是十二月一號有首例,潛伏期為7天,確診後的人漸漸恢複,在四個星期時,95.5%的人痊愈。死亡率在確診後三星期設在3.5%,四個星期後在4.5%

傳染係數在十二月和大部分的一月設在3.2,後麵依時降低,二月份嚴控後R03.2急降到0.4。由於潛伏期我設得較實際情況簡單,又遠小於最大潛伏時間,所以在二月底及往後取一個比較小又漸減的值。

傳染率深究下來也是很複雜的。這個簡化版是一個病人在7天內傳給他人的數目,在病人得到確診後被隔離到病區,不再傳播給未感染者。

 

這些前提下,最終感染數在十二萬多,最終死亡數在近五千多人。最高留院人數在靠近二月底,達近十萬人。最高死亡在近期,二月中下旬,同時也是痊愈人員大增期。預計二月底將看到疫情的拐點,走向日新增病例和死亡人數的雙下降和痊愈人數的急增。

 

但是,這個模式比起直實情況顯然低估了。

 

模式各項數據在二月初比官方報告的都少很多。原因之一是初始條件十二月一號僅有一例。官報有一例,但實際上肯定不止這個數。你看見了一隻狼,這意味著林子裏有一群狼。如果初始病例往前挪,十二月一號有同代的兩例的話,所有模擬數據要乘二!如果糾正到二月初的死亡人數,那是個好幾倍的關係。這個比較讓人出冷汗。

 

另一個看出模式低估了是:圖三紅線的每天死亡都在200人以下。這個是不現實的,殯儀館的情況看出比這個要高得多。

 

所以,一個簡單的糾正是要乘以一個倍數,去和殯儀館的相吻合。讀者可以自己估計一個,因為官方數據的不透明性,每個人又有自己的想法。我的估計是各個數字乘以三倍。那麽死亡數最終可能達一萬五千人。我設立的最終死亡率是4.5%。這個讀者自行感受並調整。

 

時間線上會有誤差。R0的設定不完美,而這個又是個指數關係。潛伏期設在7天。

但總體上會看到一些好的趨勢,二月底或者三月初會見到拐點。這個是因為二月份的嚴控,R0急速降低。如果有更多的藥物出來,後期數據走向還可能好一點。

 

此外我調控過R0來做敏感度試驗。越早控製越重要。對於如此烈性的病毒,時間就是生命。提前一個星期嚴控對最終結果有極大的影響。

 

所有跟帖: 

我記得傳染係數在3.77左右可能到4,當然都是估算的。不過從現有和需要預備的床位基本能知道大致範圍。 -itistrue- 給 itistrue 發送悄悄話 itistrue 的博客首頁 (0 bytes) () 02/16/2020 postreply 16:27:32

初期的我曾設成3點8,結果會嚇人好多。我往下調了,一個原因是想reproduce官方的數據,至少是它的trend。 -簡體- 給 簡體 發送悄悄話 簡體 的博客首頁 (0 bytes) () 02/16/2020 postreply 16:30:08

可以從床位數反推大致知道能收住多少,還有不住院的有輕症,或者遺漏的。。。 -itistrue- 給 itistrue 發送悄悄話 itistrue 的博客首頁 (0 bytes) () 02/16/2020 postreply 16:32:59

你說的很有道理。我也想過。那些輕症自愈的有多少包括在內都是疑問。官方數據應用同時通報預報的才好,就像我們台風預警一樣。 -簡體- 給 簡體 發送悄悄話 簡體 的博客首頁 (0 bytes) () 02/16/2020 postreply 16:43:27

打仗你管他戰場多大幹啥呢?你看上多少兵力嘛 -顛三倒四了- 給 顛三倒四了 發送悄悄話 (295 bytes) () 02/16/2020 postreply 16:46:54

這樣簡單呀。兵力是可以隨時變動補充的,但床位數變動是受限製的。 -itistrue- 給 itistrue 發送悄悄話 itistrue 的博客首頁 (0 bytes) () 02/16/2020 postreply 16:51:15

地鋪算不算 躺椅算不算 -顛三倒四了- 給 顛三倒四了 發送悄悄話 (73 bytes) () 02/16/2020 postreply 16:55:03

這個就是變動係數,可以估算的。因為也是受限製的。 -itistrue- 給 itistrue 發送悄悄話 itistrue 的博客首頁 (0 bytes) () 02/16/2020 postreply 16:56:13

誰能保證病毒不變異 -fuz- 給 fuz 發送悄悄話 fuz 的博客首頁 (0 bytes) () 02/16/2020 postreply 16:32:31

到目前,官方沒說有變異。如果變弱就好了。 -簡體- 給 簡體 發送悄悄話 簡體 的博客首頁 (0 bytes) () 02/16/2020 postreply 16:34:29

官方已經說了每傳一代測基因序列都有一點點不一樣. 這就是變異. -janeice65- 給 janeice65 發送悄悄話 (0 bytes) () 02/16/2020 postreply 17:22:40

哦。我的模型很簡單。投入值在圖頂端。r0動態,但死亡率設成一定的了。隻供看個大概吧。full model 要花時間的。 -簡體- 給 簡體 發送悄悄話 簡體 的博客首頁 (0 bytes) () 02/16/2020 postreply 17:28:13

我作為一個科技盲,是不是這樣理解, -要管3721- 給 要管3721 發送悄悄話 要管3721 的博客首頁 (215 bytes) () 02/16/2020 postreply 17:55:31

變異得太快疫苗會來不及:( -簡體- 給 簡體 發送悄悄話 簡體 的博客首頁 (348 bytes) () 02/16/2020 postreply 18:04:02

可是,如果是生化武器,這種變弱似乎隻是願望。 -要管3721- 給 要管3721 發送悄悄話 要管3721 的博客首頁 (602 bytes) () 02/16/2020 postreply 18:18:39

妹妹 -簡體- 給 簡體 發送悄悄話 簡體 的博客首頁 (155 bytes) () 02/16/2020 postreply 18:40:58

For viruses that mutate fast, we may have to --worldance-- 給 -worldance- 發送悄悄話 (347 bytes) () 02/16/2020 postreply 20:14:24

可是,厲害國不讓大家去呀。 -要管3721- 給 要管3721 發送悄悄話 要管3721 的博客首頁 (56 bytes) () 02/16/2020 postreply 20:17:57

昨天我在JT 說這個死亡人數,人家告訴我把官方加一個零,加兩個零隨我的便,現在看來最後的數字加一個零都可能是低估,也許發聲的都是 -綠蟻采菊- 給 綠蟻采菊 發送悄悄話 綠蟻采菊 的博客首頁 (260 bytes) () 02/16/2020 postreply 16:33:27

正常生病死人 本地醫生很輕鬆就對付了 -顛三倒四了- 給 顛三倒四了 發送悄悄話 (73 bytes) () 02/16/2020 postreply 16:37:16

小顛你總是三言兩語道出最明白的道理~我…… -簡體- 給 簡體 發送悄悄話 簡體 的博客首頁 (216 bytes) () 02/16/2020 postreply 16:40:28

他一點都不顛三倒四,裝的。為你的認真點讚一個。 -Bali- 給 Bali 發送悄悄話 Bali 的博客首頁 (0 bytes) () 02/16/2020 postreply 17:05:03

今天高級了 喝的威士忌了 -顛三倒四了- 給 顛三倒四了 發送悄悄話 (0 bytes) () 02/16/2020 postreply 17:06:17

嗯,倒過來不再顛三倒四樂 -早晨從中午開始- 給 早晨從中午開始 發送悄悄話 (0 bytes) () 02/16/2020 postreply 17:11:06

你現在是早晨呢 還是中午呢? -顛三倒四了- 給 顛三倒四了 發送悄悄話 (0 bytes) () 02/16/2020 postreply 17:12:03

恍惚中,我也去喝點,這些日子真鬱悶啊 -早晨從中午開始- 給 早晨從中午開始 發送悄悄話 (0 bytes) () 02/16/2020 postreply 17:19:43

喝高了負負得正顛回去了 -Walede- 給 Walede 發送悄悄話 (0 bytes) () 02/16/2020 postreply 17:11:30

哎嘛 為了跟得上灌水節奏 浪費不少二鍋頭了 -顛三倒四了- 給 顛三倒四了 發送悄悄話 (0 bytes) () 02/16/2020 postreply 17:12:58

謝謝巴厘妹妹:)就是想自己看拐點給逼出來的。 -簡體- 給 簡體 發送悄悄話 簡體 的博客首頁 (80 bytes) () 02/16/2020 postreply 17:19:23

你就是閑著愛動腦子 我們不動腦子的 就看幾大航空公司停飛的日期就好 請耐心等待四月底 -顛三倒四了- 給 顛三倒四了 發送悄悄話 (0 bytes) () 02/16/2020 postreply 17:21:31

嗯哪~抱著希望,學會等待 -簡體- 給 簡體 發送悄悄話 簡體 的博客首頁 (0 bytes) () 02/16/2020 postreply 17:33:24

早破萬了 -Walede- 給 Walede 發送悄悄話 (0 bytes) () 02/16/2020 postreply 16:58:26

按這個模式現在是高峰,我的模式很簡單,也參考了官方的數據。繼續看吧。我希望最終死亡率隨著藥物的引入,和更好的醫護,病毒的自身削弱 -簡體- 給 簡體 發送悄悄話 簡體 的博客首頁 (0 bytes) () 02/16/2020 postreply 16:38:20

而降低。 -簡體- 給 簡體 發送悄悄話 簡體 的博客首頁 (0 bytes) () 02/16/2020 postreply 16:38:45

這就是我今早發帖的起點。因為看到有人質疑你的說法。 -pollyli- 給 pollyli 發送悄悄話 pollyli 的博客首頁 (0 bytes) () 02/16/2020 postreply 18:11:02

厲害。我一直想的是,確診人數=治愈人數+死亡人數嗎? -夏茗- 給 夏茗 發送悄悄話 夏茗 的博客首頁 (0 bytes) () 02/16/2020 postreply 17:12:11

是的,這個平衡檢查過了。 -簡體- 給 簡體 發送悄悄話 簡體 的博客首頁 (0 bytes) () 02/16/2020 postreply 17:17:12

沒治愈的在哪裏 -yukimama- 給 yukimama 發送悄悄話 (0 bytes) () 02/16/2020 postreply 17:23:06

第三張圖的紅線,是每天待床的人數。最後的治愈和死亡等於所有的確診人數。 -簡體- 給 簡體 發送悄悄話 簡體 的博客首頁 (0 bytes) () 02/16/2020 postreply 17:25:42

從一美國nerd 網站估算死亡人數,武漢火葬爐總共有49個, 就算二小時一個人, 24小時一天再乘以多少天就是死亡人數 -janeice65- 給 janeice65 發送悄悄話 (0 bytes) () 02/16/2020 postreply 17:27:26

我也想過。但那個數太大了,用來做calibrate的話 -簡體- 給 簡體 發送悄悄話 簡體 的博客首頁 (264 bytes) () 02/16/2020 postreply 17:32:29

50分鍾1人。以前正常運行是從早上6點前到中午12點。 -看客2010- 給 看客2010 發送悄悄話 看客2010 的博客首頁 (0 bytes) () 02/16/2020 postreply 17:42:53

那就一天上千人。跟我的有10倍關係。 -簡體- 給 簡體 發送悄悄話 簡體 的博客首頁 (116 bytes) () 02/16/2020 postreply 17:59:16

還沒算上世衛組織去武漢前二天燒醫療垃級導致SO2 大幅升高的估值. -janeice65- 給 janeice65 發送悄悄話 (0 bytes) () 02/16/2020 postreply 18:09:42

是啊,算多少爐子隻是保守數據。 -看客2010- 給 看客2010 發送悄悄話 看客2010 的博客首頁 (0 bytes) () 02/16/2020 postreply 18:18:49

這個用來定量,不確定性比較大,雖然私下認為是個好的指標性信息。但是, -簡體- 給 簡體 發送悄悄話 簡體 的博客首頁 (218 bytes) () 02/16/2020 postreply 18:36:32

看到的最高數字是,武漢火葬場總容量782人/天 -味苦- 給 味苦 發送悄悄話 味苦 的博客首頁 (0 bytes) () 02/16/2020 postreply 19:09:00

加跟帖:

當前帖子已經過期歸檔,不能加跟帖!