摘要
目的 通过深入分析2002至2021年20年间职业性尘肺病死亡率的变化趋势，探究其有关联性的社会经济影响因素，旨在为制定有效的预防措施和公共健康政策提供科学依据，以降低职业性尘肺病的死亡率。 方法 通过Joinpoint模型分析中国2002—2021年间总人口、不同性别、城乡间居民尘肺病粗死亡率及标化死亡率，计算年度变化百分比（APC）和平均年度变化百分比（AAPC）。通过年龄—时期—队列模型，分析三种效应对城乡间不同人群职业性尘肺病死亡率变化的影响。通过对数线图综合描述城镇化率、能源消费总量、第二产业GDP、工伤保险参保人数、文盲人口占比，五种社会经济因素的变化趋势。通过绘制散点图并进行斯皮尔曼秩相关分析五种社会经济因素与职业性尘肺病死亡率的相关性，通过主成分分析中提取主成分，采用贝塔分析方法探究各因素与职业性尘肺死亡率的关联性。 结果 2002—2021年粗死亡率变化呈上升趋势（APC=AAPC=5.5%，95%CI：4.1~7.0， P＜0.001）；标准化死亡率在2004—2021年呈上升趋势（APC=5.6%，95%CI：3.6~7.7， P＜0.001），差异均有统计学意义。不同年份男性粗死亡率与标准化死亡率均高于女性。2002—2021年男性粗死亡率变化呈上升趋势（APC=AAPC=5.8%，95%CI：4.4~7.3，P＜0.001）；男性标准化死亡率在2004—2021年呈上升趋势（APC=6.0%，95%CI：4.1~8.0，P＜0.001）差异均有统计学意义。2002—2021年女性粗死亡率与标准化死亡率变化呈下降趋势（P＞0.05）差异无统计学意义。城市标准化死亡率整体呈下降趋势（APC=-0.2%，P＞0.05），其中，农村标准化死亡率2004—2021年呈上升趋势（APC=7.9%，P＜0.001）。年龄效应显示城乡尘肺病死亡率均随年龄增高而增大。时期效应随时间增加，且在城乡间有不同表现，城市死亡风险呈下降趋势，农村死亡风险呈上升趋势。队列效应显示，城市居民出生越晚死亡风险越低，农村居民出生队列呈波动趋势。对社会经济影响因素综合分析，发现城镇化率、能源消费总量、第二产业GDP、工伤保险参保人数呈上升趋势，文盲人口占比逐年减少。通过斯皮尔曼秩相关分析，除文盲人口占比呈负相关外外，均呈正相关关系（P＜0.05）差异有统计学意义。主成分分析提取两主成分，第一主成分可表示为经济发展，第二主成分与教育程度，方差解释率99.49%。通过贝塔回归方法构建尘肺病死亡率（1/100万）与主成分F1、F2的拟合模型，得出R2=0.81，拟合良好。Phi值较大，表明精度良好。对尘肺病死亡率（1/100万）影响最大的指标依次是：文盲人口占比、第二产业GDP、城镇化率、工伤保险参保人数、能源消费总量，系数分别为0.180、0.178、0.167、0.082、0.072均大于0为正向影响。 结论 2002—2021年职业性尘肺病总率呈上升趋势。男性死亡率高于女性，农村死亡率高于城市，男性、农村尘肺病死亡率每年增加。农村尘肺病死亡率随年龄、时期、队列效应的增加而产生的影响均高于城市。农村为尘肺病死亡率高发的重点地区、男性为尘肺病死亡的高风险人群。文盲率对职业性尘肺病死亡率影响最大，教育程度的提高可以降低尘肺病的死亡率。其次是第二产业GDP与城镇化率对尘肺病死亡率呈正向影响。能源消费总量、工伤保险参保人数，对尘肺病死亡率呈正向影响但影响相对较弱。
NETL