Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
2.80030819572237 0.546110929167911 -1.08882457056234 1.85116612407754 1.81431171737364 1.10336738027983 -1.1245171435247 1.71646640041695 0.793858512491152 -1.07059339259452 1.37785351123131 -0.79301277417615 -0.213022653692917 0.325551388183876 -3.07152677312816 -1.6593526905903 -0.977898156347797 2.77746716194583 0.659684775709612 -0.399824712085136 -1.94483744059989 2.63902152884315 0.830351272801263 0.239452118133438 -1.98023780133566 -3.14893729642629 1.52344843344618 1.52324864513037 3.07560476768955 1.68333908953913 -0.0411419750371334 -1.42576307236563 0.603469563697563 -2.02817409689518 2.36128930601495 -4.84900520812879 -0.296289875821015 -0.273406341522266 -2.41797545929271 -0.735276518673968 -0.758208236878704 0.719357808400705 3.82579134038112 0.751147004565213 -0.686160672385733 2.7839041091014 -0.353356298124068 -1.19565368260915 1.51622971016679 0.687899021959271 -1.03746150319703 0.232258589016444 0.954677970598071 0.253726601750312 0.0906111040418593 -0.337464765144114 -0.146645835344418 -2.14580429676713 -2.64273856198557 -1.69885035349378 0.0886900288640063 2.94885012445766 2.74086814209461 -2.78292060200689 2.11354456172898 1.5019821601947 0.928451976610637 0.251386956762454 1.92686717410181 -0.637152825120997 -3.00404371957499 1.53170453971748 -0.46034424964922 -0.435142320623821 0.241189107006538 -0.439858435973453 -1.66304983417195 1.54593083942837 0.0543382827538832 -2.76444539571847 -1.5519700760652 1.99785210877644 1.02671169946072 0.965014171087942 0.833645548583962 -1.78443597034853 -1.35462195476152 0.0905629201358726 1.2672859043733 1.3367497955578 1.00922794934105 2.97970594008853 -0.0198870686489594 1.37558503989304 1.62130116503375 -0.380035665253275 0.140633397218948 1.088271459054 -4.14041609113554 -1.36191185169048 -3.501894103396 1.82567198282555 1.00217218900356 -1.70508751233354 -1.54985320913675 1.52749755170791 -1.11909042570084 0.782613205637535 0.394327211581645 1.03045019427192 -1.98850357932067 -2.3911175541089 -2.49477174792859 1.43479131886814 -0.23693199243609 1.54252687780381 0.753916630559178 -1.12999945988543 0.434385612031916 0.497954292194177 -0.415657064566058 1.25212054898028 0.538408325374887 0.363650201264021 -0.525700021635337 -1.63043333745253 0.651092396769283 0.605424752110409 1.27467788471578 -3.66701110704255 -1.09455091524368 -2.01546316552593 -0.129017313491364 1.80505691986507 -2.93931661204688 0.266590824347748 1.08312554226064 -0.0657619162545335 1.90954638692403 -1.52044597892013 0.720915677264753 -1.30223881899939 1.12662982290344 -0.651780253148442 0.873322331496494 -0.545999882625613 -0.386205647008669 2.73997287136118 -0.246626845056549 2.49376474567805 -2.78811070614103 -0.470542681359586 1.71855478201569 1.18881581030952 2.56216432022585 0.443671331515407 -1.46466491183079 -2.05654394411988 1.31101180554756 -1.38802306664201 -1.42282760393946 -0.770545943306666 0.204099941303661 -3.73653059676392 2.41609709604484 0.459211471769727 -1.95214272484412 -0.467671553906553 -2.54700443680883 0.530260345025432 2.59956805328044 -0.048523881273836 -1.31831688661928 1.0921477226447 1.49992331623734 -1.21692016754482 -0.350713444316563 -0.709451653003092 0.660733517233603
# 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) par2 <- as.numeric(par2) if (par3 == '0') bw <- NULL if (par3 != '0') bw <- as.numeric(par3) if (par1 < 10) par1 = 10 if (par1 > 5000) par1 = 5000 library(modeest) 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 s.mode <- mlv(s[i], method='mfv')$M s.kernelmode <- mlv(s[i], method='kernel', bw=bw)$M c(s.mean, s.median, s.midrange, s.mode, s.kernelmode) } x<-na.omit(x) (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='plot7.png') plot(r$t[,4],type='p',ylab='simulated values',main='Simulation of Mode') grid() dev.off() bitmap(file='plot8.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() bitmap(file='plot9.png') densityplot(~r$t[,4],col='black',main='Density Plot',xlab='mode') dev.off() bitmap(file='plot10.png') densityplot(~r$t[,5],col='black',main='Density Plot',xlab='mode of kernel dens.') 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='plot11.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',10,TRUE) a<-table.row.end(a) 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' } 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,par2)) a<-table.element(a,signif(p05,par2)) a<-table.element(a,signif(q1,par2)) a<-table.element(a,signif(r$t0[1],par2)) a<-table.element(a,signif(q3,par2)) a<-table.element(a,signif(p95,par2)) a<-table.element(a,signif(p99,par2)) a<-table.element( a,signif( sqrt(var(r$t[,1])),par2 ) ) a<-table.element(a,signif(q3-q1,par2)) 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,par2)) a<-table.element(a,signif(p05,par2)) a<-table.element(a,signif(q1,par2)) a<-table.element(a,signif(r$t0[2],par2)) a<-table.element(a,signif(q3,par2)) a<-table.element(a,signif(p95,par2)) a<-table.element(a,signif(p99,par2)) a<-table.element(a,signif(sqrt(var(r$t[,2])),par2)) a<-table.element(a,signif(q3-q1,par2)) 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,par2)) a<-table.element(a,signif(p05,par2)) a<-table.element(a,signif(q1,par2)) a<-table.element(a,signif(r$t0[3],par2)) a<-table.element(a,signif(q3,par2)) a<-table.element(a,signif(p95,par2)) a<-table.element(a,signif(p99,par2)) a<-table.element(a,signif(sqrt(var(r$t[,3])),par2)) a<-table.element(a,signif(q3-q1,par2)) 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,par2)) a<-table.element(a,signif(p05,par2)) a<-table.element(a,signif(q1,par2)) a<-table.element(a,signif(r$t0[4],par2)) a<-table.element(a,signif(q3,par2)) a<-table.element(a,signif(p95,par2)) a<-table.element(a,signif(p99,par2)) a<-table.element(a,signif(sqrt(var(r$t[,4])),par2)) a<-table.element(a,signif(q3-q1,par2)) 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,par2)) a<-table.element(a,signif(p05,par2)) a<-table.element(a,signif(q1,par2)) a<-table.element(a,signif(r$t0[5],par2)) a<-table.element(a,signif(q3,par2)) a<-table.element(a,signif(p95,par2)) a<-table.element(a,signif(p99,par2)) a<-table.element(a,signif(sqrt(var(r$t[,5])),par2)) a<-table.element(a,signif(q3-q1,par2)) 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