Description Usage Arguments Details Value Note Author(s) References See Also Examples
Construct an optimal two-stage or three-stage designs with a time-to-event endpoint evaluated at a pre-specified time (e.g., 6 months) comparing treatment versus either a historical control rate with possible stopping for futility (single-arm), or an active control arm with possible stopping for both futility and superiority (two-arm), after the end of Stage I utilizing time to event data. The design minimizes either the expected duration of accrual (EDA), expected sample size (ES), or the expected total study length (ETSL).
1 2 3 4 5 |
B.init |
A vector of user-specified time points (B1, ..., Bb) that determine a set of time intervals with uniform accrual. |
m.init |
The projected
number of patients that can be accrued within the time intervals
determined by |
alpha |
Type I error. |
beta |
Type II error. |
param |
Events should be defined as poor outcomes (e.g. death, progression). Computations
and reporting are based on the proportion without an event at a
pre-specified time, |
x |
Pre-specified time for the event-free endpoint (e.g., 1 year). |
target |
The expected duration of
accrual (EDA) is minimized with |
sf |
Spending function for |
num.arm |
Number of arms: a single-arm design with |
r |
Proportion of patients randomized to the treatment arm when |
num.stage |
Number of stages: a two-stage design with |
pause |
The pause in accrual following the scheduled times for interim
analyses. Data collected during the pause on the previously accrued
patients are included in the interim analysis conducted at the end of the
pause. Accrual resumes after the pause without interuption as if no pause had
occurred. Default is |
control |
An optional list of control settings. See
|
... |
No additional optional parameters are currently implemented. |
OptimDes
finds an two-stage or three-stage design with a time to event endpoint
evaluated at a pre-specified time with potential stopping after the first stage.
For single arm designs, it implements the Case and Morgan (2003) and Huang, Talukder and Thomas (2010) generalizaton of the Simon (1989) two-stage design for comparing a treatment to a known standard rate with possible stopping for futility at the interim.
For randomized two-arm comparative designs, it allows an early stopping for both
futility and superiority. The spending function for superiority can be chosen with
argument sf
.
The design minimizes either the expected duration of accrual (EDA), expected sample size (ES), or the expected total study length (ETSL).
The design calculations assume Weibull distributions for the event-free
endpoint in the treatment group, and for the (assumed known, "Null") control
distribution. The function weibPmatch
can be used to select
Weibull parameters that yield a target event-free rate at a
specified time. Estimation is based on the Kaplan-Meier or
Nelson-Aalen estimators evaluated at a target time (e.g., 1 year).
The full treatment and control distributions and the accrual
distribution affect power (and alpha level in some settings), see Case and Morgan (2003)).
Accrual rates are specified by the user. These rates can differ across time intervals specified by the user (this generalizes the results in Case and Morgan).
A package vignette as user manual can be found in
the /doc subdirectory of the OptInterim
package. It can be accessed from the
HTML help page for the package.
A list with components:
target |
The optimizaton target ("EDA" or "ETSL" or "ES"). |
sf |
The alpha spending function ("futility" or "OF" or "Pocock"). |
test |
A vector giving the type I error |
design |
A vector giving the number of study arms |
accrual |
A list containing the input vectors |
result |
A 5-element vector containing the expected duration of accrual (EDA), the expected total study length (ETSL), the expected sample size of the optimal design (ES), and the probability of stopping under the Null (pstopNULL) and Alternative hypotheses (pstopAlt). |
n |
A two (or three)-element vector containing the sample size for the interim analysis and the sample size if all stages are completed. |
stageTime |
A 3 (or 4)-element vector giving the times for the interim and final analyses, and maximum duration of accrual. Interim times are given for the beginning of any pause before the analysis occurs. |
boundary |
A vector giving the rejection cutpoints (see |
se |
A vector of length 4 (or 6) with the asymptotic standard errors at the iterim and final analysis under the null hypothesis, followed by the corresponding SEs under the alternative hypothesis. These SEs must be divided by the square root of sample size. |
u |
A two (or three)-element vector giving means of interim test statistics under H1. See detailed description. It is also used to compute conditional power. |
exposure |
The expected total exposure of patients at the time
of the planned interim analysis (including any accrual pause).
Patient exposure is truncated by both
the interim analysis time (including any pause) and the
target surival time (i.e., no
exposure after |
all.info |
A data frame containing the results for all of the evaluated sample sizes. |
single.stage |
A six-element vector giving the sample
size |
The algorithm will search for the optimal n
between the
sample size for a single-stage design and the
user specified maximum sample size sum(m.init)
.
When the length of B.init
or m.init
is 1, the accrual
rate is constant as in Lin et al. (1996), Case and Morgan
(2003).
Bo Huang <bo.huang@pfizer.com> and Neal Thomas <neal.thomas@pfizer.com>
Huang B., Talukder E. and Thomas N. (2010). Optimal two-stage Phase II designs with long-term endpoints. Statistics in Biopharmaceutical Research, 2, 51–61.
Case M. D. and Morgan T. M. (2003). Design of Phase II cancer trials evaluating survival probabilities. BMC Medical Research Methodology, 3, 7.
Lin D. Y., Shen L., Ying Z. and Breslow N. E. (1996). Group seqential designs for monitoring survival probabilities. Biometrics, 52, 1033–1042.
Simon R. (1989). Optimal two-stage designs for phase II clinical trials. Controlled Clinical Trials, 10, 1–10.
np.OptimDes
, print.OptimDes
,
plot.OptimDes
, weibPmatch
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 | ## Not run:
B.init <- c(1, 2, 3, 4, 5)
m.init <- c(15, 20, 25, 20, 15)
alpha <- 0.05
beta <- 0.1
param <- c(1, 1.09, 2, 1.40)
x <- 1
# H0: S0=0.40 H1: S1=0.60
object12 <- OptimDes(B.init,m.init,alpha,beta,param,x,target="EDA",
sf="futility",num.arm=1,num.stage=2,control=OptimDesControl(n.int=c(1,5),trace=TRUE))
print(object12)
m.init <- 4*c(15, 20, 25, 20, 15)
object2 <- OptimDes(B.init,m.init,alpha,beta,param,x,target="EDA",sf="futility",num.arm=2)
print(object2)
object23O <- OptimDes(B.init,m.init,alpha,beta,param,x,target="ETSL",sf="OF",
num.arm=2,num.stage=3,control=OptimDesControl(trace=TRUE,aboveMin=c(1.05,1.10)))
print(object3)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.