康育慧 崔詠梅 曹文君

【提 要】 目的 了解我國人間布魯氏菌病(人布病)發(fā)病的長期趨勢和季節(jié)性特征,為布病防控提供科學依據。方法 收集2011年1月-2017年12月國家衛(wèi)計委發(fā)布的人布病逐月發(fā)病數據,首先利用時間序列法分解出人布病月發(fā)病數的季節(jié)性變動因素和長期趨勢因素,并對長期趨勢進行曲線擬合,然后構建基于季節(jié)性變動因素的圓形分布法并分析人布病的季節(jié)性特征。結果 2011-2017年我國人布病月發(fā)病數的長期趨勢因素的最佳擬合方程為:T=3988.01-85.983t+4.117t2-0.041t3(R2=0.956),發(fā)病數先緩慢上升后下降;圓形分布顯示我國人布病發(fā)病具有季節(jié)性集中趨勢(P<0.01),發(fā)病高峰日為6月7日,流行高峰期為3月3日-9月11日。結論 2011-2017年我國人布病發(fā)病例數呈緩慢上升后下降的趨勢,發(fā)病呈季節(jié)性,春夏季高發(fā),應在冬季做好布病防控工作。

布魯氏菌病是由布魯氏菌引起的人獸共患性傳染病,人感染布魯氏菌后在臨床上表現為長期發(fā)熱、多汗、關節(jié)疼痛及肝脾腫大等。 2000年以后,人間布魯氏菌病(人布病)在我國成為報告發(fā)病數上升速度最快的傳染病之一[1],成為影響人民群眾身體健康的重要公共衛(wèi)生問題。分析疾病的時間分布特征,提前做好疾病的預防、控制等準備工作,對降低疾病的發(fā)病,保障居民健康具有重要意義。目前,關于人布病的時間分布研究主要探討了我國部分地區(qū)或省市的疾病的季節(jié)性特征[2-5],缺少全國范圍人布病季節(jié)性的統(tǒng)計學分析;同時,缺乏對人布病連續(xù)多年發(fā)病變化規(guī)律的長期趨勢的系統(tǒng)研究。本研究采用時間序列分解法和基于季節(jié)性變動因素的圓形分布法，對2011-2017年我國人布病發(fā)病的長期趨勢和季節(jié)性進行探討,旨在了解其發(fā)病規(guī)律,為我國人布病的防控決策提供科學依據。

資料與方法

1.資料來源

資料源于國家衛(wèi)生和計劃生育委員會公布的法定報告?zhèn)魅静∫咔樵聢蟊怼?/p>

2.方法

(1)數據準備[6]

對原始數據進行校正,去除因各月天數不同及人口自然增長率可能引起的偏差。校正公式為:

(2)時間序列分解法[6-8]

傳染病時間序列通常受到長期趨勢(T)、季節(jié)性變動(S)、循環(huán)性變動(C)和隨機性變動(R)四個因素的影響。為測量人布病發(fā)病的長期趨勢和季節(jié)性變動,本研究采用乘法模型(Xi=T×C×S×R)對該病逐月發(fā)病數據資料進行分解。具體步驟如下:

①計算移動平均數MAi,消除季節(jié)性變動和隨機性變動得到循環(huán)性變動和長期趨勢序列(TC)。

②計算SRi,消除循環(huán)性變動和長期趨勢因素得到包含季節(jié)性因素和隨機性因素的序列(SR)

SRi=Xi/MAi=Xi/(T×C)=T×C×S×R/(T×C)=S×R

④對MAi序列值(即TC)進行曲線擬合獲得長期趨勢T的方程。

(3)圓形分布法[2-5，9]

將分解得出的季節(jié)性變動因素利用圓形分布法分析人布病發(fā)病的季節(jié)性特征。具體步驟如下:

①計算月中位角αi:將一年看成一個圓是360°,則一天約為0.986°,以每月的中間一天為月中值折算為角度即月中位角。

x=∑Sicosαi/∑Si,

y=∑Sisinαi/∑Si,

其中Si為月季節(jié)變動因素,αi為月中位角,r為集中趨勢,s為角標準差。

Z=nr2，n為校正的發(fā)病總數。

3.統(tǒng)計學處理

采用excel 2007對資料進行整理與分析;利用SPSS18.0作曲線擬合;雷氏檢驗的檢驗水準取0.05。

結果與分析

1.數據校正結果

2011-2017年我國年平均人口自然增長率約為5‰,依據公式(1+p月)12=(1+p年)求得月平均人口自然增長率為0.416‰。五年內我國共報告布病發(fā)病數344883例,經校正后為339253.3例。

2.時間序列分解結果

(1)計算MAi序列值并利用SPSS18.0軟件對其進行曲線擬合,獲得長期趨勢T的方程:

