杭州市人感染H7N9禽流感流行強度影響因素的兩水平模型分析.docxVIP

1、杭州市人感染H7N9禽流感流行強度影響因素的兩水平模型分析基金項目：杭州市甄大科技創(chuàng)新專項(20131813A07),杭州市衛(wèi)生科技重點項目(20I3Z10)1. 杭州市疾病預防控制中心(310021)2. 杭州市氣象局通信作者：丁華.E-mail：dh.dhua999黃春萍'丁，華'楊旭輝'劉牧文|王妹'孫晝'考慶君'吳彥'王兵I宋姝娟I葉春2金成2【提要】目的為杭州市制定人感染H7N9禽流感的防控策略提供參考。方法收集杭州市2013-2016年各地區(qū)人感染H7N9禽流感病例數(shù)、當?shù)卣块T出臺的有關活禽交易市場的政策、外環(huán)境監(jiān)測情況以

2、及氣象因素，利用ML-wiN2.25軟件擬合兩水平方差成分模型，探討不同地區(qū)影響人感染H7N9禽流感流行強度的因素c結果關閉活禽交易的政策實施、外環(huán)境監(jiān)測陽性率以及最低氣溫、當月降雨量和當月日照時數(shù)都成為影響人感染H7N9禽流感流行強度的因素。結論關閉活禽交易的政策實施在很大程度上降低了主城區(qū)和蕭山區(qū)疫情的發(fā)生水平。外環(huán)境監(jiān)測的陽性率與疫情的發(fā)生水平成正比，提示在高發(fā)季節(jié)加強外環(huán)境監(jiān)測有利于控制疫情發(fā)生規(guī)模?！娟P鍵詞】人感染H7N9禽流感流行強度影響因素兩水平統(tǒng)計模型人感染H7N9禽流感作為新發(fā)傳染病，于2013年2月首次出現(xiàn)在我國長三角地區(qū)Hi,杭州市2013年3月30日報告了首例人感染H7

3、N9禽流感病例，之后每年均有疫情發(fā)生，截止2016年底共報告88例人感染H7N9禽流感病例，死亡28例,病死率為31.82%。為有效防控人感染H7N9禽流感疫情,特開展杭州市人感染H7N9禽流感流行強度的影響因素分析，以期為早期預警和及時采取防控措施提供參考。對象與方法1. 疫情資料杭州市2013-2016年分地區(qū)分月人感染H7N9禽流感病例數(shù)來源于中國疾病預防控制信息系統(tǒng)。杭州市2013-2016年外環(huán)境H7N9樣本陽性率來源于杭州市疾病預防控制中心H7N9外環(huán)境監(jiān)測。2. 氣象資料杭州市2013-2016年分地區(qū)分月最低氣溫(P)、月度降雨量(毫米)及月度日照時數(shù)(小時)均來源于杭州市氣象

4、局。3. 研究方法將影響因素歸納為年度影響因素和月度影響因素兩個水平,以年度為水平2單位，月度為水平1單位,病例數(shù)為應變量,解釋變量包括當?shù)厥欠癯雠_關閉活禽交易市場的政策、外環(huán)境監(jiān)測陽性率、月度最低氣溫、月度降雨量、月度日照時數(shù)等，利用MLwiN2.25軟件擬合兩水平方差成分模型。變量篩選采用Wald檢驗，檢驗水準。定為0.050兩水平方差成分模型基本結構為:Zy=%+Pixa+eoa式中i=1,2,nt表示水平1單位;J=1,2，nj表示水平2單位;.為截距;但為解釋變量的回歸系數(shù);e回為隨機誤差項，即水平1單位的隨機誤差。兩水平方差成分模型建立在兩個假設的基礎上：(1) E(ev)=0,v

5、ar(ev)=兄(2) cov(u0,e0v)=0因此,兩水平方差成分模型也可以表示為坊=(角+&與)+(uOJ+eOij)其中場和&描述模型的固定效應(fixedeffects)；和廿。描述模型的隨機效應(randomeffects)，分別是水平2單位的方差成分和水平1單位的方差成分。結果1. 人感染H7N9禽流感頻數(shù)分布杭州市自2013年報告首例H7N9禽流感病例以來，每年均有病例報告，發(fā)病可分為兩個階段:第一階段為2013-2014年，該階段發(fā)病水平較高，發(fā)病地區(qū)以主城區(qū)和郊區(qū)為主;第二階段為2015年之后，該階段發(fā)病水平明顯下降，發(fā)病地區(qū)轉(zhuǎn)變?yōu)橐越紖^(qū)和郊縣為主(表1)。

