手足口病发病趋势预测及传播动力学模拟研究.docxVIP

手足口病发病趋势预测及传播动力学模拟研究 目的深入研究和了解手足口病的流行病学特征和病毒分布特征 及危险因素 , 预测和掌握该病的流行趋势及规律 , 对于更好的防控手 足口病 , 保护学龄前儿童身心健康 , 完善防控手段和制定防控策略具 有重要意义。方法 1. 利用 GIS 软件绘制县级地区发病率分布情况图。 利用 Geocoding 软件对所有报告病例进行经纬度定位 , 再导入 GIS 软件开展核密度分析 , 绘制核密度和点密度图。 2. 开展个体危险因素和群体危险因素现场调查 , 利用 Epidata 软件建立数据库 , 再利用 SAS统计软件进行多因素 Logistic 逐步回归分析 , 并计算相关危险因素的 OR值及 95%可信区间 , 筛选个体及群体危险因素。 3. 依据 2008 年～ 2016 年手足口病各月发病率数据 , 利用 ARIMA模型进行建模和预测 , 产出线性预测序列和残差序列。依据残差数据序列 , 利用 NARNN模型 进行建模和预测 , 进一步产出残差预测序列。 再将两者结果进行综合 , 最终形成发病率预测结果序列。 4. 根据 2008～2016 年各月发病率特征, 按照发病率高低将手足口病月发病率（ 1/10 万）划分为 3 个状态（发病率 lt;7.5;7.5 ≤发病率 lt;10; 发病率≥ 10）。再建立马尔科 夫链模型计算预测区间转移概率 , 得出 2017 年各月发病率所处状态 预测结果。 5. 利用传播动力学方法建立 SIR 模型 , 结合全局敏感性和 不确定性分析考察各因素间影响程度 , 利用 Matlab 软件实现疫情模 拟计算过程 , 从而判断疫情走势。结果 1. 流行特征方面 : 宜昌市 2008～2016 年间共报告手足口病病例 37988 例, 报告发病率总体呈锯 齿状 , 偶数年发病率较高而奇数年较低。全年发病呈 4～6 月和 11～ 月双高峰。长江水系流经乡镇和非流经乡镇的月累计报告发病数 差异有统计学意义（ Xsup2/sup=117.93,Plt;0.0001 ）。报告病 例点密度分析集中趋势同各县级累计报告病例数占总病例数的百分 比结果一致 , 排名前三位的地区为西陵区、夷陵区和宜都市。人群不 同性别（ Xsup2/sup=30.90,Plt;0.0001 ）、类别 Xsup2/sup=52.99,Plt;0.0001 ）和年龄 Xsup2/sup=4.06,P=0.0440 ）分布差异具有统计学意义。九年间检测阳性 3162 例, 其中 EV71病毒感染 2226 例, 占实验室检测病例数的 70.40%;CoxA16病毒感染 454 例, 占 14.36%;其他肠道病毒感染 482 例, 占 15.24%。EV71病毒检测阳性病例核密度分布主要集中于宜昌市城区和夷陵区晓溪塔街道。 其他肠道病毒检测阳性病例主要集中于秭归县茅坪镇和长阳县龙舟坪镇。 CoxA16病毒分布范围较 EV71病毒最广 , 主要集中于宜昌市城区和兴山县古夫镇。 2. 个体危险因素分析方面 : 居住在低层（ 1～3 层）为患病者的保护因素（ OR=0.698）。居住地周边有公共厕所是患病者的危险因素（ OR=1.294）。通过传统媒体了解疫苗信息（ OR=2.785）, 手足口病知识知晓情况评分为 0 分是不患病者的危险因素（ OR=2.535）。人工喂养方式（ OR=1.837）,身体状况良好（OR=3.573）是不患病者的危险因素。 参加课外学习班 , 能独立排便者 , 有家族史者 , 有过敏史者 , 共用水杯者 , 啃咬手指者 , 最 10 天去过长途汽车或者火车站者以及最近 10 天去过商场超市者是手足口病不患病者的保护因素 , 其 OR值分别为 0.648 、0.473 、0.315 、 0.568 、0.670 、0.566 、0.492 和 0.802 。3. 群体危险因素分析方面 : 教职工人数与幼儿数配比高于 1:5 标准 , 教职工人数与幼儿数配比在 1:5 ～1:7 标准范围内 , 幼儿园配备医务室以及幼儿园有专车接送学生 者是手足口病发病率较低地区的保护因素 , 其 OR值分别为 0.079 、 0.193 、0.175 和 0.101 。公共厕所有灭蝇蛆措施和粪池加盖者是手足 口病发病率较低地区的危险因素 , 其 OR值分别为 14.346 和 10.629 。 预测和模拟方面 : 根据拟合优度（BIC=3.022）和 LB 统计量（P=0.334） 选定最优模型为 ARIMA（1,0,5 ）×（ 1,1,0 ）sub1/sub2。ARIMA 模型预测结果同实际值之间存在一定差异 , 但对发病高峰等基本趋势 预测相符。根据均方误差、相关系数（ R=