ARIMA模型在湖南省常德市手足口病发病预测中应用.docVIP

ARIMA模型在湖南省常德市手足口病发病预测中的应用 　　【摘要】目的：构建ARIMA模型预测湖南省常德市手足口病发病趋势，为早期防控提供科学依据。方法：分析处理2010-2015年常德市手足口病发病率时间序列，构建ARIMA模型，对2 016-2 017年常德市手足口病发病趋势进行预测。结果：通过拟合优度检验和残差序列白噪声检验，得出本研究最优模型为ARIMA（1，0，1）（0，1，1），正态化BIC=4.235，Ljung-Box=12. 019（ P=O. 678），2016年观测值均位于预测值95%可信区间内，与预测值的相对误差为0.72%-35. 09%，预测2017年常德市手足口病月发病率在6月出现最大值，为24. 03/10万。结论：该ARIMA模型能够较好拟合常德市手足口病发病趋势，预测效果良好。 　　【关键词】ARIMA模型预测；湖南省常德市；手足口病 　　手足口病（hand-foot-mouth disease，HFMD）是一种由多种肠道病毒（主要为EV71和COX A16）引起的儿童常见传染病，以发热和手、足、口等部位出现皮疹和疱疹为主要特征；该病好发于五岁以下婴幼儿，容易在幼儿园及托幼场所引起疫情爆发，是我国法定的丙类传染病。手足口病的发生具有趋势性、季节性、短期波动等特点，本文通过收集我国传染病信息管理系统中湖南省常德市2010-2016年手足口病疫情信息，建立时间序列分析中的自回归综合移动平均模型（ARIMA模型），通过选择最佳模型对2017年常德市手足口病发病情况进行预测，为常德市疾病预防控制中心早期防控手足口病提供科学依据。 　　1 资料与方法 　　1.1 资料来源 　　疫情数据来自常德市疾病预防控制中心2010-2016年常德市手足口病疫情资料数据库。 　　1.2 研究方法 　　通过对2010-2015年常德市手足口病发病情况进行分析并构建ARIMA模型，选择最优模型比较常德市2016年手足口发病率预测值与实际观测值，评价所建立模型的预测效果，对2017年常德市手足口病发病率进行预测。 　　1.2.1 ARIMA模型 　　ARIMA模型又称B-J方法，是一种时间序列预测方法，是指将非平稳时间序列转化为平稳时间序列，然后将因变量仅对它的滞后值以及随机误差项的现值和滞后值进行回归所建立的模型。对存在季节性和周期性变化的时间序列建立ARIMA模型，能够分析时间序列的随机性、平滑性和季节性，在分析的基础上选择适当的模型，可以对具有季节性变动趋势的传染病发病情况进行预测。 　　ARIMA模型有三个参数：p，d和q，它们分别表示模型的自回归阶数、差分的阶数和移动平均阶数。如果时间序列存在季节性趋势，可采用乘积季节性模型[ARIMA （p，d，q） （P， D， Q）。 参数P表示季节性自回归阶数，O表示季节性的移动平均阶数，D表示季节差分的阶数，S表示季节周期。 　　ARIMA模型的建模过程：（1）序列特征及平稳化。应用时间序列图、自相关系数函数图（ACF图）、偏自相关系数函数图（PACF图）分析序列，ARIMA建模的前提是平稳时间序列，对于非平稳的序列要经过差分或对数变换转变为平稳序列。（2）模型的识别和定阶。采用Box-Jenkins定阶法和BIC准则相结合。 （3）利用非线性最小二乘法估计模型参数。好的模型序列残差是随机的，为白噪声序列，常用Box-Ljung Q统计量来检验。（4）判断模型的拟合优度。使用R2、标准化的BIC、平均绝对误差百分比（MAPE）比较模型的拟合效果，当R2较高、BIC值和MAPE值较小时，该模型为相对最优模型。（5）评估及应用。将最优模型预测值与观测值比较，计算相对误差来评估预测效果，对该时间序列未来发展趋势进行分析预测。 　　1.3 统计分析 　　采用Excel建立2010-2016年湖南省常德市手足口病发病情况数据库，使用SPSS20.0软件中的时间序列模型进行资料分析，检验水准a=0.05。 　　2 结果 　　2.1 手足口病发病率 　　对湖南省常德市2010年1月至2015年12月手足口发病率绘制时间序列图，见图1，可见湖南省常德市手足口发病率存在明显的季节性分布规律，以年为流行周期，每年存在两个流行高峰，大高峰出现在每年的5-7月，小高峰出现在每年的9-11月，每年1-2月为发病低谷。最高峰值出现在2014年5月，达到44.57/10万。由于该时间序列存在明显的季节性趋势，季节周期为12个月，为非平稳序列，为了排除方差波动的干扰，实现序列平稳化，需要对序列进行差分处理。季节性ARIMA模型有一般差分和季?性差分，由于该研究发病率序列季节性明显，需对序列进行一阶季节性差分，对一阶季节性差分后的新时间序列做序列图，见图2，经一阶季