Page 106 - 《中国药房》2022年11期

表2 不同曲线的最优拟合分布及参数值

PFS 曲线类型 最优拟合分布 生存函数 分布参数 AIC BIC
D-Rd组OS曲线 Exponential 0.006 956 699.87 703.79
D-Rd组PFS曲线 Exponential S（t）＝e －λt 0.010 869 874.07 877.98
Rd组OS曲线 Exponential 0.010 002 820.82 824.73
PD
Rd组PFS曲线 Exponential 0.020 193 922.95 926.86
S：生存率；t：时间；λ：指数分布的参数值
Death 1.5 成本和效用值数据
本研究因基于我国卫生体系角度评估D-Rd方案对
图1 3种健康状态的分区生存模型结构
比Rd方案治疗TNE-NDMM的经济性，故仅考虑直接医
1.4 临床数据 疗成本，主要包括药品成本、疾病管理成本、实验室检测
当模拟时间处于试验随访期内时，PFS 健康状态的 成本和不良反应处理成本。药品成本数据来源于药智
人数比例直接从 PFS 曲线获得，PD 状态的人数比例等 网（https：//www.yaozh.com）2021 年全国各省市药品挂
于总生存（overall survival，OS）曲线生存率减去 PFS 曲 网与中标价格的中位值。疾病管理成本包括诊断费、静
线生存率，死亡人数比例等于1减去OS曲线生存率。当 脉输液费、护理费和床位费，数据来源于已发表文
模拟时间超出随访期时，使用参数法 对生存曲线进行 献 。根据《中国多发性骨髓瘤诊治指南（2020 年修
[13]
[14]
外推获得生存函数。本研究首先通过 Get Data Graph 订）》 ，实验室检测成本包括血液检查、尿液检查、骨髓
[11]
[8]
Digitizer 2.26 软件从 MAIA 试验 的 K-M 曲线中提取原 检查、影像学检查和其他检查的成本，数据来源于已发表
始数据，再通过 Stata MP 16.0 软件重构 2 组患者的 K-M 文献 。不良反应处理成本来源于已发表的文献 [14，16－17] 。
[15]
曲线，并对重构曲线进行视觉检查。分别选取Exponen- 临床上，当发生严重不良反应的时候，医生可能会停药
tial 分布、Weibull 分布、Gompertz 分布、Lognormal 分布、 或换药，因此本研究假设不良反应处理成本只考虑 1
[18]
Loglogistic 分布、Generalized gamma 分布对生存曲线的 次。模型中患者的效用值来源于Lu等 的研究，这是根
参数进行拟合，根据赤池信息准则（Akaike information 据其他国家或地区 NDMM 患者的效用值，应用其他国
criterion，AIC）和贝叶斯信息准则（Bayesian information 家与中国患者的年龄和性别相匹配的普通人群相关的
criterion，BIC）结合视觉检查进行拟合优度检验（AIC和 生活质量比率，估算得到的中国 NDMM 患者的效用
BIC越小，表明模型的拟合优度越好）。2组曲线的参数 值。研究中所有成本均采用 Gamma 分布，效用值采用
分布拟合结果见表 1。根据检验结果，将最优拟合分布 Beta分布，并对成本和效用值采用半周期校正。假设患
的参数值代入生存函数中，求得各周期各状态下的患者 者体质量为65 kg，模型参数及分布见表3。
比例。不同曲线的最优拟合分布及参数值见表 2。因 1.6 敏感性分析
1～2级药物不良反应一般无须处理，且D-Rd和Rd治疗 为了检验模型的稳健性，本研究进行了单因素敏感
方案的不良反应发生率较高，故本研究仅纳入MAIA试 性分析（deterministic sensitivity analysis，DSA）和概率敏
验中 3 级及以上且发生率≥10％的不良反应，包括贫 感性分析（probabilistic sensitivity analysis，PSA）。采用
血、白细胞减少症、中性粒细胞减少症、淋巴细胞减少 DSA 来评估模型参数的不确定性对增量成本-效果比
症、白内障、低钾血症和肺炎，相应不良反应发生率从 （incremental cost-effectiveness ratio，ICER）的影响程度，
结果用旋风图呈现。模型参数变化范围按基线值±20％
MAIA 试验中提取获得。
[19]
表1 2组治疗方案OS和PFS曲线的参数分布拟合结果 计算 ，根据《中国药物经济学评价指南 2020》，贴现率
变化范围设定为0～8％ ，各参数的下限值和上限值见
[12]
D-Rd组 Rd组
曲线类型 参数分布 表3。采用PSA来评估多个不确定参数分布变化对模型
AIC BIC AIC BIC
OS曲线 Exponential 699.87 703.79 820.82 824.73 输入的影响，利用二阶蒙特卡洛模拟对各个参数的分布
Weibull 700.15 707.97 817.54 825.35 随机抽样1 000次，结果以成本-效果散点图和成本-效果
Gompertz 701.86 709.68 816.86 824.67
Lognormal 713.61 721.43 836.47 844.29 可接受曲线呈现。
Loglogistic 702.71 710.53 820.85 828.66 2 结果
Generalized gamma 700.48 712.21 818.70 830.42
PFS曲线 Exponential 874.07 877.98 922.95 926.86 2.1 基础分析结果
Weibull 875.39 883.20 922.26 930.08 基础分析结果显示，相较于Rd组，D-Rd组患者可以
Gompertz 875.17 882.99 924.59 932.41 获得更多质量调整生命年（quality-adjusted life year，
Lognormal 878.71 886.53 920.38 928.20
Loglogistic 875.55 883.37 919.53 927.35 QALY），但相应成本也更高，详见表 4。在 20 年研究时
Generalized gamma 876.91 888.63 919.58 931.32 限内，与Rd组相比，D-Rd组可多获得1.28个QALYs，但

·1376 · China Pharmacy 2022 Vol. 33 No. 11 中国药房 2022年第33卷第11期