R/efficacy.R

In shrosehu/Continual-reassessment: Functions in this package help decide stopping rules in Phase 1 and Phase 2 clinical trials

#' Efficacy Stoping bound
#'
#' This function generats an efficacy stoping bound giving the parameters p and b.
#'
#' @param p0  null response rate
#' @param p1 target response rate
#' @param N maximum sample size
#' @param n1 the sample size to start interim monitoring
#' @param alpha alpha
#' @param beta beta
#' @importFrom foreach foreach %do%
#' @export
#' @examples efficacy(N=44, n1=1, p0=0.6, p1=0.8,alpha=0.05, beta=0.2)

efficacy <- function(N, n1, p0, p1,alpha,beta){
  e_tail <- rm.tails(N, p=p0, seq(0.5, 0.999, 0.001))

  ##design operating characteristics
  design_opers <- sapply(e_tail, function(x){
    e <- find.bound(tail=x, N, p=p0)

    last  <- min(N, max(which(diff(e) == 0)) + 1) #  last chance to cross e
    f<- rep(0,N)

    fun0 <- efcont(p=p0,e, f, n1=n1)
    err1 <- fun0$effstop[N]
    es0 <- ess(1-fun0$accum)
    if(err1< alpha){
    pwr <- efcont(p=p1, e, 1, n1=n1)$effstop[N]
    }else{
     pwr <- es0 <- NA
    }
    c(err1, pwr, es0, x)
  })

  alpha_vec <- design_opers[1,]
  pwr_vec <- design_opers[2,]
  es_vec <- design_opers[3,]

  out <- data.frame("alpha"=alpha_vec,
                    "power"=pwr_vec,
                    "EN"=es_vec,
                    "e_tail"=design_opers[4,])

  out <- out[which(pwr_vec>1-beta),]

  if(is.null(dim(out))) warning('no proper design')

  out
}