2.1 Model structure
To compare the cost-effectiveness of pembrolizumab versus chemotherapy in patients with MSI-H–dMMR stage IV CRC, a PSM was conducted to simulate the process of dMMR/MSI-H mCRC depending on the clinical data from KEYNOTE 177. Three distinct health states: PFS, progressive disease (PD) and death (Fig. 1) were included in the model. The initial state is assumed to be PFS, and death is the absorbing state. The population was a cohort with the same characteristics and treatments as those in the KEYNOTE 177 trial. Unlike Markov model, the PSM use of the trials' Kaplan-Meier (K-M) curves to directly divide patients into different health states without putting forward assumptions for the transition of patients between different health states. Therefore, the estimation of patients’ proportion in each health state was acquired straightly from the cumulative survival probabilities in OS and PFS curves by parametric functions fitting and extrapolation. The cycle length was 6-week. The analysis was conducted from the Chinese health care system perspective with a life-time horizon so as to ensure there were less than 1% survivors. Costs, life years (LYs), quality-adjusted life years (QALYs), and incremental cost-effectiveness ratios (ICERs) were calculated in each treatment group. If the ICER is below $33,600 threshold (three times GDP per capita of China in 2020, ¥210,000.00), the treatment is generally considered to be cost-effective. In line with Chinese pharmacoeconomic guidelines[22], both costs and benefits are discounted at 5% (range: 0%-8%) per year.
2.2 Efficacy Estimates
We used the Engauge Digitizer (version 12.1, http://digitizer.sourceforge.net) to collect the data points from the K-M curves (PFS and OS curves) of the two arms and followed the method of Guyot et al.[23] to reconstruct estimates of underlying individual patient data (IPD) over the clinical trial time. In terms of the IPD out of the clinical trial time, standard parametric models fitting and extrapolation were used to estimate the long-term survival probabilities by using R software (version 4.1.0, https://www.r-project.org). Specifically, six parametric functions were considered, including exponential, Weibull, Gompertz, log-logistic, log-normal, and generalized gamma distributions. And then, in order to evaluate the goodness-of-fit of each parametric survival model, multiple methods were applied such as visual inspection, Akaike information criterion (AIC) test and Bayesian information criteria (BIC) test proposed by NICE DSU technical support document 14 (TSD14)[24]. Lower AIC and BIC values indicate better fit of the selected parametric model. (eTable 1–4 are given in Online Resource) Meanwhile, superposed graphs of the K-M curves from the trial and the estimated curves based on the relative better fitting parametric survival models were presented in Fig. 2 so as to intuitively inspect the survival prediction.
The generalized gamma model was chosen as the best fit model for the OS curve of both arms and PFS of pembrolizumab arm and log-normal model for the PFS curve of chemotherapy arm. Considerations were as follows: 1) the lowest or relative lower AIC and BIC values among all survival models; 2) the best fit with the observed curves based on visual inspection.
2.3 Utility Estimates
The health utility score reflects the level of physical, mental, and social functioning associated with a disease correlative health state that varies from 0 to 1, with 0 representing the worst health state/death and 1 representing the best. Average health utility of the PFS state was 0.7825, which was based on quality-of-life data collected in the KEYNOTE 177 trial.[25] In the trial, EQ-5D-3L index[26] mean utility scores were 0.77 in the pembrolizumab group and 0.75 in the chemotherapy group at baseline. At the end of week 18, mean scores were 0.84 and 0.77 respectively. So, in the simulation, we assigned a utility of 0.7825 for PFS, which is the mean of above values, with 0.75 and 0.84 as the boundaries of the range used in sensitivity analyses. Average health utility of the PD state was 0.64 (95% CI [0.576–0.704]) derived from the previously published literature[27]. The disutility values of grade 3–4 adverse events (AEs) were considered in our analysis,[28, 29] but only one-time assessment was carried out during the first cycle for simplification given the trivial influence of AE disutilities. QALYs loss caused by AEs were assessed by the product of the incidence of AEs and the corresponding disutility value. And ± 20% were the boundaries of the range in sensitivity analyses.
2.4 Cost Inputs
Only direct medical costs were considered, including cost of the drug utilization, cost of drug administration, cost of follow up, cost of main AEs management and cost of treatments for progression (including active treatments and supportive care), and stated in 2021 United States dollars (USD) by exchange rate: 1 CYN = 0.16 USD. Drug prices were estimated from the local bid-winning price.[30]
Drug administration costs were calculated as a function of administration cost per attendance and administration frequency (number of attendance per cycle). Unit administration costs were derived from local charge. The costs of follow up, including carcinoembryonic antigen level test, ultrasound of abdomen, outpatient specialist clinic, were from a published study[31], and were converted to 2020 USD using Medical Care component of Chinese Consumer Price Index (CPI). The frequency of follow up were based on the guideline of Chinese society of clinical oncology (CSCO) for colorectal cancer[12]. The costs of main AEs management were derived from the previously published literatures[32, 33] (adjusted by CPI) and were calculated only once in the first cycle. In terms of the treatments for progression, active therapeutic schedules included CAPIRI, CAPOX, FOLFOX (± bevacizumab or cetuximab), FOLFIRI (± bevacizumab or cetuximab), nivolumab, pembrolizumab.[34, 35] The prices of the drugs in subsequent treatments were also from local bid-winning price[30], and the fee of supportive care was assumed to be 0. Additionally, costs were discounted at an annual rate of 5%[22].
2.5 Clinical Inputs
According to the global phase III clinical trial - KEYNOTE-177[17], the patients in the pembrolizumab group (P group, N = 153) received pembrolizumab at a dose of 200 mg every 3 weeks intravenously (IV) for up to 35 treatments (approximately 2 years) and the patients in the chemotherapy group (C group, N = 154) received FOLFOX (oxaliplatin 85 mg/m2 IV on Day 1, leucovorin 400 mg/m2 IV on Day 1, 5-FU 400 mg/m2 IV bolus on Day 1 and then 1200 mg/m2/day IV over 2 days for total dose of 2400 mg/m2 in each 2-week cycle) or FOLFIRI (irinotecan 180 mg/m2 IV on Day 1, leucovorin 400 mg/m2 IV on Day 1, 5-FU 400 mg/m2 IV bolus on Day 1 and then 1200 mg/m2/day IV over 2 days for total dose of 2400 mg/m2 in each 2-week cycle), with or without either cetuximab (cetuximab 400 mg/m2 IV over 2 hours then 250 mg/m2 over 1 hour weekly in each 2-week cycle) or bevacizumab (5 mg/kg IV on Day 1 of each 2-week cycle). Treatment was continued until disease progression, development of unacceptable toxic effects, illness, or a decision by the physician or patient to withdraw from the trial. In C group, the percentage of patients in mFOLFOX, mFOLFOX + bevacizumab, mFOLFOX + cetuximab, FOLFIRI, FOLFIRI + bevacizumab and FOLFIRI + cetuximab were 7.69%, 44.76%, 3.50%, 11.19%, 25.17% and 7.69%, respectively. After disease progression, patients randomly assigned to the C group could cross over to pembrolizumab (to receive a maximum of 35 treatments) and patients in the P group can continue pembrolizumab (to receive a maximum of 17 treatments). The subsequent treatments were mainly composed of pembrolizumab (P group 10%, C group 42.75%), other PD-1/PD-L1 ICIs (P group 7.5%, C group 28.24%), chemotherapy (P group 43.75%, C group 15.27%), VEGF inhibitor (P group 27.5%, C group 9.92%), and EGFR (P group 11.25%, C group 3.82%) inhibitor. Referring to the other studies and guideline[34, 35], we assumed that CAPIRI, CAPOX, FOLFOX and FOLFIRI were used as the standard second-line chemotherapy, nivolumab (up to 2 years) and pembrolizumab as the ICIs, cetuximab as the EGFR inhibitor and bevacizumab as the VEGF inhibitor, which are commonly used in China. We included grade 3 to 4 AEs in the model that had obvious clinical impact and significantly different rates between the arms of the KEYNOTE 177 trial, which were diarrhea, anemia, hypokalemia, and neutropenia. The incidence of neutrophil count decrease multiplied by 0.5 is included in neutropenia. For dosage calculation, assuming that values of body surface area (BSA) and weight were 1.80 m2 and 65kg, respectively.[36]
2.6 Sensitivity Analysis
A series of sensitivity analyses were conducted to test the robustness of the model and address uncertainty in the estimation of model parameters by using Microsoft Excel (version 16.51). One-way deterministic sensitivity analyses (DSA) were used to evaluate the impact of uncertainty of a single input variable on the ICER. The range of drug prices depends on the local charge or ± 20% of the baseline values. Other parameters were adjusted within the reported 95% confidence intervals (CI) or assuming reasonable ranges of the base case values (± 20%) if 95% CIs were unavailable, in accordance with established approach[37]. In addition, a separate scenario considering patient assistance program (PAP) was evaluated in the DSA.
In the probabilistic sensitivity analysis (PSA), a Monte Carlo simulation of 500 iterations was generated by randomly sampling the key model parameters from the pre-specified distributions simultaneously. We used gamma distribution for the cost parameters and beta distribution for utility and probability parameters. Based on the data from 500 iterations, a cost-effectiveness acceptability curve (CEAC) was created to represent the likelihood that pembrolizumab would be considered cost-effective compared with chemotherapy on the basis of a willingness to pay (WTP) threshold of $33,600 per QALY in China.
Baseline values, ranges, and distributions for sensitivity analysis were summarized in Table 1
Table 1
Baseline values, ranges, and distributions for sensitivity analysis
Variable | Baseline value | Range | α | β | mean | SE | Reference for baseline value | Distribution (parameters) |
Minimum | Maximum |
Cycle | 6-week | - | - | - | - | - | - | - | - |
Horizon | Life-time | - | - | - | - | - | - | - | - |
WTP, $/QALY | 33600 | - | - | - | - | - | - | [21] | - |
Discount rate | 0.05 | 0 | 0.08 | - | - | - | - | [21] | - |
Mean body surface area, m2 | 1.80 | - | - | - | - | - | - | - | - |
Patients’ weight, kg | 65 | - | - | - | - | - | - | - | - |
Pembrolizumab group _ AEs incidence |
Diarrhea | 0.060 | 0.048 | 0.072 | 90.218 | 1413.409 | 0.060 | 0.006 | [16] | Beta |
Anemia | 0.050 | 0.040 | 0.060 | 91.188 | 1732.572 | 0.050 | 0.005 | Beta |
Hypokalemia | 0.010 | 0.008 | 0.012 | 95.070 | 9411.890 | 0.010 | 0.001 | Beta |
Neutropenia | 0.000 | 0.000 | 0.000 | - | - | 0.000 | 0.000 | Beta |
Chemotherapy group _ AEs incidence |
Diarrhea | 0.110 | 0.088 | 0.132 | 85.366 | 690.685 | 0.110 | 0.011 | [16] Assumed the incidence of neutrophil count decrease multiplied by 0.5 is included in neutropenia | Beta |
Anemia | 0.100 | 0.080 | 0.120 | 86.336 | 777.024 | 0.100 | 0.010 | Beta |
Hypokalemia | 0.060 | 0.048 | 0.072 | 90.218 | 1413.409 | 0.060 | 0.006 | Beta |
Neutropenia | 0.235 | 0.188 | 0.282 | 73.236 | 238.405 | 0.235 | 0.024 | Beta |
Pembrolizumab group second-line therapy proportion |
Pembrolizumab | 0.100 | 0.080 | 0.120 | 86.336 | 777.024 | 0.100 | 0.010 | KEYNOTE 177, assumed only the second-line treatments of a proportion of more than 2% were considered in our analysis. | Beta |
other PD-1/PD-L1 Checkpoint Inhibitor | 0.075 | 0.060 | 0.090 | 88.762 | 1094.731 | 0.075 | 0.008 | Beta |
Chemotherapy | 0.438 | 0.350 | 0.525 | 53.585 | 68.895 | 0.438 | 0.045 | Beta |
VEGF inhibitor | 0.275 | 0.220 | 0.330 | 69.354 | 182.842 | 0.275 | 0.028 | Beta |
EGFR Inhibitor | 0.113 | 0.090 | 0.135 | 85.123 | 671.526 | 0.113 | 0.011 | Beta |
Chemotherapy group second-line therapy proportion |
Pembrolizumab | 0.427 | 0.342 | 0.513 | 54.557 | 73.068 | 0.427 | 0.044 | KEYNOTE 177, assumed only the second-line treatments of a proportion of more than 2% were considered in our analysis. | Beta |
other PD-1/PD-L1 Checkpoint Inhibitor | 0.282 | 0.226 | 0.339 | 68.632 | 174.362 | 0.282 | 0.029 | Beta |
Chemotherapy | 0.153 | 0.122 | 0.183 | 81.225 | 450.797 | 0.153 | 0.016 | Beta |
VEGF inhibitor | 0.099 | 0.079 | 0.119 | 86.410 | 784.338 | 0.099 | 0.010 | Beta |
EGFR Inhibitor | 0.038 | 0.031 | 0.046 | 92.336 | 2326.872 | 0.038 | 0.004 | Beta |
Health preferences |
Utility of PFS | 0.7825 | 0.75 | 0.84 | 245.001 | 63.176 | 0.795 | 0.023 | [24] | Beta |
Utility of PD | 0.64 | 0.576 | 0.704 | 137.658 | 77.432 | 0.640 | 0.033 | [26] | Beta |
Disutility due to AEs (grade ≥ 3) |
Diarrhea | -0.090 | -0.072 | -0.108 | 14.470 | 146.308 | -0.090 | 0.023 | [27][28] | Beta |
Anemia | -0.085 | -0.068 | -0.102 | 14.555 | 156.680 | -0.085 | 0.021 | Beta |
Hypokalemia | -0.080 | -0.064 | -0.096 | 14.640 | 168.360 | -0.080 | 0.020 | Beta |
Neutropenia | -0.0607 | -0.049 | -0.073 | 14.968 | 231.623 | -0.061 | 0.015 | Beta |
Drug cost, $/mg |
Pembrolizumab | 179.18 | - | - | - | - | - | - | 2021 Local charge https://www.yaozh.com | Fix |
Oxaliplatin | 3.4 | 2.72 | 4.08 | 16.000 | 0.213 | 3.400 | 0.850 | Gamma |
Leucovorin | 0.25 | 0.2 | 0.3 | 16.000 | 0.016 | 0.250 | 0.063 | Gamma |
5-FU | 0.29 | 0.232 | 0.348 | 16.000 | 0.018 | 0.290 | 0.073 | Gamma |
Bevacizumab | 15 | 11.88 | 15 | 16.000 | 0.840 | 13.440 | 3.360 | Gamma |
Cetuximab | 12.04 | - | - | - | - | - | - | Fix |
Irinotecan | 17.73 | 14.05 | 24.9 | 16.000 | 1.217 | 19.475 | 4.869 | Gamma |
Nivolumab | 96.2 | - | - | - | - | - | - | Fix |
Capecitabine | 0.04 | 0.01 | 0.04 | 16.000 | 0.002 | 0.025 | 0.006 | Gamma |
AEs cost, $ |
Diarrhea | 392.44 | 313.95 | 470.93 | 16.000 | 24.527 | 392.439 | 98.110 | [32] | Gamma |
Anemia | 724.64 | 579.71 | 869.57 | 16.000 | 45.290 | 724.638 | 181.160 | [31] | Gamma |
Hypokalemia | 157.28 | 125.82 | 188.73 | 16.000 | 9.830 | 157.275 | 39.319 | [31], assumed same as fatigue | Gamma |
Neutropenia | 628.96 | 503.17 | 754.76 | 16.000 | 39.310 | 628.965 | 157.241 | [31] | Gamma |
Administration, $/attendance | 310.16 | 248.13 | 372.19 | 16.000 | 19.385 | 310.160 | 77.540 | Local charge | Gamma |
Cost of BSC, $/cycle | 0.00 | 0.00 | 0.00 | - | - | - | - | Assumed to be 0 | - |
Follow up, $/cycle | 31.33 | 25.06 | 37.60 | 16.000 | 1.958 | 31.330 | 7.833 | [11][30] | Gamma |
PD-1, programmed death 1; PD-L1, programmed death ligand 1; EGFR, epidermal growth factor receptor; VEGF, vascular endothelial growth factor; PD, progressive disease; PFS, progression-free survival; AE, adverse event; BSC, best supportive care. |