應(yīng)用趨勢(shì)季節(jié)模型分析我國(guó)狂犬病發(fā)病動(dòng)態(tài)

李艷琳武玉欣

應(yīng)用趨勢(shì)季節(jié)模型分析我國(guó)狂犬病發(fā)病動(dòng)態(tài)

李艷琳1武玉欣2△

目的分析我國(guó)2009-2013年狂犬病報(bào)告病例數(shù)據(jù)，應(yīng)用趨勢(shì)季節(jié)模型預(yù)測(cè)2014年全國(guó)狂犬病的發(fā)病情況。方法采用趨勢(shì)季節(jié)模型對(duì)數(shù)據(jù)進(jìn)行擬合分析并建立模型，具體計(jì)算借助SPSS15.0和excel軟件完成。結(jié)果我國(guó)大陸地區(qū)2009-2013年狂犬病的報(bào)告發(fā)病數(shù)呈逐年下降趨勢(shì)，研究建立的趨勢(shì)季節(jié)模型為：yt′＝（1766.339-20.168t）st（方差分析結(jié)果F＝33.734，P＜0.001），模型有統(tǒng)計(jì)學(xué)意義，模型預(yù)測(cè)2014年狂犬病發(fā)病例數(shù)，預(yù)測(cè)相對(duì)誤差平均為16.97%。結(jié)論趨勢(shì)季節(jié)模型可用于我國(guó)狂犬病發(fā)病情況的短期預(yù)測(cè)。

狂犬病 疫情 趨勢(shì)季節(jié)模型

狂犬病（rabies）是由狂犬病病毒（RV）引起的乙類傳染病，又稱恐水癥。一旦感染該病，病死率近達(dá)100%，是迄今為止人類病死率最高的急性傳染病［1］。狂犬病是世界性疾病，目前有100多個(gè)國(guó)家和地區(qū)存在該病。每年，全球近7萬(wàn)人死于狂犬病，平均每10分鐘就有一人因該病死亡。95%以上的狂犬病人間病例系感染犬咬傷而致［2］。我國(guó)是狂犬病的高發(fā)地區(qū)，感染狂犬病死亡人數(shù)居世界第二，僅次于印度［3］。因此，應(yīng)及時(shí)掌握我國(guó)狂犬病的流行特征及其變化規(guī)律，切實(shí)加強(qiáng)狂犬病防控研究工作，最終達(dá)到控制直至消滅狂犬病的目標(biāo)。本文收集2009-2013年全國(guó)范圍狂犬病疫情資料，擬合建立其趨勢(shì)季節(jié)模型，進(jìn)行動(dòng)態(tài)分析研究和預(yù)測(cè)，旨在了解狂犬病近年的流行規(guī)律，為疾病的防控提供參考。

資料與方法

1.資料來(lái)源

收集整理2009-2014年《中華人民共和國(guó)衛(wèi)生部公報(bào)》及《中華人民共和國(guó)衛(wèi)生和計(jì)劃生育委員會(huì)公報(bào)》轉(zhuǎn)載的歷年我國(guó)狂犬病的逐月發(fā)病疫情資料［4］。

2.方法

狂犬病的發(fā)病受自然因素、社會(huì)因素、環(huán)境因素等因素的影響，呈現(xiàn)出季節(jié)性和趨勢(shì)性變化。趨勢(shì)季節(jié)模型的建立既考慮了事件發(fā)生的趨勢(shì)性，也考慮了季節(jié)性的影響，擬合的結(jié)果較為接近實(shí)際［5］。本研究應(yīng)用趨勢(shì)季節(jié)模型對(duì)狂犬病的月發(fā)病數(shù)資料進(jìn)行模型擬合。數(shù)據(jù)的分析計(jì)算利用SPSS15.0和ExcelXP完成。

結(jié) 果

1.2009-2013年我國(guó)狂犬病疫情的一般情況

5年間，我國(guó)大陸地區(qū)共報(bào)告狂犬病病例9040例，歷年逐月報(bào)告發(fā)病人數(shù)詳見(jiàn)圖1和表1。全國(guó)狂犬病發(fā)病人數(shù)既有逐年遞減的長(zhǎng)期趨勢(shì)，每年的發(fā)病又有季節(jié)性特點(diǎn)。

圖1 2009-2013年全國(guó)狂犬病逐月報(bào)告發(fā)病人數(shù)

表1 2009-2013年全國(guó)狂犬病逐月報(bào)告發(fā)病數(shù)統(tǒng)計(jì)

2.趨勢(shì)季節(jié)模型分析過(guò)程及結(jié)果

（1）發(fā)病人數(shù)的趨勢(shì)季節(jié)模型的建立

