R/samonalpha0Tables.R

In samon: Sensitivity Analysis for Missing Data

alpha0Tables <- function( alpha0TableResultsObj, trt1lab = "Treatment 1", trt2lab = "Treatment 2" ) {

    Obsmeans <- alpha0TableResultsObj$Obsmeans

    ltrt1lab <- nchar(trt1lab)
    ltrt2lab <- nchar(trt2lab)

    lastp    <- round( (15 + 39 + ltrt1lab) / 2, 1)
    lastp2   <- round( ( 9 + 39 + ltrt2lab) / 2, 1)
    
    fmt      <- sprintf("%%%2.0fs %%%2.0fs",lastp,lastp2)
    
    dash    <- paste(rep("-",72), sep="", collapse="")
    headerT <- "                      Treatment 1                  Treatment 2        "
    headerT <- sprintf(fmt,trt1lab,trt2lab)
    header  <- "    Time        N       Mean         SD        N       Mean         SD"

    Text1 <- "Number of values available at each time-point, mean and standard"
    Text2 <- "deviation of outcome."
    cat("\n\n")
    cat(paste(Text1,  "\n"))
    cat(paste(Text2,  "\n"))
    cat(paste(dash,   "\n"))
    cat(paste(headerT,"\n"))
    cat(paste(header, "\n"))
    cat(paste(dash,   "\n"))

    for ( i in 1:nrow(Obsmeans) )  {
       line <-  sprintf("%8.0f %8.0f %10.2f %10.2f %8.0f %10.2f %10.2f\n", Obsmeans[i,1], Obsmeans[i,2], Obsmeans[i,3], Obsmeans[i,4], Obsmeans[i,5], Obsmeans[i,6], Obsmeans[i,7])
       cat(line)
    }    
    cat(paste(dash,  "\n\n\n\n"))


    # Table 2 obs
    ObsResults <- alpha0TableResultsObj$ObsResults
    ObsTtest   <- alpha0TableResultsObj$ObsTtest

    dash    <- paste(rep("-",75), sep="", collapse="")
    header  <- "Treatment                 N       Mean         SD    t-statistic    p-value"
    Text1 <- "At final time-point (observed)"
    cat(paste(Text1,  "\n"))
    cat(paste(dash,   "\n"))
    cat(paste(header, "\n"))
    cat(paste(dash,   "\n"))
    
    line    <- sprintf("%-16s %10.0f %10.2f %10.2f\n",trt1lab, ObsResults[1,"NComplete"],ObsResults[1,"meanComplete"],ObsResults[1,"sdComplete"])
    cat(line) 
    line    <- sprintf("%-16s %10.0f %10.2f %10.2f\n",trt2lab, ObsResults[2,"NComplete"],ObsResults[2,"meanComplete"],ObsResults[2,"sdComplete"])
    cat(line)

    cat("\n")
    line    <- sprintf("Difference (2-1)  %20.2f  %24.3f %10.3f\n", ObsTtest[1,"Difference"],ObsTtest[1,"ttestStatistic"],ObsTtest[1,"Pvalue"])
    cat(line)
    
    cat(paste(dash,   "\n\n\n"))



    # Table 3 LOCF
    LOCFResults <- alpha0TableResultsObj$LOCF
    LOCFTtest   <- alpha0TableResultsObj$LOCFTtest

    dash    <- paste(rep("-",75), sep="", collapse="")
    header  <- "Treatment                 N       Mean         SD    t-statistic    p-value"
    Text1 <- "At final time-point (LOCF)"
    cat(paste(Text1,  "\n"))
    cat(paste(dash,   "\n"))
    cat(paste(header, "\n"))
    cat(paste(dash,   "\n"))
    
    line    <- sprintf("%-16s %10.0f %10.2f %10.2f\n",trt1lab, LOCFResults[1,"N"],LOCFResults[1,"mean"],LOCFResults[1,"sd"])
    cat(line) 
    line    <- sprintf("%-16s %10.0f %10.2f %10.2f\n",trt2lab, LOCFResults[2,"N"],LOCFResults[2,"mean"],LOCFResults[2,"sd"])
    cat(line)

    cat("\n")
    line    <- sprintf("Difference (2-1)  %20.2f  %24.3f %10.3f\n", LOCFTtest[1,"Difference"],LOCFTtest[1,"ttestStatistic"],LOCFTtest[1,"Pvalue"])
    cat(line)
    
    cat(paste(dash,   "\n\n\n"))


    # Table 4 IF results
    IFResults <- alpha0TableResultsObj$IFResults

    CIlevel  <- 100*alpha0TableResultsObj$CIlevel
    if ( round(CIlevel) == CIlevel ) {
      headerT1 <- sprintf("                           IF            Bootstrap %2d%%           Jackknife %2d%%", CIlevel, CIlevel)
    } else {
      headerT1 <- sprintf("                           IF          Bootstrap %4.1f%%         Jackknife %4.1f%%", CIlevel, CIlevel)
    }
        
    dash     <- paste(rep("-",78), sep="", collapse="")
    headerT0 <- "                                                                Bootstrap  and"
    headerT  <- "                                         Bootstrap                Jackknife"
    header   <- "Treatment               Estimate        LCL        UCL         LCL         LCL"
    Text1 <- "At final time-point alpha = 0."
    cat(paste(Text1,   "\n"))
    cat(paste(dash,    "\n"))
    cat(paste(headerT0,"\n"))
    cat(paste(headerT1,"\n"))
    cat(paste(header,  "\n"))
    cat(paste(dash,    "\n"))
    
    line    <- sprintf("%-16s %15.3f %10.3f %10.3f %11.3f %11.3f\n",trt1lab, IFResults[1,"Estimate"],IFResults[1,"LCLBootstrap"],IFResults[1,"UCLBootstrap"],IFResults[1,"LCLJK"],IFResults[1,"UCLJK"])
    cat(line) 
    line    <- sprintf("%-16s %15.3f %10.3f %10.3f %11.3f %11.3f\n",trt2lab, IFResults[2,"Estimate"],IFResults[2,"LCLBootstrap"],IFResults[2,"UCLBootstrap"],IFResults[2,"LCLJK"],IFResults[2,"UCLJK"])
    cat(line)

    cat("\n")
    line    <- sprintf("Difference (2-1)  %14.3f %10.3f %10.3f %11.3f %11.3f\n", IFResults[3,"Estimate"],IFResults[3,"LCLBootstrap"],IFResults[3,"UCLBootstrap"],IFResults[3,"LCLJK"],IFResults[3,"UCLJK"])
    cat(line)
    
    cat(paste(dash,   "\n\n\n"))



    # Table 5 Optimal sigmaH and sigmaF
    optRes <- alpha0TableResultsObj$optRes

    dash    <- paste(rep("-",39), sep="", collapse="")
    header  <- "Treatment            SigmaH     SigmaF"
    Text1 <- "Smoothing parameters:"
    cat(paste(Text1,  "\n"))
    cat(paste(dash,   "\n"))
    cat(paste(header, "\n"))
    cat(paste(dash,   "\n"))
    
    line    <- sprintf("%-16s %10.3f %10.3f\n",trt1lab, optRes[1,"SigmaH"],optRes[1,"SigmaF"])
    cat(line) 
    line    <- sprintf("%-16s %10.3f %10.3f\n",trt2lab, optRes[2,"SigmaH"],optRes[2,"SigmaF"])
    cat(line)
    
    cat(paste(dash,   "\n\n"))
}