Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
-0.688542 0.402673 -0.658137 0.471084 0.341863 0.425477 -0.513713 0.326661 0.5699 0.250649 0.433078 0.425477 0.524293 -0.513713 -0.498511 -0.521314 0.433078 -0.536517 0.395072 -0.452903 0.471084 0.455882 -0.528916 0.478686 0.448281 0.463483 0.448281 -0.604928 0.440679 -0.490909 -0.536517 -0.612529 0.516692 0.478686 0.455882 0.471084 0.402673 -0.536517 0.417876 0.372268 -0.490909 0.326661 0.37987 0.448281 0.509091 0.425477 -0.559321 0.37987 0.455882 -0.490909 0.410275 0.501489 0.410275 0.387471 0.486287 0.410275 0.463483 -0.498511 0.433078 -0.589725 0.410275 0.417876 -0.536517 -0.498511 0.395072 -0.612529 0.402673 -0.513713 0.395072 0.433078 -0.582124 -0.62013 -0.582124 -0.612529 0.357066 0.433078 -0.475707 0.334262 -0.559321 0.37987 -0.566922 -0.551719 0.471084 -0.566922 -0.513713 -0.521314 0.471084 -0.62013 0.417876 0.433078 -0.475707 -0.551719 -0.513713 0.402673 0.425477 0.448281 -0.574523 0.387471 -0.528916 0.425477 0.448281 0.417876 0.493888 0.455882 -0.490909 0.448281 0.478686 0.440679 -0.589725 -0.62013 -0.513713 -0.566922 -0.597327 -0.544118 0.433078 0.31906 0.357066 0.440679 0.455882 -0.635333 0.357066 0.349465 -0.589725 0.440679 -0.574523 -0.711345 0.509091 0.417876 0.448281 0.417876 0.425477 0.471084 0.410275 0.387471 -0.551719 0.478686 0.463483 0.395072 0.395072 0.402673 0.471084 -0.528916 0.31906 0.357066 0.463483 0.486287 0.509091 -0.62013 -0.551719 -0.703744 0.326661 -0.62013 -0.612529 -0.62013 0.364667 -0.468106 0.387471 0.387471 -0.62013 0.402673 0.37987 0.395072 -0.528916 -0.635333 -0.566922 -0.589725 0.455882 -0.627732 0.37987 -0.703744 0.433078 -0.513713 0.440679 -0.604928 -0.551719 0.501489 -0.597327 0.455882 0.372268 0.387471 0.493888 0.349465 -0.597327 0.410275 -0.544118 -0.627732 -0.627732 0.448281 -0.521314 0.402673 -0.597327 -0.574523 -0.544118 -0.566922 0.478686 0.402673 0.501489 0.372268 -0.566922 -0.650535 0.288655 -0.521314 -0.597327 -0.582124 -0.582124 0.0758201 0.311458 -0.551719 -0.536517 0.531894 -0.544118 0.471084 0.296256 0.433078 -0.513713 0.493888 0.433078 0.471084 -0.506112 -0.475707 -0.498511 0.334262 0.387471 -0.688542 0.463483 0.455882 -0.574523 -0.475707 -0.475707 0.37987 0.478686 -0.536517 -0.589725 -0.589725 0.471084 0.493888 -0.559321 -0.566922 -0.597327 0.539496 0.486287 0.440679 0.395072 -0.513713 0.493888 0.174636 0.455882 -0.544118 -0.4301 0.417876 -0.627732 0.387471 0.37987 0.524293 -0.635333 0.387471 -0.551719 0.448281 0.486287 -0.544118 0.440679 -0.612529 0.425477 -0.566922 0.493888 -0.452903 -0.521314 0.349465 0.303857 -0.506112 0.334262 0.417876 -0.574523 0.364667 -0.68094 -0.460504 0.440679 0.463483 -0.506112 0.501489 -0.536517 -0.544118 0.425477 -0.650535 -0.528916 0.425477 0.455882 0.296256 -0.544118 0.448281 0.493888 0.440679 0.425477 -0.589725 0.493888 -0.772155 -0.551719 0.455882 -0.551719 0.410275 -0.559321 -0.536517 0.448281 -0.566922 0.463483 0.524293 -0.62013 -0.597327 -0.498511 0.493888 -0.513713 0.402673 0.455882 0.425477 0.440679 -0.513713 0.326661 -0.604928 0.334262 -0.589725 -0.475707 -0.574523 -0.483308 -0.544118 -0.551719 0.509091 0.410275 0.243047 -0.490909 -0.559321 -0.544118 -0.62013 -0.513713 0.478686 0.326661 0.425477 0.410275 -0.475707 0.357066 -0.551719 0.357066 0.417876 -0.551719 0.501489 -0.559321 -0.650535 -0.513713 -0.566922 -0.498511 -0.483308 0.425477 -0.498511 -0.635333 0.395072 0.471084 0.478686 0.402673 -0.566922 0.433078 0.410275 -0.612529 -0.582124 -0.506112 -0.62013 0.471084 -0.544118 -0.589725 -0.513713 -0.650535 -0.483308 -0.536517 0.417876 0.524293 0.387471 -0.566922 -0.483308 0.448281 -0.475707 0.448281 0.387471 0.425477 0.417876 0.455882 -0.521314 -0.688542 -0.521314 -0.528916 -0.460504 0.493888 -0.582124 -0.544118 0.509091 -0.498511 -0.468106 0.486287 -0.612529 0.509091 0.349465 -0.665738 0.402673 -0.544118 0.524293 -0.551719 0.440679 0.501489 0.37987 -0.513713 0.463483 0.37987 -0.582124 0.410275 0.440679 0.372268 0.471084 -0.468106 0.334262 -0.490909 -0.566922 0.417876 -0.635333 0.440679 0.478686 -0.627732 -0.506112 0.364667 0.402673 0.493888 0.341863 -0.597327 0.433078 0.417876 0.493888 0.471084 -0.566922 -0.566922 0.410275 -0.574523 0.463483 0.433078 -0.566922 0.524293 -0.483308 0.478686 -0.574523 -0.559321 0.478686 0.372268 0.5699 0.440679 -0.4301 -0.627732 0.410275 -0.498511 0.349465 -0.521314 0.372268 0.455882 0.341863 0.372268 0.440679 0.516692 0.539496 0.440679 -0.498511 0.471084 -0.521314 -0.498511 0.531894 0.509091 0.478686 0.547097 0.326661 0.509091 0.387471 -0.521314 0.516692 0.417876 0.440679 -0.521314 0.531894 0.235446 0.471084 0.486287 -0.498511 0.463483 0.433078 -0.566922 -0.513713 0.493888 0.372268 0.31906 0.501489 -0.536517 0.311458 -0.582124 0.516692 0.471084
# simulations
Significant digits
Bandwidth
(?)
Quantiles
P1 P5 Q1 Q3 P95 P99
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) } (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