Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
600 425 398 582 458 455 621 635 589 220 351 379 683 524 536 598 581 632 645 722 689 645 354 486 423 479 684 601 608 463 602 485 563 645 486 435 479 579 563 202 389 467 466 706 546 689 531 528 579 684 651 637 548 496 582 467 693 615 708 648 899 852 745 689 582 674 684 542 489 472 398 486 549 766 654 628 689 648 578 536 548 496 475 687 642 584 596 609 678 694 485 489 537 706 489 598
# 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
1 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation