visit (version 2.2)

vtSingleTrial: 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, ...)

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

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

• 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

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