奉化区缺血性脑卒中发病季节的SARIMA模型预测

目的利用季节性自回归滑动平均模型法(SARIMA)对奉化区缺血性脑卒中发病情况进行预测,为做好脑卒中防控工作提供依据。方法采用R 3.3.2软件对奉化区2009—2014年缺血性脑卒中的发病数据建立SARIMA模型,对原始序列做一阶差分与一阶季节差分的预处理,预处理后的序列进行KPSS平稳性检验。利用auto.arima()函数结合模型参数估计来选择最优模型,再利用sarima()函数进行参数估计与残差诊断,最后应用建立的模型对2015年奉化区缺血性脑卒中发病进行预测验证。结果最优模型为SARIMA(0,1,1)(0,1,1)12,模型的参数估计值ma1=-0.824,sma1=-0.999,2个参数经检验差异均有统计学意义(P0.05)。sarima()函数自带的残差诊断图显示模型残差为白噪声序列。2015年各月预测结果误差率最高为-25.00%,最低为0.00%,平均误差率为10.75%,2015年每月实际发病数均在预测值的95%可信区间内。结论 SARIMA(0,1,1)(0,1,1)12模型拟合效果较好,可用于奉化区缺血性脑卒中发病情况的短期预测。