Page 72 - 《中国药房》2025年20期

1.3 模型构建 表2 生存曲线的拟合结果
1.3.1 动态Markov模型 生存曲线 分布 AIC BIC 最优拟合分布参数
本研究运用TreeAge Pro Healthcare 2022软件，参照 PemEV组OS曲线 Exponential a 168.708 171.251 rate＝0.010
Gamma 170.680 175.767
泛亚洲 aUC 患者的疾病发展过程构建动态 Markov 模
Gompertz 170.505 175.591
型。该模型涵盖 PFS、PD 及死亡 3 种互斥的健康状态， Weibull 170.695 175.782
患者在每个周期内只能处于1种健康状态；所有患者最 Log-logistic 170.451 175.537
初进入模型时均处于 PFS 状态，且 PD 状态不能逆转为 Log-normal 170.057 175.143
PemEV组PFS曲线 Exponential 306.819 309.362
PFS状态，而死亡状态的患者无法再从该状态转出。据模
Gamma 308.792 313.878
型预测，20年后PemEV组和GP组中均有超过90%的患 Gompertz 306.194 311.281
者处于死亡状态，因此本研究将研究时限设为20年，模型 Weibull 308.806 313.892
的循环周期参照给药方案设为21 d。模型产出指标为成 Log-logistic 306.208 311.294
Log-normal a 303.469 308.556 meanlog＝3.172，sdlog＝1.502
本及质量调整生命年（quality-adjusted life years，QALYs）， GP组OS曲线 Exponential 282.405 284.812
评价指标为 PemEV 相较于 GP 的 ICER，WTP 阈值采用 Gamma 281.553 286.366
3 倍我国 2024 年人均国内生产总值 [18―19] ，即 287 247 Gompertz 284.386 289.199
元/QALY，以年贴现率5%对成本和健康产出进行贴现 。 Weibull 282.418 287.231
[20]
Log-logistic 279.906 284.719
1.3.2 生存分析 Log-normal a 277.474 282.287 meanlog＝2.945，sdlog＝1.140
采用 GetData Graph Digitizer 2.24 软件从 EV-302 试 GP组PFS曲线 Exponential 331.399 333.806
验泛亚洲亚组的原始生存曲线中提取数据，获得各时间 Gamma 326.389 331.202
Gompertz 333.041 337.855
[21]
点的累积生存概率。参考 Guyot 等 披露的方法，应用
Weibull 328.771 333.585
R 4.2.2 软件“survHE”程序包重构生存曲线。分别采用 Log-logistic 321.939 326.753
Exponential、Gamma、Gompertz、Weibull、Log-logistic 及 Log-normal a 318.797 323.610 meanlog＝1.810，sdlog＝0.919
Log-normal等常见分布对重构后的生存曲线进行拟合， a：最优拟合分布；rate：率参数，为Exponential分布的参数；mean‐
选取赤池信息准则（Akaike’s information criterion，AIC） log：对数均值，为Log-normal分布的参数；sdlog：对数标准差，为Log-
和贝叶斯信息准则（Bayesian information criterion，BIC） normal分布的参数。
最小的分布作为最优拟合分布。最终，除 PemEV 组 OS 围。在不良事件方面，为简化模型，本研究仅纳入3级及
曲线选取 Exponential 分布外，其余曲线均选取 Log- 以上且发生率≥20%的TRAEs处理成本，且将其视为一
normal分布。结果见表2。 次性成本计入模型首周期。上述成本数据均来源于已
参照周挺等 提出的方法，利用前述拟合参数计算 公开发表的文献 [23―24] ，并依据我国医疗保健类居民消费
[22]
动态Markov模型中各状态的转移概率。其中，第1周期 价格指数 将其调整为2024年12月的水平。
[19]
PFS状态转移至死亡状态的概率采用OS曲线的累积死 EV-302 试验暂未披露泛亚洲亚组患者准确的状态
亡函数，后续周期PFS状态转移至死亡状态的概率则假 效用值，鉴于KEYNOTE-052试验与EV-302试验的目标
[19]
设为2024年中国人口的自然死亡率 。 人群均为既往未接受过相关系统治疗的成年aUC患者，
1.3.3 模型参数
且治疗方案相似（试验组均含免疫检查点抑制剂，对照
本研究基于我国卫生体系的视角，测算aUC患者使 [25] [26]
组均为单纯化疗） ，故本研究参考 Patterson 等 关于
用 PemEV 方案或 GP 方案的直接医疗成本，包括药品、
KEYNOTE-052 试验健康相关生命质量数据的研究，将
随访、BSC、临终关怀和 TRAEs 处理的成本。维恩妥尤
PFS 状态的效用值设为 0.842，PD 状态的效用值设为
单抗、帕博利珠单抗因处于专利保护期内，仅有单一价
0.800。
格，故无需对其费用进行处理；吉西他滨、顺铂和卡铂的
[14]
价格来源于药智网2024年的药品中标信息 ，以其中位 本研究假设中国 aUC 患者的平均体重为 65 kg，平
[27―28]
2
中标价格作为基线值。 均体表面积为1.72 m ，平均肌酐清除率为70 mL/min 。
参数值及其分布见表3。
aUC 随访包括血液检验、尿液检验（含尿细胞学）、
影像学检查（含CT尿路造影）及膀胱镜检查等 [1，16] ，其成 1.4 敏感性分析
本数据参考本地医院收费标准确定。自患者进入 PFS 通过调整单个或多个参数值，观察参数调整之后
[18]
状态起开始计算，至死亡结束，aUC 随访频率为每 12 ICER的变动情况来评价基础研究结果是否稳健 。
[16]
周1次 。 1.4.1 单因素敏感性分析
BSC 涵盖营养治疗及并发症管理，本研究假设 PD 逐一调整每个参数的取值，其上下限为参数的
[29]
后未接受后线治疗或后线治疗已停止的患者接受 BSC 95%CI 或基线值的±20%（当 95%CI 无法获得时） ，借
直至死亡。除此以外，本研究还假设患者在生命尽头会 此探究单个参数在其预设范围内变化对结果的具体影
接受临终关怀，因此也将临终关怀的成本纳入考虑范 响，分析结果采用飓风图展示。

· 2550 · China Pharmacy 2025 Vol. 36 No. 20 中国药房 2025年第36卷第20期