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

這是個很簡單的模型。

 

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

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

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

 

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

 

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

 

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

 

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

 

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

 

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

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

 

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