Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
17681170.13 19858875.65 16997477.78 16969881.42 18908869.11 15692144.17 15159951.13 15806842.33 16007123.26 16059123.25 16189383.57 12522497.94 14733828.17 15686348.94 13779681.49 14423755.53 15290621.44 14308336.46 13855616.12 14384486.00 15638580.70 19711553.27 20359793.22 16141449.72 20056915.06 20605542.58 19325754.03 20547653.75 19211178.55 19009453.59 18746813.27 16471529.53 18957217.20 20515191.95 18374420.60 16192909.22 18147463.68 19301440.71 18344657.85 17183629.01 19629987.60 17167191.42 17428458.32 16016524.60 18466459.42 18406552.23 18174068.44 14851949.20 16260733.42 18329610.38 18003781.65 15903762.33 19554176.37 16554237.93 16198892.67 16571771.60 17535166.38 16198106.13 17487530.67 13768040.14
# 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) 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') bitmap(file='plot7.png') b <- boxplot(z,notch=TRUE,ylab='simulated values',main='Bootstrap Simulation - Central Tendency') grid() dev.off() b 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.end(a) table.save(a,file='mytable.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='mytable1.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