sampsize <- c(100,300,500,1500) prop <- c(0.10,0.25,0.50) rval <- 2 myfile <- "OUTPUT" cat("WIDTH OF 95% CONFIDENCE INTERVAL FOR A GIVEN PROP/SAMPLE SIZE\n", file=myfile,append=F) cat(paste(date()),"\n\n",file=myfile,append=T) cat("N=",file=myfile,append=T) for (i in 1:length(prop)) { cat("\t",prop[i],sep="",file=myfile,append=T) } cat("\n\n",file=myfile,append=T) for (i in 1:length(sampsize)) { width <- 2*1.96*sqrt(prop*(1-prop)/sampsize[i]) cat(sampsize[i],": ",sep="",file=myfile,append=T) for (j in 1:length(width)) { cat(round(width[j],rval),"\t",sep="", file=myfile,append=T) } cat("\n",file=myfile,append=T) }