simulate_tepi: Simulate operating characteristics using TEPI

In phase12designs: Comprehensive Tools for Running Model-Assisted Phase I/II Trial Simulations

Simulate operating characteristics using TEPI

Description

This function runs simulations of the TEPI design by evaluating operating characteristics over a range of cohort sizes. For each dose level within the user-specified range, it performs multiple trials and saves the results to a corresponding file.

Usage

simulate_tepi(
  ndose,
  ssizerange,
  target_t,
  lower_e,
  cohortsize = 3,
  startdose = 1,
  effint_l = c(0, lower_e, lower_e + 0.2, lower_e + 0.4),
  effint_u = c(lower_e, lower_e + 0.2, lower_e + 0.4, 1),
  toxint_l = c(0, 0.15, target_t, target_t + 0.05),
  toxint_u = c(0.15, target_t, target_t + 0.05, 1),
  psafe = 0.95,
  pfutility = 0.9,
  ntrial = 10000,
  utilitytype = 1,
  u1,
  u2,
  prob = NULL,
  save_dir = tempdir(),
  save_folder = "tepi_simulations",
  save_file = "tepi2_simulation.csv"
)

Arguments

ndose	Integer. Number of dose levels. (Required)
ssizerange	Integer vector. Range of number of cohorts to simulate. (Required)
target_t	Numeric. Target toxicity probability. (Required)
lower_e	Numeric. Minimum acceptable efficacy probability. (Required)
cohortsize	Integer. Size of a cohort. (Default is 3)
startdose	Integer. Starting dose level. (Default is 1)
effint_l	Lower efficacy bounds for dose assignment decision table. (Default is c(0,lower_e,lower_e+0.2,lower_e+0.4))
effint_u	Lower efficacy bounds for dose assignment decision table. (Default is c(lower_e,lower_e+0.2,lower_e+0.4,1))
toxint_l	Lower toxicity bounds for dose assignment decision table. (Default is c(0,0.15,target_t,target_t+0.05))
toxint_u	Lower toxicity bounds for dose assignment decision table. (Default is c(0.15,target_t,target_t+0.05,1))
psafe	Numeric. Early stopping cutoff for toxicity. (Default is 0.95)
pfutility	Numeric. Early stopping cutoff for efficacy. (Default is 0.90)
ntrial	Integer. Number of random trial replications. (Default is 10000)
utilitytype	Integer. Type of utility structure. (Default is 1) If set to 1: Use preset weights (w11 = 0.6, w00 = 0.4) If set to 2: Use (w11 = 1, w00 = 0) Other: Use user-specified values from u1 and u2.
u1	Numeric. Utility parameter w_11. (0-100)
u2	Numeric. Utility parameter w_00. (0-100)
prob	Fixed probability vectors. If not specified, a random scenario is used by default. Use this parameter to provide fixed probability vectors as a list with the following named elements: pE: Numeric vector of efficacy probabilities for each dose level. pT: Numeric vector of toxicity probabilities for each dose level. obd: Integer indicating the index of the true Optimal Biological Dose (OBD). mtd: Integer indicating the index of the true Maximum Tolerated Dose (MTD). For example: prob <- list( pE = c(0.4, 0.5, 0.6, 0.6, 0.6), pT = c(0.1, 0.2, 0.3, 0.4, 0.4), obd = 3, mtd = 2 )
save_dir	Directory to save output folders. Default is tempdir().
save_folder	Folder name. (Default is "boin12_simulations")
save_file	File name. (Default is "boin12_simulation.csv")

Value

No return value, called for side effects

Examples

simulate_tepi(
  ndose = 5,
  ssizerange = 1:2,
  target_t = 0.3,
  lower_e = 0.4,
  ntrial = 10,
)

phase12designs documentation built on Sept. 2, 2025, 9:09 a.m.