Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
-2.27165721706283 0.771506678877866 2.26979779371262 3.82681347374055 -2.20085664727154 -1.59463906900306 3.72834278293661 -1.53936572130171 -1.41932980610789 1.48223295024959 1.55303352004144 -0.0100201622858731 1.31903313911459 1.06685185696677 -0.763910329598855 -0.271657217063394 0.996051287174924 4.28532501580312 3.36219703505577 0.432997604847614 0.771506678877782 1.17136055007009 3.09448853081745 1.65754221314476 1.28532501580312 1.83285147603993 1.70070610908593 -1.53936572130171 0.58067019389211 0.439069054308411 -0.567002395152386 -0.588567619542262 -0.446966479958563 0.58067019389211 -1.27165721706339 -2.39773113455659 -2.63173151548344 -1.41932980610789 1.55303352004144 2.14372387621942 1.84837869813043 -1.56700239515239 2.53143484849014 0.291396465263916 -0.151621301869572 -3.82863944992991 -2.05318405822705 0.214524446011269 0.919179267922277 -1.1239846280189 -0.834710899390704 -0.00394871282507572 -2.82863944992991 0.389833708906438 -1.53933227414029 2.38379570660707 0.432997604847614 0.58067019389211 0.657542213144756 2.29139646526392 0.783649577799376 0.876015371981101 -0.0963479541682254 -0.637802964944235 0.58067019389211 1.68517888699543 1.90365204583177 3.65754221314476 -3.41932980610789 0.58067019389211 -2.73624020858676 -0.933148143033227 1.31296168965379 1.31296168965379 0.876015371981101 3.92525071738307 -0.348529236316041 1.50986962410026 -2.41932980610789 0.876015371981101 0.58067019389211 0.58067019389211 -0.567002395152386 0.214524446011269 1.99605128717492 0.482232950249588 0.728342782936606 0.876015371981101 0.762050906248076 -1.8131122278394 -0.0255473843763764 0.58067019389211 -0.271657217063394 -2.1239846280189 1.70070610908593 0.389833708906438 2.62990553929408 0.285325015803119 -0.151621301869572 -0.9823834884352 1.65754221314476 2.53143484849014 0.805214802189252 1.21452444601127 -1.53329427184092 1.43299760484761 0.58067019389211 1.95288739123375 0.488304399710385 0.58067019389211 -0.0963479541682254 2.50986962410026 -1.63780296494424 -2.59463906900306 1.79914335272846 -1.90551146918255 1.26375979141324 -1.63780296494424 0.728342782936606 -1.34245778685524 0.0944885308174464 -2.41932980610789 -0.687038310346209 -1.05318405822705 -1.02554738437638 0.553033520041437 1.91917926792228 0.891542594071604 -2.78547555398873 2.14372387621942 -3.89944001972176 3.11608720236875 -2.78547555398873 -1.41932980610789 -0.567002395152386 0.0729233064275707 0.362197035055765 -2.56093094569159 -0.883912797631253 -2.22242187166142 3.11608720236875 1.33456036120509 -0.173219973420872 1.21452444601127 -3.03158538667575 1.79914335272846 -1.82863944992991 1.0236879610256 0.368268484516562 -2.39773113455659 -0.757838880138058 2.33456036120509 4.82681347374055 2.00212273663572 -2.58856761954226 -0.0255473843763764 1.53146829565156 0.891542594071604 1.50986962410026 -0.440895030497766 0.411432380457738 -0.413258356647093 -0.736240208586757 0.165289100609296 1.58674164335291 -2.29322244145327 -0.637802964944235 -3.9763120389744 -1.53329427184092 1.46670572815908 1.60834031490421 0.432997604847614 -1.2655857676026 -2.9547133674231 -3.56700239515239 0.509869624100261 0.00212273663572132 -1.07474928261692 -0.299293890914067 -0.927076693572429 0.728342782936606 -0.320859115303943 1.67910743753463 -0.271657217063394 -7.56700239515239 0.876015371981101 2.36219703505577 -0.250058545512094 -0.954713367423102 0.608340314904208 -0.828639449929907 0.728342782936606 2.28535846296454 1.80521480218925 -1.15162130186957 0.0236879610255973 2.62990553929408 0.728342782936606 1.51594107356106 0.876015371981101 1.46670572815908 -0.151621301869572 -2.8131122278394 -3.1239846280189 2.11608720236875 0.214524446011269 1.36219703505576 0.432997604847614 -2.85627612378058 0.700706109085933 -2.53329427184091 -3.56700239515239 0.805214802189252 2.85445014759123 1.17136055007009 -0.151621301869572 1.91917926792228 -0.00394871282507572 1.14979532568022 -0.41932980610789 -1.82863944992991 -0.299293890914067 0.319033139114589 -1.51776704975041 0.368268484516562 -3.76387688243743 -0.742311658047555 0.925250717383075 -1.48405892643894 -0.927076693572429 1.53146829565156 -3.2655857676026 4.38379570660707 1.60830686774278 -1.29929389091407 -2.56093094569159 -7.17318652625945 -2.1239846280189 1.31903313911459 -1.65936818933411 -0.440895030497766 -1.75176743067726 0.728376230098031 1.82681347374055 0.848378698130428 -0.29322244145327 0.777578128338579 -3.37009446070592 -0.145549852408774 -1.07474928261692 -1.25005854551209 -1.9547133674231 0.679140884696057 2.77757812833858 -1.70860353473608 -0.342457786855244 1.46670572815908 2.00212273663572 -1.78547555398873 -4.97627859181298 0.362197035055765 -4.41932980610789 -0.0686778331561277 0.629905539294083
Chart options
Title:
Y-axis minimum
Y-axis maximum
R Code
geomean <- function(x) { return(exp(mean(log(x)))) } harmean <- function(x) { return(1/mean(1/x)) } quamean <- function(x) { return(sqrt(mean(x*x))) } winmean <- function(x) { x <-sort(x[!is.na(x)]) n<-length(x) denom <- 3 nodenom <- n/denom if (nodenom>40) denom <- n/40 sqrtn = sqrt(n) roundnodenom = floor(nodenom) win <- array(NA,dim=c(roundnodenom,2)) for (j in 1:roundnodenom) { win[j,1] <- (j*x[j+1]+sum(x[(j+1):(n-j)])+j*x[n-j])/n win[j,2] <- sd(c(rep(x[j+1],j),x[(j+1):(n-j)],rep(x[n-j],j)))/sqrtn } return(win) } trimean <- function(x) { x <-sort(x[!is.na(x)]) n<-length(x) denom <- 3 nodenom <- n/denom if (nodenom>40) denom <- n/40 sqrtn = sqrt(n) roundnodenom = floor(nodenom) tri <- array(NA,dim=c(roundnodenom,2)) for (j in 1:roundnodenom) { tri[j,1] <- mean(x,trim=j/n) tri[j,2] <- sd(x[(j+1):(n-j)]) / sqrt(n-j*2) } return(tri) } midrange <- function(x) { return((max(x)+min(x))/2) } q1 <- function(data,n,p,i,f) { np <- n*p; i <<- floor(np) f <<- np - i qvalue <- (1-f)*data[i] + f*data[i+1] } q2 <- function(data,n,p,i,f) { np <- (n+1)*p i <<- floor(np) f <<- np - i qvalue <- (1-f)*data[i] + f*data[i+1] } q3 <- function(data,n,p,i,f) { np <- n*p i <<- floor(np) f <<- np - i if (f==0) { qvalue <- data[i] } else { qvalue <- data[i+1] } } q4 <- function(data,n,p,i,f) { np <- n*p i <<- floor(np) f <<- np - i if (f==0) { qvalue <- (data[i]+data[i+1])/2 } else { qvalue <- data[i+1] } } q5 <- function(data,n,p,i,f) { np <- (n-1)*p i <<- floor(np) f <<- np - i if (f==0) { qvalue <- data[i+1] } else { qvalue <- data[i+1] + f*(data[i+2]-data[i+1]) } } q6 <- function(data,n,p,i,f) { np <- n*p+0.5 i <<- floor(np) f <<- np - i qvalue <- data[i] } q7 <- function(data,n,p,i,f) { np <- (n+1)*p i <<- floor(np) f <<- np - i if (f==0) { qvalue <- data[i] } else { qvalue <- f*data[i] + (1-f)*data[i+1] } } q8 <- function(data,n,p,i,f) { np <- (n+1)*p i <<- floor(np) f <<- np - i if (f==0) { qvalue <- data[i] } else { if (f == 0.5) { qvalue <- (data[i]+data[i+1])/2 } else { if (f < 0.5) { qvalue <- data[i] } else { qvalue <- data[i+1] } } } } midmean <- function(x,def) { x <-sort(x[!is.na(x)]) n<-length(x) if (def==1) { qvalue1 <- q1(x,n,0.25,i,f) qvalue3 <- q1(x,n,0.75,i,f) } if (def==2) { qvalue1 <- q2(x,n,0.25,i,f) qvalue3 <- q2(x,n,0.75,i,f) } if (def==3) { qvalue1 <- q3(x,n,0.25,i,f) qvalue3 <- q3(x,n,0.75,i,f) } if (def==4) { qvalue1 <- q4(x,n,0.25,i,f) qvalue3 <- q4(x,n,0.75,i,f) } if (def==5) { qvalue1 <- q5(x,n,0.25,i,f) qvalue3 <- q5(x,n,0.75,i,f) } if (def==6) { qvalue1 <- q6(x,n,0.25,i,f) qvalue3 <- q6(x,n,0.75,i,f) } if (def==7) { qvalue1 <- q7(x,n,0.25,i,f) qvalue3 <- q7(x,n,0.75,i,f) } if (def==8) { qvalue1 <- q8(x,n,0.25,i,f) qvalue3 <- q8(x,n,0.75,i,f) } midm <- 0 myn <- 0 roundno4 <- round(n/4) round3no4 <- round(3*n/4) for (i in 1:n) { if ((x[i]>=qvalue1) & (x[i]<=qvalue3)){ midm = midm + x[i] myn = myn + 1 } } midm = midm / myn return(midm) } (arm <- mean(x)) sqrtn <- sqrt(length(x)) (armse <- sd(x) / sqrtn) (armose <- arm / armse) (geo <- geomean(x)) (har <- harmean(x)) (qua <- quamean(x)) (win <- winmean(x)) (tri <- trimean(x)) (midr <- midrange(x)) midm <- array(NA,dim=8) for (j in 1:8) midm[j] <- midmean(x,j) midm bitmap(file='test1.png') lb <- win[,1] - 2*win[,2] ub <- win[,1] + 2*win[,2] if ((ylimmin == '') | (ylimmax == '')) plot(win[,1],type='b',main=main, xlab='j', pch=19, ylab='Winsorized Mean(j/n)', ylim=c(min(lb),max(ub))) else plot(win[,1],type='l',main=main, xlab='j', pch=19, ylab='Winsorized Mean(j/n)', ylim=c(ylimmin,ylimmax)) lines(ub,lty=3) lines(lb,lty=3) grid() dev.off() bitmap(file='test2.png') lb <- tri[,1] - 2*tri[,2] ub <- tri[,1] + 2*tri[,2] if ((ylimmin == '') | (ylimmax == '')) plot(tri[,1],type='b',main=main, xlab='j', pch=19, ylab='Trimmed Mean(j/n)', ylim=c(min(lb),max(ub))) else plot(tri[,1],type='l',main=main, xlab='j', pch=19, ylab='Trimmed Mean(j/n)', ylim=c(ylimmin,ylimmax)) lines(ub,lty=3) lines(lb,lty=3) grid() dev.off() load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Central Tendency - Ungrouped Data',4,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Measure',header=TRUE) a<-table.element(a,'Value',header=TRUE) a<-table.element(a,'S.E.',header=TRUE) a<-table.element(a,'Value/S.E.',header=TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Arithmetic Mean',header=TRUE) a<-table.element(a,signif(arm,6)) a<-table.element(a, signif(armse,6)) a<-table.element(a,signif(armose,6)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, 'Geometric Mean',header=TRUE) a<-table.element(a,signif(geo,6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, 'Harmonic Mean',header=TRUE) a<-table.element(a,signif(har,6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, 'Quadratic Mean',header=TRUE) a<-table.element(a,signif(qua,6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) for (j in 1:length(win[,1])) { a<-table.row.start(a) mylabel <- paste('Winsorized Mean (',j) mylabel <- paste(mylabel,'/') mylabel <- paste(mylabel,length(win[,1])) mylabel <- paste(mylabel,')') a<-table.element(a, mylabel,header=TRUE) a<-table.element(a,signif(win[j,1],6)) a<-table.element(a,signif(win[j,2],6)) a<-table.element(a,signif(win[j,1]/win[j,2],6)) a<-table.row.end(a) } for (j in 1:length(tri[,1])) { a<-table.row.start(a) mylabel <- paste('Trimmed Mean (',j) mylabel <- paste(mylabel,'/') mylabel <- paste(mylabel,length(tri[,1])) mylabel <- paste(mylabel,')') a<-table.element(a, mylabel,header=TRUE) a<-table.element(a,signif(tri[j,1],6)) a<-table.element(a,signif(tri[j,2],6)) a<-table.element(a,signif(tri[j,1]/tri[j,2],6)) a<-table.row.end(a) } a<-table.row.start(a) a<-table.element(a, 'Median',header=TRUE) a<-table.element(a,signif(median(x),6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a, 'Midrange',header=TRUE) a<-table.element(a,signif(midr,6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) mymid <- 'Midmean' mylabel <- paste(mymid,'Weighted Average at Xnp',sep=' - ') a<-table.element(a,mylabel,header=TRUE) a<-table.element(a,signif(midm[1],6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) mymid <- 'Midmean' mylabel <- paste(mymid,'Weighted Average at X(n+1)p',sep=' - ') a<-table.element(a,mylabel,header=TRUE) a<-table.element(a,signif(midm[2],6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) mymid <- 'Midmean' mylabel <- paste(mymid,'Empirical Distribution Function',sep=' - ') a<-table.element(a,mylabel,header=TRUE) a<-table.element(a,signif(midm[3],6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) mymid <- 'Midmean' mylabel <- paste(mymid,'Empirical Distribution Function - Averaging',sep=' - ') a<-table.element(a,mylabel,header=TRUE) a<-table.element(a,signif(midm[4],6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) mymid <- 'Midmean' mylabel <- paste(mymid,'Empirical Distribution Function - Interpolation',sep=' - ') a<-table.element(a,mylabel,header=TRUE) a<-table.element(a,signif(midm[5],6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) mymid <- 'Midmean' mylabel <- paste(mymid,'Closest Observation',sep=' - ') a<-table.element(a,mylabel,header=TRUE) a<-table.element(a,signif(midm[6],6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) mymid <- 'Midmean' mylabel <- paste(mymid,'True Basic - Statistics Graphics Toolkit',sep=' - ') a<-table.element(a,mylabel,header=TRUE) a<-table.element(a,signif(midm[7],6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) mymid <- 'Midmean' mylabel <- paste(mymid,'MS Excel (old versions)',sep=' - ') a<-table.element(a,mylabel,header=TRUE) a<-table.element(a,signif(midm[8],6)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Number of observations',header=TRUE) a<-table.element(a,signif(length(x),6)) a<-table.element(a,'') a<-table.element(a,'') 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