※ 引述《ttv (天晴)》之铭言:
: 专家忧 染登革热 未发病估9万人
: 嘘 RLH: 这些专家很会危言耸听 从SARS到后来的变种禽流感都夸大不实
我不觉得这个数字在唬烂,且听我道来。
我刚看到这篇新闻的时候也觉得数字太夸大,但是这几天新闻的数字跳的实在太快,我
到疾管署去查了一下资料,发现疾管署从1998年起就开始统计每周的登革热增加病例。
台湾过去的案例绝大多数发生在高屏地区,我特别研究了南台湾的历史数据,归纳出了
一些结论:
一、每年的首起案例差不多都是在第 20 周左右。
二、差不多第 30 周起,每周的新增发病数会进入直线上升阶段。
三、大概到第 40 周起每周新增病例的增幅会开始减缓,高峰通常落在第 43-45 这几周
。
四、接着新病例数会快速下降,在第 52 周后就几乎不太增加,但是会拖点小尾巴到隔年
前几周去,再些零星的新增案例。
过去五年几乎都是类似的模式,唯一的例外是 2013 年(高峰延至第48周,之后病例快速
下滑)。
*如下图 (facebook)
https://www.facebook.com/photo.php?fbid=10205257297774862
*或这里 (imgur)
这种对称型的钟形分布,我直觉想到四边形固定边长总和下最大面积为正方形的数学题目
。化学上就是 A+B -> 2B 的连锁反应,给定一个初始浓度 [A],A 碰到 B 就自己变成 B
,在密闭系统下最后会全部变成 B,反应速率为 k*[A][B]。
因此可以尝试用这个数学模型解释登革热的传播:最早的一只带原的蚊子咬了第一个人,
成为第一个(或第一批)带原者,接着这些人又被其他的蚊子咬,中标的蚊子又去咬别的
人让他们中标。所以当"带原者B" 去碰到"非带原者A"时,会变成两个带原者B,中间的媒
介是蚊子。最后系统内所有人都变成僵尸 B。
接着我把过去几年的例子用这个式子套了一下,发现可以 fit 出很吻合的曲线(如图)
,而且疫情如果愈严重,也就是病患数愈多的话,曲线与实际数值愈吻合(参考 2014 年
的例子,R^2 = 0.98)。
这边稍微解释一下,[A]+[B] 的最大上限应该是区域内的人口总数,但是因为蚊子一生当
中飞行的距离有限,也不是每个地方都有合适的水域提供蚊子世代繁衍,所以两者有个关
系是:(总人口数)*(一个小于1的常数)=(最终患者数),我把这个小于1的常数称为防疫系
数。防疫系数愈高,代表防疫做得愈差,可能得病的人数也愈高。(另外这个式子也可以
用于预测一个地区的僵尸数量,因为僵尸的传播也是类似的方式。)
所以重头戏来了,到底这次台南市的登革热疫情结束了没?
从图上看起来是没有,而且可怕的还在后面,从这四五周的上升趋势来看,最后发病的总
数可能真的会超过十万,甚至是二十万。登革热死亡率差不多有1-5%,假设用低标1%再打
五折 =0.5% 来看,十万人会死 500 人,威力远远超越八仙尘暴。
登革热一般潜伏期约一到二周,蚊子的生长周期约两周,也就是从前一个带原者被咬后到
下一个被咬者的发病日大概有2-4周的迟滞,现在是第37周,在进入最大高峰的44周还有
2-3个传播世代。以现在这个发病的规模来看根本不是地方政府有办法处里的了,我建议
政府应该比照 SARS 模式,立即把台湾能动员的化学兵都投入防疫,趁著还有点机会多救
一些人吧。