南京市疾病預(yù)防控制中心(210003) 許陽(yáng)婷

【提 要】 目的 分析南京市1951-2019年細(xì)菌性痢疾流行變化特征，為細(xì)菌性痢疾的防控提供依據(jù)。方法 運(yùn)用移動(dòng)平均、joinpoint分析和小波分析方法，對(duì)南京市1951-2019年細(xì)菌性痢疾年發(fā)病率時(shí)間變化進(jìn)行分析。結(jié)果 細(xì)菌性痢疾發(fā)病率在南京市整體呈下降趨勢(shì)，5年移動(dòng)平均曲線呈現(xiàn)先臺(tái)階式后滑道樣線形下降。joinpoint分析顯示平均年度變化百分比為-5.60，細(xì)菌性痢疾流行過(guò)程中有3個(gè)轉(zhuǎn)折點(diǎn)，分別在1965年從升到降、1968年由降到略升和1984年下降。Morlet小波法時(shí)頻分析小波系數(shù)實(shí)部和模方顯示12和31為2個(gè)中心尺度，呈現(xiàn)的周期為20年和38年。小波方差顯示有2個(gè)主周期，第一主周期在31尺度。從細(xì)菌性痢疾的周期、趨勢(shì)分析結(jié)果看，目前南京市進(jìn)入細(xì)菌性痢疾高發(fā)周期中。結(jié)論 依據(jù)細(xì)菌性痢疾的變化特征，當(dāng)前應(yīng)該加強(qiáng)對(duì)細(xì)菌性痢疾的監(jiān)測(cè)，及時(shí)做好疫情處置。

細(xì)菌性痢疾(bacillary dysentery)是由志賀氏疾病桿菌引起的一種急性腸道傳染病，人群普遍易感，病后免疫力持續(xù)時(shí)間短暫，不同型別菌株間無(wú)交叉免疫，缺乏有效的疫苗免疫。其發(fā)病率居于我國(guó)法定乙類報(bào)告?zhèn)魅静〉那?位[1]。因而有流行范圍廣、傳播速度快、對(duì)健康危害大等特點(diǎn)。南京地處中國(guó)東部、長(zhǎng)江下游、北亞熱帶濕潤(rùn)氣候，春秋短、冬夏長(zhǎng)，四季分明，雨水充沛，相對(duì)濕度76%[2]。此地理氣候條件適宜志賀菌屬的繁殖與傳播。有研究顯示，近年南京市有些地區(qū)細(xì)菌性痢疾年均發(fā)病率維持在10/10萬(wàn)以上[3]。本文通過(guò)對(duì)南京市1951-2019年69年來(lái)細(xì)菌性痢疾發(fā)病率趨勢(shì)變化和周期等特征的分析，為南京市細(xì)菌性痢疾防控、合理資源配置提供理論支撐。

資料與方法

1.資料來(lái)源

1951-2000年細(xì)菌性痢疾發(fā)病率資料來(lái)源于南京市衛(wèi)生防疫站年報(bào)，2001-2003年細(xì)菌性痢疾資料來(lái)源于南京市疾病預(yù)防控制中心年報(bào)，2004-2019年病例資料和人口學(xué)資料來(lái)源于“中國(guó)疾病監(jiān)測(cè)信息系統(tǒng)”的病例報(bào)告系統(tǒng)和基本信息系統(tǒng)，發(fā)病數(shù)按發(fā)病時(shí)間導(dǎo)出。

2.方法

主要用小波分析和joinpoint方法對(duì)細(xì)菌性痢疾發(fā)病特征、趨勢(shì)、周期、突變點(diǎn)等變化進(jìn)行詳細(xì)的分析；

(1)移動(dòng)平均

指定時(shí)間段，對(duì)時(shí)間序列數(shù)據(jù)進(jìn)行移動(dòng)計(jì)算平均值。即x的n移動(dòng)平均值：

(2)joinpoint模型[4]

通過(guò)模型擬合將一個(gè)長(zhǎng)期趨勢(shì)線分成若干有統(tǒng)計(jì)學(xué)意義的趨勢(shì)區(qū)段，各段用連續(xù)的線形進(jìn)行描述[5]。y為年份，Ry為該年份發(fā)病率，表達(dá)式為：

log(Ry)=b0+b1y

