时间序列自回归综合移动平均模型在新城疫疫情预测中的应用.doc

PAGE PAGE 4 时间序列自回归综合移动平均模型在新城疫疫情预测中的应用 李静1基金项目：*作者：李静（1982-），女，山东青岛人，博士研究生，主要研究方向：兽医信息学，E-mail：qingzhu2872@**通讯作者：王洪斌（1958-），男，内蒙古赤峰人，教授，博士生导师，主要研究方向：兽医信息学、小动物腹腔镜及动物麻醉等，Tel：0451E-mail：hbwang@，许丹宁2，肖建华1，马海鹍1， 基金项目： *作者：李静（1982-），女，山东青岛人，博士研究生，主要研究方向：兽医信息学，E-mail：qingzhu2872@ **通讯作者：王洪斌（1958-），男，内蒙古赤峰人，教授，博士生导师，主要研究方向：兽医信息学、小动物腹腔镜及动物麻醉等，Tel：0451E-mail：hbwang@ 东北农业大学动物医学学院 哈尔滨，150030；2. 仲恺农业工程学院生命科学学院 广州，510225） 摘要：目的 验证自回归综合移动平均（ARIMA）模型预测新城疫疫情的发病可行性，并为家禽养殖中新城疫的防治提供数据支持。方法 选择我国南方A、B两省作为研究对象，选取《兽医公报》中2000年1月~2007年12月8年间的月发病数以及A、B两省的家禽当年出栏量计算得出的新城疫月发病率作为时间序列模型的数据源。利用SPSS软件进行时间序列模型的构建，通过2008年的发病率进行模型检验。并据此对模型预测新城疫发病情况分析的可行性、建模步骤及准确性验证进行了探讨。结果 结果显示ARIMA模型计算出的预测值与实际值拟合较好，可用于对未来的新城疫发病率进行预测，为新城疫的防控工作提供可靠的参考依据。 关键词：新城疫；时间序列模型；月发病率；预测 新城疫（newcastle disease，ND）是由新城疫病毒引起禽的一种急性、热性、败血性和高度的接触性传染病。以高热、呼吸困难、下痢、神经紊乱、黏膜和浆膜出血为特征。具有很高的发病率和病死率，是危害养禽业的一种主要传染病。世界动物卫生组织（OIE）将其列为A类疫病，在我国将其列为一类传染病。对新城疫发病疫情的预测能够最大限度的降低该病对家禽养殖业的影响。疫病的发病率、死亡率等是随时间变化的，而且有着明显的时间先后顺序，这种一时间顺序排列起来的一系列观测值称为时间序列。时间序列分析是专门用于分析这种时间序列资料的统计模型，它考虑的不是变量间的因果关系，而是重点考察变量在时间方面的变化规律，并为之建立数学模型。 1 材料与方法 1.1数据来源 新城疫月发病数为2000.1-2008.12年间A、B两省新城疫数据（来自农业部出版的《兽医公报》），家禽出栏数来自《中国农业统计年鉴》。数据分别用A序列、B序列表示。 ND月发病率= ND月发病数 ×K 当年年底出栏量 注：（K为单位系数，本研究中月发病率的单位为1/万只） 1.2 时间序列自回归综合移动平均模型简介 1.2.1 ARIMA模型全称为自回归综合移动平均模型(Autoregressive Integrated Moving Average Model，简记ARIMA)，是由博克思(Box)和詹金斯(Jenkins)于70年代初提出的一著名时间序列预测方法，所以又称为box-jenkins模型、博克思-詹金斯法。其中ARIMA（p，d，q）称为差分自回归移动平均模型，AR是自回归，p为自回归项；MA为移动平均，q为移动平均项数，d为时间序列成为平稳时所做的差分次数。 1.2.2 将预测对象随时间推移而形成的数据序列视为一个随机序列，用一定的数学模型来近似描述这个序列。这个模型一旦被识别后就可以从时间序列的过去值及现在值来预测未来值。 1.2.3 ARIMA模型预测的基本程序[1 4步建模法：①序列的平稳化：首先从散点图中确定序列是否平稳，若不平稳则要对序列进行平稳化。②模型识别：利用自相关图等，把握模型的大致方向，为目标模型定阶，并选1-2个特定的模型拟合所分析的时间序列数据。③模型的参数估计与检验：参数估计是对识别阶段提供的粗模型参数估计并假设检验，用以判断模型是否恰当，并通过调试比较，确定最佳模型；如不佳则返回第一阶段，重新选定模型。④预测：利用2008年的新城疫发病数据检验2000-2007年间的新城疫月发病率为基础数据拟合的模型预测效果。 1.3 分析方法 采用SPSS（Statistical Product and Service Solutions）13.0统计分析软件进行模型的建立。运用其中的相关模块进行数据处理和分析，建立月发病率数据库，其中2000.1-2007.12年8年间的新城疫月发病率用于建立模型，2008年数据用于验证模型的预测效果。 2 建模与模型验证 2