根据2020年第7次人口普查数据,我国≥60岁老年人口已达2.64亿,占总人口18.7%
[1]。随着人口老龄化加剧,老年人肺功能水平随年龄增长自然下降,不仅与各类呼吸系统疾病密切相关,还会增加心血管疾病、阿尔兹海默症及脑卒中等多种疾病的患病风险,甚至与全因死亡风险呈负相关,肺功能水平越低,全因死亡率越高,严重威胁老年居民健康
[2-5]。
肺功能检查(PFT)是评估肺功能水平的重要手段,已被证实与改善患者预后、降低全因死亡率相关
[6]。中华医学会明确建议,将PFT纳入≥60岁人群的定期体检项目,对肺功能减退的老年人群,每3~6个月进行1次监测筛查
[7]。但临床实践中,PFT对受试者配合度要求较高,部分老年患者常因自身疾病、认知障碍等原因而导致数据缺失,给肺功能评估带来了挑战
[6, 8]。
机器学习作为人工智能核心分支,凭借参数优化与特征工程优势,已在医学领域展现出良好的应用前景
[9]。现有研究虽在特定呼吸系统疾病的诊断预测方面取得一定进展,但针对老年患者的肺功能评估研究仍较为匮乏
[10, 11]。基于此,本研究以用力肺活量(FVC)和1s用力呼气量(FEV
1)为目标,整合易获取的临床指标,利用机器学习算法构建接受度高、普适性广的肺功能评估模型,为老年患者肺功能异常的早期识别、及时干预及个体化诊疗决策提供参考。
1 资料和方法
1.1 研究对象
选取2020年9月1日~2024年12月31日在广东省人民医院进行PFT的7720例老年患者作为研究对象。根据中华医学会呼吸病学分会发布的《肺功能检查指南(第二部分)—肺量计检查》,肺量计对FVC和FEV
1的标准测量量程为0.5~8 L
[12]。纳入标准:年龄≥60岁;行PFT前未使用支气管扩张剂。排除标准:FVC或FEV
1检测数据缺失者;FVC或FEV
1超出仪器标准测量范围者。本研究已通过广东省人民医院伦理委员会审批(伦理批号:KY2024-1206-01)。
1.2 研究方法
1.2.1 肺功能检查
由经验丰富的技术人员严格遵循美国胸科学会和欧洲呼吸学会的相关标准规范操作
[13],采用Cosmed K4b2肺功能仪进行检测。呼吸高危人群的定义标准为:使用支气管扩张剂前FVC占正常预计值的百分比(FVC%)或FEV
1占正常预计值的百分比(FEV
1%)低于80%
[14]。
1.2.2 数据收集及预处理
从电子病历中收集研究对象的数据:基本信息:年龄、性别、婚姻状况、体质量指数(BMI)、吸烟史、饮酒史;疾病相关因素:住院次数、糖尿病、脑卒中、高血压、冠状动脉粥样硬化(CAD)、肺动脉高压以及呼吸系统疾病(其中包括肺部感染、肺炎、支气管炎、哮喘、慢性阻塞性肺病、支气管扩张、间质性肺炎、肺结核、肺气肿、肺不张和肺癌);实验室检查:中性粒细胞比例(NEUT%)、乳酸脱氢酶(LDH)、凝血酶原时间(PT)、凝血酶时间(TT)、嗜酸性粒细胞百分比(EOS%)、大血小板比率(P-LCR)、氯(Cl)、活化部分凝血活酶时间(APTT)、红细胞分布宽度变异系数(RDW-CV)、红细胞分布宽度标准差(RDW-SD)、红细胞压积(HCT)、红细胞计数(RBC)、葡萄糖(Glu)、血小板计数(PLT)、血浆纤维蛋白原(FIB)、血红蛋白(HGB)、钙(Ca)、钠(Na)、钾(K)。对缺失比例≥30%的临床变量进行删除,缺失比例<30%的变量使用多重插补法进行填补(
表1)。
1.2.3 变量筛选与模型构建
首先通过单因素分析初步筛选预测变量,然后采用LASSO回归进一步筛选变量,通过十折交叉验证确定最优惩罚参数,以减少多重共线性。按照8∶2的比例将数据集随机划分为训练集(6176例)和测试集(1544例)。使用线性回归(LR)、随机森林(RF)、极端梯度提升(XGBoost)、梯度提升机(GBM)和支持向量回归(SVR)构建筛查模型。采用Pearson相关系数(r)、组内相关系数(ICC)、均方误差(MSE)、均方根误差(RMSE)和平均绝对误差(MAE)评估模型的预测性能并筛选出最优模型。同时比较完整数据集与多重插补后数据集的模型性能差异,以评估缺失数据处理方法对模型预测结果的影响。使用沙普利加性解释(SHAP)方法分析预测变量对最优模型的影响。此外,通过确定最佳截断值,比较模型预测值与真实测量值的受试者工作特征曲线(ROC)下面积(AUC)、敏感度和特异度,评估模型对呼吸高危人群的分类效能,并与临床标准截断值进行对比。
1.2.4 统计学分析
根据FVC和FEV1的四分位数,将研究对象分为4个水平组进行基线特征分析。其中,连续变量以均数±标准差描述,组间比较采用t检验或方差分析,分类变量使用频数(百分比)描述,使用卡方检验或Fisher确切概率法比较组间差异。P<0.05为差异有统计学意义,所有统计分析均使用R4.5.1进行。
2 结果
2.1 研究对象的基线特征及单因素分析
本研究共纳入7720例研究对象,其中,男性4197例(54.4%),女性3523例(45.6%);年龄60~96(68.56±6.27)岁。肺功能指标方面,FVC为0.5~5.19 L,2.45±0.69 L;FEV
1为0.5~3.93 L,1.99±0.54 L。单因素分析结果显示,FVC中共计28个变量的差异有统计学意义,FEV
1中27个变量的差异有统计学意义(表
2、
3)。
2.2 预测变量的LASSO回归筛选
LASSO回归结果显示,FVC的
和
分别为0.0013和0.0191,分别筛选出27个和14个变量;FEV
1的
和
分别为0.0004和0.0099,分别筛选出29个和16个变量。为平衡模型拟合度与复杂度,避免过拟合,本研究基于
筛选变量,最终FVC的预测变量为年龄、BMI、性别、呼吸系统疾病、肺动脉高压、脑卒中、PT、Na、Cl、APTT、RDW-CV、FIB、HGB和K,FEV
1的预测变量为年龄、性别、吸烟史、呼吸系统疾病、肺动脉高压、脑卒中、PT、EOS%、P-LCR、Cl、APTT、RDW-CV、RDW-SD、PLT、FIB和HGB(
图1)。
2.3 5种机器学习模型比较
在插补后数据集中构建5种机器学习模型,结果显示XGBoost对FVC和FEV
1的预测性能最佳,r值分别为0.733和0.711,ICC值为0.694和0.668,MSE值为0.210L
2和0.136L
2,RMSE值为0.459L和0.369L,MAE值为0.364L和0.291L(
表4、
图2)。
2.4 完整数据集与插补后数据集模型性能对比
结果显示,完整数据集下所有模型的r和ICC均略有降低,MSE、RMSE、MAE等误差指标较插补后数据集有所升高,提示多重插补对模型预测性能有所改善。各模型的预测性能差异幅度较小,且XGBoost在两种数据集中均保持了接近最优的表现,表明本研究基于插补后数据集构建的模型结论具有一定稳健性(
表5)。
2.5 XGBoost的SHAP解释
采用SHAP方法对最优XGBoost模型进行可解释性分析。变量重要性排序显示(
图3A、B),性别、年龄是影响FVC与FEV
1预测的共同核心因素。在FVC预测中,贡献度前10位的变量依次为性别、年龄、PT、HGB、RDW-CV、FIB、BMI、Cl、K、Na;在FEV
1预测中,前10位变量则为性别、年龄、HGB、FIB、PT、吸烟史、RDW-CV、Cl、APTT、呼吸系统疾病,提示FVC与FEV
1的影响因素既有共性,也存在特异性差异。蜂群图直观展示了各变量的影响方向与强度(
图3C、D),其中男性均表现为显著正向影响,女性则表现为显著负向影响;年龄呈明确的负向影响,随着年龄升高,对FVC与FEV
1预测值的负向影响持续增强。单样本SHAP力图进一步展示了变量对个体预测值的作用方向与幅度(
图4)。
2.6 XGBoost模型对呼吸高危人群的分类性能
基于XGBoost模型分析得出,FVC%与FEV
1%的最佳截断值分别为82.93%、85.46%,对应AUC为0.700、0.706。进一步以临床标准截断值进行对照评估,结果显示,FVC%和FEV
1%的AUC均为0.688。整体对比可见,相较于临床截断值,模型最优截断值下的区分效能更佳,整体分类表现更稳定(
表6、
图5)。
3 讨论
PFT是临床评估肺功能的重要手段,但其检测对受试者配合度要求较高,老年患者常因认知、体力、疾病等因素难以完成规范检查,易导致数据缺失,进而阻碍肺功能评估
[6, 8]。为此,本研究基于临床常规易获取的指标,使用机器学习算法构建接受度高、普适性好的肺功能评估模型,为临床早期干预与诊疗决策提供参考依据。
在5种机器学习算法中,XGBoost模型表现出最优预测效能。XGBoost模型基于决策树的集成学习框架,具有良好的泛化能力与拟合性能,已在疾病诊断、分型及预后预测等复杂临床任务中展现出优异的应用价值
[15-17]。为提升模型的可解释性,本研究采用SHAP算法进行特征分析。其中,男性对FVC及FEV
1均呈正向影响,这可能与男性在肺部体积、气道管径及胸廓结构上所具备的天然解剖生理优势密切相关
[18]。年龄对肺功能的影响呈现负向趋势,高龄个体的肺功能预测值更低,其机制可能与肺功能的生理性衰退有关,包括肺组织弹性下降、肺泡表面积减少及呼吸肌力减弱等,这些变化共同削弱肺部的防御与代偿能力
[19, 20];也可能与环境暴露的累积效应有关,随着年龄增长,个体接触烟草、职业粉尘及空气污染物等有害因素的总量增加,导致气道与肺实质的慢性进行性损伤
[20, 21]。同时,凝血指标(PT、APTT、FIB)也有重要的预测价值。从病理机制来看,血液高凝状态由肺微小血管炎性损伤、血管内皮受损、血小板活化及慢性缺氧等多种因素诱发,而这种高凝状态会进一步影响肺气体交换,加重肺功能损伤
[22, 23]。多项研究证实,哮喘患者的凝血功能指标与肺功能指标存在显著相关性
[24, 25]。此外,FEV
1水平与呼吸系统疾病密切相关。现有研究已表明,哮喘、慢性阻塞性肺疾病、肺炎等呼吸系统疾病均可显著影响FEV
1水平
[26, 27]。上述结果反映了不同肺功能指标的病理生理特征,可为临床开展个体化肺功能评估提供可解释性依据。
本研究进一步评估了XGBoost模型对呼吸高危人群的分类效能。结果显示模型具备一定的高危人群识别能力,但整体效能处于中等水平,尚未达到高精度筛查与鉴别要求。此外,模型最佳截断值(FVC%:82.93%,FEV1%:85.46%)均高于临床80%的标准阈值,存在一定的阈值偏移。核心原因可能在于呼吸高危人群肺功能状态受多因素共同调控,年龄、性别、基础疾病、生活习惯及环境等变量间存在复杂交互关系,增加了模型精准区分的难度。同时,肺功能评估存在显著个体异质性,不同个体生理基础与疾病进展差异较大,这也是临床肺功能预测模型普遍面临的挑战。
本研究构建的模型具有一定优势:以FVC和FEV1为预测目标,利用常规易获取的临床资料构建模型,无需患者积极配合即可完成肺功能预测,贴合临床实际需求;样本量充足,为模型构建提供了坚实的数据基础,有助于获得稳健的统计推断结果。
然而,本研究也存在一定的局限性:研究采用横断面及回顾性设计,虽总体样本量较大,但仅在单中心数据集进行验证,后续需开展多中心、多人群外部验证,以全面评估模型的泛化能力;模型对FVC、FEV1的预测及高危人群分类效能均处于中等水平,未来可通过扩大样本量、增加预测变量及融合影像特征等方式,进一步提升预测准确性。
综上,本研究构建了一个肺功能评估模型,为老年患者肺功能异常的早期识别、及时干预及个体化诊疗决策提供参考。本研究结果也为后续模型优化提供了方向,为进一步完善肺功能评估体系、探索更精准的高危人群识别策略奠定了基础。
广东省科技计划项目重点领域研发计划(2025B0101120008)