Skip to contents

This function runs the Partitioned Survival Model (PSM) that estimates the cost-effectiveness of 'Isotretinoin' vs 'Dinutuximab β' in the treatment of Neuroblastoma in children.

Usage

run_psm(models_fit = NeuroblastomaPSM::parametric_models, l_params)

Arguments

models_fit

A named list of fitted survival models.

l_params

A list of model parameters including the list of parameters in the list l_psm_parameters in addition to:

  • time_horizon: The time horizon for the model in years.

  • cycle_length: The length of a model cycle measured in years.

  • disc_rate_costs: The annual discount rate for incurred costs.

  • disc_rate_qalys: The annual discount rate for accrued QALYs.

Value

A vector containing two scalars: discounted costs and discounted QALYs for each treatment.

Examples

if (FALSE) {
# Load the fitted Gompertz model parameters
models_fit <- NeuroblastomaPSM::parametric_models

# Define model parameters
params <- c(
  time_horizon = 10,
  cycle_length = 1/12,
  disc_rate_costs = 0.035,
  disc_rate_qalys = 0.015,
  NeuroblastomaPSM::l_psm_parameters
)

# Run the Partitioned Survival Model
v_psm_results <- NeuroblastomaPSM::run_psm(
  models_fit = NeuroblastomaPSM::parametric_models,
  l_params = params
)

v_psm_results

# Run the Partitioned Survival Model
v_psm_results <- NeuroblastomaPSM::run_psm(
  l_params = params
)

v_psm_results
}