Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
38.1 32.4 34.5 20.7 21.5 23.1 29.7 36.6 36.1 20.6 20.4 30.1 38.7 41.4 37 36 37 38 23 26.7 27.5 21.7 22.9 26.2 36.5 41.8 21.5 19.2 25 28.9 23.2 31.5 36.2 38.2 26.4 20.9 21.5 30.2 33.4 32.6 22.2 21.7 30 35.7 32.8 39.3 25.5 23 19.9 21.3 20.8 21.7 23.8 29 23.7 21.3 28.5 33.6 34.6 34.2 27 24.2 19.9 19.7 21.5 30.6 30 19 19.6 20.6 23.6 17.9 17.3 21.4 24.1 20.9 30.1 32.6 21.3 19.5 19.9 21 25.4 17.5 20.4 26.8 25.8 20.9 19.4 25.8 26.3 29.6 30.3 23.6 28.4 20.7 24.1 27.3 23.2 18.3 24.6 27.4 20.4 18.1 25.2 19.8 21 23.7 19.6 18.1 20.8 26 18.4 22 14.4 19.9 22.6 13.7 15.9 21.2 23.7 24 17.2 23.2 25.2 17.2 16 15.6 13.4 16 16.8 14.6 19.4 21 19.5 18.5 13.3 13.7 14.3 14.1 11.4 13.6 16.6 17.6 14.6 17.2 14.4 16.4 17.3 17.6 17.2 17.7 14.2 16.6 15.7 13.7 14.7 13.1 12.9 15.4 11.9 15.2 15.3 16.5 16.1 11.7 11.2 11.5 10.8 16.1 14.8 13.6 13.8 9.7 10.7 11 15.3 15.3 17 16 16.3 15.7 14.5 10.8 10.5 13.4 12.2 13.2 13 12.4 13.1 9.8 10.5 13.4 11 13.1 15 16.7 16.1 18.2 15.7 17.7 15.9 15.1 15.2 14.7 13.3 14.5 11.1 13.1 13.7 14.6 12.9 12.8 15.2 14.5 17.2 14.5 14.4 11 13.1 13.6 14.6 12.7 13.6 12.7 15.5 17.4 15.2 14.2 17.7 19.2 12.5 14.2 15.3 15.7 17 19 13.1 13.2 13.2 15.7 14.1 15.6 15.5 15.9 15.1 16 19.4 21.5 23.7 18.7 23.8 18 16.2 18.5 20.6 18.3 22.5 26.9 19.4 15.9 20.5 21.2 19.5 14.7 17.6 15.8 17.7 14.3 16.8 18.6 21.9 21.4 20.8 14 17 23 26.4 19.6 22.7 26.9 14.7 15.2 19.8 26.9 20.2 14.3 14.8 18.5 21.7 21.4 21.8 18.2 15.8 15.3 18.5 19.2 28.5 32.2 21.8 22.1 20.7 17 24.7 26.2 29 21.6 17.1 16.9 19.1 24.7 25.4 19.8 18.2 16.3 17 17.7 15.5 14.7 15.8 19.9 20.4 23.3 20.2 28.8 31.2 17.4 18.5 26.8 34.3 30.1 20.5 20.5 19.8 27 21 33 22.6 28.3 21.1 19 17.3 27 30.2 24.8 17.9 17.9 20.7 30.9 36.2 21 20.2 21.3 24.2 21 20.7 17.8 19.6 22.6 20.5 24.1 22.2
# simulations
blockwidth of bootstrap
Significant digits
Quantiles
P1 P5 Q1 Q3 P95 P99
P0.5 P2.5 Q1 Q3 P97.5 P99.5
P10 P20 Q1 Q3 P80 P90
bandwidth
Chart options
R Code
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) s.midrange <- (max(s) + min(s)) / 2 c(s.mean, s.median, s.midrange) } (r <- tsboot(x, boot.stat, R=par1, l=12, sim='fixed')) bitmap(file='plot1.png') plot(r$t[,1],type='p',ylab='simulated values',main='Simulation of Mean') grid() dev.off() bitmap(file='plot2.png') plot(r$t[,2],type='p',ylab='simulated values',main='Simulation of Median') grid() dev.off() bitmap(file='plot3.png') plot(r$t[,3],type='p',ylab='simulated values',main='Simulation of Midrange') grid() dev.off() bitmap(file='plot4.png') densityplot(~r$t[,1],col='black',main='Density Plot',xlab='mean') dev.off() bitmap(file='plot5.png') densityplot(~r$t[,2],col='black',main='Density Plot',xlab='median') dev.off() bitmap(file='plot6.png') densityplot(~r$t[,3],col='black',main='Density Plot',xlab='midrange') dev.off() z <- data.frame(cbind(r$t[,1],r$t[,2],r$t[,3])) colnames(z) <- list('mean','median','midrange') bitmap(file='plot7.png') boxplot(z,notch=TRUE,ylab='simulated values',main='Bootstrap Simulation - Central Tendency') grid() dev.off() load(file='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.row.start(a) a<-table.element(a,'midrange',header=TRUE) q1 <- quantile(r$t[,3],0.25)[[1]] q3 <- quantile(r$t[,3],0.75)[[1]] a<-table.element(a,q1) a<-table.element(a,r$t0[3]) a<-table.element(a,q3) a<-table.element(a,sqrt(var(r$t[,3]))) a<-table.element(a,q3-q1) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab')
Compute
Summary of computational transaction
Raw Input
view raw input (R code)
Raw Output
view raw output of R engine
Computing time
0 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation