R/print.boot.two.R

In wBoot: Bootstrap Methods

print.boot.two <-
function(x, ...)
{
test <- !is.null(x$Null)
hist(x$Boot.values,breaks=20,xlab=paste("bootstrap",x$Statistic),
 main=paste("Histogram of bootstrap ",x$Statistic,"s",sep=""))
abline(v=x$Observed,col="2")
abline(v=x$Mean,col="3")
abline(v=c(x$Confidence.limits),col="4")
if (test) abline(v=x$Null,col="5")
leg.text <- if (test) expression(Observed,Mean.boots,Confidence.interval,Null.value)
 else expression(Observed,Mean.boots,Confidence.interval)
legend("topright",leg.text,col=2:5,lwd=2,cex=.6)
cat("\n\n",x$Header,"\n\n")
if (x$Stacked || !is.null(x$Variable))
 {if (identical("proportion",x$Parameter))
   print(data.frame(SUMMARY="STATISTICS",Variable=x$Variable,
         Pop.1=x$Pop.1,Pop.2=x$Pop.2,n.1=x$n.1,
         x.1=x$Observed.1*x$n.1,n.2=x$n.2,x.2=x$Observed.2*x$n.2,
         Statistic=x$Statistic,Observed=x$Observed),row.names=FALSE) 
  else 
   print(data.frame(SUMMARY="STATISTICS",Variable=x$Variable,
         Pop.1=x$Pop.1,Pop.2=x$Pop.2,n.1=x$n.1,
         n.2=x$n.2,Statistic=x$Statistic,Observed=x$Observed),row.names=FALSE)}
else
 {if (identical("proportion",x$Parameter))
   print(data.frame(SUMMARY="STATISTICS",Pop.1=x$Pop.1,
         Pop.2=x$Pop.2,n.1=x$n.1, x.1=x$Observed.1*x$n.1,
         n.2=x$n.2,x.2=x$Observed.2*x$n.2,Statistic=x$Statistic,
         Observed=x$Observed),row.names=FALSE) 
  else
   print(data.frame(SUMMARY="STATISTICS",Pop.1=x$Pop.1,
         Pop.2=x$Pop.2,n.1=x$n.1,n.2=x$n.2,Statistic=x$Statistic,
         Observed=x$Observed),row.names=FALSE)} 
cat("\n")
print(data.frame(BOOTSTRAP="SUMMARY",Replications=x$Replications,Mean=x$Mean,
 SE=x$SE,Bias=x$Bias,Percent.bias=x$Percent.bias),row.names=FALSE)
cat("\n")
if (test) print(data.frame(HYPOTHESIS="TEST",Null=x$Null,
  Alternative=x$Alternative,P.value=x$P.value),row.names=FALSE)
if (test) cat("\n")
print(data.frame(CONFIDENCE="INTERVAL",Level=x$Level,Type=x$Type,
 Confidence.interval=x$Confidence.interval),row.names=FALSE)
cat("\n\n")
}