b0為截距，b1為回歸系數(shù)，反映各時(shí)間段內(nèi)發(fā)病率的趨勢(shì)變化[5]。用t檢驗(yàn)判斷總體回歸系數(shù)是否為0，即趨勢(shì)變化是否有統(tǒng)計(jì)學(xué)意義。平均年度變化百分比(average annual percent change，AAPC)[5]和年度變化百分比(annual percent change，APC)[5]表達(dá)式如下：

(3)小波分析

小波分析(wavelet analysis)亦稱多分辨率分析(multiresolution analysis)[6]，用Morlet小波進(jìn)行時(shí)-頻小波分析

其中，a為頻率參數(shù)(亦稱尺度因子)，b為時(shí)間平移參數(shù)(亦稱平移因子)；Δt為取樣時(shí)間間隔；n為樣本數(shù)[10]。

(4)小波方差分析

將小波系數(shù)的平方值在b域上積分，就可得到小波方差，公式如下：

var(a)=[ωf(a，b)]2

其隨尺度a的變化過(guò)程，即為小波方差圖[11]。

3. 統(tǒng)計(jì)分析

應(yīng)用Excel 2007對(duì)數(shù)據(jù)進(jìn)行整理、制表及分析，用Joinpoint Trend Analysis Software 4.8.0.1和MATLAB 2008b軟件進(jìn)行發(fā)病率趨勢(shì)、變化點(diǎn)和周期分析；小波方差檢驗(yàn)水準(zhǔn)為α=0.01，Joinpoint檢驗(yàn)水準(zhǔn)=0.05。

表1 南京市1951-2019年細(xì)菌性痢疾發(fā)病率年內(nèi)分布情況

結(jié) 果

1.細(xì)菌性痢疾發(fā)病率年內(nèi)各月發(fā)病情況

年內(nèi)細(xì)菌性痢疾發(fā)病呈現(xiàn)夏秋季單峰，季節(jié)聚集性明顯，占年內(nèi)發(fā)病構(gòu)成的70.74%，冬春季為低發(fā)病期。

2.細(xì)菌性痢疾發(fā)病率年際發(fā)病情況

各年發(fā)病率呈下降趨勢(shì)，發(fā)病率5年移動(dòng)平均曲線呈現(xiàn)先臺(tái)階下降，后滑梯形下降。1951-1966年為高發(fā)病年，發(fā)病率在高位波動(dòng)，年平均發(fā)病率為378.95/10萬(wàn)；1970-1984年為中發(fā)病年，年平均發(fā)病率為262.00/10萬(wàn)；隨后一路下降，2014年以來(lái)發(fā)病率維持在5/10萬(wàn)上下波動(dòng)，見(jiàn)圖1。joinpoint分析顯示1951-2019年發(fā)病率AAPC為-5.60，呈下降趨勢(shì)，出現(xiàn)3個(gè)轉(zhuǎn)折點(diǎn)，分別在1965年、1968年和1984年，其中1965-1968年發(fā)病率年度變化無(wú)統(tǒng)計(jì)學(xué)意義(t=-1.20，P=0.30>0.05)，見(jiàn)圖2和表2。

圖1 南京市1951-2019年細(xì)菌性痢疾發(fā)病情況

圖2 南京市1951-2019年細(xì)菌性痢疾發(fā)病率joinpoint變化圖

表2 南京市1951-2019年細(xì)菌性痢疾發(fā)病率趨勢(shì)分析

3.細(xì)菌性痢疾年發(fā)病率序列的時(shí)頻分析

對(duì)南京市1951-2019年發(fā)病率時(shí)間序列進(jìn)行 Morlet 小波變換后的小波系數(shù)實(shí)部等值線圖，見(jiàn)圖3，圖中小波系數(shù)0等值線為突變點(diǎn)；小波系數(shù)正值等值線為高發(fā)病年份；0等值線外空白區(qū)是小波系數(shù)負(fù)值區(qū)為低發(fā)病年份。由圖3可以看出，有二個(gè)周期震蕩，分別在8～15和25～32時(shí)間尺度，尤其在25～32時(shí)間尺度上，發(fā)病率有3個(gè)偏多周期和3個(gè)偏少周期。模方圖中(圖4)在相同兩個(gè)時(shí)間尺度上也表現(xiàn)出了周期現(xiàn)象，特別在25～32時(shí)間尺度上周期表現(xiàn)的非常明顯，出現(xiàn)的時(shí)間大約為1951-2010年。兩個(gè)周期中心分別在12和31時(shí)間尺度上。從主周期不同的時(shí)間尺度變化趨勢(shì)可見(jiàn)(圖5)，在12時(shí)間尺度上，呈現(xiàn)3個(gè)周期變化特征，周期約為20年左右；在31尺度上，呈現(xiàn)兩個(gè)周期變化特征，周期約為38年左右。兩個(gè)主周期時(shí)間尺度都顯示當(dāng)前發(fā)病處在高發(fā)病區(qū)域。

