R/plotPEstopVScE.R

In wiesenfa/BDP2: Bayesian Adaptive Designs for Phase II Trials with Binary Endpoint

plotPEstopVScE=function(n,vn.int,pF,cF,pE,cEvec,p0=NULL,p1=NULL,shape1F,shape2F,shape1E=NULL,shape2E=NULL,col0="red",col1="green",cex.legend=1, show=TRUE,...){
  args=as.list(match.call(expand.dots = FALSE))
  if (is.null(shape1E)) shape1E=shape1F
  if (is.null(shape2E)) shape2E=shape2F

  vn.int=vn.int[vn.int<n]
  args$vn.int=vn.int
  args$shape1E=shape1E
  args$shape2E=shape2E

  pEstop_p0 <-pFstopEstop_vec(n=n, p=p0, pE=pE, pF=pF, vn.int=vn.int, cEvec, cF,
                              shape1F=shape1F,shape2F=shape2F,shape1E=shape1E,shape2E=shape2E)
  pEstop_p1 <-pFstopEstop_vec(n=n, p=p1, pE=pE, pF=pF, vn.int=vn.int, cEvec, cF,
                              shape1F=shape1F,shape2F=shape2F,shape1E=shape1E,shape2E=shape2E)

  res=list(args=args,pEstop_p0=pEstop_p0,pEstop_p1=pEstop_p1,
                 x.p0=as.numeric(colnames(pEstop_p0$summary)),y.p0=pEstop_p0$summary[length(vn.int)+1,],
                 x.p1=as.numeric(colnames(pEstop_p1$summary)),y.p1=pEstop_p1$summary[length(vn.int)+1,]
            )

  if (show){
    plot(as.numeric(colnames(pEstop_p0$summary)),pEstop_p0$summary[length(vn.int)+1,],xlab=expression(c[E]),
         ylab="Cumulative Probability of Efficacy at Final",
    #     sub=paste("Analyses at",paste(c(vn.int,n),collapse=", ") ),
          sub=paste0("Interim analyses at (",paste(vn.int,collapse=", "),"), final analysis at ",n,", pF=",pF,", cF=",cF,", pE =",pE),
         ylim=c(0,1),pch=20,col=1,xlim=c(min(cEvec),max(cEvec)), type="n",las=1,...) #final
    lines(as.numeric(colnames(pEstop_p0$summary)),pEstop_p0$summary[length(vn.int)+1,],
          col=col0,lwd=2)
    lines(as.numeric(colnames(pEstop_p1$summary)),pEstop_p1$summary[length(vn.int)+1,],
          col=col1,lwd=2)
    legend("bottomleft",legend=c(paste("ptrue=",p0),paste("ptrue=",p1)),lty=1,col = c(col0,col1),cex=par()$cex*cex.legend)
  }
 
  class(res)="cE_vs_pEstop"
  invisible(res)
}


plot.cE_vs_pEstop=function(x,...){
    args=as.list(match.call(expand.dots = FALSE))
    do.call("plot",c(list(x=as.numeric(colnames(x$pEstop_p0$summary)),y=x$pEstop_p0$summary[length(x$args$vn.int)+1,],xlab=expression(c[E]),
         ylab="Cumulative Probability of Efficacy at Final",
    #     sub=paste("Analyses at",paste(c(vn.int,n),collapse=", ") ),
          sub=paste0("Interim analyses at (",paste(x$args$vn.int,collapse=", "),"), final analysis at ",x$args$n,", pF=",x$args$pF,", cF=",x$args$cF,", pE =",x$args$pE),
         ylim=c(0,1),pch=20,col=1,xlim=c(min(x$args$cEvec),max(x$args$cEvec)), type="n",las=1),x$args$...,args$...))

    lines(as.numeric(colnames(x$pEstop_p0$summary)),x$pEstop_p0$summary[length(x$args$vn.int)+1,],
          col=x$args$col0,lwd=2)
    lines(as.numeric(colnames(x$pEstop_p1$summary)),x$pEstop_p1$summary[length(x$args$vn.int)+1,],
          col=x$args$col1,lwd=2)
    legend("bottomleft",legend=c(paste("ptrue=",x$args$p0),paste("ptrue=",x$args$p1)),lty=1,col = c(x$args$col0,x$args$col1),cex=par()$cex*x$args$cex.legend)
  
}