vtSingleTrial: Simulate a single trial

View source: R/visit_simu.R

vtSingleTrialR Documentation

Simulate a single trial

Description

Simulation function for simulating a single trial

Usage

vtSingleTrial(trueps, size.cohort = 3, size.level = NULL,
  etas = c(0.1, 0.3), dec.cut = 0.65, prob.mdl = c("NONPARA",
  "NONPARA+", "PARA", "PARA+"), priors = NULL, ...)

Arguments

trueps

True \theta's. A VTTRUEPS object made from vtScenario

size.cohort

Size of each cohort

size.level

Maximum number of patients for each dose level

etas

Vector of length 2 representing (p_L, p_U). p_L: lower bound of DLT risk, below which the current dose is considered absolutely safe; p_U: upper bound of DLT risk above which the current dose is considered too toxic

dec.cut

Thresholds C_1,C_2,C_3. If the vector length is shorter than 3, it is repeated to have 3 elements. See visit for details.

prob.mdl

Option of the probability models:

  • NONPARA: non-parametric+ model

  • NONPARA+: non-parametric model

  • PARA: partially parametric model

  • PARA+: partially parametric+ model

Default value is NONPARA. See visit for details.

priors

A class VTPRIOR object created by vtPriorPar for PARA and PARA+ model.

...

Optional arguments for vtPost

Value

  • dose: Optimal dose level

  • n.patients: Number of patients for each dose level and each cohort

  • ptox: Posterior mean of DLT risk rate after each interim analysis

  • pres: Posterior mean of immune response rate after each interim analysis

  • region: Identified region in the decision map after each interim analysis

  • prob: Posterior mean of \theta's after each interim analysis

  • smps: Observed data after each cohort

Examples

rst.sce <- vtScenario(tox = c(0.05, 0.05, 0.08),
                      res = c(0.2, 0.3, 0.5),
                      rho = 1)
rst.simu  <- vtSingleTrial(trueps = rst.sce, size.cohort=3, size.level=12,
                           prob.mdl="NONPARA");


visit documentation built on Aug. 9, 2023, 5:08 p.m.