R/print.sDesign.R

In binGroup: Evaluation and Experimental Design for Binomial Group Testing

"print.sDesign" <-
function(x,...)

{
args<-list(...)
if(is.null(args$digits))
 {digits<-4}
else{digits<-args$digits}

 if(x$alternative=="less")
  {alt.hyp <- paste("true proportion is less than",x$p.hyp )
   ptrue <- paste("assumed true proportion","=", x$p.hyp - x$delta)}

 if(x$alternative=="greater")
  {alt.hyp <- paste("true proportion is greater than", x$p.hyp )
   ptrue <- paste("assumed true proportion","=", x$p.hyp + x$delta)}

 if(x$alternative=="two.sided")
  {alt.hyp <- paste("true proportion is not equal to ",x$p.hyp )
   ptrue <- paste("assumed true proportion","=", x$p.hyp - x$delta,"or", x$p.hyp + x$delta)}


 if(x$power.reached==TRUE && x$bias.reached==FALSE)
  {cat("Power was reached without violating bias restriction","\n")
   cat("for s","=",x$sout,"with power","=",signif(x$powerout, digits),"\n")
   cat("and bias","=",signif(x$biasout, digits),"\n")

   cat("alternative hypothesis:", alt.hyp ,"\n")
   cat( ptrue, "\n")
  }

 if(x$power.reached==FALSE && x$bias.reached==FALSE && x$powerout!=0)
  {cat("Power was not reached in the range of s","=",min(x$sit),",",max(x$sit),"\n")
   cat("Maximal power was reached for s","=",x$sout,"with power","=",signif(x$powerout, digits),"\n")
   cat("and bias","=",x$biasout,"\n")

   cat("alternative hypothesis:", alt.hyp ,"\n")
   cat( ptrue, "\n")
  }

 if(x$power.reached==FALSE && x$bias.reached==TRUE)
  {cat("Power can not be reached without violating bias restriction","\n")
   cat("Maximal power without violating biasrest","=", x$biasrest,"\n")
   cat("was reached for s","=",x$sout,"\n")
   cat("with power","=",signif(x$powerout, digits),"and bias","=",signif(x$biasout, digits),"\n")

   cat("alternative hypothesis:", alt.hyp ,"\n")
   cat( ptrue, "\n")
  }

 if(x$power.reached==FALSE && x$bias.reached==FALSE && x$powerout==0)
  {cat("Power can not be reached in the range of s","=",min(x$sit),",",max(x$sit),"\n")
   cat("null hypothesis can not be rejected for any group size in this range\n")

   cat("alternative hypothesis:", alt.hyp ,"\n")
   cat( ptrue, "\n")
  }

invisible(x)

}