※ 引述《ttv (天晴)》之銘言:
: 專家憂 染登革熱 未發病估9萬人
: 噓 RLH: 這些專家很會危言聳聽 從SARS到後來的變種禽流感都誇大不實
我不覺得這個數字在唬爛,且聽我道來。
我剛看到這篇新聞的時候也覺得數字太誇大,但是這幾天新聞的數字跳的實在太快,我
到疾管署去查了一下資料,發現疾管署從1998年起就開始統計每週的登革熱增加病例。
台灣過去的案例絕大多數發生在高屏地區,我特別研究了南臺灣的歷史數據,歸納出了
一些結論:
一、每年的首起案例差不多都是在第 20 週左右。
二、差不多第 30 週起,每週的新增發病數會進入直線上升階段。
三、大概到第 40 週起每週新增病例的增幅會開始減緩,高峰通常落在第 43-45 這幾週
。
四、接著新病例數會快速下降,在第 52 週後就幾乎不太增加,但是會拖點小尾巴到隔年
前幾週去,再些零星的新增案例。
過去五年幾乎都是類似的模式,唯一的例外是 2013 年(高峰延至第48週,之後病例快速
下滑)。
*如下圖 (facebook)
https://www.facebook.com/photo.php?fbid=10205257297774862
*或這裡 (imgur)
http://imgur.com/BHTIyNs
這種對稱型的鐘形分布,我直覺想到四邊形固定邊長總和下最大面積為正方形的數學題目
。化學上就是 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 模式,立即把台灣能動員的化學兵都投入防疫,趁著還有點機會多救
一些人吧。