Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
-7.34658641160771 -3.83393270509369 0.601079165691683 10.1708123315412 -19.3839869794366 76.9922425740188 5.53070349660019 -40.6642925799496 7.8923727519013 17.2902637570394 12.6760697817719 19.6464948843585 -22.8607017650403 24.2145566602915 -102.787405880591 -22.146815361688 19.5158080760987 -27.5133483228973 32.8225102202477 -23.6235769736771 9.83643314161895 4.60698806707138 70.5427891086114 50.22081902765 -0.610172080927745 109.512600045334 38.3844237614283 -17.1329178076636 9.8914651620462 8.93506884166259 -50.6089193933982 6.24830624992741 -15.4577854355838 17.4048085530475 -39.2380533193169 -41.6349321334386 5.69033922999012 31.4039155628394 20.6034212187875 4.54632743294802 -32.8317960857245 5.29772239563854 18.8811727936388 21.0417879652906 -34.8082213857578 32.3312855186524 -31.5422133188891 -6.96524384146907 -51.1887131351437 -19.6210349753855 16.3293278006868 -54.3249050913562 11.0791037716048 51.4231162089253 -17.4375959401962 -28.8303943807569 23.0112582539784 -6.31347225533266 -53.84077421835 -12.249523516243 13.1904099411244 8.39642203599641 -42.7811420681512 59.8084957547755 36.2483062499274 -29.8633361331903 38.3824605416578 76.0113792797905 45.9664338814935 -1.17612597762016 21.4146723322159 -26.7205970597838 -36.2752826830959 -1.78601535890883 -12.2530318706035 2.56265100567706 1.16142391932059 9.31058249798772 -27.9481685557842 -35.9069347340566 -16.4696314381812 7.28921972470831 -19.3587966655329 1.19328162042654 22.2468389243388 0.588096592448664 -19.1512627690298 -32.8384189427771 -19.9199397354674 28.5582080595397 -27.2895444904127 20.3171435437495 39.2813121694082 -9.73130430948931 -0.942375677488173 9.60548926754502 18.1714456892322 -4.93840200425281 -7.56232751201247 12.8399739909504 -30.4867554511851 -3.63900096482928 -17.8985321242085 30.4073651037041 -5.7272065943191 -9.33170040228832 -5.54115129593736 -14.7961174338129 -55.8251391362126 -7.45031070418756 7.46339499060845 3.41546646113856 -77.9853494545382 -70.9245602064598 25.0086309727327 29.2440128930041 10.2707157618025 -13.3739967262292 60.4465269458381 27.8033040934051 28.783673988851 -65.8301244767567 18.4628713819454 -89.4987213214184 -19.623252462839 -11.2502101051705 9.82813335318713 -20.7816766235179 40.6664595968392 -28.5945915143623 -10.6826220633439 26.0439497578208 -9.82613598785401 -7.37666572389216 -50.7149932897269 -40.8176172241267 -1.70087235956341 -99.8017389500321 14.1375276385751 -11.3700983552581 -11.5851921361629 43.8639647275284 32.0648275157806 1.17526571711362 15.5921920420822 5.30726906493567 0.309518089321935 -25.6242515771811 0.855347084597148 45.82596028482 34.6050928277554 -22.2974114549406 25.724966525478 38.2553520465631 -20.4191182276973 5.13598732690237 36.9150764336932 -22.4415575475537 17.9049156953637 -14.7084067380213 20.6994781446868 4.38462512248504 35.8516667812385 -124.97712092823 64.3234036524464 -31.3808056049636 22.9818734549562 9.22512141838141 7.16267812199563 -6.14015849546619 10.5281938816587 19.0310230212635 40.8664924012808 38.6731382905774 38.6731382905774 28.9974598982497 -6.40376507741122 23.1365593461027 56.3239790494825 -81.4622605682183 -6.58429493929693 26.0776190049686 19.861114006955 5.20928321135139 57.0212696008941 9.45682352090344 7.45682352090344 12.4568235209034 16.306500717659 -34.7239638045267 -27.969678691805 33.4925597898478 -10.2838683863397 22.2460928279007 -13.1050628287336 -18.4535558221218 9.35706193361561 3.31777871462932 -9.6207490651828 18.7281864038017 89.783673988851 -23.6846828893207 13.0052033197453 21.7361067500842 -3.67581318460958 -28.9608756213021 -76.9932529907464 -32.6912391121651 49.37539265306 17.3415441214018 -49.6868480326369 25.4616632245656 -21.8286948794475 -12.784343749941 14.6582683760197 -21.787055413815 9.07381431646217 48.8846677481931 -33.938148589021 10.7993040682602 13.2288911666416 -33.8784514839943 -23.4321103218372 7.38081929087432 39.5195070781681 -41.8079987923958 9.44010451918122 1.73576760864494 -11.368712656493 14.4103106601748 36.7926918156455 11.4415062241488 3.58294561440226 42.2508217717385 -29.784849680975 -79.348398518174 -47.2453527162769 22.0010904512135 -55.9351724842194 -1.49485326878918 -25.0235049393696 -28.1146659080171 15.240449060843 64.3234036524464 -1.8963136481531 29.7923853622617 31.3960625611379 -13.8095321274412 -7.75535772813713 10.663145044613 12.0411114253833 -54.9164187819607 5.60624912728185 -11.4532615656104 36.1205139843701 14.201605571223 -21.5211895057792 -47.5902501757404 11.8037424584266 25.0071469888737 -31.3808056049636 -9.64807162553253 7.74461357566997 15.240449060843 8.14127045873426 -17.663610396475 -23.8619452129315 22.111100502846 24.4276514774178 -35.2637422077927 25.4358785498882 -11.1428696590303 9.01713310814046 12.6078628048154 19.1899681873745 -22.1745442232816 21.1920012076042 -67.4087619189757
# 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) par3 <- as.numeric(par3) if (par1 < 10) par1 = 10 if (par1 > 5000) par1 = 5000 if (par2 < 3) par2 = 3 if (par2 > length(x)) par2 = length(x) library(modeest) library(lattice) library(boot) boot.stat <- function(s) { s.mean <- mean(s) s.median <- median(s) s.midrange <- (max(s) + min(s)) / 2 s.mode <- mlv(s,method='mfv')$M s.kernelmode <- mlv(s, method='kernel')$M c(s.mean, s.median, s.midrange, s.mode, s.kernelmode) } (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='plot7a.png') plot(r$t[,4],type='p',ylab='simulated values',main='Simulation of Mode') grid() dev.off() bitmap(file='plot8a.png') plot(r$t[,5],type='p',ylab='simulated values',main='Simulation of Mode of Kernel Density') 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],r$t[,4],r$t[,5]) ) colnames(z) <- list('mean','median','midrange','mode','mode.k.dens') bitmap(file='plot7.png') boxplot(z,notch=TRUE,ylab='simulated values',main='Bootstrap Simulation - Central Tendency') grid() dev.off() if (par4 == 'P1 P5 Q1 Q3 P95 P99') { myq.1 <- 0.01 myq.2 <- 0.05 myq.3 <- 0.95 myq.4 <- 0.99 myl.1 <- 'P1' myl.2 <- 'P5' myl.3 <- 'P95' myl.4 <- 'P99' } if (par4 == 'P0.5 P2.5 Q1 Q3 P97.5 P99.5') { myq.1 <- 0.005 myq.2 <- 0.025 myq.3 <- 0.975 myq.4 <- 0.995 myl.1 <- 'P0.5' myl.2 <- 'P2.5' myl.3 <- 'P97.5' myl.4 <- 'P99.5' } if (par4 == 'P10 P20 Q1 Q3 P80 P90') { myq.1 <- 0.10 myq.2 <- 0.20 myq.3 <- 0.80 myq.4 <- 0.90 myl.1 <- 'P10' myl.2 <- 'P20' myl.3 <- 'P80' myl.4 <- 'P90' } load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Estimation Results of Blocked Bootstrap',10,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'statistic',header=TRUE) a<-table.element(a,myl.1,header=TRUE) a<-table.element(a,myl.2,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,myl.3,header=TRUE) a<-table.element(a,myl.4,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]] p01 <- quantile(r$t[,1],myq.1)[[1]] p05 <- quantile(r$t[,1],myq.2)[[1]] p95 <- quantile(r$t[,1],myq.3)[[1]] p99 <- quantile(r$t[,1],myq.4)[[1]] a<-table.element(a,signif(p01,par3)) a<-table.element(a,signif(p05,par3)) a<-table.element(a,signif(q1,par3)) a<-table.element(a,signif(r$t0[1],par3)) a<-table.element(a,signif(q3,par3)) a<-table.element(a,signif(p95,par3)) a<-table.element(a,signif(p99,par3)) a<-table.element( a,signif( sqrt(var(r$t[,1])),par3 ) ) a<-table.element(a,signif(q3-q1,par3)) 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]] p01 <- quantile(r$t[,2],myq.1)[[1]] p05 <- quantile(r$t[,2],myq.2)[[1]] p95 <- quantile(r$t[,2],myq.3)[[1]] p99 <- quantile(r$t[,2],myq.4)[[1]] a<-table.element(a,signif(p01,par3)) a<-table.element(a,signif(p05,par3)) a<-table.element(a,signif(q1,par3)) a<-table.element(a,signif(r$t0[2],par3)) a<-table.element(a,signif(q3,par3)) a<-table.element(a,signif(p95,par3)) a<-table.element(a,signif(p99,par3)) a<-table.element(a,signif(sqrt(var(r$t[,2])),par3)) a<-table.element(a,signif(q3-q1,par3)) 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]] p01 <- quantile(r$t[,3],myq.1)[[1]] p05 <- quantile(r$t[,3],myq.2)[[1]] p95 <- quantile(r$t[,3],myq.3)[[1]] p99 <- quantile(r$t[,3],myq.4)[[1]] a<-table.element(a,signif(p01,par3)) a<-table.element(a,signif(p05,par3)) a<-table.element(a,signif(q1,par3)) a<-table.element(a,signif(r$t0[3],par3)) a<-table.element(a,signif(q3,par3)) a<-table.element(a,signif(p95,par3)) a<-table.element(a,signif(p99,par3)) a<-table.element(a,signif(sqrt(var(r$t[,3])),par3)) a<-table.element(a,signif(q3-q1,par3)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'mode',header=TRUE) q1 <- quantile(r$t[,4],0.25)[[1]] q3 <- quantile(r$t[,4],0.75)[[1]] p01 <- quantile(r$t[,4],myq.1)[[1]] p05 <- quantile(r$t[,4],myq.2)[[1]] p95 <- quantile(r$t[,4],myq.3)[[1]] p99 <- quantile(r$t[,4],myq.4)[[1]] a<-table.element(a,signif(p01,par3)) a<-table.element(a,signif(p05,par3)) a<-table.element(a,signif(q1,par3)) a<-table.element(a,signif(r$t0[4],par3)) a<-table.element(a,signif(q3,par3)) a<-table.element(a,signif(p95,par3)) a<-table.element(a,signif(p99,par3)) a<-table.element(a,signif(sqrt(var(r$t[,4])),par3)) a<-table.element(a,signif(q3-q1,par3)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'mode k.dens',header=TRUE) q1 <- quantile(r$t[,5],0.25)[[1]] q3 <- quantile(r$t[,5],0.75)[[1]] p01 <- quantile(r$t[,5],myq.1)[[1]] p05 <- quantile(r$t[,5],myq.2)[[1]] p95 <- quantile(r$t[,5],myq.3)[[1]] p99 <- quantile(r$t[,5],myq.4)[[1]] a<-table.element(a,signif(p01,par3)) a<-table.element(a,signif(p05,par3)) a<-table.element(a,signif(q1,par3)) a<-table.element(a,signif(r$t0[5],par3)) a<-table.element(a,signif(q3,par3)) a<-table.element(a,signif(p95,par3)) a<-table.element(a,signif(p99,par3)) a<-table.element(a,signif(sqrt(var(r$t[,5])),par3)) a<-table.element(a,signif(q3-q1,par3)) 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