6、從分月發(fā)病圖(圖1)上可見，人感染H7N9禽流感病例發(fā)病集中于冬春季節(jié)。表1杭州市各地區(qū)2013-2016年H7N9禽流感發(fā)病情況地區(qū)名稱地區(qū)分類發(fā)病數(shù)2013年2014年2015年2016年上城區(qū)主城區(qū)8300下城區(qū)主城區(qū)4000江干區(qū)主城區(qū)5220拱墅區(qū)主城區(qū)1111西湖區(qū)主城區(qū)3011濱江區(qū)主城區(qū)0300蕭山區(qū)郊區(qū)51123余杭區(qū)郊區(qū)2626富陽區(qū)郊區(qū)，0222淳安縣郊縣0001建德市郊縣0100桐廬縣郊縣0111臨安市郊縣1012合計293012172. 人感染H7N9禽流感影響因素分析r以年度為水平2單位，月度為水平1單位,病例數(shù)為應變量（case）,解釋變量包括當?shù)厥欠駥嵤╆P閉活禽

7、交易市場的政策（close）、外環(huán)境監(jiān)測陽性率（rate）、月度最低氣溫（temp）、月度降雨量（rain）、月度日照時數(shù)（sun）等（表2）。由于杭州市主城區(qū)關閉活禽交易市場的時間統(tǒng)一，故將上城區(qū)、下城區(qū)、西湖區(qū)、拱墅區(qū)、江干區(qū)、濱江區(qū)等6個區(qū)合并為主城區(qū)進行擬合,桐廬縣、淳安縣、建德市、臨安市等4個郊縣因病例數(shù)偏少而合并為郊縣地區(qū)進行擬合，因此獲得杭州市主城區(qū)、蕭山區(qū)、余杭區(qū)、富陽區(qū)、郊縣地區(qū)共5類地區(qū)的人感染H7N9禽流感影響因素分析。表2變量賦值說明指標名稱變量名定義及賦值應變量caseNO的整數(shù)年度year水平2單位月度month水平1單位解釋變地是否實施關閉活禽交易市場的政策clo

8、se。為實施;1為未實施外環(huán)境監(jiān)測陽性率rate實際值月度最低氣溫temp實際值（P）月度降雨筮rainN0的小數(shù)（毫米）月度日照時數(shù)sunN0的小數(shù)實際值（小時）（1）主城區(qū)由主城區(qū)模型可見，不同月份的人感染H7N9禽流感疫情發(fā)生情況存在統(tǒng)計學差異（片。=2.610P=0.0027）,實施關閉活禽交易市場的政策對主城區(qū)人感染H7N9禽流感疫情起到了至關重要的作用彼云=106.704,=0.0001),其次為外環(huán)境監(jiān)測陽性率越高發(fā)病水平越高(*嘉=5.932,=0.0149)ocase,=Boicons+17.448(1.689)close,+3.892(1.598)rate,/30i=0.4

9、23(0.425)+%,.ewN(0,Q)：Q=2.610(0.870)-2xloglikelihood(IGLSDeviance)=68.352(18of48casesinuse)(2) 蕭山區(qū)由蕭山區(qū)模型可見，不同月份的人感染H7N9禽流感疫情發(fā)生情況存在統(tǒng)計學差異(廿。=0.U2P=0.0108),實施關閉活禽交易市場的政策對蕭山區(qū)人感染H7N9禽流感疫情起到了至關重要的作用(兄心=177.785,Pclose=0.0001),其次為外環(huán)境監(jiān)測陽性率越高發(fā)病水平越高以嘉=17.003,%.=0.0001),另外，當月最低氣溫越低發(fā)病水平越高(ZU=6.373=0.0116),當月降雨量越

10、多發(fā)病水平越高(*嘉=5.745,=0.0116)0case.=pQijcons+8.455(0.634)close+4.981(1.208)ratev-0.026(0.010)temp。+0.004(0.001)rain0*Poij=-1-373(0.927)+uOi+eQlj%N(0,Q)：Q=3.187(2.291)e0JN(0,Q)：Qg0.112(0.044)-2xlogHkeHhood(IGLSDeviance)=29.475(17of48casesinuse)(3) 余杭區(qū)由余杭區(qū)模型可見，不同月份的人感染H7N9禽流感疫情發(fā)生情況存在統(tǒng)計學差異(<=0.610,P=0.0

11、001),實施關閉活禽交易市場的政策對余杭區(qū)人感染H7N9禽流感疫情起到了至關重要的作用舊*=25.864,Ry=0.0001),其次為當月最低氣溫越低發(fā)病水平越高(xWp=5.246,Pwmp=0.0220)。casey=Bovcons+3.921(0.771)case。一0.027(0.012)temp,yI3oij=-0.521(0.912)+uOi+eOiJgN(0,Q“)：Q“=3.043(2.181)e。N(0,Q)：Q=.0.610(0.130)-2xloglikelihood(IGLSDeviance)=128.896-(48of48casesinuse)(4) 富陽區(qū).由富陽

