Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
-0.0447135254057231 -0.0681917505477611 0.197918269607264 0.364977510288303 151.337.996.933.282 -0.359461693682944 0.457275231630557 -0.576200433829858 -0.35806144578699 126.000.911.969.899 -134.266.139.996.902 -0.35328957148349 0.161395641520833 -0.857823002973426 -0.555144912300313 -0.727027913044969 -0.369042261654297 -0.0578971497923536 -0.769036818583236 -0.853703554579139 0.700854566671402 -0.676058766811168 0.868578036660101 -0.108845889757913 -157.887.204.236.239 -112.268.950.514.506 0.396923851905675 0.0148495697982022 0.140259353790972 0.36986509106495 -0.279776101924479 0.304546129479045 0.892466302712004 0.0893800831529647 0.135269394216119 -120.736.758.947.528 -0.230372958786352 -0.0877971490505129 -0.333232967805547 0.601966275119466 0.489790331921813 -0.300400372795234 0.101351037939037 0.596346950586006 -0.47683886806596 -0.418379377480445 -0.117113003681212 -0.14718075434479 0.532293994236984 -0.976528360213378 0.331179682756783 0.869619341824002 -0.452578488121401 -0.421389494638729 -0.0683600581132898 0.326102644105208 0.848949464326408 0.455827855734715 0.823473300306631 0.933736135720925 123.343.911.513.023 0.408339888885549 0.287976014213179 -0.211203400023215 -0.267468886934195 -111.645.429.924.937 -0.0317247398908312 0.547972212918 0.120802505224662 -0.868321571498795 -0.311155680966178 -0.380689102119736 -0.305618643653553 -0.409940938228127 -0.00967793776262078 0.495426344923938 -0.704347377621235 -0.063469879208485 -0.51389881333356 0.58951146076434 -0.349423192929551 0.641305316981306 0.0343658031275834 -0.052752170661352 -0.880138087903201 -0.110617811300259 0.725017567873533 -0.505210289918464 101.506.155.749.305 0.432165488939686 -0.605966485103475 -0.999657590804615 -0.273898494723388 0.276038635783788 100.336.898.137.961 0.0129978236300762 -0.560227111581836 -0.549704221644895 -0.264710628773928 0.279628858544481 0.678638849882516 0.661069134033666 -0.704450399639958 -0.201260460274066 0.350125398121403 0.512074040089319 0.855206906548568 0.194069541456819 0.723929039883118 118.131.067.430.535 0.142881527617021 0.017925799932583 -0.497956664933812 -0.299719686679989 0.135788266531892 -0.170331262015063 -103.309.848.477.483 -0.173819270986012 -0.963024089117913 0.698379739279324 -0.367089110910225 -0.263848476436583 -0.418819019840626 -112.067.684.130.117 0.050058520577361 0.608903578527243 0.134690293435698 0.330041347787194 0.127838692569412 0.581319132516513 -0.0387079435373926 -0.870097194938493 -0.462113796849748 -0.775456514875374 1.395.869.467.623 -0.373271731449615 -0.268080176510708 108.656.232.048.462 -0.325879955921928 0.306618547884575 -0.552774192453413 0.928735518612731 0.0930554068876903 0.793725485236845 -0.0655196055780801 0.319389737259029 -0.466224915112025 -0.262872803122903 0.0101377300968685 0.166709542388814 -0.279488697025563 -0.416437753013679 -0.221155653982297 -0.180257931962942 -0.698512272287848 -0.0309661506328094 -0.218027303452013 -0.373584442457897 0.0529514273142018 0.16272689702866 0.827781138510815 -0.724325872700113 -0.470004756981956 104.570.491.743.386 -0.192863384714999 -0.57118024241685 0.531291140789378 -0.305718227359383 0.338050048474844 0.476242122330853 -0.813898476355091 -0.0347526256916884 0.3072752726873 0.298265453953271 -0.356205547514185 -0.356604247457695 0.1619143745199 0.157622107243565 0.301928819166098 -0.561981939800843 -0.054344435964539 -0.242625453728938 -0.0291188126122691 0.228470615616782 -0.510675391591998 112.387.561.916.988 -112.598.799.910.622 0.290323760058784 0.190186117613264 0.0853119358914465 -0.70803898396401 0.0821216895898958 -0.307060856617796 0.525053351110503 -0.601608778527328 0.53818181553833 -0.306512070731858 0.654542044064766 -0.537374859392012 -0.186945154315006 -0.0415286501896012 -0.0882630415428376 -0.256492309570936 -0.413059227096436 -0.270151685278362 -0.435615700218065 0.546649355950973 0.323884533274175 0.661932000134738 0.402774353219003 -0.428179650529287 -0.185363732952725 -0.146275810216839 0.177901302311314 -0.370503387950714 0.205238677066994 -0.0894107718995766 -0.0207802658241447 -0.00240155924277742 0.027237259338887 -0.304012117517099 150.787.083.654.866 0.259049823633338 -0.546357645795879 0.581469268528633 0.330213798462246 -0.983366580775191 -0.736262986208729 -0.338744406461511 0.76053719760488 -0.261709569948618 -0.476048053943579 -0.110367750063661 169.080.994.427.217 0.149871963217241 -0.894443959356019 0.0854707859389636 -0.117916015475136 -0.215892867114885 -0.394092028282386 0.0924263195771785 0.058591669911839 0.253245871086277 0.370568028643805 -0.386523945873864 0.44711088562219 0.623154288482276 -0.184203126128357 0.705478214339145 -0.317030724633186 -0.927475247660286 -0.0393758868189332 100.274.756.958.968 0.812619294578401 0.298293415694001 0.152551035780635 -0.150639291567251 0.110741928054607 0.552786542794576 0.043495317643748 0.363209762554804 -0.022797008470584 0.504138678148394 0.138779289305952 -0.0853035308417428 -0.496484911315823 -0.0852251911408867 -0.247830972313322 -0.120849214187901 -0.451770225462982 0.612150462832856 0.350625301752204 -0.359290488889893 -0.484436529031112 0.339430852270173 -0.0417944819768093 -0.0100068969567184 -0.403912575239055 0.151831207953865 -0.229701003816159 -0.217361643531165 -0.332188937105677 0.264417994027107 0.176869376138016 -0.176677800263361 -0.135670367154218 -0.827944200296059 -0.0722734399204414 -0.117607845832689 0.314635053724339 -0.154386181086592 0.163002765901339 -0.244739215669379 -0.204481562449311 0.0360259309576101 0.00417168114425691 0.266111866388402 -0.651003114649863 0.589994534676515 0.312322955620105 0.552669928748546 -0.0962287489408798 -0.467083613141848 -0.198196785314342 0.395575343631614 0.175735234004858 0.315653622056923 -0.161149078505793 0.882040919470625 0.0637848764428182 0.858633851994241 0.822761597216032 17.726.688.669.086 -0.566981098222408 0.153261675817818 -0.279232930441738 0.049572050471916 -116.973.330.350.582 -0.0821094023296669 0.0681021355632607 -0.201544972117859 0.0764370043960761 -0.52636014292932 -0.0804139922908793 -0.39060373403853 -0.366201530340365 -0.236325627413792 -0.0197748907207889 -0.400593230734952 0.273305313443463 0.571670833577912 0.355107022563698 -0.598053247795535 0.0680124314430453 0.0825157705042455 -0.236688909377822 -0.720993035857933 0.445855527021645 -0.278184399175499 -0.813695706192672 0.0382223634028053 0.287981462997034 -0.397401294088841 0.453441800201533 -0.336993631715903 0.0350016901299891 -0.12907406621208 -0.92948724921503 0.112538313206602 -0.266214839380531 0.318919353798501 -0.235823743806231 0.312056037491925 -0.607348899940247 0.821439552875221 -0.330769416907408 -0.0368470896547259 -0.22341151284797 -0.0162393530813855 0.494328753928939
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,hyperlink('http://www.xycoon.com/arithmetic_mean.htm', 'Arithmetic Mean', 'click to view the definition of the Arithmetic Mean'),header=TRUE) a<-table.element(a,arm) a<-table.element(a,hyperlink('http://www.xycoon.com/arithmetic_mean_standard_error.htm', armse, 'click to view the definition of the Standard Error of the Arithmetic Mean')) a<-table.element(a,armose) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,hyperlink('http://www.xycoon.com/geometric_mean.htm', 'Geometric Mean', 'click to view the definition of the Geometric Mean'),header=TRUE) a<-table.element(a,geo) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,hyperlink('http://www.xycoon.com/harmonic_mean.htm', 'Harmonic Mean', 'click to view the definition of the Harmonic Mean'),header=TRUE) a<-table.element(a,har) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,hyperlink('http://www.xycoon.com/quadratic_mean.htm', 'Quadratic Mean', 'click to view the definition of the Quadratic Mean'),header=TRUE) a<-table.element(a,qua) 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,hyperlink('http://www.xycoon.com/winsorized_mean.htm', mylabel, 'click to view the definition of the Winsorized Mean'),header=TRUE) a<-table.element(a,win[j,1]) a<-table.element(a,win[j,2]) a<-table.element(a,win[j,1]/win[j,2]) 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,hyperlink('http://www.xycoon.com/arithmetic_mean.htm', mylabel, 'click to view the definition of the Trimmed Mean'),header=TRUE) a<-table.element(a,tri[j,1]) a<-table.element(a,tri[j,2]) a<-table.element(a,tri[j,1]/tri[j,2]) a<-table.row.end(a) } a<-table.row.start(a) a<-table.element(a,hyperlink('http://www.xycoon.com/median_1.htm', 'Median', 'click to view the definition of the Median'),header=TRUE) a<-table.element(a,median(x)) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,hyperlink('http://www.xycoon.com/midrange.htm', 'Midrange', 'click to view the definition of the Midrange'),header=TRUE) a<-table.element(a,midr) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) mymid <- hyperlink('http://www.xycoon.com/midmean.htm', 'Midmean', 'click to view the definition of the Midmean') mylabel <- paste(mymid,hyperlink('http://www.xycoon.com/method_1.htm','Weighted Average at Xnp',''),sep=' - ') a<-table.element(a,mylabel,header=TRUE) a<-table.element(a,midm[1]) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) mymid <- hyperlink('http://www.xycoon.com/midmean.htm', 'Midmean', 'click to view the definition of the Midmean') mylabel <- paste(mymid,hyperlink('http://www.xycoon.com/method_2.htm','Weighted Average at X(n+1)p',''),sep=' - ') a<-table.element(a,mylabel,header=TRUE) a<-table.element(a,midm[2]) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) mymid <- hyperlink('http://www.xycoon.com/midmean.htm', 'Midmean', 'click to view the definition of the Midmean') mylabel <- paste(mymid,hyperlink('http://www.xycoon.com/method_3.htm','Empirical Distribution Function',''),sep=' - ') a<-table.element(a,mylabel,header=TRUE) a<-table.element(a,midm[3]) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) mymid <- hyperlink('http://www.xycoon.com/midmean.htm', 'Midmean', 'click to view the definition of the Midmean') mylabel <- paste(mymid,hyperlink('http://www.xycoon.com/method_4.htm','Empirical Distribution Function - Averaging',''),sep=' - ') a<-table.element(a,mylabel,header=TRUE) a<-table.element(a,midm[4]) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) mymid <- hyperlink('http://www.xycoon.com/midmean.htm', 'Midmean', 'click to view the definition of the Midmean') mylabel <- paste(mymid,hyperlink('http://www.xycoon.com/method_5.htm','Empirical Distribution Function - Interpolation',''),sep=' - ') a<-table.element(a,mylabel,header=TRUE) a<-table.element(a,midm[5]) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) mymid <- hyperlink('http://www.xycoon.com/midmean.htm', 'Midmean', 'click to view the definition of the Midmean') mylabel <- paste(mymid,hyperlink('http://www.xycoon.com/method_6.htm','Closest Observation',''),sep=' - ') a<-table.element(a,mylabel,header=TRUE) a<-table.element(a,midm[6]) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) mymid <- hyperlink('http://www.xycoon.com/midmean.htm', 'Midmean', 'click to view the definition of the Midmean') mylabel <- paste(mymid,hyperlink('http://www.xycoon.com/method_7.htm','True Basic - Statistics Graphics Toolkit',''),sep=' - ') a<-table.element(a,mylabel,header=TRUE) a<-table.element(a,midm[7]) a<-table.element(a,'') a<-table.element(a,'') a<-table.row.end(a) a<-table.row.start(a) mymid <- hyperlink('http://www.xycoon.com/midmean.htm', 'Midmean', 'click to view the definition of the Midmean') mylabel <- paste(mymid,hyperlink('http://www.xycoon.com/method_8.htm','MS Excel (old versions)',''),sep=' - ') a<-table.element(a,mylabel,header=TRUE) a<-table.element(a,midm[8]) 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,length(x)) 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
1 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation