Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
20 25 15 15 25 25 25 21 30 25 20 40 13 30 25 20 25 20 25 20 20 15 15 12 20 5 20 15 25 22 20 22 25 20 20 35 30 25 20 20 20 25 25 15 20 35 25 25 30 23 10 22 25 25 22 30 20 25 25 22 25 25 25 22 25 12 18 20 20 22 30 25 22 20 50 30 25 20 30 22 25 30 22 25 22 22 25 25 25 20 22 15 20 30 20 25 30 35 22 12 30 15 10 30 9 25 20 20 35 25 35 30 12 25 15 25 25 20 20 6 15 40 20 40 25 25 20 15 15 22 24 22 20 25 25 25 35 40 20 22 22 20 25 25 18 25 20 25 30 20 22 35 22 25 25 25 25 22 23 35 15 25 18 22 25 25 28 30 20 25 25 30 22 30 10 10 25 20 22 25 25 15 22 25 25 28 22 30 25 20 25 25 20 30 20 30 50 19 20 28 20 25 35 25 25 15 16 20 20 25 30 20 25 25 25 20 20 25 25 30 22 20 25 25 18 18 20 25 25 30 25 20 25 20 20 20 22 18 22 20 15 25 25 20 25 15 22 25 25 15 12 25 30 22 15 22 25 12 18 30 25 25 40 24 25 15 25 20 25 25 25 20 30 20 25 30 22 25 25 25 50 19 50 25 35 20 20 20 20 20 25 25 25 20 20 20 20 25 18 25 22 22 30 30 8 20 25 30 50 22 20 10 25 25 25 25 18 25 20 25 30 18 20 25 22 22 20 20 25 20 20 20 20 25 20 10 20 25 30 25 50 30 30 50 15 25 25 22 20 22 30 25 18 22 22 30 40 25 20 10 20 9 15 20 15 20 30 12 15 12 20 15 12 25 20 25 25 25 30 20 25 15 15 22 10 15 10 20 25 20 20 38 20 20 20 40 25 25 30 25 10 20 25 12 15 25 20 22 22 20 25 25 25 15 40 20 20 16 25 15 20 25 20 30 50 20 25 20 30 30 25 25 12 25 25 25 20 20 20 15 20 25 15 25 50 30 20 20 25 12 15 20 20 35 22 15 18 30 22 12 12 20 20 15 25 15 20 20 25 18 30 20 25 25 25 20 20 25 20 22 15 15 22 20 10 25 20 20 15 12 20 5 20 15 15 25 25 25 15 25 22 25 20 18 22 25 35 25 25 25 35 30 22 30 50 15 25 24 20 25 25 25 12 15 22 25 25 25 25 15 20 20 15 35 30 20 22 65 20 25 22 20 25 25 20 25 15 20 12 15 10 25 15 30 35 25 25 25 25 25 40 40 25 25 20 25 25 22 25 30 25 25 30 25 25 30 25 25 20 22 22 20 25 22 25 22 40 25 25 25 22 20 35 20 35 25 22 25 25 25 25 25 40 25 30 25 20 25 25 30 22 22 20 15 15 25 25 20 20 15 25 15 20 22 25 15 15 18 5 15 25 18 40 25 25 20 30 20 25 25 25 22 22 25 25 30 25 25 25 25 20 20 25 25 25 25 20 30 25 22 30 20 20 30 25 25 30 20 25 25 24 25 30 18 15 22 22 25 22 22 25 15 20 22 18 35 20 20 20 25 25 30 15 25 22 26 25 20 25 25 25 22 25 25 20 22 30 15 30 25 20 25 25 35 22 20 25 20 20 18 20 22 25 10 20 25 20 20 30 25 20 15 20 25 10 20 25 22 22 25 25 15 25 20 10 25 16 25 35 25 15 25 25 30 25 10 22 20 25 20 20 25 22 18 30 19 25 20 25 20 25 20 22 12 30 12 22 25 25 25 25 30 30 10 22 22 25 20 22 20 25 20 15 25 20 25 20 30 15 40 25 20 22 22 30 20 40 20 25 20 25 20 50 50 25 25 40 30 22 30 20 25 25 30 25 25 20 18 18 28 25 22 15 40 40 12 12 18 12 25 26 18 25 22 15 25 15 15 15 25 15 12 22 20 20 25 20 12 9 15 12 15 25 20 20 15 15 30 21 25 22 22 50 15 25 15 25 22 18 50 20 50 20 20 30 25 20 22 25 50 40 25 25 25 25 30 40 25 30 20
# simulations
Significant digits
Bandwidth
(?)
Quantiles
P1 P5 Q1 Q3 P95 P99
P0.5 P2.5 Q1 Q3 P97.5 P99.5
P10 P20 Q1 Q3 P80 P90
Chart options
R Code
par1 <- as.numeric(par1) if (par1 < 10) par1 = 10 if (par1 > 5000) par1 = 5000 library(lattice) library(boot) boot.stat <- function(s,i) { s.mean <- mean(s[i]) s.median <- median(s[i]) s.midrange <- (max(s[i]) + min(s[i])) / 2 c(s.mean, s.median, s.midrange) } (r <- boot(x,boot.stat, R=par1, stype='i')) 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 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