Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
-6.53896045840745 -3.06308753452553 -1.32849034271104 12.7672516462708 -14.9743201582554 78.336943559759 3.92904820929143 -42.253386069434 13.1799431488537 22.0214909199079 16.2761904490883 16.8333441639479 -23.798907164877 24.0622554750583 -100.792081443483 -23.9387987784751 18.8029092011636 -22.402665234084 33.1899573065305 -26.6586786137629 9.49594583928348 7.41136738485014 73.6652455125498 46.5242209628917 2.4422206732717 109.918118397906 44.7374480719684 -19.4945832579308 8.3738069847269 5.35907814921374 -49.4609976475698 5.39965379972281 -11.8364675650002 13.6743740830477 -36.4058371904869 -44.9931456879041 7.30660022450749 27.5543976706867 21.2763860130414 7.08420115987175 -30.4821084759031 3.29753569931074 20.824058887433 19.0194478956298 -37.4385272759867 31.7086923845053 -33.6868495853794 -10.6969935514748 -47.6222918850906 -16.8912320745019 12.4312807086323 -57.0021828073131 6.34247853071848 51.2890990144928 -17.5785513400908 -28.6441318872035 20.4827418834763 -2.99928109345461 -54.7254894264013 -7.2017666187367 17.2218648484855 14.5889742271337 -46.7181178739351 55.5093354213615 35.3996537997228 -35.6582826661818 33.7164524886967 69.6433386856724 46.8594614657069 -1.44846675507214 26.7235779260356 -27.1846013672942 -35.5553685672257 -2.734553951598 -12.9266851423848 3.00299533808705 1.77708369499397 8.06402958110021 -28.9595571212668 -29.7849997308445 -15.7307815606041 6.39558429580351 -20.5935294019399 -0.379216755003139 22.7479873854765 4.13242089158884 -24.1769953661005 -35.8716172056208 -20.1947453333895 31.7400576737947 -25.1598390377148 18.3878338366768 42.2446740350395 -11.1798597918593 0.244674035039473 11.8420190681896 22.4376101474271 -6.45586571332611 -6.734553951598 10.6198051374991 -34.9238743790089 -7.65126138083469 -18.0003353469881 30.1106707707285 -5.35654189641016 -9.91383196188583 -2.28236296103412 -12.2302624936632 -59.082687753892 -6.86115825765499 5.67045540375544 2.56147272401332 -75.6335680883034 -74.6711284884266 23.7213914353847 24.1091161213783 16.4595565311831 -13.3835384976332 66.1636372859341 22.7611778830374 29.9076239062552 -65.806288897219 16.7573666341074 -88.2490542166487 -21.1421524044952 -10.8266997646441 9.94518430637842 -21.2385952686829 41.4550385256574 -26.9126339994821 -13.5365247895247 25.350989660276 -9.27615566880611 -7.52881783952696 -55.7669709180233 -37.6749397373565 -3.67112848842657 -100.473236200408 12.6942893668162 -13.2942901046896 -13.4828623145492 39.0883441080158 31.6860030142997 -0.438179080106497 17.1530973722577 4.72532476419147 2.68415555397467 -25.0603879918084 1.77250639940191 43.576328788517 32.9522055917256 -20.534999617201 28.7509457833513 38.0862088560365 -22.658793159882 3.56729166910799 37.1091161213783 -20.1843009524068 23.5579022406013 -16.2428916208819 23.6925059395032 5.31472073259509 35.6239638043875 -127.159490841835 64.9711941777751 -34.2895485292547 25.8491840624071 4.83016962811697 2.3024892955133 -7.08773867030167 13.829569193526 16.9418660616771 48.2834748612232 37.9745124224764 37.9745124224764 25.5307093764826 -8.79880231394266 22.3641010957643 54.7329307653852 -79.8289637558119 -7.8729236978651 24.5604055204227 23.7827771792706 6.57411739050205 55.3659780927542 12.2033475875511 10.2033475875511 15.2033475875511 19.5577066766482 -28.9483677004105 -25.5831223090569 29.9391146801201 -14.5965118911349 21.8840515001458 -9.07230584837406 -16.0358271075721 14.7410063042658 0.214871667868718 -11.6960147450345 20.5516380958622 90.9076239062552 -19.2489347987313 11.2076185368135 17.1684860312024 -3.63051475790319 -31.6073319850381 -81.7004566045245 -32.9958138116397 48.3908025132016 16.3841125286377 -49.5459395818405 23.7339850189187 -25.4016356454171 -16.0587343729138 19.0486484559133 -19.0815948577776 10.4969661066793 44.7330069114209 -32.6595891584261 15.9254553478232 13.2076185368135 -36.2614307118318 -17.0254875775237 10.9628437922466 36.7665231714338 -40.230298037348 6.03703297987713 -0.981305943016233 -10.8136878554729 14.2690042806994 35.1183093857636 5.87500956857675 2.41687769094315 40.0151261529997 -30.6683560747974 -78.6034133057458 -45.5697715849148 21.2869498119416 -56.3139969857003 -3.18131710873059 -25.6515136195455 -28.5393366221305 16.3373955911881 64.9711941777751 2.77123678004024 31.2977599681625 32.1498865579188 -12.0668088425645 -7.61395321942225 9.03001169452998 14.2450222309197 -57.2746752358085 7.2390498817699 -12.8246666780483 38.9745124224764 16.2798144369203 -19.1075491031359 -47.7473560491422 13.0049433289684 26.3827160955026 -34.2895485292547 -8.07454465217676 6.56378779813592 16.3373955911881 7.11146746821514 -25.5247898955711 -23.1992626399727 23.8506268676657 26.2420651321557 -34.3255837570543 25.0276181875614 -13.170430806474 6.92358984606289 14.2035431515042 18.2820784280881 -21.6517091834986 22.769701962652 -61.6807388446955
# 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