simulate_pite: Simulate operating characteristics using PRINTE

View source: R/simulate_pite.R

simulate_piteR Documentation

Simulate operating characteristics using PRINTE

Description

This function runs simulations of the PRINTE 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_pite(
  ndose,
  ssizerange,
  target_t,
  target_e,
  lower_e,
  cohortsize = 3,
  startdose = 1,
  eps1 = 0.05,
  eps2 = 0.05,
  psafe = 0.95,
  pfutility = 0.9,
  ntrial = 10000,
  utilitytype = 1,
  u1,
  u2,
  prob = NULL,
  save_dir = tempdir(),
  save_folder = "pite_simulations",
  save_file = "pite_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)

target_e

Numeric. Target efficacy 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)

eps1

Numerical. Width of the subrectangle.

eps2

Numerical. Width of the subreactangle.

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

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
)
simulate_pite(
  ndose = 5,
  ssizerange = 1:2,
  target_t = 0.3,
  target_e = 0.5,
  lower_e = 0.4,
  ntrial = 10,
  prob = prob,
)

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