2011-2017年我國(guó)人間布魯氏菌病發(fā)病的長(zhǎng)期趨勢(shì)和季節(jié)性研究*

康育慧 崔詠梅 曹文君

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

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

資料與方法

1.資料來(lái)源

資料源于國(guó)家衛(wèi)生和計(jì)劃生育委員會(huì)公布的法定報(bào)告?zhèn)魅静∫咔樵聢?bào)表。

2.方法

(1)數(shù)據(jù)準(zhǔn)備[6]

對(duì)原始數(shù)據(jù)進(jìn)行校正,去除因各月天數(shù)不同及人口自然增長(zhǎng)率可能引起的偏差。校正公式為:

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

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

①計(jì)算移動(dòng)平均數(shù)MAi,消除季節(jié)性變動(dòng)和隨機(jī)性變動(dòng)得到循環(huán)性變動(dòng)和長(zhǎng)期趨勢(shì)序列(TC)。

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

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

④對(duì)MAi序列值(即TC)進(jìn)行曲線擬合獲得長(zhǎng)期趨勢(shì)T的方程。

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

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

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

x=∑Sicosαi/∑Si,

y=∑Sisinαi/∑Si,

其中Si為月季節(jié)變動(dòng)因素,αi為月中位角,r為集中趨勢(shì),s為角標(biāo)準(zhǔn)差。

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

3.統(tǒng)計(jì)學(xué)處理

采用excel 2007對(duì)資料進(jìn)行整理與分析;利用SPSS18.0作曲線擬合;雷氏檢驗(yàn)的檢驗(yàn)水準(zhǔn)取0.05。

結(jié)果與分析

1.數(shù)據(jù)校正結(jié)果

2011-2017年我國(guó)年平均人口自然增長(zhǎng)率約為5‰,依據(jù)公式(1+p月)12=(1+p年)求得月平均人口自然增長(zhǎng)率為0.416‰。五年內(nèi)我國(guó)共報(bào)告布病發(fā)病數(shù)344883例,經(jīng)校正后為339253.3例。

2.時(shí)間序列分解結(jié)果

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

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

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

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

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

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

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

月份月中值α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合計(jì)--12.000--1.314-2.803

3.圓形分布結(jié)果

討 論

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

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

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