T=3988.01-85.983t+4.117t2-0.041t3(R2=0.886)

趨勢方程曲線顯示:2011-2017年我國人布病發(fā)病數呈先緩慢上升后下降的趨勢(圖1)。

圖1 我國人布病季節(jié)指數、TC序列及曲線擬合的TC序列對比

(2)計算出規(guī)范化的季節(jié)因素S(表1),由曲線S可見:我國布病發(fā)病呈單峰型,每年6月份病例數最多,之后迅速下降,11月份時出現一次微小上升波動,其后繼續(xù)下降(圖1)。

圖2 2011-2017年人布病校正發(fā)病數與模擬校正發(fā)病數對比

圖3 2011-2017年人布病原始發(fā)病數與模擬原始發(fā)病數對比

月份月中值αiSisinαicosαiSisinaαiSicosαi115.515.290.5190.2640.9650.1370.501 245.044.380.6860.6700.7150.4590.490375.073.481.1000.959 0.2841.0550.3124105.0103.561.2970.972-0.2341.260-0.303 5135.5133.641.5060.724-0.6901.090-1.039 6166.0163.731.5880.280-0.9600.445-1.5247196.5193.811.440-0.239-0.971-0.344-1.3988227.5224.381.193-0.699-0.715-0.834-0.8539258.0254.470.814-0.964-0.268-0.785-0.218 10288.5284.500.552-0.968 0.251-0.535 0.13911319.0314.630.648-0.712 0.703-0.461 0.45512349.5344.710.658-0.264 0.965-0.174 0.635合計--12.000--1.314-2.803

3.圓形分布結果

討 論

對疾病連續(xù)數年的觀察,分析疾病發(fā)病或死亡的長期趨勢,不僅能體現疾病的動態(tài)變化情況,還可以反映疾病的致病因素和防治對策的綜合作用效果。本研究應用時間序列分解法探索2011-2017年我國人布病發(fā)病的長期趨勢。經模擬發(fā)病數與實際發(fā)病情況比較發(fā)現該模型擬合效果良好,表明該分解法較科學地反映了客觀實際。研究結果表明，經過人口自然增長率校正后,我國人布病發(fā)病數在2011-2014年呈緩慢上升趨勢,2015年開始緩慢下降,2016年和2017年仍保持下降趨勢。前四年發(fā)病人數的升高可能與進入21 世紀布病在世界范圍內肆虐有關[1];2012年5月國務院辦公廳提出了全國布病防治的2015和2020年的階段性目標,為控制人布病的流行開展了大量防治工作,經過三年努力到2015年我國人布病防治工作初見成效,而且根據長期趨勢擬合曲線推測我國人布病在短期內仍會繼續(xù)下降。

在研究疾病的季節(jié)性特征時,如果不剔除其他非季節(jié)性因素影響而直接使用原始病例數會導致誤差[8,10],時間序列分解雖可剔除非季節(jié)因素,得到反映疾病季節(jié)性強弱的季節(jié)性變動因素,但是缺少證明其具有顯著性意義的統(tǒng)計指標[6,8]。本研究首先利用時間序列分解法剔除其他非季節(jié)性因素,得到人布病月發(fā)病數的季節(jié)性變動因素,然后應用基于季節(jié)性變動因素的圓形分布法分析該病發(fā)病的季節(jié)性,最后對季節(jié)性集中趨勢進行統(tǒng)計學檢驗,增加了結果的準確性[8]和客觀性。基于季節(jié)性變動因素的圓形分布結果顯示，2011-2017年七年間我國人布病發(fā)病有一定季節(jié)性,春夏高發(fā),發(fā)病高峰日為6月7日,流行高峰期為3月3日-9月11日。與李錫太等[2 - 4]研究的北京、山東等華北地區(qū)相比,全國范圍內的人布病發(fā)病高峰日有所推后,流行高峰期跨度略長;與趙雪松[5]研究的內蒙古通遼地區(qū)比較,高峰日一致,但全國范圍內的人布病流行高峰期縮短。這可能與全國各地布病流行菌種、菌型、氣候及各牧區(qū)管理與政府重視程度等諸多因素有關[11]。

世界上已經有多個國家成功根除了布病,2011-2017年我國人布病數雖然在經歷上升后于近兩年有所下降,但形勢依然嚴峻。布病發(fā)病具有一定季節(jié)性,因此有關疾控部門應在布病高發(fā)的春夏時節(jié)強化對動物及其產品的生產、加工、流通等環(huán)節(jié)的監(jiān)督檢查,以降低發(fā)病率;同時應于冬季加強對農牧民,尤其是布病聚集地區(qū)養(yǎng)殖人員防控知識的宣傳和培訓,增強防范意識,以便有效地控制疾病的流行和暴發(fā)。