Page 28 - 中国全科医学2022-16期
P. 28

http://www.chinagp.net   E-mail:zgqkyx@chinagp.net.cn  ·1945·


                               2
                                                     2
           12 阶统计量分别为 χ =2.403(P=0.879)和 χ =5.151              和 RMSE 值均小于 NNAR(1,1),即 ARIMA 模型更优,
           (P=0.953),提示模型拟合效果良好。在训练集上,                         见表 2。
           ARIMA(0,1,1)MAPE、MAE、RMSE 分别为 0.810%、               2.4 COPD DALYs 率预测模型构建 ARIMA(0,1,
           0.730、0.965;NNAR(1,1) 模 型 的 MAPE、MAE 和              2)(AIC=225.59,BIC=230.62) 为 AIC 与 BIC 最小的
           RMSE 分别为 1.033%、0.921、1.107。从模型的拟合和                 COPD DALYs 率预测模型。对模型进行 Ljung-Box 检验,
                                                                                            2
           预测情况来看,ARIMA(0,1,1)和 NNAR(1,1)                      延迟 6 阶、12 阶统计量分别为 χ =0.726(P=0.994)和
                                                                 2
           预测值的动态趋势与实际情况基本一致(图 3~4);                           χ =3.534(P=0.991),提示模型拟合效果良好。在训
           由 图 3 可 知,2017—2019 年 实 际 COPD 死 亡 率 均 在            练集上,ARIMA(0,1,2)MAPE、MAE、RMSE 值均
           ARIMA(0,1,1)预测值 80%CI 内。无论是在训练集                     低 于 NNAR(1,1)(0.622% 比 0.823%,11.305 比
           还是在测试集上,ARIMA(0,1,1)的 MAPE、MAE                      14.982,15.321 比 18.240)。从模型的拟合和预测情况

                                                                               5
                                         表 1 1990—2019 年中国 COPD 疾病负担变化情况(1/10 )
                                        Table 1 Changes in the burden of COPD in China,1990—2019
              年份                患病率                          死亡率                            DALYs 率
             (年)       男         女        总体        男         女         总体         男          女         总体
              1990    2 115.35  2 588.10  2 344.40  103.35   106.94     105.09    2 240.82  2 170.08   2 206.55
              1991    2 156.32   2 636.89   2 388.88   104.12   107.54   105.78   2 240.63   2 174.07   2 208.42
              1992    21 96.38   2 682.11   2 431.19   104.88   108.12   106.44   2 236.94   2 175.45   2 207.22
              1993    2 235.23   2 724.94   2 471.77   105.14   108.27   106.65   2 223.03   2 168.12   2 196.50
              1994    2 276.57   2 768.57   2 514.10   105.03   106.44   105.71   2 203.45   2 132.30   2 169.10
              1995    2 316.70  2 812.08  2 555.82  105.28   103.22     104.28    2 191.19  2 066.37   2 130.94
              1996    2 359.38   2 852.67   2 597.55   104.69   101.64   103.22   2 162.80   2 026.29   2 096.89
              1997    2 401.21   2 882.54   2 633.73   103.38   98.87   101.20    2 122.79   1 967.25   2 047.65
              1998    2 441.64   2 907.81   2 667.03   102.94   95.77   99.47     2 103.08   1 909.73   2 009.60
              1999    2 484.93   2 933.72   2 702.14   102.71   95.45   99.20     2 091.82   1 894.37   1 996.25
              2000    2 523.99  2 958.72  2 734.63  103.28   95.12      99.33     2 097.41  1 876.95   1 990.59
              2001    2 562.02   2 987.09   2 768.22   102.25   92.74   97.64     2 065.48   1 829.58   1 951.05
              2002    2 593.82   3 010.99   2 796.41   101.07   91.32   96.33     2 036.33   1 804.99   1 923.99
              2003    2 619.92   3 031.32   2 819.89   99.54   90.01    94.90     2 003.20   1 781.53   1 895.45
              2004    2 645.39   3 049.97   2 842.22   99.69   88.46    94.23     1 994.70   1 747.64   1 874.51
              2005    2 665.43  3 064.01  2 859.48  97.60    84.97      91.46     1 811.85  1 938.76   1 678.11
              2006    2 682.11   3 066.22   2 869.24   91.62   80.05    85.98     1 822.25   1 589.21   1 708.72
              2007    2 688.32   3 049.42   2 864.34   88.50   76.28    82.54     1 752.14   1 516.16   1 637.11
              2008    2 686.63   3 028.75   2 853.48   88.01   74.19    81.27     1 726.05   1 466.42   1 599.43
              2009    2 687.15   3 022.96   2 851.01   88.17   72.82    80.68     1 708.05   1 429.76   1 572.25
              2010    2 690.39  3 047.51  2 864.75  88.12    71.29      79.90     1 689.60  1 395.42   1 545.97
              2011    2 703.01   3 111.97   2 902.79   86.26   69.03    77.84     1 650.52   1 358.96   1 508.09
              2012    2 713.18   3 196.18   2 949.25   84.46   65.46    75.18     1 614.66   1 309.26   1 465.39
              2013    2 719.78   3 291.78   2 999.51   82.37   63.22    73.01     1 571.04   1 276.73   1 427.11
              2014    2 729.63   3 390.67   3 053.09   80.95   61.76    71.56     1 539.28   1 256.75   1 401.04
              2015    2 727.91  3 469.86  3 091.15  80.03    60.17      70.30     1 514.74  1 236.33   1 378.44
              2016    2 672.59   3 430.31   3 043.71   80.41   60.28    70.56     1 507.87   1 228.31   1 370.94
              2017    2 639.44   3 390.06   3 007.24   79.80   60.42    70.31     1 490.14   1 220.84   1 358.18
              2018    2 691.60   3 478.29   3 077.23   80.17   61.57    71.05     1 494.11   1 240.29   1 369.69
              2019    2 768.81  3 597.85  3 175.37  82.25    63.27      72.94     1 523.84  1 272.75   1 400.71
           变化率(%)     30.89     39.01    35.45     -20.42    -40.84     -30.59    -32.00     -41.35     -36.52
           AAPC(95%CI)  0.92    1.13      1.04     -0.83     -1.83      -1.29      -1.37     -1.86      -1.56
                                                                    a
                                                                              a
                                                                                                   a
                                                                                         a
                                               a
                                     a
                            a
                                                         a
             (%)    (0.85,0.98) (1.06,1.19) (0.97,1.11) (-1.41,0.51) (-2.17,-1.50) (-1.67,-0.97) (-1.63,-1.11) (-2.11,-1.62) (-1.86,-1.26) a
                                                         a
              注:DALYs= 伤残调整寿命年,AAPC= 平均年度变化百分比; 表示 P<0.001;COPD= 慢性阻塞性肺疾病
   23   24   25   26   27   28   29   30   31   32   33