摘要
目的了解辽宁省-年人间和家畜布鲁氏菌病（布病）疫情情况，建立时间序列模型，探索畜牧业因素对人间疫情影响，为布病预测防控提供依据。方法收集辽宁省-年布病月发病数，以发病数建立时间序列构建自回归滑动平均混合（ARIMA）模型，对年1-月进行短期预测，结合实际观测值评价预测效果；收集辽宁省同期羊存栏量和疫情，分析与人间发病存在的关联。结果构建模型ARIMA（1，1，1）（1，1，1）12，月发病数拟合值与实际观测值时间序列基本一致，预测效果良好，但精确度较低。羊年存栏量与人间布病发病数存在正关联性。结论ARIMA模型预测辽宁省布病发病情况基本可行，加强人畜布病的监测工作，以揭示人与家畜流行病学关系。前言
布鲁氏菌病（布病）是由布鲁氏菌属所引起的传染-变态反应性人畜共患病。人间布病主要通过接触病畜或被病畜污染的食物、产品等致使人感染发病。辽宁省自1年人间布病出现上升趋势，各县（区）均有新发病例报告，1年发病率为0.32/10万，年上升至.5/10万[1-2]。近年来时间序列模型被广泛应用到各类传染病预测研究[3-5]，对布病疫情进行科学预测以及合理评价防控干预措施效果方面可能会有积极作用。本研究通过探索自回归滑动平均混合（ARIMA）模型在辽宁省布病发病预测中的应用，并对畜牧业因素与布病疫情的关联性进行初步探索。
1
材料与方法
1.1疫情监测资料来源人间布病监测资料来自辽宁省疾病预防控制中心全国法定传染病监测信息报告管理系统，收集－年辽宁省人间布病月发病数，人口资料及羊年末存栏量来自于辽宁省统计年鉴公布数据，家畜布病疫情资料来自农业部出版的《疫情公报》。1.2ARIMA模型构建与预测ARIMA模型是Box-Jenkins方法中重要的时间序列分析预测模型，基本思想是将预测对象随时间的推移所形成的数据集视为一个随机序列，它是依赖于时间t的随机变量。试图解决的问题包括[-7]：（1）分析时间序列的随机性、平稳性和季节性；（2）在对时间序列进行分析的基础上选择合适的模型进行预测。本研究通过SPSS22.0软件构建模型，将年1月至年12月的月发病数据纳入模型构建，并对年前个月的月发病数进行预测，进而与真实监测值进行比较来判断预测效果。2
结果
2.1模型识别首先绘制月发病数的时间序列图，由图1可见，我国布病发病报告数呈明显上升趋势，存在周期性和季节性规律，高峰出现在每年5－9月之间，周期为12个月，不平稳序列，对原始数据进行转换从而获得平稳序列。经过探索，对数据进行对数转换后，再进行一次差分和一次季节差分，即d=1，D=1，布病的月发病数序列图趋于平稳，利用Eviews8.0软件对序列进行ADF检验，得到P0.05，即序列消除了明显的增长趋势后，围绕一个常数上下波动。模型形式可以确定为ARIMA（p，1，q）（P，1，Q）12。对数转换及d=1，D=1差分后的ACF图见图2。由图2可以发现，从1阶后出现迅速衰减，视为拖尾，且偏相关回归图也是自1阶后迅速衰减，因此初步确定Q（q）为0或1，P（p）为0或1。
2.2模型参数的估计和模型诊断根据阶数从小到大依次组合测试，对模型进行拟合优度检验（残差序列检验），Ljung-BoxQ值检验结果应为P0.05，提示模型的残差序列为白噪声，可以进行模型预测。经过筛选，得到4种白噪声序列的模型，通过比较R2，确定最优模型为ARIMA（1，1，1）（1，1，1）12，4种模型比较见表1。通过最优模型进行参数估计，见表2。绘制模型的残差序列ACF和PACF图对模型进行诊断，本模型残差序列ACF和PACF图显示残差基本都落在95%区间内，P0.05（图3）。因此本模型可以用来预测辽宁省人间布病月发病报告数。2.3预测应用
运用ARIMA（1，1，1）（1，1，1）12对辽宁省－年人间布病月发病报告数进行拟合预测（图4），结果显示观测值和拟合值都在预测值的95%CI内，但绝对误差和相对误差整体较大，预测效果良好，但精确度不高。年观测值与拟合值比较结果见表3，结果显示，观测值和拟合值的整体趋势基本一致，拟合值整体略高于观测值。2.4人间布病疫情与羊只年末存栏量、畜间疫情的关联性分析辽宁省－年人间布病疫情年报告由例持续增加到例，发病率由1.08/10万增加到.8/10万，且与全国的疫情形势基本一致（发病率由1.55/10万增加到4.34/10万）。辽宁省－年羊只年末存栏量由74万只增加到万只，见图5。与人间布病疫情关联性分析结果显示，秩相关系数为0.（P0.05），可以认为羊只年末存栏量与布病年发病报告数之间存在正相关性。辽宁省家畜间布病疫情在年报告发病48只，9年报告发病较多，当年报告布病家畜只，捕杀只，其后2、4和年无报告发病。3
讨论
本研究通过对－年辽宁省人间布病月发病数建立ARIMA模型，获得人间布病月发病数拟合值与实际观测值两条曲线基本吻合，由此对年1－月的月发病数进行预测，所得的拟合值略高于实际观测值，预测精度下降。其原因可能为：一是来自于模型本身的局限性，尽管ARIMA模型在很大程度上能够提取时间序列信息，但是从原理上决定了其依旧是一个线性模型，而传染病时间序列数据往往是非线性甚至是随机的，此外影响发病的因素包括宿主因素、环境因素、社会干预等多方面，而目前无法收集到这些因素的月度信息资料[8-9]；二是近年疾控部门与动物疫控部门执行《国家布鲁氏菌病防治计划（－年）》，加强联防联控，进行免疫家畜，起到了一定的防控效果。综上可能是对年的预测出现偏差的原因，但实际观测值并未超出拟合值95%CI范围，尚可以认为预测效果较好，对疫情预测仍具有一定意义。
全国人间布病发病近10年呈明显的上升趋势，并且不同地区的季节性差异比较明显[10]。在以往将ARIMA模型应用到人间布病研究的文献中，白永飞等[11]对山西省布病监测资料创建ARIMA（1，0，1）（1，1，0）12模型，娄鹏威等[12]对新疆布病监测资料创建ARIMA（0，1，2）（1，1，0）12模型，均取得了较好的预测效果，二者与本研究所采纳的模型不同，提示不同地区不同监测资料背景下需要根据不同的资料序列特征来进行趋势调整和季节调整，进一步建立适合于当前序列的模型来进行预测应用。此外，有文献显示具有周期性的数据在应用时间序列模型时需要足够的序列期数，同时短期预精度较高，不适用于长期预测[13]。
从流行病学分析，布病作为人畜共患病，人感染布病的传染源是家畜，与家畜疫情及相关畜牧业因素应存在必然关联。由于辽宁省人间布病多为羊种布鲁氏菌引起[14]，随着羊只存栏量的增多，人间疫情也会加重，从某个角度看符合流行病学的规律。但是，近年从农业部国家的家畜布病疫情信息系统，看到与人间的疫情明显差异，提示布病的控制需要加强与畜间疫情的联防联控的空间，疾控系统需要认真报告、核实人间病例的同时，加强漏报率管理。动物疫情控制机构也应加强家畜布病的监测，科学地反映家畜疫情信息，进一步强化检疫工作，及时处置淘汰病畜，以揭示辽宁省家畜与人间疫情真实流行病学关系，有效控制传染源，防止人、畜间疫情的传播扩散。
作者
王子江,孙英伟,毛玲玲,于维君,雷露,孙广玖,姚文清
第一作者简介
王子江，男，辽宁省沈阳市人，主管医师，硕士在读，主要从事传染病流行病学研究作者贡献:王子江：负责文章撰写与模型的建立孙英伟、毛玲玲：主要负责布病的实验室检测工作于维君：负责人间疫情数据的收集雷露：负责协助完成布病实验室检测孙广玖：负责布病疫情的调查姚文清：负责文章的修改工作扫码查看文章预览时标签不可点收录于话题#个上一篇下一篇