圖3 南京市1951-2019年細(xì)菌性痢疾發(fā)病率小波系數(shù)實(shí)部等值線圖

圖4 南京市1951-2019年細(xì)菌性痢疾發(fā)病率小波系數(shù)模平方等值線圖

圖5 南京市1951-2019年細(xì)菌性痢疾發(fā)病率小波系數(shù)實(shí)部多時(shí)間尺度變化趨勢(shì)

4.細(xì)菌性痢疾發(fā)病率序列的主要周期分析

在小波模方差曲線出現(xiàn)了2個(gè)較明顯的峰值，見(jiàn)圖6，方差值由大到小依次對(duì)應(yīng)的尺度為31和12。其中，曲線的最大峰值對(duì)應(yīng)的尺度為31，說(shuō)明整體發(fā)病率變化31時(shí)間尺度上周期震蕩最明顯，為第一主周期尺度；結(jié)合模方圖(圖4)可見(jiàn)，31尺度上的周期震蕩在整個(gè)時(shí)間序列上表現(xiàn)非常明顯。

圖6 南京市1951-2019年細(xì)菌性痢疾發(fā)病率小波方差

討 論

南京市1951-2019年細(xì)菌性痢疾的發(fā)病率從500/10萬(wàn)左右到5/10萬(wàn)左右，經(jīng)歷了三個(gè)轉(zhuǎn)折變化點(diǎn)，分別為1965年的下降，1968年的上升和1984年的下降過(guò)程。發(fā)病情況與國(guó)民經(jīng)濟(jì)的增長(zhǎng)、居民收入的提高等社會(huì)經(jīng)濟(jì)發(fā)展因素，基本衛(wèi)生服務(wù)的水平提高及環(huán)境衛(wèi)生條件改善密切相關(guān)。這與劉東立等人的報(bào)道一致[12-13]。

流行病學(xué)常用時(shí)間序列來(lái)研究疾病發(fā)生發(fā)展的過(guò)程。時(shí)間序列是由時(shí)域和頻域組成。目前在研究疾病發(fā)生發(fā)展過(guò)程中常用的分析方法多是要求時(shí)間序列為平穩(wěn)序列，平穩(wěn)是對(duì)于頻域提出的要求，疾病的發(fā)生發(fā)展隨時(shí)間的變化往往受病原體、環(huán)境、人群和政策等多因素的綜合影響，呈現(xiàn)的疾病時(shí)間序列多為不平穩(wěn)序列。這種不平穩(wěn)除表現(xiàn)在趨勢(shì)、周期性上，還表現(xiàn)在隨機(jī)性和突變性上，形成“多尺度”結(jié)構(gòu)，具有多層次演變規(guī)律[14]。對(duì)非平穩(wěn)時(shí)間序列的研究，通常需要某一頻段對(duì)應(yīng)的時(shí)間信息，或某一時(shí)段的頻域信息。然而，單獨(dú)時(shí)域分析和頻域分析對(duì)此均無(wú)能為力[14-16]。20世紀(jì)80年代初，由Morlet提出的具有時(shí)-頻多分辨功能的小波分析，為研究時(shí)間序列問(wèn)題提供了可能，將隱藏在時(shí)間序列中的多種變化周期，通過(guò)不同時(shí)間尺度將時(shí)間序列的變化趨勢(shì)充分地反映出來(lái)，并能對(duì)未來(lái)發(fā)展趨勢(shì)進(jìn)行定性評(píng)估[17-19]。小波系數(shù)實(shí)部等值線圖是在時(shí)域中呈現(xiàn)出不同時(shí)間尺度的年發(fā)病率時(shí)間序列周期變化情況，可以用于推斷年發(fā)病率的未來(lái)變化趨勢(shì)。模方相當(dāng)于小波能量譜，能量譜展示出不同時(shí)間尺度的周期變化。從小波方差分析主周期趨勢(shì)圖中可以得到不同的時(shí)間尺度下[19]，發(fā)病率存在的平均周期及高-低發(fā)病率情況的變化特征。