Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
4000 3600 1600 2500 3800 3500 4100 2000 3300 4000 4100 4200 1100 1900 4800 4000 3250 1200 3600 3400 3600 3300 1300 2500 3800 3800 3000 1600 2100 3900 4000 3900 2300 3600 3600 4200 4400 1900 2700 3600 4100 5250 1200 3000 4000 5200 2100 3400 5200 6900 2300 4100 4300 4500 4700 2200 2900 4400 5100 2300 3500 4100 4700 5400 2500 3500 3900 4100 4500 1800 4900 4500 3700 3100 1700 4000 4200 4300 3600 2600 4600 4000 2900 1800 3000 4200 4400 4900 2300 2400 3600 3900 2800 3800 4400 4200 5400 2300 3500 3800 4400 4400 2500 3700 3500 3300 3700 2200 3200 3100 4800 4400 2500 3100 4700 4300 4600 3100 3200 4300 4000 4500 2600 3600 4700 3800 4100 1700 3200 4000 3700 3600 2500 3800 4500 3500 4100 2300 3400 4800 4900 2600 3300 3900 3800 5000 1700 2700 4000 3300 4000 2500 3600 3800 5900 3600 3300 4000 4200 4300 2900 3900 5200 4000 4700 2200 2400 3700 3700 3000 2100 3500 6300 5200 7800 2700 3400 5000 5500 4500 2600 3900 3800 4400 4300 3300 3600 5300 3700 3800 4000 5200 5100 5100 5200 3800 6700 5200 4900 5700 2700 6300 5200 4500 6400 4000 3500 4700 4700 5000 2300 5200 5100 5400 4400 2600 4700 6100 6100 3500 4900 5500 5300 6300 2800 5900 5500 5200 6400 3100 5100 5400 5400 5100 2400 5500 5300 5000 4500 2000 3800 4900 5000 4000 1900 5200 5400 3200 800 3800 4300 2700
# 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