Page 81 - 《中国药房》2026年8期

mental cost-effectiveness ratio，ICER）来评价方案的经济 表1 生存曲线拟合分布的AIC、BIC值
性。参考《中国药物经济学评价指南（2025 中英双语 生存曲线 组别 参数 Weibull Log-logistic Log-normal Gompertz Exponential Gen-gamma Gamma
[13]
版）》 ，本研究中的成本和效用值均采用 4.5% 的折现 PFS曲线 联合组 AIC 544.403 7 543.883 2 542.196 5 549.247 2 556.285 9 544.187 8 543.566 6
BIC 550.943 5 550.418 9 548.732 2 555.782 9 559.553 8 553.991 3 550.102 3
率，意愿支付（willingness-to-pay，WTP）阈值设为2025年
对照组 AIC 811.039 3 811.181 5 815.932 4 814.621 5 822.145 2 812.899 5 810.915 5
我国 2 倍人均国内生产总值（gross domestic product， BIC 817.512 2 817.654 4 822.405 3 821.094 4 825.381 6 822.608 8 817.388 4
GDP），即199 330元/QALY。 OS曲线 联合组 AIC 264.566 0 264.249 3 262.879 6 268.324 8 275.378 7 264.030 9 263.912 5
BIC 271.101 7 270.785 0 269.415 3 274.860 5 278.646 5 273.834 5 270.448 3
对照组 AIC 173.340 1 173.234 7 174.123 4 174.814 1 184.441 5 175.315 7 173.408 5
PFS PD BIC 178.530 4 178.424 9 179.313 6 180.004 3 187.036 6 183.101 1 178.598 7
表2 最优拟合分布及参数
死亡
组别 生存曲线 最优拟合分布 参数值
图1 分区生存模型 联合组 PFS曲线 Log-normal meanlog＝2.908；sdlog＝1.026
OS曲线 Log-normal meanlog＝3.709；sdlog＝0.903
1.3 生存曲线的拟合和外推 对照组 PFS曲线 Gamma shape＝1.643；rate＝0.070
本研究拟基于 CULMATE-1 研究的总生存（overall OS曲线 Log-logistic shape＝2.741；scale＝25.957
survival，OS）曲线与 PFS 曲线，获取不同健康状态下的 meanlog：平均值；sdlog：标准差；shape：形状参数；rate：率参数；
scale：尺度参数。
患 者 分 布 比 例 ，进 而 开 展 成 本 - 效 用 分 析 。 由 于
CULMATE-1 研究仅公布了随访期内 OS 曲线与 PFS 曲 成本、最佳支持治疗成本及临终关怀治疗成本等。为简
线对应的患者比例数据，对于超出随访期的长期生存数 化模型，本研究根据《美国国家癌症研究所不良事件通
据，则需借助参数法对生存曲线进行重构、拟合与外推。 用术语标准 5.0 版》，仅纳入 CULMATE-1 研究报告的 3
对此，本研究首先使用 Web Plot Digitizer 软件从原始生 级及以上且发生率≥5%的不良反应，包括贫血、中性粒
存曲线取点；之后，利用 R 语言 4.4.3 软件重构两组患者 细胞减少症、白细胞减少症以及低血钾等。不良反应发
[10]
的个体水平数据，并采用相应的参数分布（包括 Expo‐ 生率均来自CULMATE-1研究 。药品价格为药智网中
nential、Gamma、Gompertz、Weibull、Log-logistic、Log- 2025 年各药品在各省份中标价的中位数。由于药智网
normal、Gen-gamma 分布）对重构的个体水平数据进行 未公开库莫西利胶囊的挂网价格，笔者通过咨询该药品
[14]
生存分析拟合 ；最后，根据赤池信息量准则（Akaike in‐ 生产企业获取库莫西利胶囊（60 mg×28粒）的官方拟定
formation criterion，AIC）和贝叶斯信息准则（Bayesian 价格为3 010.13元，此价格与北京、重庆、广东、上海等地
information criterion，BIC），选择 AIC、BIC 值最小的分 区药房实际销售价格一致。其余所需成本数据均来源
[15]
布，并结合视觉检验选取最终的最优拟合分布 。本研 于已发表的相关文献 [11，16―20] ，并通过中国居民消费价格
究最终选择 Log-normal 分布来拟合联合组患者的 PFS 指数调整至2025年价格水平，确保成本核算的时效性与
曲线和 OS 曲线，选择 Gamma、Log-logistic 分布来拟合 准确性。
对照组的PFS曲线和OS曲线。两组患者的生存曲线拟 由于CULMATE-1研究无针对我国HR+/HER2－乳
合情况及AIC、BIC值分别见图2、表1，最优拟合分布及 腺癌患者的健康效用值，故本研究的效用值参考其他文
参数见表2。 献 [21―22] ，分 别 取 PFS 和 PD 状 态 的 效 用 值 为 0.837 和
100 100 0.443。此外，本研究还考虑了不良反应导致的效用损
80 80 失 [18―20，23] 。模型参数设置见表3。
生存率/% 60 库莫西利联合氟维司群：原始 生存率/% 60 库莫西利联合氟维司群：原始 1.5 敏感性分析
Gamma
Gamma
40
40
Gompertz
Gompertz
Exponential
20 Log-normal 20 Log-normal
Exponential
Gen-gamma
Gen-gamma
Log-logistic Log-logistic 本研究围绕各项参数开展单因素敏感性分析和概
0 Weibull 0 Weibull
0 10 20 30 0 10 20 30 率敏感性分析，以验证模型的稳定性。药品价格的取值
时间/月 时间/月
A.联合组PFS拟合曲线 B.联合组OS拟合曲线 范围以 2025 年药智网公布的中标最低价与最高价为参
100 100 照；其余参数的取值范围优先依据其 95% 置信区间设
80 80 [22]
生存率/% 60 氟维司群：原始 生存率/% 60 氟维司群：原始 定；对于未明确 95% 置信区间的参数，参考张森等 的
Gamma
Gamma
40
Gompertz
40
Gompertz
Exponential
Exponential
20 Log-normal 20 Log-normal 研究，将其变动范围设为基准值的±20%。通过改变参
Gen-gamma
Log-logistic Gene-gamma
Log-logistic
0 Weibull 0 Weibull 数的变化临界值进行单因素敏感性分析，以探究参数变
0 10 20 30 0 10 20 30
时间/月 时间/月 化对ICER的影响程度。概率敏感性分析通过二阶蒙特
C.对照组PFS拟合曲线 D.对照组OS拟合曲线 卡罗模拟1 000次来进行，结果以成本-效果散点图和成
图2 两组患者生存曲线拟合情况
本-效果可接受曲线来展示。概率性敏感性分析中，假
1.4 模型参数 设本研究涉及的成本参数均服从 Gamma 分布，效用值
本研究从中国卫生体系角度出发，仅纳入直接医疗 和不良反应发生率参数均服从Beta分布，其余参数同样
成本，包括二线治疗药物成本、随访成本、不良反应处理 服从Gamma分布（见表3）。
中国药房 2026年第37卷第8期 China Pharmacy 2026 Vol. 37 No. 8 · 1035 ·