# For 90 percent power (pow), a call to gsearlyModel provides a feasible design
fp <- c(0.0000,0.0010,0.0250)
tn <- c(0.4800,0.7200,0.9750)
rctdesign <- gsearlyModel(rmodel="fix", trecruit=36, s=3, tfu=c(3,6,12),
tinterims=c(16,31), pow=0.8, cmodel="exponential", sd=20,
rho=0.75, theta=8, fp=fp, tn=tn)
fixedSampsize(rctdesign)
Run the code above in your browser using DataLab