powered by
Calculates the probability to declare efficacy at the end of study given the Z value at the current look
curtail(lower.boundary,upper.boundary,look,t,noncen,current=lower.boundary[look])
numeric
lower boundaries
upper boundaries
current look number
time of looks
noncentrality parameter
current Z statistic
t<-c(0.33,0.67,1) f<- function(t) 0.025*t^4 g<-function(t) 0.20*t^3 a<-f(t) b<-g(t) noncen<-pnorm(0.975)+pnorm(0.8) curtail(b,a,1,t,noncen)
Run the code above in your browser using DataLab