首先將2009-2013年60個(gè)月份順次編號(hào)（t），2005年1月t＝1，2005年2月t＝2，……，2009年12月t＝60。擬合月報(bào)告發(fā)病人數(shù)（yt）和月份t之間的直線回歸方程：

該回歸方程有統(tǒng)計(jì)學(xué)意義（方差分析結(jié)果F＝33.734，P＜0.000）。依方程求2009-2013年各月的估計(jì)值。各月實(shí)際報(bào)告發(fā)病人數(shù)與估計(jì)發(fā)病人數(shù)的比值即為該月的季節(jié)指數(shù)。相同月份季節(jié)指數(shù)的平均值即調(diào)整季節(jié)指數(shù)st。所建模型1～12月的調(diào)整季節(jié)指數(shù)見(jiàn)表2第（4）欄。

利用以上分析建立的直線回歸方程（1）和各月的調(diào)整季節(jié)指數(shù)，建立趨勢(shì)季節(jié)模型如下：

其中，yt′預(yù)測(cè)月份的估計(jì)發(fā)病人數(shù)，t為預(yù)測(cè)月份的序號(hào)，st為預(yù)測(cè)月份對(duì)應(yīng)的調(diào)整季節(jié)指數(shù)。

（2）2014年狂犬病發(fā)病人數(shù)的預(yù)測(cè)

2014年1～12月對(duì)應(yīng)的月份序號(hào)為61～72，將各月份的序號(hào)和對(duì)應(yīng)的調(diào)整季節(jié)指數(shù)帶入方程（2），即可求得2014年各月的發(fā)病人數(shù)，見(jiàn)表2第（5）欄。

表2 2014年我國(guó)大陸地區(qū)狂犬病逐月發(fā)病人數(shù)預(yù)測(cè)結(jié)果

2014年1～12月份全國(guó)狂犬病報(bào)告的發(fā)病人數(shù)見(jiàn)表2，預(yù)測(cè)相對(duì)誤差為1.15%～40.63%，平均誤差僅為16.97%。

討 論

趨勢(shì)季節(jié)模型計(jì)算相對(duì)簡(jiǎn)單、易于理解，模型兼顧了資料的趨勢(shì)性和季節(jié)性的特點(diǎn)。對(duì)既有一定的長(zhǎng)期趨勢(shì)，又同時(shí)具有明顯季節(jié)周期性的疾病流行趨勢(shì)的短期預(yù)測(cè)效果較好。2009-2013年全國(guó)狂犬病報(bào)告發(fā)病數(shù)據(jù)資料顯示，呈逐年下降趨勢(shì)，年內(nèi)各月份的發(fā)病具有明顯的季節(jié)性特點(diǎn)。因此，可應(yīng)用趨勢(shì)季節(jié)模型擬合其發(fā)病模型，并進(jìn)行回代和預(yù)測(cè)。

根據(jù)建立模型預(yù)測(cè)2014年月發(fā)病數(shù)的相對(duì)誤差17%左右，預(yù)測(cè)效果較為理想。由于影響疾病發(fā)生的因素非常復(fù)雜，包括自然因素、環(huán)境因素、社會(huì)因素等，而該模型僅將時(shí)間因素作為所有影響因素的一個(gè)綜合因子，引入模型，具體的影響因素在模型中無(wú)體現(xiàn)，因此趨勢(shì)季節(jié)模型可作為短期預(yù)測(cè)模型，其長(zhǎng)期預(yù)測(cè)效果并不理想。具體應(yīng)用該模型時(shí)，要根據(jù)實(shí)際疫情的變化，判斷是否還適合模型擬合，并及時(shí)根據(jù)最新數(shù)據(jù)更新模型參數(shù)。

1.唐丹丹，朱駿，曹向英等.狂犬病研究進(jìn)展.畜牧獸醫(yī)科技信息，2014（11）：4-6.

2.Bernard Vallat.No more Deaths from Rabies.OIE Bulletin，2014（3）：1-2.

3.蒙建作.狂犬病的流行因素分析與防制策略.醫(yī)學(xué)動(dòng)物防制，2013，29（4）：386-390.

4.中華人民共和國(guó)衛(wèi)生部.2009-2014年中華人民共和國(guó)衛(wèi)生部公報(bào). http：／／www.nhfpc.gov.cn／mohbgt／pwsbgb／wsbgb_index.shtml.

5.王瑞平.兩種趨勢(shì)季節(jié)模型預(yù)測(cè)法的應(yīng)用比較.中國(guó)衛(wèi)生統(tǒng)計(jì)，2011，28（1）：77-79.

（責(zé)任編輯：張 悅）

1.沈陽(yáng)市和平區(qū)疾病預(yù)防控制中心（110003）

2.中國(guó)醫(yī)科大學(xué)（110122）

△通信作者：武玉欣，E-mail：yxwu＠m(xù)ail.cmu.edu.cn