R/LR.stats.SM.R

In CLAST: Exact Confidence Limits after a Sequential Trial

#' @export

LR.stats.SM <-
function(obj,alpha=0.05,type="upper"){
    #  Calculates all possible upper limits for p from an object
    #  that contains all possible data sets from a group sequential trial.
    lims=NULL
    for(h in 1:length(obj$S)){
        if(type=="upper"){
            lims=c(lims,LR.upper(obj$S[h],
              sum(obj$design[,1][1:obj$M[h]]),a=alpha))}
        if(type=="lower"){
            lims=c(lims,LR.lower(obj$S[h],
              sum(obj$design[,1][1:obj$M[h]]),a=alpha))}
    }
    out=obj
    out$lims=signif(lims,6)
    out
}