Page 73 - 2023-05-中国全科医学

·586· http: //www.chinagp.net E-mail: zgqkyx@chinagp.net.cn February 2023, Vol.26 No.5

2016 年全国肺结核发病率空间分布 2017 年全国肺结核发病率空间分布 2018 年全国肺结核发病率空间分布

rate rate rate
0.583~1.583 1.742~2.793 1.782~2.585
1.584~2.583 2.794~3.583 2.586~3.275
2.584~4.500 3.584~4.192 3.276~4.257
4.501~6.750 4.193~5.260
6.751~9.083 5.261~7.322 4.258~5.436
5.437~7.088
9.084~13.083 7.323~12.898 7.089~13.888
13.084~26.000 12.899~16.882 13.889~25.412
1∶48000000 1∶48000000 1∶48000000
图 1 全国肺结核发病率空间分布图
Figure 1 Spatial distribution of pulmonary tuberculosis incidence in China during 2016—2018
2.2 模型对比结果 经过多重共线性和空间自相关检 表 1 OLS、GWR、GTWR 模型比较的结果
2
2
验（全局 Moran's I= 0.376），剔除气温这一具有强共线 Table 1 Values of AICc，R and adjusted R of OLS，GWR and GTWR
models
性的变量（VIF Tmax =48.01，VIF Tmin =48.34）后，分别建立
模型 AICc R 2 Adjusted R 2
OLS、GWR 和 GTWR 模型，评估并比较模型优度，结 OLS 5 871.92 0.226 0.220
2
2
果如表 1 所示。GTWR 模型的 R 值与 Adjusted R 值均 GWR 4 746.57 0.573 0.572
比 OLS 和 GWR 模型要高，同时 GTWR 模型的 AICc 值 GTWR 3 806.78 0.902 0.901
均比 OLS 和 GWR 模型要小，表明 GTWR 模型能更好 注：AICc=赤池信息量，OLS=普通最小二乘，GWR=地理加权回归，
地解释自变量对肺结核发病情况的影响，更能解释具有 GTWR= 时空地理加权回归
时空特征的数据。
表 2 GTWR 模型的拟合系数
2.3 拟合系数的时空特性 使用均数、最小值、最大 Table 2 Fitting coefficients of meteorological and air quality factors
值、四分位数间距来描述 GTWR 模型的拟合系数，结 included in GTWR model
果如表 2 所示。基于 GTWR 模型的拟合系数结果绘制 变量 均数 最小值 最大值 第一四分位数 中位数 第三四分位数
每个变量的核密度图，结果如图 2 所示。拟合系数的核 风速 -0.454 25 -4.006 1.823 -0.874 00 -0.458 50 -0.074 50
密度图结果显示：自变量风速呈现多峰分布，主峰约 湿度 0.013 85 -0.316 0.100 0.000 01 0.016 00 0.041 75
为 -0.5；自变量湿度呈现多峰分布，主峰约为 0.01， PM 2.5 0.014 04 -0.054 0.378 -0.010 00 -0.002 00 0.015 75
且各峰系数均大于 0；自变量 PM 2.5 呈现左偏峰分布， PM 10 0.002 42 -0.081 0.060 -0.005 00 0.006 00 0.014 00
主峰约为 -0.01；自变量 PM 10 呈现单峰分布，主峰约为 SO 2 0.028 45 -0.393 0.634 -0.016 00 0.004 00 0.050 00
0.005；自变量 SO 2 呈现多峰分布，主峰约为 0；自变量 NO 2 -0.062 42 -0.345 0.102 -0.106 00 -0.041 00 -0.014 00
NO 2 呈现多峰分布，主峰约为 -0.03；自变量 CO 呈现 CO 0.101 42 -25.231 7.941 -0.617 00 0.174 00 1.356 75
右偏峰分布，主峰约为 -0.5；自变量 O 3 呈现多峰分布， O3 -0.002 43 -0.113 0.053 -0.007 00 -0.001 00 0.001 00
主峰约为 0。
基于 GTWR 模型的拟合系数，分别绘制各个变量 部和西部的拟合系数较低、东南部和中东部的拟合系数
的时空分布图（图 3）。拟合系数的时空分布图结果显示： 较高的分布格局，其中 PM 10 对浙江、上海、福建等地
自变量风速呈现中部和东南部的拟合系数较低、东北和 区的肺结核发病率具有显著的促进作用；自变量 SO 2 呈
西部的拟合系数较高的分布格局，其中风速对青海、甘 现中部和北部的拟合系数较低、西部和南部的拟合系数
肃、湖南、江西等地区的肺结核发病率具有显著的抑制 较高的分布格局，其中 SO 2 对西藏、广东、新疆等地区
作用，对西藏、新疆、黑龙江等地区的肺结核发病率具 的肺结核发病率具有显著的促进作用；自变量 NO 2 呈现
有显著的促进作用；自变量湿度呈现西部和北部的拟合 西北部和西南部的拟合系数较低，东北部、中部的拟合
系数较低、东南和南部的拟合系数较高的分布格局，其 系数较高的分布格局，其中 NO 2 对黑龙江、吉林等地区
中湿度对海南、广东、广西等地区的肺结核发病率具有 的肺结核发病率具有显著的促进作用；自变量 CO 呈现
显著的促进作用；自变量 PM 2.5 呈现东北部和中东部的 西部和西南部的拟合系数较低、东南部和中西部的拟合
拟合系数较低、西南部和西部的拟合系数较高的分布格 系数较高的分布格局，其中 CO 对福建、青海、甘肃等
局，其中 PM 2.5 对广西、云南、新疆、西藏等地区的肺 地区的肺结核发病率具有显著的促进作用；自变量 O 3
结核发病率具有显著的促进作用；自变量 PM 10 呈现中 呈现西北部和南部的拟合系数较低、东北部和西南部的