simulate_efftox: Simulate operating characteristics using EffTox

View source: R/simulate_efftox.R

simulate_efftoxR Documentation

Simulate operating characteristics using EffTox

Description

This function runs simulations of the EffTox 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_efftox(
  ndose,
  ssizerange,
  target_t,
  lower_e,
  startdose = 1,
  ntrial = 10000,
  utilitytype = 1,
  prob = NULL,
  save_dir = tempdir(),
  save_folder = "efftox_simulations",
  save_file = "efftox_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)

startdose

Integer. Starting dose level. (Default is 1)

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)

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 of 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

prob <- list(
  pE = c(0.4, 0.5),
  pT = c(0.1, 0.2),
  obd = 2,
  mtd = 2
)
simulate_efftox(
  ndose = 2,
  ssizerange = 1,
  target_t = 0.3,
  lower_e = 0.4,
  ntrial = 2,
  prob = prob,
)

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