Page 34 - 2023-04-中国全科医学
P. 34

2023年2月   第26卷   第4期                                 http: //www.chinagp.net   E-mail: zgqkyx@chinagp.net.cn  ·419·

           炎发病率极低，本研究仅对伤寒与副伤寒、病毒性肝炎、                           标（LISA）聚集性地图（受限于篇幅，仅绘制 2008、
           痢疾的发病例数和发病率数据进行分析。于2022年3月，                         2013、2018 年 3 个时间节点的 LISA 聚集性地图），空
           从由中国疾病预防控制中心负责建设的对外数据服务门                            间关联模式可分为 4 种类型：“高 - 高聚集区”“低 -
           户，即公共卫生科学数据中心              ［11］ 提取 2008 年 1月至       低聚集区”“高 - 低聚集区”“低 - 高聚集区”                 ［16］ 。
           2018 年 12 月中国 31 个省份（不含香港特别行政区、                     1.2.3 时空扫描统计分析 时空扫描统计分析用于识别
           澳门特别行政区和台湾地区）甲乙类肠道传染病〔伤寒                            疾病的高发地区，分析疾病的时间聚集性                  ［17］ 。以年度
           与副伤寒、病毒性肝炎（甲型、戊型和未分型肝炎）、                            为周期对 2008—2018 年中国各省份甲乙类肠道传染病
           痢疾〕的月发病例数和月发病率数据。同期人口数据来                            发病数据进行时空扫描分析。基于离散型 Poisson 分布
           源于《中国统计年鉴》（2009—2019 年）             ［12］ 。全国行       模型  ［18］ ，构建动态的圆柱体扫描窗口，其中圆柱体的
           政区矢量地图及经纬度坐标通过国家基础地理信息中心                            底面表示地理区域，圆柱体的高度表示对应的时间区间。
           获取  ［13］ 。                                          对每一个扫描窗口，根据实际发病例数和人口数计算出
           1.2 研究方法                                            预期发病例数后，利用扫描窗口内、外的实际发病例
           1.2.1 描述性分析 以发病例数、发病率等为主要指标                         数和预期发病例数计算对数似然比（LLR）和相对危险
           对甲乙类肠道传染病的流行病学特征进行描述性分析，                            度（RR），判断不同时间和地理区域疾病发病的聚集
                        2
           采用趋势性 χ 检验分析 2008—2018 年中国甲乙类肠                      程度  ［17］ 。
                                                                            C
                                                                          {
                                                                                      ］
           道传染病发病率的变化趋势，并呈现 2008—2018 年中                             LLR=log （—） c ［ （C-c） （Ｃ -c） }
                                                                                 （C-n）
                                                                            n
           国甲乙类肠道传染病发病率的季节性和周期性变化。绘                                其中 C 为总发病例数，c 为扫描窗口内的实际发病
           制中国甲乙类肠道传染病发病率地图（受限于篇幅，仅                            例数，n 为扫描窗口内的预期发病例数。RR 是指扫描
           绘制 2008、2013、2018 年 3 个时间节点的发病率地图），                 窗口内和扫描窗口外的发病风险（实际发病例数 / 预期
           发病率分级采用 ArcGIS 自然间断点分级法。检验水准                        发病例数）之比。基于数据统计特征，结合既往研究经
           α=0.05。                                             验进行相关参数的设置，将 Kulldorff 圆形扫描统计量中
           1.2.2 空间自相关分析 按年度对 2008—2018 年中国                    最大空间扫描窗口设定为风险人口数的 25%，采用蒙特
           各省份甲乙类肠道传染病发病率进行空间自相关分析。                            卡罗模拟的方法检验 LLR 的 P 值，蒙特卡罗模拟次数
           空间自相关分析主要用于研究一个地区空间单元上的                             为 999 次，选取 LLR 最大且有统计学意义（P<0.01）的
           某一属性值是否与邻近地区空间单元上同一属性值存                             扫描窗口作为Ⅰ类聚集区域。
           在相关性，是度量空间单元属性值聚集离散程度的方                             1.3 统计学方法 采用 Excel 2019 软件和 SPSS 26.0 软
           法 ［14］ ，包括全局空间自相关分析和局部空间自相关分                        件进行描述性统计分析。使用 ArcGIS 10.8 软件绘制发
           析。其中全局空间自相关主要反映属性值在整个区域的                            病率地图和 LISA 聚集性地图。采用 SaTScan v10.0.2 软
           空间分布特征。通过计算全局 Moran's I 值（取值范围为                     件进行时空扫描统计分析。
           -1~1），并对全局 Moran's I 值进行 Z 检验，判断属性                  2 结果
           值的空间分布属于分散分布、随机分布、集聚分布中的                            2.1 2008—2018 年中国甲乙类肠道传染病流行特征
           何种类型。P>0.05 时，属性值的空间分布呈随机分布；                        2.1.1 2008—2018 年中国甲乙类肠道传染病发病情况
           P<0.05 时，属性值的空间分布呈非随机分布。P<0.05                      及变化趋势 2008—2018 年全国累计报告甲乙类肠道
           且全局 Moran's I 值 >0，表示属性值呈现空间正相关，                    传染病发病 3 220 480 例，其中伤寒与副伤寒 142 490
           全局 Moran's I 值越大，属性值在空间分布上呈现出的                      例（4.43%），病毒性肝炎 991 336 例（30.78%），痢
           集聚效应越明显；P<0.05 且全局 Moran's I 值 <0，表示                疾 2 086 654 例（64.79%）；甲乙类肠道传染病年均发
                                                                            5
           属性值呈现空间负相关，全局 Moran's I 值越小，属性                      病率为 7.25/10 ，发病率整体呈下降趋势（P<0.05）；
           值在空间分布上呈现出的分散效应越明显                  ［15］ 。          伤寒与副伤寒、病毒性肝炎、痢疾发病率均呈下降趋势
                         n  n                                  （P<0.05），伤寒与副伤寒发病率由 2008 年的 1.18/10            5
                         ∑ ∑W ij （x i - ）（x j - ）
                                                                                   5
             Moran's I=  i=1  j=1 n  n  n                      降至 2018 年的 0.78/10 ，病毒性肝炎发病率由 2008 年
                           ∑ ∑W ij ∑ （x i - ） 2                        5                    5
                            i  j  i                            的 9.81/10 降至 2018 年的 4.34/10 ，痢疾发病率由 2008
                                                                                               5
                                                                          5
               其中 n 表示所要研究的对象中空间的区域数，x i 为                     年的 23.65/10 降至 2018 年的 6.56/10 ，见表 1。
           第 i 个区域的属性值，x j 为第 j 个区域的属性值，              表示       2.1.2 2008—2018 年中国甲乙类肠道传染病发病率的
           所有区域属性值的平均值，W ij 表示空间权重矩阵 W 中                       季节性和周期性变化 2008—2018 年全国每月均有甲
           的第（i，j）个元素。局部空间自相关可反映属性值在                           乙类肠道传染病病例出现，甲乙类肠道传染病流行强度
           局部区域的空间分布特征            ［14］ 。绘制空间关联局部指              呈逐渐减弱趋势。伤寒与副伤寒、痢疾发病具有较强的
   29   30   31   32   33   34   35   36   37   38   39