12、區(qū)模型可見，不同月份的人感染H7N9禽流感疫情發(fā)生情況存在統(tǒng)計學差異(B。=0.213,P=0.0144),外環(huán)境監(jiān)測陽性率越高發(fā)病水平越高(=11.032=0.0009)0case.=Bacons+2.475(0.745)rate,pQi=0.140(0.157)+eQi%N(0,OQ:Q,=0.213(0.087)-2xloglikelihood(IGLSDeviance)=15.523(12of48casesinuse)(5) 郊縣地區(qū)由郊縣地區(qū)模型可見，不同月份的人感染H7N9禽流感疫情發(fā)生情況存在統(tǒng)計學差異(片。=0.043,=0.0001),當月最低氣溫越低發(fā)病水平越iWUU=5.

13、657,Plerop=0.0174)ocase,=角cons-0.004(0.002)temp,缶=0.076(0.019)+e0<e0/N(0,Q,):Q=0.043(0.004)-2xloglikelihood(IGLSDeviance)=-57.493(192of192casesinuse)討論國內(nèi)多項研究表明，活禽市場的關閉對人感染H7N9禽流感疫情控制有顯著作用，且長時間的休市對H7N9病毒控制的作用優(yōu)于短期休市人群發(fā)病與外環(huán)境H7N9病毒監(jiān)測結果相吻合。因此，本研究將杭州市各區(qū)縣關閉活禽交易的政策實施、外環(huán)境監(jiān)測陽性率納入到影響因素分析中進行建模。研究結果顯示:關閉活禽交易的

14、政策實施、外環(huán)境監(jiān)測陽性率以及最低氣溫、當月降雨量和當月日照時數(shù)都成為影響人感染H7N9禽流感發(fā)生水平的因素。通過人為干預政策的實施,可以在很大程度上改變當?shù)厝烁腥綡7N9禽流感發(fā)生水平。例如從上述蕭山區(qū)兩水平統(tǒng)計模型可以推算出：關閉活禽交易政策的實施使得該地區(qū)每月平均減少8個病例的發(fā)生，雖然2013-2016年，簫山區(qū)共發(fā)生病例21例，占全市病例數(shù)的23.86%,居于各區(qū)縣首位，但永久性關閉當?shù)鼗钋萁灰渍叩膶嵤O大地改變了當?shù)匾咔榈陌l(fā)生水平，政策實施后的2015年和2016年蕭山區(qū)僅分別報告2例和3例人感染H7N9禽流感病例。本研究中多個模型證實了外環(huán)境監(jiān)測陽性率與當?shù)厝烁腥綡7N9禽流感

15、病例發(fā)生水平呈正相關,提示加強外環(huán)境監(jiān)測有利于當?shù)匾咔榈脑缙陬A警。與國內(nèi)同類研究類似”創(chuàng)，本研究還發(fā)現(xiàn)氣象因素與人感染H7N9禽流感病例的發(fā)生水平有關，杭州市分月發(fā)病圖顯示,人感染H7N9禽流感病例發(fā)病集中于冬春季節(jié),說明病毒適宜低溫環(huán)境，因此，本研究在進行氣象因素選擇時，依據(jù)病毒的生物學特性以月度最低氣溫作為可能的影響因素納入建模,分析結果表明最低氣溫越低、降雨量越大，人感染H7N9禽流感病例的發(fā)生水平越高，這符合H7N9病毒在低溫、潮濕環(huán)境中更容易生存的生物學特性。本研究在對不同地區(qū)的影響因素分析中，均獲得因月份不同而疫情發(fā)生水平有所不同的結果，證實了該疫情的發(fā)生存在季節(jié)性高發(fā)的特點。本研

16、究的局限性在于各地區(qū)開展外環(huán)境監(jiān)測不均衡，導致某些地區(qū)有的月份外環(huán)境監(jiān)測數(shù)據(jù)缺失，造成模型擬合時樣本含量的減少，降低了模型擬合的效果。參考文獻1 賈瀟岳，柱沛玲,徐鎮(zhèn)喜，等.人畬流感H7N9病例的流行病學特征分析.汕頭大學醫(yī)學院學報,2016.29(4)：229-230.2 任瑞琦，周蕾，向妮娟，等.中國內(nèi)地人感染H7N9禽流感疫情流行病學特征分析.中華流行病學雜志,2014,35(12):1362-1365.3 楊琨,李曉松，主編.醫(yī)學和公共衛(wèi)生研究常用多水平統(tǒng)計模型.北京大學醫(yī)學出版社.2007.4 許可,鮑倡俊.人感染H7N9禽流感流行病學研究進展.江蘇預防醫(yī)學,2015,26(1):43-47.5 肖琦，張鈞，陳立凌，等.蘇州市人感染H7N9禽流感流行病學特征及外環(huán)境監(jiān)測結果.職業(yè)與健康,2017,33(14):1942-1944.6