Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
1.72856105193432 0.537017873831517 0.867980543338276 -0.469786520936907 -1.65809537332321 0.146636501699633 -0.0672386349137408 -0.13498404516995 2.12223240279242 -1.19320887060888 1.27933928266697 2.96694884090103 1.39545186365274 -2.1451992140415 -2.43398275751913 0.632864700711983 -0.682707813636235 2.07131204525417 -0.104531746575379 -1.43248560723516 2.06729526607686 -0.762199524917308 -0.656220032961174 -0.272891878011051 -1.38725178006058 2.73448056175203 -0.406928566252687 -1.06644693820818 -0.267090064307712 0.603099469782987 -2.24501673460205 0.704074348654696 -0.14323212663897 1.96584936240914 -0.644282435815208 0.51371072453424 -0.219904883989661 1.53433608960596 0.873410873202486 -0.0810741216198124 2.22709173221706 0.172181871131915 0.156343873949484 -0.119778140438289 0.201567470778364 0.72333345086232 1.94101724859539 0.96842651759471 2.62676169681312 -1.00490408093144 1.4317321639036 0.363326953488502 -1.19455057992329 0.874908023486463 -1.33600415902196 0.777564046322009 1.7059917118051 -1.51728851106142 0.280648212437612 0.628560037292622 -0.380970785994089 -1.10180574201106 0.488976525397815 0.548138796768258 2.34695896401334 1.8680589844253 2.42664197141552 0.772644235160232 -0.148150496596256 -0.66314751310936 1.18319180419306 -1.85385290900922 -0.192985025699194 -1.88497624723516 -1.60767366973841 -1.96420262347225 -1.83163580988418 0.356353819920789 -0.139424675968572 1.03002834236818 -0.87383034679475 1.0942591515073 -0.75970903681286 -1.86572130831276 0.0307532695368973 1.45943548478181 -0.968566254388339 -0.27359733158947 -0.779667001008091 -0.450534010320096 -3.78219026865631 0.570020715704707 -0.908921286671703 -0.550266728957996 0.0788218936010009 1.29430043581788 -0.406246225019453 -0.83598776925008 -0.719487754941648 -0.207533768023798 3.51491483206611 0.741532568484479 0.0363195667803673 0.636057003639965 -1.76614953455778 -1.79476412153827 -1.94638170405188 0.0661695996320116 0.664714600509923 0.39974322111863 -3.75942479132569 1.63115007447768 2.20831840196341 -0.268443445172315 0.783843599184011 1.27822692217559 -0.354551216060357 -0.134550526350053 -3.22470676015799 1.20211034507058 -0.368260965732864 -0.248380427983121 1.57769314333774 -0.198083020227395 -0.505520871287781 1.04816849249361 -1.52550138468108 0.153544508257518 1.72694554270275 -1.19023140010116 -0.742410307645107 0.485072601254582 1.20775671780436 -0.183102393486175 -1.34917720203908 1.35985348068112 -0.0818385426949573 1.10174444630446 0.614187843931161 -2.76102656252592 0.960157328132108 0.0691798496834958 -2.56236140147681 0.135687170538491 -1.08769931577129 0.385580902989905 -1.02234900782415 -0.6986228107012 -1.49310553498315 1.53209806266973 -0.0385009625534758 -1.3395232933726 -1.08832226555382 -2.10030958695191 -1.11463205606125 0.667893597484437 -0.105826159943256 0.547492310686565 -0.285955141601865 1.26088518311837 -1.99925770819766 2.40402841094769 -0.789845751580985 3.53418064863639 0.308019428734996 -1.87700317613243 -0.383228286520615 -1.07287170169405 1.87214923937838 0.125541199509366 0.574284797540187 0.337349699639604 -0.653471487322515 -0.494244507381467 -0.359198646483944 0.902932440228196 1.4871314356976 -1.69322526038076 -1.45487145328323
# 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