用数学模型看病毒源问题

—— by 万精油 2020-03-22

中国政府发言人在推特上说“美国军人把病毒带到武汉”。这样的话可以说是很严重的指控。用”might be”似乎是想淡化这种严重性,但实际上却说明是空穴来风,没有根据(否则绝不会以这样的方式出现)。英文推特下的评论基本都是骂他的,中文(微博)评论大部分观点却是“川大嘴都可以乱说,为什么我们不可以”。对此观点,我的一个读者评论道:“川普胡说八道,有美国主流媒体批判,而中国官方胡说八道,喉舌集体转发,网民欢呼叫好,这是媒体和数亿人的问题,不是官方几个人的问题,区别很大”。我认为这个评论点到了问题本质。

空穴来风的指控,给美国反华派提供借口。川普喊出了“中国病毒”,一些电视台也开始用“中国病毒”,“武汉病毒”了。更有少数极端分子开始以此为由歧视华裔。这要归功于罩大嘴逞一时口舌之快。当然,罩大嘴是政府发言人,代表政府说话。只是苦了众多华裔。所谓神仙打仗,凡人遭殃。

抛开政治因素不说,我们单纯从数字上来看看这个“美国军人把病毒带到武汉”的可能性。

病毒传播是指数模式,其快慢程度由一些参数决定。目前看到的两个模式是:

1. R0=3, G=4; 也就是说平均一个感染者传递给3个人。4天一代。这个模型是著名网络科普者李永乐所用的模型,他说他的依据是现有数据,被广为接受。

2. 加倍期 d = 3天。也就是说每三天,感染人数翻倍。中国人根据中国数据发的文章说加倍期是一天半(武汉),两天(湖北),三天(新疆)等等。英文文章(包括福布斯,纽约人,新英格兰医学)说3天,4天,6天。我看到的最全的文章里列出了世界各地的数字,在无控制的初期,加倍时间不到3天。有一定控制后,平均6天。所以,在无控制的状态下,我用3天算是比较保守的。

这两种模型实际上是等价的,只是参数不同而已。

军运会是10月18号开始的。到一月23号封城,有97天。这97天是完全无控制,而且还有万人餐之类的几万人参加的活动帮助扩散(增大R0)。另外,作为九州通衢的武汉,这些天是水陆空全开,加上春运前期,可以说是最大效率扩散。在这种完全无控制的情况下,两种模型扩散97天的结果是什么?

附图是我用MATLAB做的。97天下来,两者都会传染给全世界所有人。图形中,横坐标是从10月18号开始的天数。纵坐标是人数。蓝线是第一种模型的增长曲线。黑线是第二种模型的增长曲线。上面的红色横线是世界现有人口数。右面的红色竖线是武汉封城的1月23号。可以看到,到1月23号,两种模型都达到世界人口数。实际上在第一种模型下,只需82天。当然,这两个模型都是简单模型,实际上当感染到人口一定比例后(比如超过一半),感染率就会降下来。但是,不要说一半,即使感染全世界人口十分之一,也是超出想象的恐怖,而不是武汉封城时的情况。真有全世界十分之一被感染的话,封不封城都无所谓了。

换一个角度。李文亮受到训诫是12月31号。那个时候离军运会开始已经74天。在第一种模型下,74天已经感染十亿人。现在公认的死亡率是1%(老年人死亡率更高),十亿人的1%是一千万。如果有1千万人死去(或正在死去),是一种什么情况?就算加上滞后期,只有几十万人死亡,这个数量的死亡还用得着有人吹哨吗?

结论:“美国军人把病毒带到武汉”数学上讲不通。QED。

注:第一种模型的公式是:

N = (R0^(t/G+1)-1)/(R0-1)

第二种模型的公式是:

N = 2^(t/d)

上面的t是天数。

附录:在微博上看到有人声称湖北是唯一一个最多只需跨两省就到达中国大陆各省的地方。而且有人给了一个地图,果然如此。中间蓝色是湖北。绿色省份一步到位。黄省两步。深绿色省份三步到位,中间只需跨两省。九州通衢不是白叫的。

Share this page

發佈留言

發佈留言必須填寫的電子郵件地址不會公開。 必填欄位標示為 *

這個網站採用 Akismet 服務減少垃圾留言。進一步了解 Akismet 如何處理網站訪客的留言資料