Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
68 70 90 95 75 69 39 80 100 80 81 55 100 75 70 100 61 83 90 60 44 90 100 82 80 88 100 83 75 55 75 80 70 100 60 95 95 90 90 80 76 65 78 70 60 65 100 65 90 90 90 80 91 100 90 80 100 88 61 90 58 90 88 90 80 85 60 95 100 61 80 90 75 85 95 78 90 80 100 30 70 67 35 100 96 74 70 80 100 100 74 59 35 80 91 75 82 100 90 50 90 60 100 80 65 100 72 98 93 45 88 82 90 100 75 100 100 80 71 71 60 100 80 85 100 61 73 80 60 100 83 70 91 45 70 95 80 78 100 100 100 82 90 87 75 100
# 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