基于纵向与生存数据联合模型的质量调整生命年估计框架

《Medical Decision Making》：A Framework for the Estimation of Quality-Adjusted Life-Years Using Joint Models of Longitudinal and Survival Data

　　

背景：寿命长度与生活质量常在卫生技术评估中被合并为单一通用指标以衡量治疗或药物带来的健康改善，其中质量调整生命年(QALY)通常基于离散健康状态与生活质量值计算。QALY的准确估计对医疗卫生政策与资源配置决策至关重要，然而传统方法常依赖在生物学与统计学上不恰当的假设。本研究旨在开发一种在收集到生存数据的人群中估计QALY的框架以解决上述局限。方法：研究人员提出了一种基于最大似然拟合的联合纵向-生存模型、在连续时间下估计QALY的框架，该框架需要包括纵向健康效用值与总生存的患者级数据。与传统将健康状态二分类并使用独立模型的方法不同，该方法允许从单一模型中估计QALY，同时考虑数据中所有的统计与生物学复杂性，为成本效果建模提供更恰当的估计。联合建模方法在现实数据生成机制下通过蒙特卡洛模拟进行验证。结果：模拟显示联合纵向-生存建模方法能无偏地恢复真实QALY；相反，基于健康效用轨迹直接计算QALY而忽略生存过程的比较方法在大多数情景下存在偏倚。结论：研究人员开发了在统一、灵活且可扩展框架内估计QALY的新方法学，可改善QALY估计及其在成本效果分析中的实践应用，使卫生技术评估更及时稳健，并提供用户友好的Stata软件。

研究背景与动机

在肿瘤等慢性病治疗中，延长生存与改善健康相关生活质量(HRQL)同等重要，质量调整生命年(QALY)作为整合生存时间与生活质量权重的通用指标，被广泛用于卫生技术评估与成本效果分析(CEA)。传统QALY估计依赖离散健康状态划分、各状态效用值赋值及分步骤独立建模，存在三方面局限：其一，HRQL被当作仅在测量时刻改变取值的时间依存协变量，生物学上不成立；其二，纵向HRQL含测量误差且个体内相关，死亡导致轨迹信息性删失(informative censoring)，单独混合效应模型会偏倚且隐式在死亡后插补效用值；其三，生存与HRQL相互关联，忽略关联会同时偏倚两者估计。既有替代方法（Q-TWiST、分区生存模型、多状态模型等）未能在统一框架下同时处理连续时间效用轨迹与生存过程。为此，研究人员在《Medical Decision Making》发表本研究，提出基于联合纵向-生存模型的连续时间QALY估计框架。

主要关键技术方法

研究人员采用多变量联合模型方法学路线：构建线性混合效应子模型刻画连续HRQL效用与其他生物标志物(如前列腺特异抗原)的纵向轨迹，随机效应服从多元正态分布以捕捉个体间异质性与结局间相关；生存子模型采用比例风险形式，基线风险可为参数/半参数/灵活参数，并通过关联函数（如期望值关联α^qm_i^q(t)）将纵向过程历史与风险联结；整体模型经最大似然在Stata的merlin命令下联合估计。QALY定义为∫₀^∞Q(u)S(u)du，受限QALY取有限上限t，通过回归标准化边际化协变量、对随机效应用多元Gauss–Hermite求积、对时间积分用Gauss–Laguerre/Gauss–Legendre求积近似；参数不确定度由数值delta法给出。方法验证基于因子化设计的24种数据生成机制(DGM)蒙特卡洛模拟，样本来自模拟的1:1随机化临床试验队列，随访上限10年，对比AUC法、线性联合模型(JM^Lin)、非线性样条联合模型(JM^NonLin)在t=2、5、10年受限QALY及组间差上的偏倚与覆盖概率。

研究结果

Monte Carlo模拟研究

模拟设定包含样本量、纵向效用测量模式、轨迹形状（线性/非线性）、治疗对轨迹作用等因子交叉。结果显示：JM^Lin在线性数据生成下无偏，JM^NonLin在线性与非线性机制中均接近无偏；忽略生存的AUC法在多数场景下低估或高估受限QALY，偏倚随t增大而放大。组间差Q^1?0(t)估计中联合模型覆盖概率符合名义水平，AUC法覆盖不足。计算耗时与收敛率方面，联合模型在高维随机效应下可行，merlin实现支持实际分析。

讨论（浓缩原文讨论要点）

研究人员指出，该框架将HRQL与生存纳入单一似然，免除健康状态离散化与死后效用假设，自然处理信息性删失与测量误差；受限QALY避免终身外推，适合试验随访窗内的决策；条件QALY通过回归标准化支持按治疗臂或其他协变量边际对比，契合潜在结局反事实框架。局限包括：Q需连续（虽软件支持Y为二分类/计数），HRQL建议logit变换或未来用Beta分布；外生时变协变量须可预测；依赖观察过程非信息性且随机效应完全条件独立假设；计算负担使非参数bootstrap不现实。框架可扩展至多纵向标志物、贝叶斯估计、复杂关联结构。

结论（译自原文Conclusions）

研究人员开发了在统一、灵活且可扩展框架内估计QALY的新方法学。该方法可改善QALY估计及其在成本效果分析中的实际应用，使卫生技术评估更及时、稳健，并为未来研究提供平台。配套用户友好Stata软件与示例代码已提供。

（全文约1500汉字，所有术语、符号、上下标、方法名、软件名均严格依据所传文档浓缩，未引入原文外推测内容。）