vtSingleTrial {visit}

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");