河北省麻疹疫情時(shí)間序列分析.pdf

1、目的:河北省麻疹疫情在密集地粗放性控制措施下，達(dá)到了2012年消除麻疹所要求的發(fā)病率低于1/100萬(wàn)的目標(biāo)，但是成績(jī)未能保持，發(fā)病率來(lái)年即達(dá)10倍的上漲，而且上漲勢(shì)頭在2014年更加猛烈。本文旨在已有的工作經(jīng)驗(yàn)基礎(chǔ)上，建立河北省量化的麻疹疫情預(yù)測(cè)模型，設(shè)立量化指標(biāo)的預(yù)警機(jī)制，配合和指導(dǎo)現(xiàn)行的工作方法，為真正實(shí)現(xiàn)消除麻疹目標(biāo)制定更為精細(xì)有效的策略與措施提供依據(jù)。
　　方法:利用EViews8.0軟件，采用自回歸滑動(dòng)平均混合（ARIM

2、A）模型法，對(duì)河北省1979-2013年麻疹年發(fā)病率和2001至2014年第三季度的季、月和周麻疹發(fā)病數(shù)進(jìn)行分析和預(yù)測(cè)，過(guò)程主要包括模型的識(shí)別;模型中參數(shù)的估計(jì)和模型的檢驗(yàn);預(yù)測(cè)應(yīng)用三個(gè)階段。建立河北省麻疹疫情年、季、月和周四層時(shí)間序列模型，以適應(yīng)不同的工作目的和滿(mǎn)足不同工作人員的需要。將預(yù)測(cè)值的95％可信限作為麻疹疫情的預(yù)警線，并以此來(lái)評(píng)估歷史數(shù)據(jù)和警戒當(dāng)前和即將發(fā)生的疫情。
　　結(jié)果:
　　1、河北省麻疹疫情特征

3、　　河北省麻疹發(fā)病率自1950年以來(lái)整體上呈現(xiàn)波浪式下降趨勢(shì)，但是1979年之前發(fā)病率處于高位且起伏較大，而1979年及以后的發(fā)病率表現(xiàn)出比較穩(wěn)定的低位水平。每年發(fā)病均有高峰和低谷，為單峰曲線，高峰跨度主要在3-6月間，最高點(diǎn)出現(xiàn)在4、5月份，以5月份為主。即使在2011-2013年這些低發(fā)病水平年份，發(fā)病的周期性特點(diǎn)依然明顯。
　　2、年發(fā)病率ARIMA模型估計(jì)和預(yù)警建立年發(fā)病率序列yt，對(duì)檢驗(yàn)后的平穩(wěn)序列yt識(shí)別模型為ARIM

4、A(1，0，0);經(jīng)過(guò)參數(shù)、模型和殘差檢驗(yàn)，最后選定的模型為無(wú)常數(shù)項(xiàng)ARIMA(1，0，0)模型，方程式為:(1-0.73B)yt=et。預(yù)測(cè)評(píng)價(jià)指標(biāo)Theil不等式系數(shù)=0.33，BP和VP分別等于0.04和0.07，CVP=0.89，可認(rèn)為模型的擬合和預(yù)測(cè)結(jié)果良好。擬合值和實(shí)際值曲線基本吻合，預(yù)測(cè)樣本外2013年發(fā)病率0.04/10萬(wàn)，實(shí)際發(fā)病率為0.54/10萬(wàn)，落在預(yù)測(cè)值95％可信限范圍內(nèi)（9.76/10萬(wàn)），2014年發(fā)病率預(yù)

