R/xtable.R

In gsDesign: Group Sequential Design

xtable.gsDesign <- function (x, caption = NULL, label=NULL, align=NULL, digits=NULL, display=NULL,
                             footnote = NULL, fnwid = "9cm", deltaname="delta", 
                             Nname="N",logdelta=FALSE, ...) 
{
  k <- x$k
  deltafutility <- gsDelta(x=x,i=1:x$k,z=x$lower$bound[1:x$k])
  deltaefficacy <- gsDelta(x=x,i=1:x$k,z=x$upper$bound[1:x$k])
  deltavals <- c(x$delta0,x$delta1)
  if (logdelta){
    deltafutility <- exp(deltafutility)
    deltaefficacy <- exp(deltaefficacy)
    deltavals <- exp(deltavals)
  }
  stat <- c("Z-value", "p (1-sided)",
            paste(deltaname,"at bound"),
            paste("P\\{Cross\\} if ", deltaname,"=", 
                  deltavals[1], sep = ""), 
            paste("P\\{Cross\\} if ", deltaname,"=", 
                  deltavals[2],sep = ""))
  st <- stat
  for (i in 2:k) stat <- c(stat, st)
  an <- array(" ", 5 * k)
  tim <- an
  enrol <- an
  fut <- an
  eff <- an
  an[5 * (0:(k - 1)) + 1] <- c(paste("IA ", 
      as.character(1:(k - 1)), ": ", as.character(round(100 * x$timing[1:(k - 1)],1)),
      "\\%", sep = ""), "Final analysis")
  an[5 * (1:(k - 1)) + 1] <- paste("\\hline", an[5 * (1:(k -1)) + 1])
  an[5 * (0:(k - 1)) + 2] <- paste(Nname,":", ceiling(x$n.I[1:k]))
  fut[5 * (0:(k - 1)) + 1] <- as.character(round(x$lower$bound,2))
  eff[5 * (0:(k - 1)) + 1] <- as.character(round(x$upper$bound,2))
  asp <- as.character(round(pnorm(-x$upper$bound), 4))
  asp[asp == "0"] <- "$< 0.0001$"
  eff[5 * (0:(k - 1)) + 2] <- asp
  bsp <- as.character(round(pnorm(-x$lower$bound), 4))
  bsp[bsp == "0"] <- " $< 0.0001$"
  fut[5 * (0:(k - 1)) + 2] <- bsp
  asp <- as.character(round(deltafutility[1:x$k], 4))
  fut[5 * (0:(k - 1)) + 3] <- asp
  bsp <- as.character(round(deltaefficacy[1:x$k], 4))
  eff[5 * (0:(k - 1)) + 3] <- bsp
  asp <- as.character(round(cumsum(x$upper$prob[, 1]), 4))
  asp[asp == "0"] <- "$< 0.0001$"
  eff[5 * (0:(k - 1)) + 4] <- asp
  bsp <- as.character(round(cumsum(x$lower$prob[, 1]), 5))
  bsp[bsp == "0"] <- "$< 0.0001$"
  fut[5 * (0:(k - 1)) + 4] <- bsp
  asp <- as.character(round(cumsum(x$upper$prob[, 2]), 4))
  asp[asp == "0"] <- "$< 0.0001$"
  eff[5 * (0:(k - 1)) + 5] <- asp
  bsp <- as.character(round(cumsum(x$lower$prob[, 2]), 4))
  bsp[bsp == "0"] <- "$< 0.0001$"
  fut[5 * (0:(k - 1)) + 5] <- bsp
  neff <- length(eff)
  if (!is.null(footnote)) 
    eff[neff] <- paste(eff[neff], "\\\\ \\hline \\multicolumn{4}{p{", 
                       fnwid, "}}{\\footnotesize", footnote, "}")
  x <- data.frame(cbind(an, stat, fut, eff))
  colnames(x) <- c("Analysis", "Value", "Futility", "Efficacy")
  xtable(x,caption = caption, label=label, align=align, digits=digits, display=display, ...)
}