x <- c(14,18,11,12,16,18,14,14,15,15,17,19,10,16,18,14,14,17,14,16,18,11,14,12,17,9,16,14,15,11,16,13,17,15,14,16,9,15,17,13,15,16,16,12,12,11,15,15,17,13,16,14,11,12,12,15,16,15,12,12,8,13,11,14,15,10,11,12,15,15,14,16,15,15,13,12,17,13,15,13,15,16,15,16,15,14,15,14,13,7,17,13,15,14,13,16,12,14,17,15,17,12,16,11,15,9,16,15,10,10,15,11,13,14,18,16,14,14,14,14,12,14,15,15,15,13,17,17,19,15,13,9,15,15,15,16,11,14,11,15,13,15,16,14,15,16,16,11,12,9,16,13,16,12,9,13,13,14,19,13,12,13) par2 = '12' par1 = '500' par1 <- as.numeric(par1) par2 <- as.numeric(par2) if (par1 < 10) par1 = 10 if (par1 > 5000) par1 = 5000 if (par2 < 3) par2 = 3 if (par2 > length(x)) par2 = length(x) library(lattice) library(boot) boot.stat <- function(s) { s.mean <- mean(s) s.median <- median(s) c(s.mean, s.median) } (r <- tsboot(x, boot.stat, R=par1, l=12, sim='fixed')) z <- data.frame(cbind(r$t[,1],r$t[,2])) colnames(z) <- list('mean','median') postscript(file="/var/www/rcomp/tmp/13q6z1321112263.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) b <- boxplot(z,notch=TRUE,ylab='simulated values',main='Bootstrap Simulation - Central Tendency') grid() dev.off() b #Note: the /var/www/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab load(file="/var/www/rcomp/createtable") a<-table.start() a<-table.row.start(a) a<-table.element(a,'Estimation Results of Blocked Bootstrap',6,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'statistic',header=TRUE) a<-table.element(a,'Q1',header=TRUE) a<-table.element(a,'Estimate',header=TRUE) a<-table.element(a,'Q3',header=TRUE) a<-table.element(a,'S.D.',header=TRUE) a<-table.element(a,'IQR',header=TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'mean',header=TRUE) q1 <- quantile(r$t[,1],0.25)[[1]] q3 <- quantile(r$t[,1],0.75)[[1]] a<-table.element(a,q1) a<-table.element(a,r$t0[1]) a<-table.element(a,q3) a<-table.element(a,sqrt(var(r$t[,1]))) a<-table.element(a,q3-q1) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'median',header=TRUE) q1 <- quantile(r$t[,2],0.25)[[1]] q3 <- quantile(r$t[,2],0.75)[[1]] a<-table.element(a,q1) a<-table.element(a,r$t0[2]) a<-table.element(a,q3) a<-table.element(a,sqrt(var(r$t[,2]))) a<-table.element(a,q3-q1) a<-table.row.end(a) a<-table.end(a) table.save(a,file="/var/www/rcomp/tmp/23e4t1321112263.tab") a<-table.start() a<-table.row.start(a) a<-table.element(a,'95% Confidence Intervals',3,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'',1,TRUE) a<-table.element(a,'Mean',1,TRUE) a<-table.element(a,'Median',1,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Lower Bound',1,TRUE) a<-table.element(a,b$conf[1,1]) a<-table.element(a,b$conf[1,2]) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Upper Bound',1,TRUE) a<-table.element(a,b$conf[2,1]) a<-table.element(a,b$conf[2,2]) a<-table.row.end(a) a<-table.end(a) table.save(a,file="/var/www/rcomp/tmp/3daxc1321112263.tab") try(system("convert tmp/13q6z1321112263.ps tmp/13q6z1321112263.png",intern=TRUE))