Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
13900 14050 14050 14050 14150 14350 14350 14400 14400 14350 14350 14650 14550 14500 14500 14500 14500 14650 14650 15050 14950 15100 15150 15200 15200 14650 14850 14900 14650 14650 14400 14650 14450 14750 14800 14850 14800 14800 15000 14950 15000 15100 15100 15150 14850 14850 14800 14800 14600 14500 14625 14625 14600 14575 14550 14600 14550 14550 14850 15075 14950 15075 14950 15500 15425 15175 15325 15600 15650 15650 15875 15700 15700 16025 15700 16025 15900 16250 16200 16275 16275 16700 17000 16850 17000 17150 17150 17575 17850 17850 17775 17050 17425 17075 17050 16150 16200 16650 16650 16300 16325 16225 16200 15900 15900 15750 15750 15525 15575 15500 15550 14900 14400 14500 14500 14450 14300 14600 14950 14800 14900 14950 14800 14775 14775 15475 15450 15650 15675 15750 16025 16300 16750 17050 16650 16050 16175 15850 15500 15675 15775 15900 15975 15925 16175 16275 16150 16100 16000 16000 15950 16050 15800 15700 15650 15600 15350 15375 15575 15555 15600 15600 15550 15350 15300 15450 15550 15600 15750 15650 15650 15800 15050 14950 14900 14900 14500 14700 14700 14450 14550 14750 14750 14700 14950 15100 15150 15200 14950 14450 14250 14350 14650 14650 14550 14550 14650 15100 15200 15050 15000 15200 14950 15200 14950 14950 14950 15150 15100 15200 15450 15600 15800 15650 15650 15450 15750 15500 15550 15500 15450 15500 15500 15600 15550 15550 15600 15550 15500 15450 15600 15500 15450 15350 15100 15250 15200 15200 15150 15200 15250 15050 15050 15050 15050 15000 15250 15300 15400 15400 15400 15400 15250 15200 14950 14950 15050 15050 15150 15405 15450 15450 15600 15450 15650 15650 15800 15950 15950 16000 15850 15850 16050 16100 16000 16100 16100 16150 16200 16300 16300 16350 16250 16200 16350 16300 16050 16450 16450 16600 16600 16200 16200 16050 16200 16050 16050 15650 15800 15850 15750 15500 15650 15550 15700 15650 15800 15800 15850 15850 15850 15850 15850 15850 15850 15850 15900 16100 16100 16050 16050 16000 16000 16150 16200 16150 16150 16050 16100 16200 16000 16000 16000 15800 15850 15800 15750 16050 16150 16150 16150 16050 15800 15850 15850 15850 15750 15500 15600 15700 15700 15750 15550 15550 15600 15650 15600 15750 15900 15850 15700 15850 15700 15600 15650 15500 15600 15650 15700 15650 15700 15650 15550 15500 15450 15250 15300 15400 15300 15450 15300 15300 15300 15300 15450 15400 15350 15500 15400 15350 15400 15500 15650 15750 15700 15700 15700 15700 15450 15450 15450 15500 15500 15550 15450 15550 15550 15500 15650 15700 15800 15800 15800 15550 15600 15550 15650 15550 15550 15550 15550 15550 15600 15700 15700 15750 15950 16200 16200 16250 16450 16300 16300 16350 16550 16550 16550 16750 16650 16500 16550 16400 16600 16650 16550 16550 16350 16650 16700 16650 16750 16700 16850 16900 17100 17000 16950 17200 16950 16950 16950 17100 17300 17300 17500 17350 17600 17450 17700 17950 18350 18200 18100 17850 17850 17650 17700 17250 17200 17000 17000 17300 17450 17800 17450 17800 17800 17250 17200 17200 17200 17350 17450 17400 17400 17500 17600 17450 17450 17550 17550 17850 17800 17800 17950 18100 18100 18200 18200 18550 18800 18600 18600 19300 19050 19300 19300 19300 19350 19200 19100 19300 18950 19400 19450 19900 19950 20050 19900 19850 19900 19600 19400 19450 19750 19750 20250 20250 19800 19900 19800 19700 19750 20100 20450 20350 20400 20050 20350 20200 20450 20600 20650 20650 20250 20300 20200 20200 20200 20250 20375 20700 20250 20100 19150 19150 19000 19200 19250 19100 18850 18600 18150 18400 18350 18600 18650 18500 18775 18625 18600 18760 18700 18900 18950 18500 18450 18350 18200 18075 18200 18150 18150 17575 17850 17950 18050 18050 18050 18250 18050 17900 17800 18125 18550 18600 18750 18850 18600 18750 18750 18450 18450 18400 18200 18300 18200 18150 18025 18200 18200 18175 17950 18025 18050 18150 18150 18300 18475 18550 18175 18200 18025 18500 18750 18700 19000 19000 18950 18850 18850 18700 18850 18850 19350 19450 19500 19400 19300 19300 19400 19400 18800 18925 18900 18900 18750 18650 18625 18700 18225 18100 18250 18200 18350 17350 17450 17825 17050 17325 17125 17600 17700 17875 17775 17775 17950 17950 18050 18100 17700 17700 17650 17800 18050 18050 17550 17550 17150 17300 17300 17550 19150 19150 19100 19300 19200 19300 19075 19300 20050 19900 19900 19400 19450 20350 20850 20850 21500 20350 19650 19850 20100 19150 18950 19000 18800 18225 17650 18650 19175 18700 18550 18300 18200 18400 18400 18700 18600 18760 18650 18325 18300 18850 18650 18650 18875 19100 19850 20250 20000 20125 20100 20400 19450 19600 19325 19200 19150 19125 19225 19850 19450 19650 19450 19300 19125 19050 19250 19250 19250 19200 19850 19900 19650 19650 20075 20050 20050 20350 19775 20000 20300 19700 19850 19675 19800 20100 20350 21000 20900 21150 22100 21900 21550 21475 22450 22750 22450 22000 22550 22850 22450 22250 22550 23275 23250 23550 24150 24650 24725 24750 24400 24800 23750 23750 23925 24275 23250 23250 23100 23600 23600 23000 22600 22825 22800 22800 22500 22200 22100 22350 22250 21950 21900 21950 21950 22350 22100 22200 22100 21450 20700 20900 21450 21450 21100 21450 21450 21500 21250 21300 21800 21750 21875 22050 22200 21950 21950 21625 21200 21350 21500 21600 21850 21700 21400 21350 21550 21650 21650 22100 21650 21650 21750 21750 21700 21750 21650 21850 21950 21950 21950 22000 21675 22000 21850 21800 21700 21700 21600 21400 21500 21500 21425 21350 21400 21050 20900 20900 21350 21250 21300 21050 21200 21200 21200 21100 21150 20850 20850 20800 20950 21100 21200 21150 21250 21250 21250 21250 21300 21350 21300 21100 21025 21150 21350 21100 21350 21300 21300 21450 21350 21300 21250 21325 21225 21175 21125 21125 21025 21225 21050 21075 21125 21075 21300 21150 21400 21800 22150 22050 22150 22000 21750 21850 21900 21800 21950 22050 21950 21750 21800 21900 21800 21550 21600 21600 21700 21950 21350 21850 22050 22550 22800 22700 22700 22725 22850 22600 22350 22475 22650 22550 22450 22500 22325 22250 22250 22350 22500 22550 22875 23400 23300 23375 23550 23400 23550 23800 23675 24000 24200 24350 24375 24450 24800 24900 25050 25075 24550 24850 25050 25650 25700 25400 24750 24900 24450 24400 24550 24450 24400 24650 25000 24500 24750 24400 24250 24475 24600 24450 24300 24375 24375 24550 24725 24825 25100 24950 25325 25325 24800 24975 25125 25125 25125 25400 25175 24650 24775 24675 24825 24775 24675 24750 24875 25400 25400 24750 24900 24825 24875 24975 25375 25600 26000 25900 25850 26075 26275 26050 26000 25825 26075 26150 26275 26475 26500 26575 26425 26275 26375 25900 25850 25625 25900 26050 26150 26275 26100 25975 25975 26125 26175 26225 26225 26200 26275 26275 26275 26750 27075 27475 27525 27125 27000 26950 27075 27150 26875 26925 27150 27150 27425 27625 27475 28075 28075 28175 28350 28350 28500 29350 30225 29575 30125 30125 31150 31350 32175 31725 31600 30800 30800 29700 30875 31275 31500 31375 31400 31650 31975 31650 31975 32575 32025 33050 32300 32100 32250 32050 31975 32100 32025 32275 32100 32275 31975 32175 32375 32300 32450 32425 30800 30850 30750 30175 30350 30125 30625 30375 30425 30325 29825 29450 29100 29450 29550 29575 29425 29050 28525 28575 28500 28875 28625 28625 28925 28925 28950 28950 29100 29700 30000 30400 30375 30425 30625 30700 30825 30800 31100 31175 31025 30975 31025 31350 31075 31125 30900 31150 31575 31575 31375 31100 30975 31200 31125 31075 31275 31175 30950 30725 30900 30700 30625 30700 30650 30525 30850 30725 31025 30975 30550 30900 31000 31000 31000 31000 31325 31000 31000 30300 30575 30575 30775 30550 30750 31025 31000 30850 30600 31150 31800 32500 32325 31800 31850 31625 31750 31650 31525 32075 32725 32900 32775 32825 33200 34100 33800 33525 33775 34000 33425 33550 33400 33300 33400 33000 33500 33550 33725 33700 33600 33550 33500 34200 34000 33600
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