识别慢性术后疼痛基线预测因素的两部分统计模型

Two-part Statistical Model for Identifying Baseline Predictors of Chronic Postsurgical Pain

作者信息Stephan G Frangakis, Xuran Meng, Mark C Bicket, Vidhya Gunaseelan, Sawsan As-Sanie, Andrew Urquhart, Yi Li, Chad M Brummett
PMID41941700
发布时间2026-07-01
DOI10.1097/ALN.0000000000006080

实验完整度

基于前瞻性单中心队列(n=3,925)的观察性分析,采用重复随机拆分训练/测试集(400次)评估模型性能,无体外或动物功能验证。

主要模型

前瞻性单中心队列(n = 3,925) 训练集(n = 3,000)与测试集(n = 925)

重点核对

数据随机拆分为训练集(n = 3,000)和测试集(n = 925),重复400次 两部分模型的第一部分使用logistic回归估计疼痛发生概率(NRS > 0 vs. NRS = 0) 两部分模型的第二部分使用gamma分布log链接的广义线性模型估计疼痛严重程度(NRS > 0) 主要结局为术后3个月手术部位平均疼痛NRS评分(0-10) 模型性能指标:R2、RMSE、MAE,使用Wilcoxon匹配配对符号秩检验比较

摘要

背景:相当比例的患者报告术后无疼痛,导致零值过多,给传统统计模型的分析带来挑战。本研究旨在检验以下假设:常用于医疗支出研究的两部分模型在预测术后疼痛方面优于传统模型,并能更好地识别这一临床重要结局的预测因素。方法:本研究分析了一个前瞻性收集的单中心慢性术后疼痛数据集(n = 3,925),将新型两部分建模框架与logistic回归和线性回归进行比较。两部分模型首先估计术后3个月发生疼痛的概率(二分类结局),然后估计疼痛患者的疼痛严重程度(1至10数字评定量表)。为获得无偏的模型性能评估,数据被随机分为训练集(n = 3,000)和测试集(n = 925)。模型在训练集上训练,在测试集上评估。该过程重复400次以计算平均性能估计。作为次要目标,研究评估了15个基线因素(包括经过验证的患者报告疼痛、功能和心理测量、合并症及手术细节)在不同建模方法中的关联。结果:与单独线性回归相比,两部分模型表现出更优的预测性能,平均R2值更高(0.075 vs. 0.050;P < 0.0001),均方根误差更低(1.466 vs. 1.485;P < 0.0001),平均绝对误差更低(1.020 vs. 1.030;P < 0.0001)。由于共享二分类成分,无法与单独logistic回归进行统计比较。两部分模型识别出七个基线术前协变量与术后疼痛独立相关,而这些被logistic或线性模型或两者遗漏:自报种族、教育水平、美国麻醉医师学会分级、整体身体疼痛、广泛性疼痛、症状严重程度和焦虑水平。在三种模型中均识别出的显著基线因素为患者性别、手术类型和手术部位疼痛。结论:两部分模型可能提供优于线性回归和logistic回归的统计方法,更好地识别慢性术后疼痛的患者和临床护理相关风险因素,主要通过区分驱动慢性疼痛发生的因素与与疼痛严重程度相关的因素。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
与传统的多变量线性和logistic回归模型相比,两部分模型能否更好地建模术后疼痛,并识别出独特的患者和/或临床护理相关独立关联因素?
核心机制
两部分模型通过分离两个过程来建模术后疼痛:第一部分使用logistic回归估计发生任何疼痛的概率,第二部分使用gamma分布log链接的广义线性模型估计疼痛患者的疼痛严重程度,从而区分影响疼痛发生与疼痛严重程度的因素。
主要证据
在n = 3,925的前瞻性队列中,通过400次随机训练/测试拆分,两部分模型在R2(0.075 vs. 0.050)、RMSE(1.466 vs. 1.485)和MAE(1.020 vs. 1.030)上均优于线性回归(P < 0.0001),并识别出7个被传统模型遗漏的显著基线协变量。
研究意义
两部分模型可能提供优于线性和logistic回归的统计方法,更好地识别慢性术后疼痛的患者和临床护理相关风险因素,为未来慢性术后疼痛分析提供更准确、可解释的协变量估计,并清晰分离影响疼痛发生与严重程度的因素。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

队列构建与数据收集

建立用于分析慢性术后疼痛预测因素的前瞻性队列

从Michigan Genomics Initiative和Analgesic Outcome Study两个单中心前瞻性观察队列中,纳入18岁及以上、2010年3月1日至2024年7月31日期间接受择期手术的患者,收集术前基线协变量和术后3个月疼痛结局。

2

基线协变量与疼痛表型评估

获取用于预测术后疼痛的候选基线因素

通过术前调查和医疗记录提取人口学数据、BMI、ASA分级、手术类型、当前阿片类药物使用;使用经过验证的问卷评估疼痛严重程度、广泛性疼痛、焦虑和抑郁症状。

3

两部分模型构建

建立能处理零膨胀半连续疼痛结局的统计模型

第一部分使用logistic回归估计发生任何术后疼痛的概率(NRS > 0 vs. NRS = 0),训练于全部数据集(n = 3,925);第二部分使用gamma分布log链接的广义线性模型估计疼痛患者的疼痛严重程度(NRS > 0),仅训练于疼痛患者子集。两部分独立训练。

4

模型验证与性能评估

评估两部分模型相对于传统线性回归的预测性能和泛化能力

将患者随机分为训练集(n = 3,000)和测试集(n = 925),重复400次;每次迭代在训练集上训练模型,在测试集上评估R2、RMSE和MAE;使用Wilcoxon匹配配对符号秩检验比较模型性能。

5

协变量显著性比较与识别

比较两部分模型与logistic、线性模型识别的显著基线协变量

在每个训练队列中,识别在logistic、线性和两部分模型中显著的协变量(P ≤ 0.05),统计每个协变量在400次独立随机训练队列中显著的次数;比较各模型在50%、75%和90%迭代中识别的协变量。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
R软件(版本4.2.2)----
PROsetta Stone转换表----

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
队列构建与纳入排除
纳入标准:18岁及以上,择期手术,2010年3月1日至2024年7月31日;排除标准:无法提供书面知情同意、非英语使用者、囚犯;完整病例分析,n = 3,925
阅读提示:Methods中Study Design and Recruitment小节及Figure 1 STROBE图
基线协变量与疼痛表型评估
15个基线因素包括人口学、BMI、ASA分级、手术类型、当前阿片类药物使用、Brief Pain Inventory平均疼痛、Michigan Body Map、症状严重程度指数、HADS或PROMIS 4a焦虑抑郁评分;术前和术后3个月评估
阅读提示:Methods中Patient Phenotypes小节及Table 1、Table 2
两部分模型构建
第一部分logistic回归估计NRS > 0概率,第二部分gamma分布log链接广义线性模型估计NRS > 0患者疼痛严重程度;排除完美预测患者37例;排除最高手术部位疼痛变量
阅读提示:Methods中Two-part Model小节及Table 6
模型验证与性能评估
随机拆分训练集n = 3,000,测试集n = 925;重复400次;性能指标R2、RMSE、MAE;Wilcoxon匹配配对符号秩检验,双侧α = 0.05
阅读提示:Methods中Model Validation小节及Table 3
协变量显著性比较与识别
在每个训练队列中识别显著协变量(P ≤ 0.05),统计400次迭代中显著次数;比较50%、75%、90%迭代阈值
阅读提示:Results中Covariate significance across model iterations部分及Table 4、Table 5