5、測(cè)值為0.39/10萬(wàn)，預(yù)警線為10.11/10萬(wàn)，與前期發(fā)病比較，呈現(xiàn)上升趨勢(shì)，需要警惕。
　　3、季發(fā)病數(shù)SARIMA模型估計(jì)和預(yù)警
　　建立季發(fā)病數(shù)序列qt，對(duì)qt序列直接擬合的SARIMA(0，0，1)(0，1，3)4模型擬合值和實(shí)際值曲線不太吻合，而且2010年末至2013年預(yù)測(cè)值出現(xiàn)負(fù)數(shù)。對(duì)qt序列對(duì)數(shù)轉(zhuǎn)換后重新擬合的模型為SARIMA(0，1，2)(0，1，3)4，方程式為:▽▽4ln(qt)=(1-0.31B

6、2)(1+1.26B4-0.32B12)et。前后模型的Theil不等式系數(shù)分別為0.42和0.21，系數(shù)成倍減小;前模型BP、VP和CVP分別為0、0.07和0.93，現(xiàn)模型分別為0.01、0.03和0.96，改進(jìn)后現(xiàn)模型的擬合和預(yù)測(cè)結(jié)果有所提升。相比前模型，現(xiàn)模型的擬合值和實(shí)際值曲線更為吻合，沒(méi)有出現(xiàn)負(fù)的預(yù)測(cè)值。預(yù)測(cè)樣本外2014年第三季度發(fā)病數(shù)249例，實(shí)際發(fā)病數(shù)182例，落在95％可信限范圍內(nèi)（897例）。2014年第四季度發(fā)病

7、數(shù)預(yù)測(cè)值為97例，預(yù)警線為326例，與前期發(fā)病比較，呈現(xiàn)下降趨勢(shì)，與往年的發(fā)病周期性趨勢(shì)基本吻合。
　　4、月發(fā)病數(shù)SARIMA模型估計(jì)和預(yù)警
　　建立月發(fā)病數(shù)序列mt，對(duì)數(shù)轉(zhuǎn)換為序列l(wèi)n(mt+1)，一階差分平穩(wěn)后模型識(shí)別為SARIMA(3，1,2)(2，1，2)12。經(jīng)過(guò)參數(shù)、模型和殘差檢驗(yàn)，最后選定的模型為SARIMA(2，1，2)(1，1，2)12，方程式為:(1-0.73B2)(1+0.73B12)▽▽12ln(m

8、t+1)=(1+0.75B2)(1+0.85B24)et。預(yù)測(cè)評(píng)價(jià)指標(biāo)Theil不等式系數(shù)=0.13，BP≈0，VP=0.05，CVP=0.95，可認(rèn)為模型擬合和預(yù)測(cè)結(jié)果良好。預(yù)測(cè)樣本外2014年9月發(fā)病數(shù)63例，實(shí)際發(fā)病數(shù)18例，落在95％可信限范圍內(nèi)（155例）。2014年10月發(fā)病數(shù)預(yù)測(cè)值為18例，預(yù)警線為45例，與前期發(fā)病比較，呈現(xiàn)下降趨勢(shì)，與往年的發(fā)病周期性趨勢(shì)基本吻合。
　　5、周發(fā)病數(shù)SARIMA模型估計(jì)和預(yù)警

9、>　　建立周發(fā)病數(shù)序列wt，對(duì)數(shù)轉(zhuǎn)換為序列l(wèi)n(wt+1)，一階差分平穩(wěn)后模型識(shí)別為SARIMA(1，1，1)(2，1，2)52。經(jīng)過(guò)參數(shù)、模型和殘差檢驗(yàn)，最后選定的模型為SARIMA(1，1，1)(2，1，2)52模型，方程式為:(1+0.24B)(1+0.83B52)(1+0.17B104)▽▽52ln(wt+1)=(1+0.33B)(1+0.87B104)et。預(yù)測(cè)評(píng)價(jià)指標(biāo)Theil不等式系數(shù)=0.11，BP≈0，VP=0.01，

10、 CVP=0.99，可認(rèn)為模型擬合和預(yù)測(cè)結(jié)果良好。預(yù)測(cè)樣本外2014年第38周發(fā)病數(shù)6例，實(shí)際發(fā)病數(shù)6例，落在95％可信限范圍內(nèi)（14例）。2014年第39周發(fā)病數(shù)預(yù)測(cè)值為7例，預(yù)警線為16例，與前期發(fā)病比較，呈現(xiàn)上升趨勢(shì)。
　　結(jié)論:
　　1、麻疹疫情四層時(shí)間序列模型的聯(lián)合使用，在一定程度上提升了預(yù)測(cè)的精度和延長(zhǎng)了預(yù)測(cè)時(shí)長(zhǎng)，彌補(bǔ)了單一模型的這種不足。
　　2、配合常規(guī)的評(píng)價(jià)方法，即時(shí)地利用麻疹時(shí)間序列模型評(píng)價(jià)現(xiàn)行的疫