吉林省2008~2016年手足口病流行特征分析及时间序列模型的应用

《吉林省2008~2016年手足口病流行特征分析及时间序列模型的应用》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

分类号：R183.4单位代码：10183研究生学号：201472D058密级：公开吉林大学硕士学位论文（专业学位）吉林省2008～2016年手足口病流行特征分析及时间序列模型的应用Epidemiccharacteristicsofhand-foot-mouthdiseaseandtheapplicationoftimeseriesmodelinJilinprovincefrom2008to2016作者姓名：栾博类别：公共卫生硕士领域（方向）：疾病预防与控制指导教师：李波教授培养单位：公共卫生学院2017年12月 吉林省2008～2016年手足口病流行特征分析及时间序列模型的应用Epidemiccharacteristicsofhand-foot-mouthdiseaseandtheapplicationoftimeseriesmodelinJilinprovincefrom2008to2016作者姓名：栾博专业名称：公共卫生硕士指导教师：李波教授学位类别：公共卫生硕士答辩日期：2017年12月2日 未经本论文作者的书面授权，依法收存和保管本论文书面版本＼电子版本的任何单位和个人，均不得对本论文的全部或部分内容进行任何形式的复制、修改、发行＼出租、改编等有碍作者著作权的商业性使用（但纯学术性使用不在此限）。否则应承担侵权的法律责任。，吉林大学硕士学位论文原创性声明本人郑重声明：所呈交学位论文，是本人在指导教师的指导下，独立进行研究工作所取得的成果２除文中己经注明引用的内容外，本论文不包含任何其他个人或集体已经发表或撰写过的作品成果。对本文的研究做出重要贡献的个人和集体，均己在文中以明确方式标明。本人完全意识到本声明的法律结果由本人承担。学位论文作者签名：曰期：々年ｕ月ｒ日ｑ 中文摘要吉林省2008～2016年手足口病流行特征分析及时间序列模型的应用目的：分析吉林省手足口病流行特征及病原学特点，构建手足口病时间序列模型并评价预测效果，从而探究手足口病的发病与气象因素之间的关系，为吉林省今后制定更加科学、有效的手足口病预防控制措施提供合理化的建议。方法：使用Excel2010软件建立吉林省手足口病数据库，采用描述性流行病学方法，对2008～2016年吉林省手足口病疫情资料进行分析。使用IBM_SPSS21.0软件对2008年～2017年5月年吉林省手足口病月发病数构建时间序列预测模型，利用2017年6～8月月发病数进行预测和评价。使用IBM_SPSS21.0软件进行Spearman相关分析探讨长春市2014～2016年手足口病发病数与平均气压（hpa）、平均温度（℃）、平均相对湿度（％）、降水量（mm）、平均风速（m/s）的相关关系、各个气象因素之间的相关关系及相关程度。结果：1.2008～2016年吉林省HFMD累计发病146538例，年平均发病率59.34/102万，发病率在不同年份间具有统计学差异(X=30469.57,P＜0.001)。2010年、2011年和2014年重症发生率(重症数/病例数)高于其他年份，不同年份之间的重症率2差异有统计学意义（X=477.21,P＜0.001）。除2013年、2015年和2016年外各年均有死亡病例报告。2.2008～2016年吉林省HFMD患者主要集中在5岁以下儿童，占79.95%，男女性别比为1.46:1。病例以散居儿童最多，占发病例数的58.68%；幼托儿童位居第二位，占发病例数的34.66%；学生位居第三位，占6.09%。3.2008～2016年呈现隔1年流行强度增强的特点。全年各月均有病例报告，I 各年发布趋势基本相同。4月开始发病数上升，高发季节为6～8月，10月后疫情下降，除2016年发病高峰出现在8月外，其他年度发病高峰均出现在7月。4.全省9个地区均有病例报告，年平均发病率前三位的地区分别为延边朝鲜族自治州（92.53/10万，2020例）、通化市（89.57/10万，2246例）以及四平市（70.19/10万，2500例）。5.2008～2016年吉林省实验室诊断手足口病病例总体以EV71阳性为主占35.04%，但各年主要流行型别不尽相同，2008～2010年及2014年以EV71阳性为主，2012年以CoxA16阳性为主，2013、2015和2016年以其他EV阳性为主。6.2008～2016年吉林省共报告HFMD重症1186例，其中男女性别比1.73:1；散居儿童占79.43%，托幼儿童占18.80%，学生占1.77%；1～3岁年龄组报告的重症数较多。病原学检测显示EV71阳性占52.02%，CoxA16阳性占6.16%，其他EV阳性占8.60%，阴性占33.22%。2008～2016年吉林省共报告HFMD死亡病例35例，其中有34例为重症病例；35例死亡病例中，19例为男性，16例女性；死亡年龄均为5岁以下，其中2～岁年龄组最高；病原学检测显示以EV71阳性为主。7.用吉林省2008～2017年5月手足口病疫情数据建模，经过逐个筛选和参数检验拟合出ARIMA（0,1,0）（0,1,1）12模型，并在此基础上预测2017年6月~8月吉林省手足口病的实际发病情况。通过对预测值与实际值进行比较，结果表明实际值包含于预测值的95%可信区间，平均相对误差23.78%，该模型具有较强的预测能力。8.长春市HFMD周发病数与周平均气温呈高度正相关，与周平均气压及平均风速呈中度及弱的负相关，与相对湿度及降水量呈弱及中度的正相关。长春市各气象因素之间也具有一定的相关性，其中平均气压与平均温度、平均降水量之间分别呈高度及中度负相关，平均温度与平均降水量之间呈中度正相关，平均相对湿度与平均降水量之间呈中度正相关，平均相对湿度与平均风速之间呈中度负相关。II 结论：1.吉林省手足口病在2008年～2016年期间具有明显的季节性，6~8月为本地手足口病的高发季节；高发地区为延边朝鲜族自治州、通化市和四平市；人群分类上以散居儿童为主，年龄5岁以下儿童为主，男性高于女性；实验室诊断病例总体以EV71阳性为主；重症及死亡病例分布特征与整体一致。2.用吉林省2008～2017年5月手足口病疫情数据拟合出ARIMA（0,1,0）（0,1,1）12模型，对2017年6～8月病例预测效果较好。3.长春市2014～2016年手足口病周发病数与平均气压（hpa）、平均温度（℃）、平均相对湿度（％）、降水量（mm）、平均风速（m/s）存在相关性。关键词：手足口病，相关性，时间序列，预测，ARIMAIII ABSTRACTEpidemiccharacteristicsofhand-foot-mouthdiseaseandtheapplicationoftimeseriesmodelinJilinprovincefrom2008to2016Objective：Toanalyzetheepidemiologicalandpathogeniccharacteristicsofhand-foot-mouthdisease（HFMD）inJilinProvince,toestablishatime-seriesanalysismodelofHFMDandassessitspredictiveeffect,andtoexaminetheassociationbetweenHFMDmorbidityandmeteorologicalfactors,soastoproviderationalizationproposalforthescientific,effectivepreventionandcontrolofHFMDinJilinProvinceinthefuture.Methods：Excel2010wasusedtoestablishthedatabaseofhand-foot-mouthdiseaseinJilinprovince,andtheepidemicsituationdataofHFMDfrom2008to2016inJilinprovincewasanalyzedwiththemethodofdescriptiveepidemiology.IBM_SPSS21.0wasusedtoconstructthetimeseriespredictionmodelforHFMDmonthlycasesfrom2008toMay,2017,andpredictionandevaluationwasmadewiththeHFMDcasesfromJunetoAugustin2017.IBM_SPSS21.0wasusedtoconductSpearmancorrelationanalysis,todiscusstherelationshipbetweenHFMDcasesfrom2014to2016inChangchunandtheaveragepressure(hpa),theaveragetemperature(℃),theaveragerelativehumidity(%),precipitation(mm),theaveragewindspeed(m/s),aswellastherelationshipbetweenvariousmeteorologicalfactorsandthedegreeofthecorrelation.Results：1.Atotalof146538casesofhand-foot-mouthdiseasewerereportedinJilinProvinceduring2008~2016,withanannualincidencerateof59.34/100000,therewas2asignificantdifferenceindifferentyears(X=30469.57,P＜0.001).TheincidenceofIV severeillness(severenumber/case)washigherin2010,2011and2014thanotheryears,thedifferenceofsevereillnessincidencebetweendifferentyearswas2statisticallysignificant（X=477.21,P＜0.001）,deathsarereportedonhardlyeveryyearexpect2013,2015and2016.2.MostHFMDcaseswerechildrennuder5yearsold,accountingfor79.05%ofthetotal,thesexratioofmale:femaleis1.46:1.Mostcaseswerescatteredchildren,accountingfor58.68%ofthenumberoftotalcases,theproportionofkindergarten’schildren,whichwasthesecondhighest,was34.66%,theproportionofstudentwas6.09%.3.From2008to2016,thestrengthofepidemicstrengthincreasedeveryotheryear.Therearecasereportsineachmonthoftheyears,andthetrendsarebasicallythesameineachyear.ThenumberofcasesbegantoriseinApril,thehighincidenceseasonwas6~8months,andtheepidemicdeclinedafterOctober,thepeakincidenceusuallyappearedinJuly,exceptfor2016,inwhichyearthepeakincidencewasinAugust.4.TherearecasereportsinallthenineprovincesinJilinProvince.ThetopthreeregionsoftheannualaverageincidenceratewereYanbianKoreanAutonomousPrefecture(92.53/10million,2020cases),Tonghua(89.57/10million,2246cases)andSipingCity(70.19/10million,2500cases).5.Thehand-foot-mouthdiseasecasesdiagnosedbylaboratoryinProvincewasmainlypositiveEV71,accountingfor35.04%,butthemajorepidemictypesarenotthesamethroughouttheyears,in2008~2010yearsand2014thereweremainlypositiveEV71,positiveCoxA16in2012,andotherpositiveEVin2013,2015and2016.6.Atotalof1186severeHFMDcaseswerereportedinJilinProvincefrom2008to2016,theproportionofmaleandfemalewas1.73:1.Scatteredchildrenaccountedfor79.43%,kindergarten’schildrenaccountedfor18.80%,studentsaccountedfor1.77%.The1~3-year-oldgroupreportedahighernumberofseverecases.PathogenicdetectionshowedthatpositiveEV71accountedfor52.02%,positiveCoxA16V accountedfor6.16%,andotherpositiveEVaccountedfor8.60%,negativevirusaccountedfor33.22%.In2008~2016years,35casesofHFMDdeathswerereportedinJilinProvince,andtherewere34severecasesamongthem.Therewere19malesand16femalesinthedeathcases,alllessthan5yearsold,andthenumberofcasesof2~3yearsoldwasthehighest,TheetiologytestshowedthatpositiveEV71wasthemajorpathogen.7.ModelwasestablishedbyHFMDdatainJilinProvincefrom2008toMay,2017,thenARIMA(1,0,1)(1,1,0)12wasfittedthoughscreeningandtexting,andtopredicttheactualdataseriesofHFMDinJilinProvincefromJunetoAugustin2017.Comparedthepredictedvalueandtheactualvalue,theresultsshowthattheactualvaluesarewithinthe95%confidenceintervalofthepredictedvalue,andtheaveragerelativeerroris23.78%,thepredictiveabilityofourmodelwasgood.8.TheweeklycasesofHFMDinChangchunwaspositivelycorrelatedwiththeaveragetemperature,showedmoderateandweaknegativecorrelationwithmeanbarometricpressureandmeanwindspeed,weaklyandmoderatelycorrelationwithrelativehumidityandprecipitation.TherewerealsosomecorrelationsamongthemeteorologicalfactorsinChangchun,meanpressureshowedhighandmoderatenegativecorrelationswithaveragetemperatureandaverageprecipitation,therewasamoderatepositivecorrelationbetweenaveragetemperatureandaverageprecipitation,alsoamoderatepositivecorrelationbetweentheaveragerelativehumidityandtheaverageprecipitation,andamoderatenegativecorrelationbetweenmeanrelativehumidityandmeanwindspeed.Conclusion：1.InJilinProvince,hand-foot-mouthdiseaseshowedanobviousseasonaldistributionfrom2008to2016,thehigh-occurrenceseasonwasintheperiodofJunetoAugust;thehighincidenceareaswereYanbianKoreanAutonomousPrefecture,TonghuaCityandSipingCity;scatteredchildrenwerethemajorpatients,mostcaseswerechildrenunder5yearsold,malesweremorethanfemales;laboratorydiagnosticshowedthatthemajorityofcaseswerepositiveEV71;thedistributioncharacteristicsVI ofsevereanddeathcaseswereinlinewiththeoverallsituation.2.TheARIMA(0,1,0)(0,1,1)12Model,fittedthoughtheepidemicdataofhand-foot-mouthdiseaseinJilinProvincefrom2008toMay2017,wassatisfactoryforthepredictionofthecasesfromJunetoAugustin2017.3.TherewasacorrelationbetweenweeklyHFMDcasesfrom2014to2016inChangchunandmeteorologicalfactors,suchasaveragepressure(hpa),theaveragetemperature(℃),theaveragerelativehumidity(%),precipitation(mm),theaveragewindspeed(m/s).Keywords:Hand-foot-mouthdisease,correlation,timeseries,prediction,ARIMAVII 目录第1章前言........................................................................................11.1手足口病.........................................................................................11.1.1手足口病概述..........................................................................11.1.2病原学.......................................................................................11.1.3临床表现..................................................................................21.1.4治疗原则..................................................................................21.1.5疫苗研发..................................................................................21.1.6流行特征..................................................................................21.1.7国内外流行现状......................................................................31.1.8手足口病研究进展..................................................................41.2时间序列分析.................................................................................51.2.1时间序列概述..........................................................................51.2.2时间序列分析基本原理..........................................................5第2章材料与方法................................................................................72.1研究内容.........................................................................................72.2研究对象.........................................................................................72.2.1资料来源..................................................................................72.2.2病例定义..................................................................................72.2.3标本采集及实验室检测..........................................................82.2.4病例的发现与报告..................................................................82.3HFMD时间序列模型构建及预测.................................................82.3.1基本理论..................................................................................82.3.2ARIMA(p,d,q)(P,D,Q)s模型构建基本步骤..........................82.4HFMD发病与气象因素的相关分析.............................................92.5统计学方法及应用.......................................................................10VIII 第3章结果......................................................................................113.1HFMD三间分布...........................................................................113.1.1HFMD人群分布....................................................................113.1.2HFMD时间分布....................................................................133.1.3HFMD地区分布....................................................................143.1.4病原学特征............................................................................153.1.5重症及死亡分析....................................................................163.2HFMDARIMA(P,D,Q)(P,D,Q)S模型构建及预测.........................173.2.1序列平稳化............................................................................173.2.2模型识别与检验....................................................................183.2.3模型预测................................................................................193.3HFMD发病与气象因素的关系探讨...........................................203.3.1气象资料描述........................................................................203.2.2相关性分析............................................................................21第4章讨论......................................................................................224.1手足口病流行特征分析...............................................................224.1.1时间分布特征分析................................................................224.1.2人群分布特征分析................................................................224.1.3病原学分布特征分析............................................................234.2ARIMA模型的建立.....................................................................244.3手足口病发病与气象因素的关系分析.......................................244.4手足口病防控措施建议...............................................................25第5章结论......................................................................................27参考文献...................................................................................................28作者简介及在学期间所取得的科研成果..............................................36致谢.......................................................................................................37IX 英文缩略词英文缩写英文全称中文全称HFMDHand-Foot-MouthDisease手足口病CoxA16CoxsackievirusA16柯萨奇病毒A组16型EV71Humanenterovirus71肠道病毒71型ARIMAAutoregressiveIntegratedMovingAverageModel差分自回归移动平均模型ACFAutocorrelationFunction自相关函数PACFPartialAutocorrelationFunction偏自相关函数ARMAAuto-regressiveandMovingAverageModel自回归移动平均模型X 第1章前言1.1手足口病1.1.1手足口病概述手足口病（Hand-Foot-MouthDisease，HFMD）是一种急性传染性疾病，由肠道病毒引起，多为coxsackievirusA16(CoxA16)和humanenterovirus71(EV71)[1]感染，感染者的年龄通常为5岁以下。HFMD具有广泛的传播途径，如咽喉分泌物（如唾液、痰和鼻涕）、病人疱疹液等。病人、隐性感染者以及无明显症状[2]的病毒携带者均可成为该病的传染源。由于HFMD具有复杂的传播途径以及[3]较强的传染性，所以如果预防控制不当，该病容易在短时间内产生大流行。为了降低手足口病的传播带来的社会负担，卫生部早在2008年5月便将HFMD列为我国丙类传染病进行严格管理，并采取了有效的控制措施，随时监测国内的HFMD疫情。1.1.2病原学手足口病病毒具有多种型别，其中人类A组肠道病毒是手足口病最重要的[4]病原体，特别是柯萨奇病毒A16和EV71。这两种病毒具有较高的遗传学同源性，感染人体后的症状也基本相似。但不同的是，CoxA16感染人体后一般症状较轻微，且能够在短期内自愈；而EV71感染起病急，病情进展迅速，累计呼吸[5]及中枢神经系统，严重者可导致患者死亡。近些年来在世界范围内，EV71感染在HFMD患者中所占比例呈升高趋势，尤其是重症及死亡病例，已成为该病[6-9]主要的病原体。肠道病毒由于其特殊的理化性质，使其能够在湿热环境下大量繁殖、传播。在4℃环境中该病毒可存活1年，在-20℃环境中长期存活。该病毒对5%来苏和75%酒精具有抵抗性，不能被去氯胆酸盐、乙醚等灭活，但可以在特殊的条件下将该病毒灭活，如干燥环境、56℃30分钟、紫外线照射、碘酒、甲醛、漂白粉[10]及高锰酸钾等氧化剂。1 1.1.3临床表现[10]HFMD的潜伏期2~10天，一般为3~5天，病程一般为7~10天。起病较急，伴发热，口腔黏膜、手足及臀部可见疱疹，疱疹内部有少量液体，周围伴炎性红晕。疱疹感染人体后，可引发咳嗽、流涕等呼吸道症状，部分患者有食欲不振等消化道症状。多数患者症状较轻，且愈后良好，但对于少数EV71感染的患者则症状较严重，可出现脑炎、脑膜炎、脑脊髓炎等症状，严重者可产生后遗症，[10]甚至出现死亡病例。1.1.4治疗原则支持疗法为主，多数患者可自愈，少部分患儿发展为重症，如不进行及时的对症治疗，病情得不到有效控制，可导致患儿死亡。因此，针对HFMD病毒的传播途径和特点，随时监测HFMD的发病情况，采取有效的预防控制措施，做到早发现、早诊断、早治疗，可以有效降低HFMD的感染，并降低人群死亡率。1.1.5疫苗研发我国对于HFMD的疫苗研究进展较为迅速，日前我国针对EV71的疫苗已经上市，对于有效降低HFMD病毒感染具有重要价值。研究发现，我国的EV71疫苗有效性较好，对于轻症病例的有效率为94.8%，重症病例已达100%。另外，[11]临床试验表明该疫苗的安全性也较好，尚无明显的不良反应。1.1.6流行特征传染源：HFMD的主要传染源通常为该病的患者和隐性感染者。其中隐性感染者通常为成人，因为成人较儿童来讲，具有较强的抵抗力，感染肠道病毒后不具有典型症状而被忽视，成为隐形传染源而使抵抗力较弱的儿童发生感染。而5岁以下儿童由于抵抗力低，感染后常表现出典型的临床症状。但根据新的研究表明，肠道EV71病毒可以广泛的存在于人群中无症状传播，研究结果显示在流行季节，5岁以下健康儿童的病毒携带率为27.33%，健康成人的病毒携带率为[12-15]36%，这一较高的比例提示手足口病具有极高的隐秘性传染性特征。传播途径：HFMD具有广泛的传播途径，粪-口途径为最主要的传播方式，此外也可经呼吸道和接触传播。接触患者的疱疹液、口鼻分泌物以及其他被污染2 的物品均有造成感染的风险。患者感染HFMD后，大量的病毒存在于疱疹液中，接触疱疹液可导致病毒传播；患者的唾液及咽喉分泌物中也含有手足口病病毒，咳嗽或打喷嚏可以将病毒播散到空气中，从而在人群间进行传播；患者的粪便中也含有大量病毒，接触后感染的几率也将加大，即使病后数周，粪便中仍可检到病毒。流行期间，人群密集场所尤其是抵抗力较低的低龄儿童集中的学校、幼儿园等托幼机构容易引起手足口病的大面积传播；手足口病病毒通过污染水和食物也能够造成人群间的感染；儿童之间由于管理不严，交叉接触食具、奶具、毛巾等用品以及玩具、床上用品等也可以造成感染。此外，医疗机构交叉感染也会引起肠道病毒的感染传播。1.1.7国内外流行现状手足口病呈全球流行趋势，最早的报道可追溯到1957年，新西兰人对其进[16]行了报道并描述了手足口病的症状。1958年，一项加拿大的研究得到了[17]CoxA16病毒，并初步确定了CoxA16就是HFMD的病原体。1959年，英国学者成功分离出CoxA16，并证实其为导致手足口的主要病原体，同时，证实了该病主要症状为口腔炎，手足部斑疹、丘疹或疱疹，“手足口病”便因此而得名，[18]并沿用了下来。1969年，被分离出的EV71改变了人们认为CoxA16是引起手[19]足口病的主要病原体这一观点，1972年美国专家首度证实EV71同为手足口[20]病的致病病毒。至此，医学界确认手足口病是以CoxA16和EV71两种病原体导致。1975年和1978年，在欧洲东部的保加利亚和匈牙利分别发生了手足口病的暴发，这两次暴发的病原体均为EV71，患者的症状主要累及中枢神经系统，最[21-22]终分别造成了保加利亚和匈牙利各40多例死亡；1994年的英国，由CoxA16[23]感染而引起了HFMD的暴发，但多数患者没有发生严重的症状。其它国家如加拿大、德国、俄罗斯、巴西等国家，也先后不同程度地出现手足口病的散发或暴发。近年来，手足口病在亚太地区也广泛流行，中国、日本、马来西亚、印度、菲律宾等国家先后报道了手足口病的流行情况。日本是资料记载中大规模暴发疫情最多的国家之一，1969～1970年暴发的疫情以CoxA16感染为主，1973和1978年暴发的疫情以EV71感染为主。1982年日本再次暴发手足口病疫情，报3 告病例数达到14万多；1997年日本的手足口病发病例数再度明显上升，病原体[24-26]以EV71、CoxA16为主。20世纪末以来，手足口病在亚太地区开始肆虐，1997年马来西亚、1999年澳大利亚、2005年越南均报道了HFMD的暴发，病[27-31]原体均为EV71。在我国，手足口病也长期处于高发态势。1981年上海首次发现手足口病后，全国范围内多个省市均报道了手足口病，范围包括华北、华中、东北和沿海地区。[32]1983年和1986年天津暴发手足口病，经检验主要病原体为CoxA16。1998年[33-34]台湾暴发手足口病，经检验主要病原体为EV71。2000年山东、2001年北京、2002年上海，也都有不同程度的HFMD暴发的报道。根据我国国内的报告，2006年除西藏以外，我国各省、自治区、直辖市均发现了HFMD病例，全国共计13637例，死亡病例6例。截止目前，手足口病已波及我国全部省自治区直辖市。为了降低手足口病的传播带来的社会负担，卫生部在2008年5月将HFMD列为我国丙类传染病进行严格管理，并采取了有效的控制措施，随时监测国内的HFMD疫情。1984年吉林省首次报道手足口病，2008年纳入法定传染病网络直报以来始终呈高发态势，各年发病率均居全省甲乙丙类传染病的前五位，2008～2016年年均发病率为59.34/10万。1.1.8手足口病研究进展HFMD患者及隐性感染者是HFMD的主要传染源。HFMD病毒在人群间传播的方式多种多样，人群间的密切接触是传播的重要方式，此外，通过空气飞沫、水和食物也可以引起HFMD病毒的感染。此前，HFMD病毒的疫苗研究已经获得了成功。同时，手足口病的防控措施以避免接触病人、全面消毒受污染的物品及场所、改善个人卫生习惯等为主。3岁以下婴幼儿是HFMD最主要的易感人[4][35]群。有研究表明手足口病的人群有种族差异。国内手足口病的季节特点略有不同，南方呈现5～6月为第一个发病高峰，9～10月是第二个发病高峰，而[36]北方多为仅6～9月份一个发病高峰，冬季发病较少。HFMD易发生在人口聚集、卫生条件差、经济欠发达的地区以及高温、低压、降水量大的地区，上述因素适宜病毒生存、繁殖及传播。研究表明，发病地区的纬度会对HFMD的流行4 情况及流行特征产生一定的影响。在纬度较低的热带及温带地区，HFMD发病率更高，每年的发病高峰也更持续。另外，许多气象因素也与HFMD的发病存[37]在一定的关联。目前，越来越多的研究应用时间序列模型探索手足口病的发病趋势，但尚未有对吉林省手足口病时间序列的研究，本研究将应用2008～2017年5月吉林省手足口病疫情数据构建手足口病时间序列模型并评价预测效果。1.2时间序列分析1.2.1时间序列概述任何事物都处在不断的发展变化之中，人类总是希望能够透过现象看到事物的本质和它发展变化的规律，这就需要探索一定的方法去解决问题。时间序列分析的思想是按照时间顺序将数据排列起来，从而研究、分析这些数据随时间变化的规律，然后建立起相应的模型，并用这个模型对未来进行预测。公元前5000年左右的古埃及，为了降低尼罗河泛滥给农业带来的损失，古埃及人将尼罗河每天的水位情况记录了下来，从而掌握了尼罗河河水的涨落规律，保障了农业发展，[38]这就是可以追寻到的时间序列分析方法的最早的雏形。进入20世纪，时间序[39]列分析得到实质性的发展，这得益于英国学者尤尔创建的平稳性自回归模型。[39]上世纪30年代，Yule建立了AR（2）和AR（4）模型，将时间序列变量之间[40]的关系清晰地呈现了出来，促进了时间序列分析的进步。此后，Walker和[41][42]Slutsky又相继建立了AR(S)和AR(h)模型。40年代末期，Wold证实了离散平稳过程由线性回归和隐周期组成，这一思想为使用模型拟合平稳序列奠定了理[43][43]论基础。1970年GeorgeEdwardPelhamBox和GwilymMeirionJenkins讨论了非平稳自回归移动平均模型(ARIMA)及建模使用方法，此理论将原有时间序列分析方法提升了一个新的高度，目前应用广泛且发展迅速，近年来在公共卫生领域取得了良好的预测效果。1.2.2时间序列分析基本原理时间序列分析以数理统计学和随机过程理论为基础，对一系列观测随机数据所遵从的统计规律进行研究并对未来值进行预测。5 时间序列分析法是一种定量预测方法，它的基本原理是：一是认为事物的发展是延续的。通过对观测数据的分析就能预测出事物未来的趋势。二是考虑事物发展的随机性。事物发展规律不全是必然的，也受偶然因素的影响，这就需要对数据进行处理。时间序列分析通常包括建立统计模型，并进行统计推断，完成自相关分析等一般的统计学分析，另外还包括预测、控制时间序列以及滤波等工作[44]。时间序列一般可预测三种变化规律：随机性、趋势性、周期性。时间序列建立模型的基本步骤是：①在观测系统中通过观测、调查、统计、抽样等方法得到一组具有一定时序的动态数据。②根据得到的动态数据得出相关图，并通过对相关图进行相关分析，得到自相关模型中的函数。③选择合适的随机模型，并进行[45]曲线拟合，即通过随机模型拟合实际的观测数据。重复上述三个步骤，直到找到最优模型进行预测。平稳的时间序列可以用ARMA模型（自回归滑动平均模型）来进行拟合，非平稳的时间序列，可先通过数据转换转化为平稳时间序列，再使用平稳时间序列的方法进行拟合。将季节因素添加到传统的ARIMA模型后，可得到SARIMA模型，后者对[56][57]于传染病的研究有着更显著的预测效果。刘涛和颜玉炳等分别通过[58]SARIMA模型，成功地预测了山东省和厦门市的HFMD发病情况。李丽丽引入季节因素应用多种模型对2014年北京市房山区手足口病发病率进行预测，结[4]果显示SARIMA(0,0,1)(0,1,0)12为最适合模型。闫莉分别用引入和不引入季节因素两种方式对吉林省长春市手足口病模型进行构建并预测，结果显示复SARIMA(2,1,2)(1,1,1)12模型预测值与实际值的平均误差小于单一SARIMA(1,1,1)(1,1,0)12模型预测平均误差。6 第2章材料与方法2.1研究内容分析吉林省手足口病流行特征及病原学特点，构建手足口病时间序列模型并评价预测效果，探究手足口病的发病与气象因素之间的关系，为吉林省手足口病防控措施的制定提供科学依据。2.2研究对象2.2.1资料来源2008～2016年吉林省HFMD疫情资料数据来源于“中国疾病预防控制信息系统”，选用现住址为吉林省HFMD临床诊断病例及实验室确诊病例为研究对象。人口资料来源于“中国疾病预防控制信息系统”及《吉林省统计年鉴》统计数据。2014～2016年吉林省长春市气象因素资料来源于吉林省气象局，主要包括平均气压（hpa）、平均温度（℃）、平均相对湿度（％）、降水量（mm）、平均风速（m/s）5个因素。2.2.2病例定义2.2.2.1临床诊断病例发病时间位于流行季节内，多数为学龄前儿童。(1)普通病例：手、足、口及臀部出现皮疹，伴或不伴发热(2)重症病例：累及神经系统，或引发循环功能障碍。经进一步实验室检查，可有血糖增高、外周血白细胞增高等表现，超声心动图及胸部X线检查可发现[10]异常。少数病例症状不典型，需结合病原学检查作出诊断。2.2.2.2实验室确诊病例符合下列条件之一即可诊断：(1)标本中检测到EV71或CVA16等人肠道病7 毒特异性核酸；(2)组织标本中分离出可导致HFMD的肠道病毒；(3)在急性期与恢复期内，血清HFMD肠道病毒特异性抗体升高4倍及以上，或血清肠道病毒[10]型抗体滴度≥1∶256。2.2.3标本采集及实验室检测医疗单位对HFMD的普通病例、重症病例及死亡病例采集疱疹液、咽拭子、粪便等标本，并将标本送至实验室，应用分子生物学技术对HFMD标本进行核[10]酸检测。2.2.4病例的发现与报告根据《中华人民共和国传染病防治法》及《传染病信息报告管理规范》的相关规定，各医疗单位对HFMD的病例24小时内通过“中国疾病预防控制信息系[10]统”进行网络报告。2.3HFMD时间序列模型构建及预测2.3.1基本理论时间序列就是把随机事件的整个变化的过程，按照时间的先后顺序进行记[49]录，记录的内容就构成了一个时间序列。对于过往的时间序列进行观察、分析，探索随机事件发展的规律，并根据这一规律对未来该随机事件的发展情况进行预[38]测，就是时间序列分析的研究内容。上世纪60年代，西方学者Box和Jenkins建立了系统、全面的时间序列分析的研究、预测和调控的方法，其中ARIMA模型成为了最常用的进行时间序列分析、预测的模型，为时间序列分析的应用和推广奠定了基础。通常来说，ARIMA模型的建模包括3个主要的步骤：模型选择、[38，参数估计和模型诊断。通过这3个步骤，可以得到研究对象的最佳预测结果49-51]。2.3.2ARIMA(p,d,q)(P,D,Q)s模型构建基本步骤2.3.2.1序列平稳化首先，时间序列分析的对象是平稳的时间序列，平稳指的是一下方面：①在8 时间序列上，各研究时点的均值、方差为一常数；②任意两个时间点的协方差相等。当研究得到了不平稳的时间序列时，在分析前可以通过以下两种方法，进行数据转化：①差分：常用的方法是季节差分和非季节差分。②数据转换：包括对数转换和平方根转换。数据转化的目的是在构建模型前，保证时间序列的平稳，使分析和预测结果更可靠。2.3.2.2模型识别确定时间序列的平稳化以后，对于序列的自相关系数和偏自相关系数进行分[52]析，确定变量p，d，q和P，D，Q的阶值。识别模型的过程中，一般先识[52]别季节性成分，再对序列的非季节成分进行分析。按照简洁性原则初步确定模型阶数。2.3.2.3参数估计与模型诊断利用差分和季节性差分的阶数确定d，D值。当d，D确定后，应用最小二[52]乘估计法或最大似然估计法，对p，q，P，Q的值进行估算。经过参数估计所得到的ARIMA模型能否应用于所进行的研究，还需要进行模型诊断。首先模型应满足的条件是：模型的残差序列为白噪声序列。白噪声序列要求ACF、PACF与Box-LjungQ之间不具有统计学差异。若时间序列模型不满足白噪声序列，说明之前在时间序列的数据提取过程中有遗漏或错误，该模型需进行校对或改进。若有若干时间序列模型均满足白噪声序列，则需要从中找到22一个相对最优的模型，常用的评价指标平稳的R、R值越大表示模型拟合越好，RMSE、MAPE、MaxAPE、MAE、MaxAE和BIC值越小，模型的拟合效果越好。2.3.2.4模型预测检验应用构建的预测模型对吉林省2017年6～8月HFMD报吿发病数进行预测，根据相对误差判定模型预测效果。2.4HFMD发病与气象因素的相关分析简要描述长春市2014～2016年各气象因素的数据分布情况。按发病日期统计2014～2016年长春市每周的HFMD病例数。周的平均气象9 数据为每周气象数据的平均值。将周HFMD发病数分别与同期平均气压（hpa）、平均温度（℃）、平均相对湿度（％）、平均降水量（mm）、平均风速（m/s）进行单因素Spearman相关分析，探讨长春市HFMD发病与各气象因素之间的关联性。应用Spearman相关分析，研究长春市各气象因素之间的相关关系，并分析其关联程度。2.5统计学方法及应用使用Excel2010软件建立吉林省2008～2016年HFMD疫情数据库，按照发病时间进行频数统计，描述全省HFMD发病数、重症数及死亡数的时间分布特征;按照不同年份分别按照年龄、性别和职业对病例进行分层分析，计算不同人群类型的发病率及其构成比;分别按照分月和分周对各年HFMD发病数的时间分布情况进行描述；描述主要的流行地区；根据确诊病例实验室检测结果分析吉林省每年的HFMD病原体主要流行型别；对重症和死亡HFMD病例进行简要的三间分布和病原体分析。通过IBM_SPSS21.0利用2008～2017年5月月发病数建立乘积季节ARIMA模型，并对该模型进行预测，在与实际病例比较后检验该模型的预测效果。使用Excel2010软件建立长春市2014～2016年气象资料数据库，简要描述各气象因素的数据分布情况，使用IBM_SPSS21.0软件进行Spearman相关分析探讨长春市2014～2016年手足口病发病数与平均气压（hpa）、平均温度（℃）、平均相对湿度（％）、降水量（mm）、平均风速（m/s）的相关关系、各个气象因素之间的相关关系及相关程度。10 第3章结果3.1HFMD三间分布2008～2016年吉林省HFMD累计发病146538例，年均发病率为59.34/102万，不同年份间HFMD发病率具有统计学差异(X=30287.02，P＜0.001)。2008～2016年呈现隔1年流行强度增强的特点。2008～2016年吉林省HFMD共报告重症1186例，死亡35例。2010年、2011年和2014年重症发生率(重症数/病例数)2高于其他年份，不同年份之间的重症率差异有统计学意义（X=477.21，P＜0.001）。除2013年、2015年和2016年外各年均有死亡病例报告。详见表3-1。表3-12008～2016年吉林省手足口病疫情情况发病率重症发生率年份发病数重症数死亡数人口总数（1/10万）（%）200898894062729994536.220.40200925258117122733994892.390.46201028838309727395576105.271.0720111991125242746265172.501.2720121201110822749409643.690.90201375722202750401027.530.2920142051532642751280074.571.5920156871802752379824.960.12201615673402753300056.920.033.1.1HFMD人群分布3.1.1.1性别、年龄分布2008～2016年吉林省HFMD患者主要集中在5岁以下儿童，占79.95%(117155/146538)，尤其2岁年龄组病例最多，占总病例数的21.99%11 （32218/146538）。详见表3-2。男性86868例，女性59670例，男女性别比为1.46:1。2008年吉林省报告的HFMD病例中，男性5979例，女性3910例，男女性别比为1.53:1；2009年，吉林省报告的HFMD病例中男性15035例，女性10223例，男女性别比为1.47:1；2010年吉林省报告的HFMD病例中，男性17165例，女性11673例，男女性别比为1.47:1；2011年吉林省报告的HFMD病例中，男性11820例，女性8091例，男女性别比为1.46:1；2012年吉林省报告的病例中，男性7257例，女性4754例，男女性别比为1.53:1；2013年吉林省报告的HFMD病例中，男性4567例，女性3005例，男女性别比为1.52:1；2014年吉林省报告的HFMD病例中男性11800例，女性8715例，男女性别比为1.35:1；2015年吉林省报告的HFMD病例中，男性4008例，女性2863例，男女性别比为1.46:1；2016年，吉林省报告的HFMD病例中，男性9237例，女性6436例，男女性别比为1.44:1。表3-22008～2016吉林省手足口病5岁以下儿童年龄构成（岁）2008200920102011201220132014201520160~345976163411647857407804237961~1188279149503594208518273566137029262~2327542663614195253216914463154636773~2305583761904039249213474082139432994~16934506438628541799871287688621413.1.1.2职业分布2008～2016年吉林省HFMD病例以散居儿童最多，占发病例数的58.68%(85984/146538)；幼托儿童位居第二位，占发病例数的34.66%(50791/146538)；学生位居第三位，占6.09%(8931/146538)。除2016年外各年的发病人群构成也均以散居儿童为主。详见表3-3。12 表3-32008～2016吉林省手足口病职业分布幼托儿童散居儿童学生其他年份发病数构成比（%）发病数构成比（%）发病数构成比（%）发病数构成比（%）2008374437.86545355.146346.41580.592009826932.741528260.5015686.211390.552010935432.441785261.9014855.151470.512011644732.381220761.3111585.82990.502012371030.89765263.715944.95550.462013199226.31523669.152993.95450.592014567627.671310063.8616287.941110.542015222132.32418560.914126.00530.772016937859.84501732.0111537.361250.80合计5079134.668598458.6889316.098320.573.1.2HFMD时间分布3.1.2.1HFMD月分布2008～2016年呈现隔1年流行强度增强的特点。全年各月均有病例报告，各年发布趋势基本相同。4月开始发病数上升，高发季节为6～8月，10月后疫情下降，除2016年发病高峰出现在8月外，其他年度发病高峰均出现在7月。详见图3-1。图3-1吉林省2008～2016年手足口病报告发病数月分布图13 3.1.2.2HFMD周分布2008～2016年，吉林省各年周病例数的最低值在0～10之间，最高值在501～2071之间，周病例平均值为312例。详见表3-4。表3-42008～2016吉林省手足口病周数据描述性统计年份均值±标准差最小值P25P50P75最大值2008190.50±199.25051642927292009477.28±524.5618017596317582010554.60±588.2734840381520712011382.77±395.85105425263614172012230.75±332.192177329313662013145.50±217.4508321507652014394.69±507.9704214956218932015129.74±157.4909622145012016301.10±371.281612452313453.1.3HFMD地区分布全省9个地区均有病例报告，年平均发病率前三位的地区分别为延边朝鲜族自治州（92.53/10万，2020例）、通化市（89.57/10万，2246例）和四平市（70.19/10万，2500例）。2008年发病率前三位的地区分布为辽源市（67.75/10万，871例）、通化市（59.83/10万，1401例）、延边朝鲜族自治州（56.58/10万，1244例）；2009年发病率前三位的地区分布为通化市（160.61/10万，3763例）、松原市（124.40/10万，3581例）、延边朝鲜族自治州（122.84/10万，2696例）；2010年发病率前三位的地区分布为延边朝鲜族自治州（150.66/10万，3303例）、四平市（133.54/10万，4521例）、白城市（132.60/10万，2768例）；2011年发病率前三位的地区分布为白城市（170.48/10万，3466例）、通化市（141.83/10万，3298例）、四平市（92.73/10万，3140例）；2012年发病率前三位的地区分布为延边朝鲜族自治州（98.65/10万，2247例）、四平市（68.04/10万，2321例）、通化市（64.19/10万，1451例）；2013年发病率前三位的地区分布为通化市（72.61/10万，1698例）、14 四平市（35.10/10万，1202例）、松原市（31.74/10万，915例）；2014年发病率前三位的地区分布为延边朝鲜族自治州（134.38/10万，3046例）、通化市（126.88/10万，2974例）、吉林市（108.22/10万，4767例）；2015年发病率前三位的地区分布为延边朝鲜族自治州（49.75/10万，1128例）、白城市（42.66/10万，864例）、通化市（30.15/10万，707例）；2016年发病率前三位的地区分布为延边朝鲜族自治州（109.91/10万，2494例）、吉林市（79.69/10万，2171例）、松原市（61.76/10万，1782例）。各地区的平均发病率和平均发病数分布见图3-2。图3-2吉林省2008～2016年手足口病平均发病地区分布图3.1.4病原学特征2008～2016年吉林省实验室诊断手足口病病例总体以EV71阳性为主占35.04%（2908/8298），但各年主要流行型别不尽相同，2008～2011年及2014年以EV71阳性为主，2012年以CoxA16阳性为主，2013、2015和2016年以其他EV阳性为主。详见图3-3。15 图3-3吉林省2008～2016年手足口病病原学特征3.1.5重症及死亡分析2008～2016年吉林省共报告HFMD重症1186例，其中男性751例，女性435例，男女性别比1.73:1；散居儿童942例（占79.43%），托幼儿童223例（占18.80%），学生21例（占1.77%）；1～3岁年龄组报告的重症数较多。病原学检测显示EV71阳性617例（占52.02%），CoxA16阳性73例（占6.16%），其他EV阳性102例（占8.60%），阴性394例（占33.22%）。详见表3-5。2008～2016年吉林省共报告HFMD死亡病例35例，其中有34例为重症病例；35例死亡病例中，19例为男性，16例女性；死亡年龄均为5岁以下，其中2～岁年龄组最高，共14例；病原学检测显示EV71阳性22例，其他阳性2例，阴性11例。表3-52008～2016吉林省手足口病重症发病情况重症发生率性别比散居儿童5岁以下儿童年份重症数（%）(男：女）发病数构成比（%）发病数构成比（%）2008400.401.67:13177.5040100.0020091170.461.66:18774.3611497.4420103091.071.89:124077.6730398.0620112521.271.80:119878.5723292.0620121080.901.84:19083.339991.672013220.291.20:11881.8222100.0020143261.591.55:126882.2128888.34201580.123.00:1675.008100.00201640.033.00:14100.004100.00合计11860.811.73:194279.43111093.5916 3.2HFMDARIMA(p,d,q)(P,D,Q)s模型构建及预测3.2.1序列平稳化2008～2017年5月吉林省手足口病的原始数据序列图如图3-4所示，为非平稳时间序列，呈现明显的季节性周期，9个年度的序列值，每隔12个月为一个周期。对发病数序列进行一次差分、一次季节差分和对数转换处理后序列数据基本平稳如图3-5所示，即d=1，D=1。图3-42008～2017年5月各月原始数据序列图图3-5一阶差分、一阶季节差分和自然对数转换后序列图17 3.2.2模型识别与检验通过自相关图和偏自相关图看出，除12阶函数外，所有函数值都落入区间内，所以p、q均取0。由于是月度序列，所以s=12。在本模型中，采取从低阶到高阶逐个尝试的方法确定参数P、Q。在模型的选择过程中，根据模型最优化的原则，选择平稳的R²=0.265最大的模型ARIMA（0,1,0）（0,1,1）12比较适合。详见图3-6，表3-6。图3-6吉林省手足口病发病数自相关和偏自相关函数图表3-6备选ARIMA模型统计量预测变模型拟合统计量Liung-BoxQ(18)*模型离群值数量数平稳的R²R²正态化BIC统计量DFSigARIMA（0，1，0）10.270.1115.0611.46170.830（0，1，1）ARIMA（0，1，0）10.130.1315.0317.22170.440（1，1，0）ARIMA（0，1，0）10.240.1115.1211.44160.780（2，1，0）*模型均经过自然对数转换根据ARIMA（0,1,0）（0,1,1）12模型作残差的自相关函数图和偏自相关函数图，得到自相关函数和偏自相关函数均在置信区间内，详见图3-7；Liung-Box18 检验差异无统计学意义（Sig=0.832），可认为残差序列为白噪声。参数检验具有统计学意义（p＜0.001）。图3-7残差的自相关和偏自相关函数图3.2.3模型预测由于ARIMA模型适用短期预测效果较好，根据前文内容分析，吉林省的手足口病高发季节为6～8月，用ARIMA（0,1,0）（0,1,1）12模型拟合数据，预测吉林省2017年6月～8月的HFMD发病情况，并对预测结果进行检验。结果表明预测值与实际值符合程度较高，实际值均在预测值的95%可信区间范围，平均相对误差23.78%，预测结果见图3-8，表3-7。19 图3-8手足口病预测发病数与实际发病数比较序列预测图表3-72017年6～8月吉林省手足口病发病数预测评价结果.95%CL相对误差月份实际值预测值下限上限（%）6月41844914610737.487月1106884168281020.098月111562777241243.773.3HFMD发病与气象因素的关系探讨3.3.1气象资料描述吉林省位于中国东北部，省内气候属温带大陆性季风气候。由于全省的气候相对复杂，且目前获得的省内气象数据有限，现以吉林省省会长春为例，研究HFMD发病与五项气象因素的相关性。2014～2016年吉林省长春市周平均气压最低值为972.03hpa，周平均气压最高值为1003.59hpa，周平均气压为987.23hpa；周平均气温最低值为-21.77℃，周平均气温最高值为26.64℃，周平均气温为6.76℃；周平均相对湿度最低值为22.71％，周平均相对湿度最高值为91.00％，平均相对湿度为60.35％；周平均降水量最低值为0，周平均降水量最高值为15.13mm，周平均降水量为1.68mm；周平均风速最低值为1.43m/s，周平均风速最高值4.24m/s，周平均风速为2.62m/s。各气象因素值随时间的变化趋势详见表3-8。20 表3-82014～2016年吉林省长春市气象因素周数据描述性统计变量均值±标准差最小值P25P50P75最大值平均气压（hpa）987.23±7.75972.03980.07987.93993.661003.59平均温度（℃）6.76±13.62-21.77-7.289.3719.5126.64平均相对湿度（％）60.35±13.2322.7151.5061.8070.1491.00降水量（mm）1.68±2.810.000.050.392.1915.13平均风速（m/s）2.62±0.641.432.172.562.994.243.3.2相关性分析吉林省长春市周发病数与周平均气压、周平均气温、平均相对湿度、平均降水量和平均风速的相关系数均有统计学意义。其中HFMD周发病数与周平均气压及平均风速呈中度及弱的负相关，与周平均气温和降水量呈高度及中度正相关，与相对湿度呈弱正相关。长春市各气象因素之间也具有一定的相关性，其中平均气压与平均温度、平均降水量之间分别呈高度及中度负相关，平均温度与平均降水量之间呈中度正相关，平均相对湿度与平均降水量、平均风速之间分别呈中度正相关和中度负相关。详见表3-9。表3-92014～2016年长春市手足口病周发病数与各气象因素之间的Spearman相关系数矩阵平均气压平均温度平均相对湿度平均降水量平均风速变量发病数（hpa）（℃）（％）（mm）（m/s）平均气压*-0.491（hpa）平均温度**0.62-0.851（℃）平均相对湿度*0.44-0.110.091（％）平均降水量****0.32-0.480.400.431（mm）平均风速**-0.35-0.04-0.11-0.55-0.121（m/s）*注：p＜0.05.21 第4章讨论手足口病是一种常见的肠道传染病，2008年纳入法定传染病按照丙类报告管理，疫情多发，具有人群普遍易感性，且具有广泛的传播途径，好发于学龄前儿童，易产生暴发流行。轻症患者可自愈，重症患者可累及心肺神经系统，甚至威胁生命。由于病原体种类繁多，可反复感染。具有明显的季节性和周期性。尚无特效的治疗药物。虽有疫苗上市，但仅限于EV71一种病毒类型。应采取综合性的防治手段防止疾病的流行，早期的监测预警能有效的提高防控工作的时效性。吉林省每年的手足口病发病率均位居法定传染病的前五位，防控形势严峻。4.1手足口病流行特征分析4.1.1时间分布特征分析吉林省2008～2016年呈现隔1年流行强度增强的特点，与其他研究结果一[53-59]致。各年度手足口病时间趋势基本相同，全年内的发病情况整体呈单峰分布，发病高峰为6月～8月，在7月达到最高峰，符合我国北方地区单个高峰、南方地区双疫情高峰的分布特点特点，与同为东北地区的黑龙江省和辽宁省分布特征[60-62]一致。出现该时间趋势主要与吉林省气候特点有关，四季分明，一年之中6-8月气温最高，湿度最大，而温度促进肠道病毒的生长繁殖。除此之外，北方冬春季节寒冷，易感染疾病的儿童外出活动减少，所以发病减少，夏季儿童活动频繁，暴露机会增加，提高患病的风险。4.1.2人群分布特征分析HFMD的发病率在性别上有差异，男性高于女性。由于男孩相比于女孩，更活泼好动，活动范围更大，接触病原体的机会大大增加。HFMD的发病人群主要为5岁以下儿童，共有117155例，占总病例数的79.95%，尤其2岁年龄组病例最多，共有32218例病例，占总病例数的21.99%。职业分布上，散居儿童最多，共85984例，占发病总例数的58.68%托幼儿童位居第二位，共50791例，占发病例数的34.66；学生位居第三位，共8931例，占发病总例数的6.09%。22 造成年龄和职业分布特征的原因可能为：（1）婴儿出生时携带母源EV71抗体，在2岁前，抗体水平逐渐下降，2岁之后再逐渐上升，到5岁时约大到50%的[63]EV71抗体阳性率水平，故5岁以下儿童由于抗体水平较低易发生手足口病感染。另外根据台湾的一项对于HFMD患者的横断面调查，婴幼儿血清EV71抗体阳性率呈现随年龄的增长而先降低后升高的趋势。该研究显示，6个月以下婴儿抗体阳性率为38~44%，7月至2岁的儿童抗体阳性率仅为0~15%，6岁以[64]上儿童血清抗体阳性率升高至50%。故2岁儿童具有最高的HFMD发病风险。[65-67]年龄分布特点与国内外报道基本一致（2）散居儿童大多由隔辈老人进行看管，不同个人和家庭的卫生习惯、文化程度和健康意识直接影响散居儿童的发病。监护人不注意个人卫生、防病知识缺乏，幼儿公共活动场所卫生条件不合格，携带病毒的监护者与幼儿之间的交叉接触等原因可能是导致散居儿童发病居高不下的因素，另外也存在一些幼儿去医院就诊时不提供托幼机构而按照散居儿童进行报告，造成病例数增加的可能。托幼儿童由于每天集中学习、游戏、就餐、睡眠易发生交叉感染，另外，部分托幼机构布局不合理、通风不畅，不定期消毒、校医责任心差、缺少晨检制度，对幼儿看护不善是导致幼托儿童高发的原因；因此要多措并举加大对手足口病防控知识的宣传力度，使散居儿童看护人员养成良好的卫生习惯，了解手足口病的基本防控措施和典型的临床表现，为疾病的预防和早期发现提供基础。针对托幼机构和学校，各相关单位和部门应强化督导和检查工作，提高管理者和教师的防控意识，严格落实晨检制度。4.1.3病原学分布特征分析2008～2016年吉林省实验室诊断手足口病病例总体以EV71阳性为主占35.04%，各年主要流行型别不尽相同，2008～2011年及2014年以EV71阳性为主，2012年以CoxA16阳性为主，2013、2015和2016年以其他EV阳性为主。从病原构成可以看出，每年均为多种病毒共存，但在HFMD高发的年份，主要的病原体均为EV71，同时在重症及死亡病例中，大部分为EV71感染，这一结[68-70]果与已有研究相符。病原体型别在不同年份发生的变化，可能与不同型别病原体基因组变化有关，需进一步对病原体基因特征进行研究。23 4.2ARIMA模型的建立时间序列是一组有着时间延续性的数据序列，可以利用数据间的延续关系来预测它未来的走势，而时间序列分析就是要发现时间序列内在的统计规律，进而推算出接下来一个时间段内的数据，通过采取有效的手段提早采取措施。时间序列分析广泛的应用于金融、气象等领域。近年来，时间序列分析也越来越多的应用在疾病发展规律探索方面，尤其是对传染病发展趋势的研究具有实际意义。本研究通过对吉林省2008年～2017年5月手足口病月发病数的分析，疫情流行特点具有以年为周期的周期性和夏季高发的季节性，所以建立了ARIMA季节乘积模型，对HFMD的发病情况进行分析和预测。ARIMA(0,1,0)(0,1,1)12模型通过了白噪声检验和参数检验，所预测2017年6～8月的发病数均在预测发病数的95%置信区间内。然而，时间序列分析所建立的预测模型为适合短期预测的工具，应根据实际情况，不断纳入新数据、建立新模型，才能保证预测的精准性。4.3手足口病发病与气象因素的关系分析很多疾病的时间趋势都能看到明显的季节周期，气象因素往往是导致疾病流［71-73］行的主要因素，可使疾病的流行趋势呈现一个或几个流行高峰。Spearman相关分析显示吉林省长春市周发病数与周平均气压呈高度负相关，与周平均气温呈高度正相关，与相对湿度呈弱正相关，与降水量呈中度正相关，与平均风速呈弱负相关。之前的研究已经证明温度、湿度均为影响HFMD流行的因素，HFMD的发病率随着环境温度、湿度的升高而增加，本研究得出的结果与该研究一致[74-80]。气象因素可以通过很多种不同的方式对疾病的发病产生影响，如温度、湿度、风速可能影响病原体的生存环境和传播途径，也影响易感人群的暴露机会和[81-83]敏感性。病原体能否存活和致病力是致病的关键。温度、湿度、阳光、酸碱[84-86]度以及盐度等条件能够直接对病原体的存活产生影响。研究表明，肠道传染[87]病的发生和温度和相对湿度相关。研究发现，湿热的环境有利于肠道病毒的繁[88-89]殖和传播。另外，易感人群的行为方式的改变可能影响手足口病发病呈现明[90]显的夏季高发的季特点。长春市地处东北平原，冬春季节寒冷干燥，夏季温暖24 多雨，孩子在夏季会从事更多的户外集体活动，从而增加接触病原体的机会，在冬季，由于寒冷，儿童更多的倾向室内活动，主要是与家庭成员之间的接触，病原体传播机会减小。这与Strikas等人研究的温带气候类型的地区，肠道病毒容[91]易在夏秋季流行的研究结果一致。多项研究结果表明发现手足口病发病与气候因素间存在滞后效应。在日本的一项研究表明，手足口病的发病与0-3周前的温[92]度和湿度呈正相关。香港的研究结果显示，HFMD的发病率会受到此前2周的空气湿度与平均气温的影响，随着湿度与气温的升高，HFMD发病率也随之[93]上升；广州的研究分析了HFMD发病率与1周前温度和相对湿度的关系，也[94][95]得到了相同的结果。郑州市的一项研究表明，HFMD周发病数与平均气温、最低气温和最高气温滞后2周呈正相关。这种气候变化与HFMD发病先后发生的现象，是气候因素影响HFMD发病的最好佐证。总之，手足口病发病与平均气压（hpa）、平均温度（℃）、平均相对湿度（％）、降水量（mm）、平均风速（m/s）的相关性提示我们高热、潮湿、低风速的季节加强对手足口病的宣传和监测工作；重点加强对托幼机构的监督管理工作。HFMD的流行不是单一因素造成的，而是包括生物、行为、环境等因素综合作[96]用而产生的的，因此，在今后的防控工作中，要注意多种因素的综合分析、多措并举，有效的制定疾病的防控措施，更好的控制疾病的发生。4.4手足口病防控措施建议预防手足口病要多措并举，多管齐下。首先政府应该高度重视，组织各行政部门，合理分工，协同合作，制定有效的防控措施，同时做好工作部署，保障措施的落实。各级疾控机构要做好疫情监测，有效利用“传染病自动预警系统”和“传染病报告信息管理系统”中的疫情数据，通过分析早期发现聚集性疫情和暴发疫情的苗头，全力开展重症及死亡病例流行病学调查和疫点疫区的处置工作；针对不同人群，通过多种形式，多种渠道广泛开展宣传教育工作。医疗机构要积极做好病人的救治工作，尤其是重症病人的早期发现与治疗至关重要。在流行季节学校和托幼机构要严格执行晨午检和消毒通风等管理制度，有疫情发生时配合疾控机构的调查工作，必要时严格按照规定封园停课。对于学校和托幼机构的防控工作，卫生监督部门应进行积极进行监督、检查。社会群众尤其是幼儿家长应广泛25 参与、积极配合加强各项防控措施的落实。对于散居儿童的预防：在儿童外出后及饭前便后要及时用肥皂、洗手液等洗手，避免儿童食用生冷食物、饮用生水，避免接触HFMD患儿；婴儿奶瓶、奶嘴在使用前后都应进行清洗；患儿的看护人员在接触患儿后均应洗手，并对污物进行处理；HFMD暴发期间，避免带儿童到人群聚集的公共场所，居室内保持卫生，勤通风，勤晒衣被；儿童出现发热、疱疹等疑似HFMD症状时应及时就诊。症状较轻的HFMD患儿，宜在家治疗，以减少院内的交叉感染，但治疗期间应避免与其他儿童接触，防止传染，监护人员要对患儿的衣物及污物妥善处理。对于学校、托幼机构等集体单位的预防：保持良好的室内通风，特别是在夏秋季HFMD的高发季节，要做到每日对儿童的个人用具、餐具、玩具等物品清洗消毒，对于室内设施如楼梯扶手桌面等进行每日擦拭消毒；工作人员进行消毒工作时应佩戴手套，并在结束清洁工作后立即对双手进行消毒；教育儿童保持良好的个人卫生，学会正确的洗手方式；对机构内所有儿童每日进行晨检，当发现儿童有可疑症状，应立即送诊，并立即对其日常用品消毒处理；当患儿数量增多时，要及时向卫生监督部门汇报。应采取多种形式广泛开展宣传教育工作，利用各种媒介大力宣传手足口病预防知识，引导正确的生活和卫生习惯。针对看护幼托儿童的老年人，由于他们获取信息的方式局限，对手足口病的认识较少，防病意识不高，要加强电视、广播、报纸等传统媒体方式对手足口病的宣传，还可以通过社区/农村宣讲，发放通俗易懂，形象生动的宣传手册，播放形象生动的宣传动画等方式广泛普及手足口病的预防知识，使看护人能够了解手足口病的临床症状，掌握最基本的预防措施，改正不良的卫生习惯。在学校、幼托机构可由校医设立健康教育课堂，教导学生养成正确的卫生习惯。医疗机构可在候诊区张贴宣传海报，使就诊的患者家长能够更多的了解手足口病的防控常识。26 第5章结论1.吉林省2008～2016年手足口病在时间分布上呈现明显的季节性，高发季节为6～8月；高发地区为延边朝鲜族自治州、通化市和四平市；人群分类上以散居儿童为主，年龄5岁以下儿童为主，男性高于女性；实验室诊断病例总体以EV71阳性为主；重症及死亡病例分布特征与整体一致。2.用吉林省2008～2017年5月手足口病疫情数据拟合出ARIMA（0,1,0）（0,1,1）12模型，对2017年6～8月病例预测效果较好。3.长春市2014～2016年手足口病周发病数与平均气压（hpa）、平均温度（℃）、平均相对湿度（％）、降水量（mm）、平均风速（m/s）存在相关性。27 参考文献[1]GOKSUGURN,GOKSUGURS.Imagesinclinicalmedicine.Hand,foot,andmouthdisease[J].TheNewEnglandandjournalofmedicine20103629(14)：e49.[2]吕晓菊.手足口病防治进展[J].华西医学,2008,23(3):236,336.[3]万朝敏,徐爱丽.手足口病的研究进展[Ｊ/ＣＤ].中华妇幼临床医学杂志(电子版),2009,5(2):106-108.[4]闫莉.长春地区手足口病地理流行病学及EV71基因亚型研究[D]．长春：吉林大学,2015．[5]TanEL,ChowVT,QuakSH,etal.Developmentofmultiplexreal-timehybridizationprobereversetranscriptasepolymerasechainreactionforspecificdetectionanddifferentiationofEnterovirus71andCoxsackievirusA16[J].DiagnosticMicrobiology&InfectiousDisease,2008,61(3):294-301.[6]MCMINNP,LINDSAYK,PERERAD,etal.Phylogeneticanalysisofenterovirus71strainsisolatedduringlinkedEpidemicsinMalaysia,Singapore,andWesternAustralia[J].JournalofVirology,2001,75(16):7732–7738.[7]CARDOSAMJ,PERERAD,BROWNBA,etal.MolecularepidemiologyofHumanenterovirus71strainsandrecentoutbreaksintheAsia-Pacificregion:ComparativeanalysisoftheVP1andVP4genes[J].EmergingInfectiousDiseases,2003,9(4):461–468.[8]SHIMIZUH,UTAMAA,ONNIMALAN,etal.MolecularepidemiologyofEnterovirus71infectionintheWesternPacificRegion[J].PediatricsInternational,2004,46(2):231–235.[9]BIBLEJM,PANTELIDISP,CHANPK,etal.Geneticevolutionofenterovirus71:epidemiologicalandpathologicalImplications[J].ReviewsinMedicalVirology,2007,17(6):371–379.[10]中华人民共和国卫生部.手足口病预防控制指南(2009版)[S].北京：中华人民共和国卫生部,2009.28 [11]黎祖秋.手足口病流行及防治研究进展[J].应用预防医学,2015,21(1):59-62[12]WITSØE1,PALACIOSG,RØNNINGENKS,etal.AsymptomaticcirculationofHEV71inNorway[J].VirusResearch,2007,123(1):19–29.[13]刘立,陆辉,田会芳.石家庄市赞皇县健康儿童肠道病毒携带情况分析[J].中国妇幼保健,2012,27(11):342-351.[14]李洋,张湘萍,包小兵.2010-2011年河南省安阳市健康人群手足口病病原体隐性感染调查及分析[J].中国卫生检验杂志,2013,1(9):123-132.[15]WITSE,PALACIOSG,CINEKO,etal.AsymptomaticcirculationofHEV71inNorway[J].VirusResearch,2007,123(1):19-29.[16]万俊峰,朱理业,刘红曹.阜阳市手足口病（EV71感染）疫情流行病学分析[J].安徽医学,2008,29(4):344–345.[17]RobinsonCR,DoaneFW,RhodesAJ.Reportofanoutbreakoffebrileillnesswithpharyngeallesionsandexanthem:Torontosummer1957;isolationofgroupACoxsackievirus[J].CanMedAssocJ,1958,79(8):615-621.[18]YamasakiK,OkunoY.Geneticdiagnosisandmolecularepidemiologicalanalysisofhand-foot-mouthdiseasewhichprevailedinOsakaPrefecturein2000[J].KansenshogakuZasshi,2001,75(11):909-915[19]沙爱龙,刘颖.手足口病的研究概况[J].生命科学仪器,2007,11(5):13-15.[20]蒋心华,潘孝彰,王岱明.现代感染病学[M].上海:上海医科大学出版社,1998:642-656.[21]文育英.健康教育在手足口患儿家属中的应用[J].护理实践与研究,2009,6(12):114-115.[22]ParkSK,ParkB,KiM,etal.Transmissionofseasonaloutbreakofchildhoodenteroviralasep-ticmeningitisandhand-foot-mouthdisease[J].JKoreanMedSci,2010,25(5):677-683.[23]徐文体,高璐,张颖,等.广东省手足口病流行特征和危险因素研究[J].中华流行病学杂志,2009,30(1):100-101.[24]余吉仙,黄劲松,吕铁峰.杭州1000例手足口病分析[J].浙江预防医学,2009,21(7):21-23.29 [25]ShinoharaM,UchidaK,ShimadaS,etal.Characterizationofenterovirustype71isolatedinSaitamaPrefecturein2000[J].KansenshogakuZasshi,2001,75(6):490-494.[26]YamazakiK,OkunoY.Geneticdiagnosisandmolecularepidemiologicalanalysesofhand,footandmouthdiseasewhichprevailedinOsakaPrefecturein2000[J].KansenshogakuZasshi,2001,75(11):909-915.[27]ChangLY,KingCC,HsuKH,eta1.Riskfactorsofenterovirus71infectionandassociatedhand,footandmouthdisease/herpanginainchildrenduringanepidemicinTaiwan[J].Pediatrics,2002,109(6):88.[28]HsuB,ChenC,WanM,eta1.DetectionandidentificationofenterovirusesfromvariousdrinkingwatersourcesinTaiwan[J].JHydrol,2009,365:134-139.[29]廖炯,黄萍.手足口病的预防与治疗进展[J].中国医药指南,2012(17):463-464.[30]于长水.天津市由coxAl6肠道病毒引起的手足口病流行[J].中华流行病学杂志,1985,6(2):66-69．[31]BendigJW,FlemingDM.Epidemiological,virological,andclinicalfeaturesofanepidemicofhand,foot,andmouthdiseaseinEnglandandWales[J].CommunDisRepCDRRev,1996,6(6):R81-86.[32]朱渭萍,储强,王勇,等.上海市浦东新区2005-2007年手足口病疫情监测分析[J].上海预防医学杂志,2008,20(8):372-375.[33]韩秀珍,李化兵.手足口病的诊断与治疗－手足口病的病原学与流行病学[J].山东医药,2008,48(17):109.[34]HerreroLJ,LeeCS,HurrelbrinkRJ,etal.Molecularepidemiologyofenterovirus71inpeninsularMalaysia,1997-2000[J].ArchVirol,2003,148(7):1369-1385.[35]SCHMIDTNJ,LENNETTEEH,HOHH.Anapparentlynewenterovirusisolatedfrompatientswithdiseaseofthecentralnervoussystem[J].JournalofInfectiousDiseases,1974,129(3):304–309.[36]孙立梅,邓爱萍,康敏,等.2009年广东省手足口病聚集性疫情流行特征分析[J].华南预防医学,2010,36(2):1-4.30 [37]刘江艺,陈明春,龚苍涛,等.泉州市2010年手足口病流行特征分析[J].海峡预防医学杂志,2011,1(3):16-20.[38]王燕.应用时间序列分析[M].中国人民大学出版社,2005:1-239.[39]YuleG.U,OnamethodofinvestigatingPeriodicitiesinDisturbedSeries,withSpecialReferencetoWolfer'sSunspotNumbers[J].PhilosophicalTransactionsoftheRoyalSocietyofLondon.SeriesA,1927,226:267-298.[40]WalkerG,OnPeriodicityinSeriesofRelatedTerms[J].ProceedingsoftheRoyalSociety,SeriesA,1991,131(818):518-532.[41]SlutskyE.TheSummationofRandonCausesAStheSourceofCyclicProcesses[J].Econometrica,1977,5(2):105-146.[42]WoldH.AStudyintheAnalysisofStationaryTimeSeries[M].Stockholm;AlmqvistandWiksell,[1938]1984:1-132.[43]BoxG.E.P,JenkinsG.M.Time-SeriesAnalysis:ForecastingandControl[M].SanFrancisco:HoldenDay,2007:1-70.[44]李守丽.时间序列模型在地级市GDP预测中的应用[D].导师：王黎明.郑州大学,2013.[45]徐先金,姜茂启.基于灰色预测模型的上海市实际利用外资预测[J].工业技术经济,2009.[46]刘涛,王显军,姜宝法.SARIMA模型预测山东省手足口病发病趋势[J].中国卫生统计,2013,1(5):697-700.[47]颜玉炳,郑惠能,刘红莲.厦门市手足口病流行特征分析及发病趋势预测方法研究[J].社区医学杂志,2014,3(21):16-20.[48]李丽丽.北京市房山区手足口病与气象因素的关系及预测模型研究[D].导师：刘起勇,中国疾控中心传染病所,2015.[49]王重建,魏晟,刘建平,等．湖北省疾病预防控制人员应急能力评价及影响因素分析［J］．中国公共卫生,2007,23(6):731–733．[50]刘刚,唐宋,孙文杰．时间序列分析法在香港结核病预测中的应用［J］．中国卫生统计,2012,29(2):226–228．31 [51]温亮,徐德忠,林明和,等．应用时间序列模型预测疟区疟疾发病率［J］．第四军医大学学报,2004,25(6):507–510．[52]杨召,叶中辉,尤爱国,等．乘积季节ARIMA模型在结核病发病率预测中应用［J］．中国公共卫生,2013,29(4):469–472．[53]申珂,邓健,郭娜娜,等.2009-2014年河北省邯郸市手足口病流行特征研究[J].中国药物与临床,2015,15(12):1728-1730.[54]曹志冬,曾大军,王全意,等.2008年北京市手足口病流行的规律与传播机制［J］.科学通报,2010,55(9):764-772．[55]许一玲,李栋.淄博市张店区手足口病流行与气象因素的关系研究［J］.热带医学杂志,2010,10(10):1237-1239．[56]张德汴.气象因子与开封市手足口病响应分析［J］.气象与环境科学,2011,34(21):135-137．[57]胡丽莉,董金荣．武威市凉州区手足口病流行与气象条件的关系［J］．甘肃科技,2011,27(17):184-185．[58]刘立,郭建花,张世勇,等.石家庄市手足口病发病与气象因素相关分析［J］.实用预防医学,2011,18(8):1389-1391．[59]刘立,韩江涛,庞志钊,等.石家庄市气象因素与手足口病流行的关系研究［J］．中国全科医学,2013,16(3):326-328．[60]任敏,张凯.2008-2010年我国手足口病发病特点分析[J].中国健康教育,2011,27(8):568-570+581.[61]陈淑红,王开利,陈露菲,等.黑龙江省2009-2010年手足口病疫情监测结果分析[J].中国公共卫生管理,2012,28(1):66-67.[62]姚文清,陈静乙,李鑫,等.辽宁省2009年手足口病流行病学分析[J].中国公共卫生,2010,26(12):1589-1590.[63]WILDER-SMITHA,FOOW,EARNESTA,etal.SeroepidemiologyofdengueintheadultpopulationofSingapore[J].TropacalMedicine&InternationalHealth,2004,9(2):305-308.[64]JEEYM,CHEONDS,KIMK,etal.GeneticanalysisoftheVP1regionofhumanenterovirus71strainsisolatedinKoreaduring2000[J].Archivesof32 Virology,2003,148(9):1735–1746.[65]XINGW,LIAOQ,VIBOUDC,et,al.Hand,foot,andmouthdiseaseinChina,2008-12:anepidemiologicalstudy[J].TheLancetInfectiousdiseases,2014,14(4):308-318.[66]TSENGFC,HUANGHC,CHICY,et,al.EpidemiologicalsurveyofenterovirusinfectionsoccurringinTaiwanbetween2000and2005:analysisofsentinelphysiciansurveillancedata[J].Journalofmedicalvirology,2007,79(12):1850-1860.[67]王江敏.2010年北京市海淀区手足口病流行病学分析[J].实用预防医学,2012,19(2):185-187.[68]MOUJ,DAWESM,LIY,et,al.Severehand,footandmouthdiseaseinShenzhen,SouthChina:Whatmattersmost?[J].Epidemiologyandinfection,2014,142(4):776-788.[69]刘白薇,王全意,李锡太,等．北京市2010年手足口病死亡病例流行病学及临床特征分析[J].现代生物医学进展,2011,11(22):4264-4266.[70]王英,何小周,赵俊伟,等．2008-2010年全国报告手足口病死亡水平及死亡病例分析[J]．疾病监测,2011,26(6):424-426.[71]FaresA．Factorsinfluencingtheseasonalpatternsofinfectiousdiseases［J］．IntJPrevMed,2013,4(2):128[72]GarskeT,FergusonNM,GhaniAC．EstimatingairtemperatureanditsinfluenceonMalariatransmissionacrossAfricaJ］．PLoSOne,2013,8(2):e56487.[73]HuangY,DengT,YuS,etal．Effectofmeteorologicalvariablesontheincidenceofhand,foot,andmouthdiseaseinchildren:atime-seriesanalysisinGuangzhou,China［J］．BMCInfectDis,2013,13:134.[74]ZhuangDF,HuWS,ＲenHY,etal．Theinfluencesoftemperatureonspatiotemporaltrendsofhand-foot-and-mouthdiseaseinmainlandChina［J］．IntJEnvironHealthRes,2014,24(1):1.[75]BoYC,SongC,WangJF,etal．Usinganautologisticregressionmodeltoidentifyspatialriskfactorsandspatialriskpatternsofhand,footandmouthdisease33 (HFMD)inMainlandChina［J］．BMCPublicHealth,2014,14:358.[76]OnozukaD,HashizumeM．Theinfluenceoftemperatureandhumidityontheincidenceofhand,foot,andmouthdiseaseinJapan［J］．SciTotalEnviron,2011,410-411:119.[77]ONOZUKAD,HASHIZUMEM.TheinfluenceoftemperatureandhumidityOntheincidenceofhand,foot,andmouthdiseaseinJapan[J].ScienceoftheTotalEnviroment,2011,41(31):119–125.[78]MAE,LAMT,WONGC,etal.Ishand,footandmouthdiseaseassociatedwithmeteorologicalparameters?[J].EpidemiolInfect,2010,138(7):1779–1788.[79]WANGY,FENGZ,YANGY,etal.Hand,FootandMouthDiseaseinChina:PatternsofSpreadandTransmissibilityduring2008–2009[J].Epidemiology,2011,22(19):781–792.[80]HUANGY,DENGT,YUS,etal.Effectofmeteorologicalvariablesontheincidenceofhand,foot,andmouthdiseaseinchildren:atimeseriesanalysisinGuangzhou,China[J].BMCInfectDis,2013,1(1):134-140.[81]PIRTLEEC,BERANGW.Virussurvivalintheenvironment[J].RevSciTech,1991,10(2):733–748.[82]SIY,WANGT,SKIDMOREAK.EnvironmentalfactorsinfluencingthespreadofthehighlypathogenicavianinfluenzaH5N1virusinwildbirdsinEurope[J].EcolSoc,2010,15(1):26-30.[83]DOWELLSF.Seasonalvariationinhostsusceptibilityandcyclesofcertaininfectiousdiseases[J].EmergInfectDis,2001,79(2):369–374.[84]BROWNJD,GOEKJIANG,POULSONR.Avianinfluenzavirusinwater:infectivityisdependentonpH,salinityandtemperature[J].VetMicrob,2009,136(3):20–26.[85]RZEZUTKAA,COOKN.Survivalofhumanentericvirusesintheenvironmentandfood[J].FEMSMicrobiolRev,2004,28(2):441–453.[86]FISCHERTK,STEINSLANDH,VALENTINER-BRANTHP.Rotavirusparticlescansurvivestorageinambienttropicaltemperaturesformorethan2months[J].JClin34 Microbiol,2002,40(8):4763–4764.[87]ONOZUKAD,HASHIZUMEM.Theinfluenceoftemperatureandhumidityontheincidenceofhand,foot,andmouthdiseaseinJapan[J].ScienceoftheTotalEnviroment,2011,41(8):119–125.[88]刘凤仁,刘渠,张起文,等．深圳市龙岗区2005～2009年手足口病暴发疫情分析［J］．中华疾病控制杂志,2010,14(7):631-634．[89]刘立,郭建花,张世勇.石家庄市手足口病发病与气象因素相关分析［J］．实用预防医学,2011,18(8):1389-1390．[90]SUMINSKIRR,POSTONWC,MARKETP,etal.Meteorologicalconditionsareassociatedwithphysicalactivitiesperformedinopen-airsettings[J].Internationaljournalofbiometeorology.2008,52(3):189-197.[91]STRIKASRA,ANDERSONLJ,PARKERRA.TemporalandgeographicpatternsofisolatesofnonpolioenterovirusintheUnitedStates,1970-1983[J].TheJournalofinfectiousdiseases,1986,152(2):346-351.[92]ONOZUKAD,HASHIZUMEM.TheinfluenceoftemperatureandhumidityOntheincidenceofhand,foot,andmouthdiseaseinJapan[J].ScienceoftheTotalEnviroment,2011,41(31):119–125.[93]MAE,LAMT,WONGC,etal.Ishand,footandmouthdiseaseassociatedwithmeteorologicalparameters?[J].EpidemiolInfect,2010,138(7):1779–1788.[84]HUANGY,DENGT,YUS,etal.Effectofmeteorologicalvariablesontheincidenceofhand,foot,andmouthdiseaseinchildren:atimeseriesanalysisinGuangzhou,China[J].BMCInfectDis,2013,1(1):134-140.[95]冯慧芬,赵秋民,段广才,等.应用时间序列分析气象因素对手足口病流行的影响［Ｊ］．郑州大学学报（医学版）,2015,50（2）：171-176.[96]周智翔,江朝强．严重急性呼吸综合征暴发的大气环境和院内感染的职业卫生因素［Ｊ］．中华劳动卫生职业病杂志,2004,22（04）：261-263.35 作者简介及在学期间所取得的科研成果作者简介栾博，女，汉族，1987年3月出生于吉林省磐石市，2010年6月毕业于郑州大学预防医学专业，获得医学学士学位，2014年9月就读于吉林大学公共卫生学院，攻读公共卫生硕士学位。在学期间所取得的科研成果发表的学术论文：黄飚，栾博，史玉婷，戚潇丹，赵庆龙.吉林省2008-2014年手足口病流行病学特征分析[J].中国公共卫生，2017，33（2）：291-294.36 致谢这篇论文能够完成首先要感谢我的导师李波教授，从选题设计到分析整理最后到修改成稿，都离不开老师的精心指导，在选题的过程中，李老师根据我所从事的手足口病防控工作精心的帮我选择了手足口病的时间序列分析这个方向，借着论文撰写的契机，使我能够有机会全面的学习时间序列分析的知识，借鉴许多前辈的宝贵经验，这些都将成为我未来工作中的宝贵财富。李老师不仅是我的专业导师，更是我的人生导师，他渊博的专业知识，严谨的工作作风，都是我学习的榜样，同时要感谢我的同事和朋友们，感谢你们在我研究过程中提出的宝贵意见，感谢田鑫师兄、高春时师弟帮助我解决了很多在数据分析过程中出现的问题，最后感谢我的家人在这两个多月的论文撰写过程中能够全力的支持我！37