丁言穩(wěn),劉 婷,蔡 奔,賈賢杰

(蚌埠醫(yī)科大學(xué) 公共衛(wèi)生學(xué)院 流行病與衛(wèi)生統(tǒng)計(jì)學(xué)教研室,安徽 蚌埠 233030)

布魯菌病(以下簡稱布病)是一種由布魯菌群侵入機(jī)體后通過淋巴途徑傳播到血液,再通過血液導(dǎo)致全身感染引起的人畜共患的傳染變態(tài)反應(yīng)性疾病[1]。該病病人負(fù)重骨關(guān)節(jié)最易受累,脊柱最常見,其次是髖、膝、骶髂、肩和踝關(guān)節(jié)等,同時(shí)可累及神經(jīng)系統(tǒng)[2]。急性期以發(fā)熱、多汗、乏力和關(guān)節(jié)痛為主要表現(xiàn),可重復(fù)感染或轉(zhuǎn)為慢性,嚴(yán)重者可致殘疾。同時(shí)布病發(fā)病癥狀缺乏特異性,臨床診斷鑒別難度大,易發(fā)生誤診誤治[2]。21世紀(jì)以來,我國布病發(fā)病率持續(xù)上升,自2006年始,布病穩(wěn)居甲乙類傳染病報(bào)告發(fā)病數(shù)前10位,布病疫情呈現(xiàn)由北向南從傳統(tǒng)牧區(qū)、半牧區(qū)向農(nóng)區(qū)擴(kuò)散,成為危及居民健康的重要公共衛(wèi)生問題之一[3]。本文在空間流行病學(xué)基礎(chǔ)上回顧性分析1950-2015年布病趨勢,通過總結(jié)布病趨勢分析,以期加深臨床醫(yī)生對(duì)布病認(rèn)知,減少誤診率,做到早診斷、早治療,改善預(yù)后,降低布病慢性化比例[4-5],減輕群眾健康的損害。現(xiàn)作報(bào)道。

1 資料與方法

1.1資料來源

數(shù)據(jù)資料來自國家人口健康共享杯服務(wù)平臺(tái)公共衛(wèi)生科學(xué)數(shù)據(jù)中心提供的《1950-2015年布病數(shù)據(jù)》。使用資料內(nèi)容包括1950-2009年全國各年、各月布病發(fā)病人數(shù)、死亡人數(shù)、發(fā)病率和死亡率;2010-2015年全國各省(自治區(qū)、直轄市)布病月發(fā)病率。

1.2方法

采用Excel2010構(gòu)建1950-2015年布病數(shù)據(jù)數(shù)據(jù)庫,并繪制60年間布病發(fā)病率、死亡率和月平均發(fā)病例數(shù)趨勢圖作直觀分析。通過Arcgis10.2軟件空間統(tǒng)計(jì)工具中的聚類分析模塊,對(duì)2000-2015年布病發(fā)病率作空間自相關(guān)分析,其中,空間自相關(guān)分析包括全局空間自相關(guān)分析和局部空間自相關(guān)分析。

1.2.1 全局空間自相關(guān)分析 該分析是對(duì)整個(gè)區(qū)域?qū)傩灾悼臻g特征的描述,表示同一個(gè)變量在不同空間位置上的相關(guān)性[6-7]。本文通過全局自相關(guān)分析描述全國各省(區(qū)、市)布病發(fā)病的總體分布情況,反映指標(biāo)為全局空間自相關(guān)系數(shù)Moran′sI,其取值一般在-1～1,>0表示正相關(guān),<0表示負(fù)相關(guān),=0表示不相關(guān);Moran′sI值越趨近于1,總體空間差異越小,值越趨近于-1,總體空間差異越大[8]。Moran′sI計(jì)算公式如下:

1.2.2 局部空間自相關(guān)性分析 反映指標(biāo)為LocalMoran′sI,能夠反映區(qū)域單元屬于高值集聚或低值集聚的空間分布,空間關(guān)聯(lián)模式可分為4種類型:高高關(guān)聯(lián)、高低關(guān)聯(lián)、低高關(guān)聯(lián)和低低關(guān)聯(lián)[8]。LocalMoran′sI計(jì)算公式如下:

式中n是參與分析的空間單元數(shù),即全國各省市數(shù);xi和xj分別表示在i省(區(qū)、市)和j省(區(qū)、市)的發(fā)病率,i≠j;wij是空間權(quán)重矩陣,表示i省(區(qū)、市)和j省(區(qū)、市)的鄰近關(guān)系[7]。

2 結(jié)果

2.11950-2009年布病發(fā)病及死亡分布

1950-2009年60年間,我國共報(bào)告布病發(fā)病927 987例,死亡547例,年平均發(fā)病率為0.201 2(1/100 000),年平均死亡率為0.000 3(1/100 000)。其中1950-1979年間布病發(fā)病率波動(dòng)幅度較大;1980-1998年間布病發(fā)病率較為平穩(wěn),無明顯的升高或降低;1999-2009年間布病發(fā)病率持續(xù)上升(見圖1)。死亡率方面,1950-1999年間我國布病死亡率相對(duì)較高;2000-2009年間雖發(fā)病率升高,但死亡率較低。其中布病發(fā)病致死較嚴(yán)重年份包括1957年(0.002)、1960年(0.002)、1964年(0.001)、1967年(0.004)、1973年(0.001)、1976年(0.001)(見圖2)。

2.22000-2015年布病發(fā)病率地區(qū)分布

