treatsel.sim {asd} | R Documentation |
ASD simulation for treatment selection
Description
Function treatsel.sim
runs simulations for a trial design that tests a number of experimental treatments against a single control treatment group in a seamless adaptive trial. Test treatments are compared to the control treatment using Dunnett's many-to-one testing procedure. An interim analysis is undertaken using an early outcome measure for each treatment (and control). A decision is made on which of the treatments to take forward, using a pre-defined selection rule. Data are simulated for the final outcome measure, and data from the interim and final analyses for the final outcome measure are combined together using either the inverse normal or Fisher combination test, and hypotheses tested at the selected level.
Usage
treatsel.sim(n=list(stage1=32,stage2=32),
effect=list(early=c(0,0,0),final=c(0,0,0)),
outcome=list(early="N",final="N"),nsim=1000,
corr=0,seed=12345678,select=0,epsilon=1,
weight=NULL,thresh=1,level=0.025,ptest=c(1),
method="invnorm",fu=FALSE,file = "")
Arguments
n |
List giving sample sizes for each treatment group at stage 1 (interim) and stage 2 (final) analyses |
effect |
List giving effect sizes for early and final outcomes |
outcome |
List giving outcome type for early and final outcomes; available options are “ |
nsim |
Number of simulations (maximum=10,000,000) |
corr |
Correlation between early and final outcomes |
seed |
Seed number |
select |
Selection rule type ( |
epsilon |
For |
weight |
Optional user set weight for combination test; default is to use those suggested by Jenkins et al. (2011) |
thresh |
For |
level |
Test level (default=0.025) |
ptest |
Vector of treatment numbers for determining power; for example, c(1,2) will count rejections of one or both hypotheses for testing treatments 1 and 2 against the control |
method |
Select combination method; available options are “ |
fu |
Logical indicating whether patients from dropped treatments (after interim selection) should be followed-up; default |
file |
File name to dump output; if unset will default to R console |
Details
A structured description of the the methodology and the simulation model is given by Friede et al. (2011) and implementation by Parsons et al. (2012).
Value
count.total |
Number of times one or more treatments are selected |
select.total |
Number of times each test treatment is selected |
reject.total |
Number of times each hypothesis is rejected |
sim.reject |
Number of times one or more of the treatments selected using |
Author(s)
Nick Parsons (nick.parsons@warwick.ac.uk)
References
Friede T, Parsons N, Stallard N, Todd S, Valdes Marquez E, Chataway J, Nicholas R. Designing a Seamless Phase II/III Clinical Trial using Early Outcomes for Treatment Selection: an Application in Multiple Sclerosis. Statistics in Medicine 2011;30:1528-1540.
Parsons N, Friede T, Todd S, Valdes Marquez E, Chataway J, Nicholas R, Stallard N. An R package for implementing simulations for seamless phase II/III clinicals trials using early outcomes for treatment selection. Computational Statistics and Data Analysis 2012;56:1150-1160.
Bretz F, Schmidli H, Koenig F, Racine A, Maurer W. Confirmatory seamless phase II/III clinical trials with hypotheses selection at interim: General concepts. Biometrical Journal 2006;48:623-634.
See Also
Examples
# two test treatment groups
# effect size = 0.3 for group 1
# for both early and final normal outcomes
# correlation = 0.3
# select one treatment only at interim
treatsel.sim(n=list(stage1=100,stage2=300),
effect=list(early=c(0,0.3,0),final=c(0,0.3,0)),
outcome=list(early="N",final="N"),
nsim=100,corr=0.3,seed=145514,select=1,
level=0.025,ptest=c(1,2),fu=FALSE,
method="invnorm",file="")
# five test treatment groups
# correlation = 0.3
# flexible selection rule, with epsilon = 1
treatsel.sim(n=list(stage1=100,stage2=300),
effect=list(early=c(0,0.3,0.2,0.1,0.3,0.05),
final=c(0,0.2,0.3,0.2,0.1,0.5)),
outcome=list(early="N",final="N"),
nsim=200,corr=0.3,seed=145514,select=4,epsilon=1,
level=0.025,ptest=c(1:5),method="invnorm")