essCRM: Calculating ESS for time to event continual reassessment...

Description Usage Arguments Value Author(s) References Examples

Description

Calculates ESS for CRM for dose-finding in phase I clinical trials, based on a binary indicator of toxicity

Usage

1
essCRM(PI,prior,betaSD,target,mRange,numMC)

Arguments

PI

A vector of the true toxicity probabilites associated with the doses.

prior

A vector of initial guesses of toxicity probabilities associated with the doses. Must be of same length as PI.

betaSD

Standard deviation of the normal prior of the model parameter.

target

The target dose limiting toxicity (DLT) rate.

m

A positive integer specified as an maximum value in which ESS is searched.

nsim

umber of simulations for numerical approximation

obswin

The observation window with respect to which the maximum tolerated dose (MTD) is defined. Default is obswin=30.

rate

Patient arrival rate: Expected number of arrivals per observation window. Example: obswin=6 and rate=3 means expecting 3 patients arrive in 6 time units. Default is rate=2.

accrual

Patient accrual scheme. Default is accrual="poisson". Alternatively use accrual="fixed" whereby inter-patient arrival is fixed.

Value

ess

Prior effective sample size

Author(s)

Jaejoon Song <jjsong2@mdanderson.org>, Satoshi Morita <smorita@kuhp.kyoto-u.ac.jp >

References

Morita, S., Thall, P. F., and Muller, P. (2010). Evaluating the impact of prior assumptions in Bayesian biostatistics. Stat Biosci, 2, 1-17.

Morita, S., Thall, P. F., and Muller, P. (2008). Determining the effective sample size of a parametric prior. Biometrics, 64, 595-602.

O'Quigley J., Pepe M., Fisher, L. (1990).Continual reassessment method: A practical design for phase I clinical trials in cancer. Biometrics, 46, 33-48.

Examples

1
2
3
4
5
6
7
8
9
## Revisiting Example 4 in Morita et al. (2010, Stat Biosci).
## We assume five standardized doses to be (.02,.06,.10,.18,.30)
## For the first toxicity scenario given in Table 5 of Morita et al. (2010),
## we have (1mg/kg, 3mg/kg, 6mg/kg, 8mg/kg, 10mg/kg) = (.02,.06,.10,.18,.30)
## Assuming a normal prior for the beta with mean zero and variance of .5,
## we can compute the ESS as the following
essCRM(PI=c(.02,.06,.10,.18,.30),
        prior=c(.02,.06,.10,.18,.30),betaSD=sqrt(.5),
        target=0.2,m=7,nsim=200)

github-js/ess documentation built on May 16, 2019, 7:11 p.m.