通過2000-2015年全國各省(市)布病發(fā)病率標(biāo)準(zhǔn)化處理結(jié)果可見,我國布病發(fā)病率地區(qū)分布總體趨勢較為平穩(wěn),其中發(fā)病率較高的為內(nèi)蒙古、黑龍江、山西、新疆等省、自治區(qū),其次為吉林、河北、河南、遼寧等省,發(fā)病率較低的為安徽、江西、湖南、廣東、廣西、貴州、云南、四川等省、自治區(qū)。

2.3布病的空間分布

2.3.1 布病全局空間自相關(guān)分析 基于Arcgis10.2統(tǒng)計(jì)軟件空間分析功能,計(jì)算2000-2015年全國布病年發(fā)病率全局Moran′sI值,結(jié)果顯示,2000-2015年我國布病發(fā)病率呈現(xiàn)空間負(fù)相關(guān)性,其中2008-2013年全局空間自相關(guān)指數(shù)差異有統(tǒng)計(jì)學(xué)意義(P<0.05),表明其發(fā)病情況的全局空間分布關(guān)系是離散模式;2000-2007年、2014-2015年我國布病發(fā)病情況的全局空間分布關(guān)系則為隨機(jī)模式(P>0.05)(見表1)。

表1 2000-2015年布病的全局空間自相關(guān)分析

2.3.2 布病局部空間自相關(guān)分析 對(duì)2000-2015年我國各省(區(qū)、市)布病發(fā)病率進(jìn)行聚類和異常值分析,結(jié)果顯示,高高空間關(guān)聯(lián)模式存在于山西(2000年)、新疆(2014-2015年)和寧夏(2014-2015年);2002-2003年、2005-2013年內(nèi)蒙古均存在高低空間關(guān)聯(lián)模式;2001-2003年、2013-2014年山西存在高低空間關(guān)聯(lián)模式;2004年西藏存在高低空間關(guān)聯(lián)模式(見表2)。

表2 2000-2015年全國各省(區(qū)、市)布病發(fā)病率聚類和異常值分析

3 討論

布病是由布魯菌侵入機(jī)體引起的一種人畜共患的傳染變態(tài)反應(yīng)性疾病[9],是《中華人民共和國傳染病防治法》規(guī)定的乙類傳染病,近年我國布病流行呈逐年上升趨勢,已成為嚴(yán)重的公共衛(wèi)生問題[10]。本研究通過時(shí)間、地區(qū)和空間三個(gè)維度探究我國布病流行趨勢,著重采用空間自相關(guān)分析法描述我國布病發(fā)病的空間分布狀況。結(jié)果顯示,我國1950-2009年60年間,共報(bào)告布病發(fā)病927 987例,死亡547例,年平均發(fā)病率為0.2012(1/100 000),年平均死亡率為0.000 3(1/100 000),1999-2009年間我國布病發(fā)病率持續(xù)上升,年平均增長率為0.232 1(1/100 000),2009年發(fā)病率高達(dá)2.696 9(1/100 000)。這與國內(nèi)相關(guān)研究結(jié)論一致,LI等[11]明確從1999年開始我國布病發(fā)病急劇上升;白永飛等[12]研究指出,國家相關(guān)干預(yù)政策與布病疫情發(fā)展有關(guān),影響布病的危險(xiǎn)因素雖持續(xù)存在,但居民就診意識(shí)和醫(yī)療機(jī)構(gòu)報(bào)告敏感性提高,導(dǎo)致布病報(bào)告發(fā)病率迅速上升。我國布病報(bào)告發(fā)病例數(shù)平均每月約370例,1-6月逐月上升,6-11月逐月減少,布病發(fā)病高峰期是3-8月。這與楊麗等[13-14]研究一致,布病病例報(bào)告主要集中在春夏季節(jié),可能與季節(jié)放牧、羔羊接生等季節(jié)性活動(dòng)有關(guān)。

本研究顯示,2000-2015年我國布病發(fā)病率地區(qū)分布總體較為穩(wěn)定,其中布病發(fā)病率較高的是內(nèi)蒙古、黑龍江、山西、新疆等省和自治區(qū)。這與國內(nèi)相關(guān)研究一致,如LIN等[15]在文章中指出我國北方布病發(fā)病率較南方高;同時(shí)有研究[16]基于2009年全國布病檢測數(shù)據(jù)指出內(nèi)蒙古、山西省位于我國布病發(fā)病率前2位。本研究結(jié)果顯示,我國布病發(fā)病呈現(xiàn)全局空間負(fù)相關(guān)性,高低空間關(guān)聯(lián)模式多存在于內(nèi)蒙古,高高空間關(guān)聯(lián)模式多存在于山西、新疆和寧夏。ZHANG等[17]對(duì)我國布病2004-2010布病發(fā)病率進(jìn)行空間分析發(fā)現(xiàn),高風(fēng)險(xiǎn)地區(qū)主要集中在內(nèi)蒙古地區(qū)及其周圍以北的一些省份,這與本文的研究結(jié)果一致。

本文采用的數(shù)據(jù)為全國布病(粗)發(fā)病率,無發(fā)病人群信息并且有部分?jǐn)?shù)據(jù)缺失,這會(huì)導(dǎo)致進(jìn)行空間自相關(guān)分析時(shí)存在一定偏倚,使計(jì)算結(jié)果與實(shí)際情況存在一定偏差。此外,本研究采用的Moran′sI統(tǒng)計(jì)量容易受空間聚集區(qū)域大小的影響[18],可能會(huì)對(duì)真實(shí)結(jié)果的計(jì)算帶來一定影響。今后的研究中有待進(jìn)一步完善布病病例數(shù)據(jù),如性別、年齡、職業(yè)等指標(biāo);同時(shí)考慮布病患病相關(guān)因素,如溫度、濕度等生態(tài)指標(biāo)對(duì)布病發(fā)病的影響。