Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
-0.0447135254057231 -0.0681917505478084 0.197918269607118 0.364977510288305 151.337.996.933.159 -0.359461693683337 0.457275231643101 -0.576200433829823 -0.358061445786517 126.000.911.969.509 -134.266.139.991.407 -0.353289571497946 0.161395641494175 -0.85782300297493 -0.555144912298359 -0.727027913044961 -0.369042261655499 -0.0578971497992717 -0.769036818576891 -0.853703554581848 0.700854566672244 -0.676058766832693 0.868578036694034 -0.10884588975301 -157.887.204.236.819 -112.268.950.514.398 0.39692385190415 0.0148495697973627 0.140259353790028 0.369865091064737 -0.279776101924426 0.304546129475365 0.892466302715145 0.0893800831442478 0.135269394231309 -120.736.758.947.938 -0.230372958786898 -0.0877971490509319 -0.33323296780618 0.60196627511901 0.48979033192109 -0.300400372794546 0.101351037939174 0.596346950582894 -0.476838868063972 -0.418379377485807 -0.117113003673828 -0.147180754345273 0.532293994236285 -0.976528360213732 0.331179682756605 0.869619341823574 -0.452578488121829 -0.42138949463861 -0.0683600581129025 0.326102644103438 0.848949464327804 0.455827855730794 0.823473300311693 0.933736135720844 123.343.911.512.994 0.408339888885158 0.287976014212968 -0.211203400023637 -0.267468886934765 -111.645.429.924.972 -0.0317247398904611 0.547972212916532 0.120802505225608 -0.868321571501203 -0.311155680962737 -0.380689102119501 -0.305618643654202 -0.409940938228271 -0.00967793776263202 0.495426344923855 -0.704347377621387 -0.0634698792084991 -0.513898813333212 0.589511460763565 -0.349423192928886 0.641305316979721 0.0343658031299827 -0.0527521706610295 -0.880138087903553 -0.110617811300362 0.725017567873349 -0.50521028991854 1.015.061.557.493 0.432165488939549 -0.605966485103154 -0.999657590805216 -0.273898494723043 0.276038635782484 100.336.898.138.124 0.0129978236302341 -0.560227111581804 -0.549704221644982 -0.264710628774072 0.279628858544352 0.678638849882419 0.661069134033749 -0.704450399639775 -0.201260460274452 0.350125398121578 0.512074040088484 0.855206906549666 0.194069541456855 0.723929039883086 118.131.067.430.518 0.142881527616885 0.0179257999324864 -0.497956664934098 -0.299719686680144 0.135788266531753 -0.170331262015486 -103.309.848.477.463 -0.173819270986569 -0.963024089117263 0.698379739279439 -0.367089110910285 -0.263848476436372 -0.418819019840622 -112.067.684.130.113 0.0500585205773671 0.608903578527248 0.134690293435901 0.330041347787177 0.127838692569612 0.581319132516197 -0.038707943536881 -0.870097194938441 -0.462113796849836 -0.775456514875537 139.586.946.762.297 -0.373271731449701 -0.268080176510487 108.656.232.048.461 -0.325879955921929 0.30661854788457 -0.552774192453423 0.92873551861249 0.0930554068879337 0.793725485236928 -0.0655196055781843 0.31938973725902 -0.466224915112164 -0.262872803122968 0.0101377300966978 0.166709542388761 -0.27948869702557 -0.416437753013764 -0.221155653982262 -0.180257931963139 -0.698512272287555 -0.030966150632744 -0.218027303452022 -0.373584442457798 0.052951427314271 0.16272689702871 0.827781138510899 -0.724325872700048 -0.470004756981797 104.570.491.743.371 -0.192863384714994 -0.571180242416833 0.531291140789517 -0.305718227359467 0.338050048474922 0.476242122330785 -0.81389847635506 -0.0347526256916431 0.307275272687164 0.298265453953272 -0.356205547514144 -0.356604247457693 0.16191437451988 0.157622107243432 0.301928819166272 -0.561981939800824 -0.0543444359644914 -0.242625453729012 -0.0291188126122587 0.228470615616776 -0.510675391591978 112.387.561.916.997 -112.598.799.910.629 0.290323760058927 0.190186117613125 0.0853119358914609 -0.708038983963899 0.0821216895899207 -0.307060856617881 0.525053351110559 -0.601608778527341 0.538181815538459 -0.306512070731953 0.654542044064886 -0.537374859392079 -0.186945154314925 -0.0415286501896635 -0.0882630415428836 -0.256492309570866 -0.41305922709641 -0.270151685278353 -0.435615700218057 0.546649355951026 0.323884533274196 0.661932000134874 0.402774353219053 -0.428179650529226 -0.185363732952745 -0.146275810216944 0.177901302311226 -0.370503387950712 0.205238677067021 -0.0894107718996266 -0.0207802658241128 -0.00240155924277281 0.0272372593388978 -0.30401211751709 150.787.083.654.869 0.259049823633271 -0.546357645795726 0.581469268528577 0.330213798462058 -0.983366580775167 -0.736262986208755 -0.338744406461673 0.760537197604814 -0.261709569948611 -0.476048053943411 -0.110367750063647 169.080.994.427.218 0.149871963217245 -0.894443959355815 0.0854707859389192 -0.11791601547532 -0.215892867114871 -0.394092028282406 0.0924263195771732 0.0585916699118239 0.253245871086333 0.370568028643833 -0.386523945873816 0.447110885622225 0.623154288482199 -0.184203126128336 0.705478214339162 -0.317030724633314 -0.927475247660233 -0.0393758868190344 100.274.756.958.954 0.812619294578388 0.298293415694118 0.152551035780658 -0.150639291567325 0.110741928054486 0.552786542794449 0.0434953176436756 0.36320976255479 -0.0227970084706981 0.504138678148357 0.138779289305853 -0.0853035308417627 -0.496484911315902 -0.0852251911409686 -0.247830972313436 -0.120849214187909 -0.451770225463026 0.612150462832891 0.350625301752197 -0.359290488889771 -0.48443652903107 0.339430852270146 -0.0417944819768526 -0.0100068969566531 -0.403912575239062 0.151831207953882 -0.229701003816198 -0.217361643531106 -0.332188937105688 0.264417994027136 0.176869376138034 -0.176677800263279 -0.13567036715417 -0.827944200296052 -0.0722734399204132 -0.117607845832746 0.314635053724412 -0.154386181086533 0.163002765901451 -0.244739215669361 -0.204481562449259 0.0360259309576058 0.00417168114426594 0.266111866388443 -0.651003114649842 0.589994534676577 0.312322955620027 0.55266992874864 -0.0962287489408765 -0.467083613141822 -0.198196785314411 0.395575343631519 0.175735234004809 0.315653622056974 -0.16114907850581 0.882040919470631 0.0637848764427299 0.858633851994301 0.822761597215925 177.266.886.690.859 -0.56698109822252 0.153261675817821 -0.27923293044209 0.0495720504717388 -116.973.330.350.604 -0.0821094023297278 0.0681021355630733 -0.201544972117822 0.0764370043961384 -0.526360142929291 -0.080413992290807 -0.390603734038536 -0.366201530340293 -0.236325627413761 -0.0197748907207431 -0.400593230734885 0.273305313443572 0.571670833577946 0.355107022563815 -0.598053247795437 0.0680124314430932 0.0825157705041395 -0.236688909377818 -0.720993035857936 0.445855527021637 -0.278184399175569 -0.813695706192547 0.0382223634028161 0.287981462997023 -0.397401294088736 0.453441800201652 -0.336993631715909 0.0350016901300973 -0.129074066212103 -0.92948724921499 0.112538313206619 -0.266214839380579 0.318919353798613 -0.235823743806211 0.312056037492047 -0.607348899940245 0.821439552875318 -0.330769416907474 -0.0368470896545833 -0.223411512848023 -0.0162393530813711 0.494328753928909
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
0 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation