杭州市人感染H7N9禽流感流行强度影响因素的两水平模型分析.docx

1、杭州市人感染H7N9禽流感流行强度影响因素的两水平模型分析基金项目：杭州市甄大科技创新专项(20131813A07),杭州市卫生科技重点项目(20I3Z10)1. 杭州市疾病预防控制中心(310021)2. 杭州市气象局通信作者：丁华.E-mail：dh.dhua999黄春萍'丁，华'杨旭辉'刘牧文|王妹'孙昼'考庆君'吴彦'王兵I宋姝娟I叶春2金成2【提要】目的为杭州市制定人感染H7N9禽流感的防控策略提供参考。方法收集杭州市2013-2016年各地区人感染H7N9禽流感病例数、当地政府部门出台的有关活禽交易市场的政策、外环境监测情况以

2、及气象因素，利用ML-wiN2.25软件拟合两水平方差成分模型，探讨不同地区影响人感染H7N9禽流感流行强度的因素c结果关闭活禽交易的政策实施、外环境监测阳性率以及最低气温、当月降雨量和当月日照时数都成为影响人感染H7N9禽流感流行强度的因素。结论关闭活禽交易的政策实施在很大程度上降低了主城区和萧山区疫情的发生水平。外环境监测的阳性率与疫情的发生水平成正比，提示在高发季节加强外环境监测有利于控制疫情发生规模。【关键词】人感染H7N9禽流感流行强度影响因素两水平统计模型人感染H7N9禽流感作为新发传染病，于2013年2月首次出现在我国长三角地区Hi,杭州市2013年3月30日报告了首例人感染H7

3、N9禽流感病例，之后每年均有疫情发生，截止2016年底共报告88例人感染H7N9禽流感病例，死亡28例,病死率为31.82%。为有效防控人感染H7N9禽流感疫情,特开展杭州市人感染H7N9禽流感流行强度的影响因素分析，以期为早期预警和及时采取防控措施提供参考。对象与方法1. 疫情资料杭州市2013-2016年分地区分月人感染H7N9禽流感病例数来源于中国疾病预防控制信息系统。杭州市2013-2016年外环境H7N9样本阳性率来源于杭州市疾病预防控制中心H7N9外环境监测。2. 气象资料杭州市2013-2016年分地区分月最低气温(P)、月度降雨量(毫米)及月度日照时数(小时)均来源于杭州市气象

4、局。3. 研究方法将影响因素归纳为年度影响因素和月度影响因素两个水平,以年度为水平2单位，月度为水平1单位,病例数为应变量,解释变量包括当地是否出台关闭活禽交易市场的政策、外环境监测阳性率、月度最低气温、月度降雨量、月度日照时数等，利用MLwiN2.25软件拟合两水平方差成分模型。变量筛选采用Wald检验，检验水准。定为0.050两水平方差成分模型基本结构为:Zy=%+Pixa+eoa式中i=1,2,nt表示水平1单位;J=1,2，nj表示水平2单位;.为截距;但为解释变量的回归系数;e回为随机误差项，即水平1单位的随机误差。两水平方差成分模型建立在两个假设的基础上：(1) E(ev)=0,v

5、ar(ev)=兄(2) cov(u0,e0v)=0因此,两水平方差成分模型也可以表示为坊=(角+&与)+(uOJ+eOij)其中场和&描述模型的固定效应(fixedeffects)；和廿。描述模型的随机效应(randomeffects)，分别是水平2单位的方差成分和水平1单位的方差成分。结果1. 人感染H7N9禽流感频数分布杭州市自2013年报告首例H7N9禽流感病例以来，每年均有病例报告，发病可分为两个阶段:第一阶段为2013-2014年，该阶段发病水平较高，发病地区以主城区和郊区为主;第二阶段为2015年之后，该阶段发病水平明显下降，发病地区转变为以郊区和郊县为主(表1)。

6、从分月发病图(图1)上可见，人感染H7N9禽流感病例发病集中于冬春季节。表1杭州市各地区2013-2016年H7N9禽流感发病情况地区名称地区分类发病数2013年2014年2015年2016年上城区主城区8300下城区主城区4000江干区主城区5220拱墅区主城区1111西湖区主城区3011滨江区主城区0300萧山区郊区51123余杭区郊区2626富阳区郊区，0222淳安县郊县0001建德市郊县0100桐庐县郊县0111临安市郊县1012合计293012172. 人感染H7N9禽流感影响因素分析r以年度为水平2单位，月度为水平1单位,病例数为应变量（case）,解释变量包括当地是否实施关闭活禽

7、交易市场的政策（close）、外环境监测阳性率（rate）、月度最低气温（temp）、月度降雨量（rain）、月度日照时数（sun）等（表2）。由于杭州市主城区关闭活禽交易市场的时间统一，故将上城区、下城区、西湖区、拱墅区、江干区、滨江区等6个区合并为主城区进行拟合,桐庐县、淳安县、建德市、临安市等4个郊县因病例数偏少而合并为郊县地区进行拟合，因此获得杭州市主城区、萧山区、余杭区、富阳区、郊县地区共5类地区的人感染H7N9禽流感影响因素分析。表2变量赋值说明指标名称变量名定义及赋值应变量caseNO的整数年度year水平2单位月度month水平1单位解释变地是否实施关闭活禽交易市场的政策clo

8、se。为实施;1为未实施外环境监测阳性率rate实际值月度最低气温temp实际值（P）月度降雨筮rainN0的小数（毫米）月度日照时数sunN0的小数实际值（小时）（1）主城区由主城区模型可见，不同月份的人感染H7N9禽流感疫情发生情况存在统计学差异（片。=2.610P=0.0027）,实施关闭活禽交易市场的政策对主城区人感染H7N9禽流感疫情起到了至关重要的作用彼云=106.704,=0.0001),其次为外环境监测阳性率越高发病水平越高(*嘉=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) 萧山区由萧山区模型可见，不同月份的人感染H7N9禽流感疫情发生情况存在统计学差异(廿。=0.U2P=0.0108),实施关闭活禽交易市场的政策对萧山区人感染H7N9禽流感疫情起到了至关重要的作用(兄心=177.785,Pclose=0.0001),其次为外环境监测阳性率越高发病水平越高以嘉=17.003,%.=0.0001),另外，当月最低气温越低发病水平越高(ZU=6.373=0.0116),当月降雨量越

10、多发病水平越高(*嘉=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) 余杭区由余杭区模型可见，不同月份的人感染H7N9禽流感疫情发生情况存在统计学差异(<=0.610,P=0.0

11、001),实施关闭活禽交易市场的政策对余杭区人感染H7N9禽流感疫情起到了至关重要的作用旧*=25.864,Ry=0.0001),其次为当月最低气温越低发病水平越高(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) 富阳区.由富阳

12、区模型可见，不同月份的人感染H7N9禽流感疫情发生情况存在统计学差异(B。=0.213,P=0.0144),外环境监测阳性率越高发病水平越高(=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) 郊县地区由郊县地区模型可见，不同月份的人感染H7N9禽流感疫情发生情况存在统计学差异(片。=0.043,=0.0001),当月最低气温越低发病水平越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)讨论国内多项研究表明，活禽市场的关闭对人感染H7N9禽流感疫情控制有显著作用，且长时间的休市对H7N9病毒控制的作用优于短期休市人群发病与外环境H7N9病毒监测结果相吻合。因此，本研究将杭州市各区县关闭活禽交易的政策实施、外环境监测阳性率纳入到影响因素分析中进行建模。研究结果显示:关闭活禽交易的

14、政策实施、外环境监测阳性率以及最低气温、当月降雨量和当月日照时数都成为影响人感染H7N9禽流感发生水平的因素。通过人为干预政策的实施,可以在很大程度上改变当地人感染H7N9禽流感发生水平。例如从上述萧山区两水平统计模型可以推算出：关闭活禽交易政策的实施使得该地区每月平均减少8个病例的发生，虽然2013-2016年，箫山区共发生病例21例，占全市病例数的23.86%,居于各区县首位，但永久性关闭当地活禽交易政策的实施极大地改变了当地疫情的发生水平，政策实施后的2015年和2016年萧山区仅分别报告2例和3例人感染H7N9禽流感病例。本研究中多个模型证实了外环境监测阳性率与当地人感染H7N9禽流感

15、病例发生水平呈正相关,提示加强外环境监测有利于当地疫情的早期预警。与国内同类研究类似”创，本研究还发现气象因素与人感染H7N9禽流感病例的发生水平有关，杭州市分月发病图显示,人感染H7N9禽流感病例发病集中于冬春季节,说明病毒适宜低温环境，因此，本研究在进行气象因素选择时，依据病毒的生物学特性以月度最低气温作为可能的影响因素纳入建模,分析结果表明最低气温越低、降雨量越大，人感染H7N9禽流感病例的发生水平越高，这符合H7N9病毒在低温、潮湿环境中更容易生存的生物学特性。本研究在对不同地区的影响因素分析中，均获得因月份不同而疫情发生水平有所不同的结果，证实了该疫情的发生存在季节性高发的特点。本研

16、究的局限性在于各地区开展外环境监测不均衡，导致某些地区有的月份外环境监测数据缺失，造成模型拟合时样本含量的减少，降低了模型拟合的效果。参考文献1 贾潇岳，柱沛玲,徐镇喜，等.人畲流感H7N9病例的流行病学特征分析.汕头大学医学院学报,2016.29(4)：229-230.2 任瑞琦，周蕾，向妮娟，等.中国内地人感染H7N9禽流感疫情流行病学特征分析.中华流行病学杂志,2014,35(12):1362-1365.3 杨琨,李晓松，主编.医学和公共卫生研究常用多水平统计模型.北京大学医学出版社.2007.4 许可,鲍倡俊.人感染H7N9禽流感流行病学研究进展.江苏预防医学,2015,26(1):43-47.5 肖琦，张钧，陈立凌，等.苏州市人感染H7N9禽流感流行病学特征及外环境监测结果.职业与健康,2017,33(14):1942-1944.6