threep3: Calculate all possible trial pathways for the standard 3+3...

Description Usage Arguments Details Value Author(s) References See Also Examples

View source: R/bcrm_0.4.6.R

Description

All possible pathways of a standard 3+3 design are calculated and assigned a probability of occurring. This facilitates the calculation of operating characteristics, using print.threep3 and plot.threep3.

Usage

1
threep3(truep, start = 1, dose = NULL)

Arguments

truep

A vector of length k (the number of doses being considered in the trial), with values equal to the true probabilities of toxicity at the dose levels.

start

Starting dose level. Defaults to 1, i.e. the lowest dose level

dose

Optional vector of length k of actual doses for presentation purposes

Details

The dose-escalation schema used here relates to that defined by Chang et al. (2006), which incorporates dose de-escalation also. Variations of this design exist in the literature (see Storer 1989, Reiner et al. 1999).

The first cohort of three patients are administered the starting dose (usually the lowest dose). The trial then proceeds as follows:

Escalation / de-escalation rules to the next dose level for subsequent cohorts proceed as follows:

If none of the rules above are satisfied then the trial stops. If the current dose level has at most one DLT observed then this is claimed to be the MTD, otherwise the dose level below is deemed to be the MTD.

If dose-escalation extends to doses outside of that defined by dose, the MTD is determined to be the largest dose in dose.

Value

threep3 returns an object of class "threep3". The function print (i.e. print.threep3) can be used to obtain operating characteristics of the design used.

An object of class "threep3" is a list with the following components:

prob

A vector with the probabilities of each design occurring. As all possible designs are calculated, this vector sums to one

ssize

A vector with the sample size of each design

mtd

A vector of dose levels giving the recommended maximum tolerated dose (MTD) at the end of the trial

exp

A vector of length k giving the average trial experimentation proportions at each dose level

dlt.no

A vector with the number of toxicities (DLTs) that occur in each trial

truep

The true probabilities of toxicity at each dose level, specified by the user

dose

The actual doses as supplied in the function arguments

n.average

The average number of patients dosed at each level

dlt.average

The average number of DLTs experienced at each dose level

all.designs

A matrix containing all possible 3+3 designs, with each row representing a different design. Columns labelled "d k" and "tox k" represent the dose level and number of toxicities for the kth cohort, respectively.

Author(s)

Graham Wheeler [email protected] (MRC Biostatistics Unit, Cambridge, UK) and

Michael Sweeting [email protected] (University of Cambridge, UK)

References

Sweeting M., Mander A., Sabin T. bcrm: Bayesian Continual Reassessment Method Designs for Phase I Dose-Finding Trials. Journal of Statistical Software (2013) 54: 1–26. http://www.jstatsoft.org/article/view/v054i13

Chang A., Ganz P., Hayes D., Kinsella T., Pass H., Schiller J., Stone R., Strecher V. Oncology: An Evidence-Based Approach. Springer (2006).

Storer B. Design and Analysis of Phase I Clinical Trials. Biometrics (1989) 45: 925–937.

Reiner E., Paoletti X., O'Quigley J. Operating characteristics of the standard phase I clinical trial design. Computational Statistics & Data Analysis (1999) 30: 303–315.

Neuenschwander B., Branson M., Gsponer T. Critical aspects of the Bayesian approach to phase I cancer trials. Statistics in Medicine (2008) 27: 2420–2439.

See Also

threep3

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
## What are the operating characteristics of a standard 3+3 design if we conside only the first 
## 12 doses of the dose-escalation cancer trial example as described in Neuenschwander et al 2008.
## Pre-defined doses
dose<-c(1,2.5,5,10,15,20,25,30,40,50,75,100)
## Pre-specified probabilities of toxicity
p.tox0<-c(0.010,0.015,0.020,0.025,0.030,0.040,0.050,0.100,0.170,0.300,0.400,0.500)

## Not run: 
design.threep3<-threep3(p.tox0,dose)
print(design.threep3)
plot(design.threep3)

## End(Not run)

bcrm documentation built on May 29, 2017, 7:59 p.m.