Finding the doses of next cohort for a phase I clinical trial based on Escalation with Overdose Control (EWOC) design considering the classic parametrization for binary response and two agents.
ewoc2(dose.a, dose.b, resp, theta, alpha, Min.Dose.A, Max.Dose.A, Min.Dose.B, Max.Dose.B,
a01, b01, a10, b10, a00, b00, a, b, delta1x, delta1y, burn, mm, delta1)# S3 method for default
ewoc2(dose.a, dose.b, resp, theta, alpha, Min.Dose.A, Max.Dose.A, Min.Dose.B, Max.Dose.B,
a01, b01, a10, b10, a00, b00, a, b, delta1x, delta1y, burn=4000, mm=2000, delta1=0.05)
a data frame containing the current doses and responses set
list of input parameters
list of prior parameters
the next recommended doses for drug A
the next recommended doses for drug B
a numeric vector of allowable doses for drug A
a numeric vector of allowable doses for drug B
a numeric vector of allowable responses, 0 or 1
a numeric value defining the proportion of expectd patients to experience a medically unacceptable, dose-limiting toxicity (DLT) if administered the MTD.
a numerical value defining the probability that dose selected by EWOC is higher than the MTD.
a numeric value defining the lower bound of the support of the MTD for drug A
a numeric value defining the upper bound of the support of the MTD for drug A
a numeric value defining the lower bound of the support of the MTD for drug B
a numeric value defining the upper bound of the support of the MTD for drug B
a numeric value for beta prior distribution associated with parameter rho01
a numeric value for beta prior distribution associated with parameter rho01
a numeric value for beta prior distribution associated with parameter rho10
a numeric value for beta prior distribution associated with parameter rho10
a numeric value for beta prior distribution associated with parameter rho00
a numeric value for beta prior distribution associated with parameter rho00
a numeric value for gamma prior distribution associated with parameter eta
a numeric value for the gamma prior distribution associated with the parameter eta
Maximum dose escalation at each step for drug A, the default is 0.2*(Max.Dose.A-Min.Dose.A if not assigned)
Maximum dose escalation at each step for drug B, the default is 0.2*(Max.Dose.B-Min.Dose.B if not assigned)
Number of iterations for adaption, see n.adapt in jags.model for detail
Number of iterations to monitor, see n.iter in code.samples for detail
Threshold for toxicity
Tighiouart M, Li Q and Rogatko A. A Bayesian adaptive design for estimating the maximuym tolerated dose curve using drug combinations in cancer phase I clinical trials. Statistics in Medicine. 2017, 36: 280-290.
test = ewoc2(dose.a=c(0,0),dose.b=c(0,0),resp=c(0,0),theta=0.33,alpha=0.25,
Min.Dose.A=0, Max.Dose.A=1, Min.Dose.B=0, Max.Dose.B=1,a01=1,b01=1,a10=1,b10=1,
a00=1,b00=1,a=0.8,b=0.0384)
print(test)
Run the code above in your browser using DataLab