Page 68 - 南京医科大学自然科学版第1期

第41卷第1期
· 62 · 南 京 医 科 大 学 学 报 2021年1月

80.00 残差ACF 残差PACF

（ 1/10万 ） 60.00 24
23
22
21
发病率 40.00 20
19
18
20.00
17
16
15
0.00 14
2010-01 2010-06 2010-11 2011-04 2011-09 2012-02 2012-07 2012-12 2013-05 2013-10 2014-03 2014-08 2015-01 2015-06 2015-11 2016-04 2016-09 2017-02 2017-07 2017-11 2018-05 2018-11 2019-04 2019-06 滞后 13
12
11
日期（年-月） 10 9
图2 2010年1月—2019年6月南通市手足口病月发病率时 8
7
间序列图 6
5
Figure 2 Time series plot of the monthly morbidity of HF⁃ 4
MO in Nantong from January 2010 to June 2019 3 2
1
1.0
-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0
0.5 残差
ACF 0.0 图 5 模型 ARIMA（1，0，1）（1，1，1） 12残差序列的自相关与
偏自相关图
-0.5
Figure 5 Autocorrelation and partial autocorrelation
-1.0
function graphs of ARIMA（1，0，1）（1，1，1） 12
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 model residual sequence
延迟数目
系数 置信度上限 置信度下限
图3 2010年1月—2019年6月南通市手足口病月发病率季 Table 3 表3 ARIMA（1，0，1）（1，1，1） 12模型参数估计
节差分自相关（ACF）函数图 Parameter estimation of ARIMA（1，0，1）（1，1，1） 12
Figure 3 Seasonal difference autocorrelation function model
graph of the monthly morbidity of HFMO in 参数 估计 标准误 t值 P值
Nantong from January 2010 to June 2019 AR 滞后1 0.52 0.10 4.98 <0.001
MA 滞后1 -0.58 0.10 -5.90 <0.001
季节性差分 1.00
1.0
MA，季节性 滞后1 0.55 022 2.54 0.010
0.5
PACF 0.0 为 7.08/10 万 、1.81/10 万 、3.74/10 万 、7.21/10 万 、

10.71/10万和11.29/10万，绘制实际值和预测值序列
-0.5
图（图7）。
-1.0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 采用两个独立样本检验法，比较 2019 年 7—
延迟数目
系数 置信度上限 置信度下限 12 月手足口病实际发病率（10.96/10万、4.30/10万、
图4 2010年1月—2019年6月南通市手足口病月发病率季 5.15/10 万、6.25/10 万、5.83/10 万和 3.82/10 万）与预
节差分偏自相关（PACF）函数图 测值的差异，经分析两者间的差异无统计学意义
Figure 4 Seasonal difference partial autocorrelation func⁃ （Z=0.48，P=0.63），预测数据与真实值比较接近，因
tion graph of the monthly morbidity of HFMO 此该模型有较好的预测效果。
in Nantong from January 2010 to June 2019

3 讨 论
异均有统计学意义（P < 0.05，表3）。
2.4.3 预测模型的效果 本研究结果显示，2010—2019 年南通市手足口
［8］
通过 ARIMA（1，0，1）（1，1，1） 12模型对 2019 年 病年平均发病率为124.36/10万，低于同时期苏州 、
［6］
［9］
7—12 月手足口病发病率进行预测，预测结果分别 无锡 和全国 平均发病水平，整体呈现隔年高发、