R version 2.12.0 (2010-10-15) Copyright (C) 2010 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-pc-linux-gnu (32-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- c(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) > ylimmax = '' > ylimmin = '' > main = 'Robustness of Central Tendency' > 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)) [1] 20452.11 > sqrtn <- sqrt(length(x)) > (armse <- sd(x) / sqrtn) [1] 148.8478 > (armose <- arm / armse) [1] 137.4029 > (geo <- geomean(x)) [1] 19833.33 > (har <- harmean(x)) [1] 19289.67 > (qua <- quamean(x)) [1] 21133.53 > (win <- winmean(x)) [,1] [,2] [1,] 20452.15 148.83819 [2,] 20452.00 148.82702 [3,] 20452.00 148.82702 [4,] 20451.68 148.77243 [5,] 20451.98 148.75275 [6,] 20451.98 148.72876 [7,] 20452.11 148.71042 [8,] 20451.49 148.66704 [9,] 20451.49 148.66704 [10,] 20451.10 148.64008 [11,] 20451.10 148.64008 [12,] 20451.57 148.62510 [13,] 20451.31 148.60762 [14,] 20451.04 148.58883 [15,] 20451.04 148.58883 [16,] 20450.73 148.50480 [17,] 20450.39 148.48214 [18,] 20450.39 148.48214 [19,] 20448.91 148.38127 [20,] 20448.13 148.25122 [21,] 20445.67 148.08659 [22,] 20444.81 148.02950 [23,] 20443.01 147.91076 [24,] 20441.61 147.81841 [25,] 20440.63 147.75455 [26,] 20439.61 147.68838 [27,] 20436.45 147.48368 [28,] 20435.90 147.34409 [29,] 20434.77 147.27170 [30,] 20434.18 147.23436 [31,] 20432.97 147.15739 [32,] 20431.72 147.07822 [33,] 20431.08 147.03750 [34,] 20431.08 147.03750 [35,] 20431.08 146.97306 [36,] 20431.78 146.95120 [37,] 20431.06 146.90570 [38,] 20428.83 146.76601 [39,] 20428.83 146.76601 [40,] 20427.27 146.59551 [41,] 20427.27 146.59551 [42,] 20428.09 146.57017 [43,] 20427.25 146.51793 [44,] 20426.39 146.46456 [45,] 20425.51 146.41007 [46,] 20425.51 146.41007 [47,] 20423.68 146.29653 [48,] 20423.68 146.29653 [49,] 20423.68 146.29653 [50,] 20423.68 146.29653 [51,] 20418.70 145.99035 [52,] 20418.70 145.80366 [53,] 20418.70 145.80366 [54,] 20416.59 145.67514 [55,] 20417.66 145.54417 [56,] 20414.38 145.34536 [57,] 20414.38 145.34536 [58,] 20415.51 145.31091 [59,] 20414.36 145.24131 [60,] 20414.36 145.13513 [61,] 20413.17 145.06340 [62,] 20413.17 145.06340 [63,] 20410.71 144.91559 [64,] 20409.46 144.84069 [65,] 20404.38 144.53766 [66,] 20403.09 144.46105 [67,] 20405.71 144.38200 [68,] 20404.38 144.30317 [69,] 20404.38 144.30317 [70,] 20403.01 144.22215 [71,] 20400.24 144.05820 [72,] 20403.05 143.97381 [73,] 20398.77 143.72187 [74,] 20397.33 143.63701 [75,] 20397.33 143.63701 [76,] 20395.84 143.55000 [77,] 20398.85 143.46032 [78,] 20398.85 143.46032 [79,] 20397.31 143.36998 [80,] 20397.31 143.36998 [81,] 20395.73 143.27750 [82,] 20395.73 143.27750 [83,] 20394.11 143.18288 [84,] 20394.11 143.18288 [85,] 20390.79 142.98952 [86,] 20390.79 142.98952 [87,] 20390.79 142.98952 [88,] 20390.79 142.98952 [89,] 20389.05 142.88852 [90,] 20389.05 142.88852 [91,] 20389.05 142.88852 [92,] 20392.64 142.78198 [93,] 20392.64 142.78198 [94,] 20392.64 142.78198 [95,] 20392.64 142.78198 [96,] 20394.52 142.56280 [97,] 20394.52 142.56280 [98,] 20394.52 142.56280 [99,] 20392.58 142.45068 [100,] 20388.68 142.22464 [101,] 20388.68 142.22464 [102,] 20388.68 142.22464 [103,] 20386.66 142.10847 [104,] 20384.63 141.99134 [105,] 20386.68 141.93122 [106,] 20386.68 141.93122 [107,] 20386.68 141.74978 [108,] 20384.57 141.62841 [109,] 20384.57 141.62841 [110,] 20384.57 141.62841 [111,] 20384.57 141.62841 [112,] 20382.39 141.50272 [113,] 20380.18 141.37606 [114,] 20384.63 141.24653 [115,] 20382.39 141.11774 [116,] 20382.39 141.11774 [117,] 20380.10 140.98687 [118,] 20380.10 140.98687 [119,] 20380.10 140.98687 [120,] 20377.76 140.65102 [121,] 20377.76 140.44791 [122,] 20377.76 140.44791 [123,] 20377.76 140.44791 [124,] 20375.34 140.30984 [125,] 20372.89 140.17082 [126,] 20372.89 140.17082 [127,] 20370.41 140.02975 [128,] 20370.41 140.02975 [129,] 20367.89 139.88663 [130,] 20357.74 139.31142 [131,] 20357.74 139.31142 [132,] 20360.32 139.01707 [133,] 20357.72 138.87048 [134,] 20357.72 138.87048 [135,] 20355.08 138.72186 [136,] 20352.43 138.57231 [137,] 20355.10 138.26827 [138,] 20347.02 137.81434 [139,] 20341.59 137.51041 [140,] 20336.12 137.20500 [141,] 20336.12 137.20500 [142,] 20336.12 137.20500 [143,] 20322.15 136.42824 [144,] 20302.46 135.34074 [145,] 20291.14 134.48214 [146,] 20293.99 134.39997 [147,] 20279.63 133.61680 [148,] 20279.63 133.61680 [149,] 20276.72 133.45860 [150,] 20265.00 132.82348 [151,] 20267.95 132.73831 [152,] 20267.95 132.49240 [153,] 20258.99 132.00883 [154,] 20228.91 130.39907 [155,] 20228.91 130.39907 [156,] 20222.82 130.07534 [157,] 20210.55 129.42616 [158,] 20210.55 129.42616 [159,] 20207.45 129.26230 [160,] 20208.07 129.24414 [161,] 20204.30 128.83983 [162,] 20175.82 127.09314 [163,] 20175.82 127.09314 [164,] 20169.41 126.76058 [165,] 20162.97 126.42684 [166,] 20159.73 126.25929 [167,] 20159.73 126.25929 [168,] 20140.04 125.24636 [169,] 20140.04 125.24636 [170,] 20116.80 124.06056 [171,] 20103.44 123.38391 [172,] 20103.44 123.38391 [173,] 20042.62 120.34930 [174,] 20029.02 119.68164 [175,] 20022.19 119.34738 [176,] 20022.19 119.34738 [177,] 20015.27 119.01030 [178,] 19977.03 117.16411 [179,] 19980.53 117.05870 [180,] 19984.04 116.95310 [181,] 19980.51 116.78388 [182,] 19973.40 116.44435 [183,] 19973.40 116.44435 [184,] 19962.62 115.93146 [185,] 19955.39 115.58901 [186,] 19951.76 115.41727 [187,] 19944.45 115.07275 [188,] 19926.09 114.21165 [189,] 19900.25 113.01144 [190,] 19889.12 112.49863 [191,] 19885.39 112.32738 [192,] 19877.89 111.98396 [193,] 19870.35 111.63993 [194,] 19855.20 110.95187 [195,] 19855.20 110.95187 [196,] 19855.20 110.95187 [197,] 19855.20 110.95187 [198,] 19859.06 110.83341 [199,] 19862.95 110.71481 [200,] 19862.95 110.71481 [201,] 19855.10 110.35987 [202,] 19855.10 110.35987 [203,] 19851.13 110.18110 [204,] 19847.15 110.00174 [205,] 19843.14 109.82181 [206,] 19843.14 109.82181 [207,] 19839.10 109.64043 [208,] 19835.04 109.45848 [209,] 19830.96 109.27596 [210,] 19830.96 109.27596 [211,] 19826.84 109.09200 [212,] 19826.84 109.09200 [213,] 19818.52 108.72152 [214,] 19818.52 108.72152 [215,] 19814.32 108.53500 [216,] 19814.32 108.53500 [217,] 19801.60 107.97211 [218,] 19801.60 107.97211 [219,] 19801.60 107.97211 [220,] 19793.87 107.56674 [221,] 19797.32 107.46050 [222,] 19792.98 107.26971 [223,] 19775.56 106.18205 [224,] 19766.81 105.80098 [225,] 19762.42 105.61009 [226,] 19758.00 105.41867 [227,] 19722.54 103.89246 [228,] 19722.54 103.89246 [229,] 19722.54 103.89246 [230,] 19722.54 103.89246 [231,] 19718.02 103.69985 [232,] 19708.96 103.31400 [233,] 19708.96 103.31400 [234,] 19681.54 102.15476 [235,] 19672.36 101.76950 [236,] 19672.36 101.76950 [237,] 19672.36 101.76950 [238,] 19667.71 101.57496 [239,] 19663.04 101.37994 [240,] 19667.73 100.89153 [241,] 19667.73 100.89153 [242,] 19658.28 100.49796 [243,] 19653.53 100.30090 [244,] 19653.53 100.30090 [245,] 19648.75 100.10259 [246,] 19639.14 99.70543 [247,] 19639.14 99.70543 [248,] 19639.14 99.70543 [249,] 19634.27 99.50498 [250,] 19634.27 99.50498 [251,] 19629.37 99.30329 [252,] 19624.45 99.10116 [253,] 19624.45 99.10116 [254,] 19624.45 99.10116 [255,] 19619.47 98.89700 [256,] 19619.47 98.89700 [257,] 19619.47 98.89700 [258,] 19614.43 98.69080 [259,] 19614.43 98.69080 [260,] 19609.35 98.48337 [261,] 19614.45 98.32341 [262,] 19614.45 98.32341 [263,] 19619.59 98.16277 [264,] 19619.59 98.16277 [265,] 19614.41 97.95145 [266,] 19614.41 97.95145 [267,] 19603.98 97.52671 [268,] 19603.98 97.52671 [269,] 19598.73 97.31330 [270,] 19598.73 97.31330 [271,] 19598.73 97.31330 [272,] 19588.10 96.88282 [273,] 19577.44 96.45223 [274,] 19577.44 96.45223 [275,] 19577.44 96.45223 [276,] 19566.66 96.01839 [277,] 19561.25 95.80124 [278,] 19555.82 95.58369 [279,] 19555.82 95.58369 [280,] 19555.82 95.58369 [281,] 19566.79 95.23861 [282,] 19555.78 94.79776 [283,] 19555.78 94.79776 [284,] 19555.78 94.79776 [285,] 19555.78 94.79776 [286,] 19550.19 94.57476 [287,] 19550.19 94.57476 [288,] 19550.19 94.57476 [289,] 19544.55 94.34980 [290,] 19533.22 93.89943 [291,] 19527.54 93.67403 [292,] 19521.83 93.44823 [293,] 19516.11 92.81558 [294,] 19510.37 92.18221 [295,] 19475.80 90.82879 [296,] 19458.45 90.15515 [297,] 19429.45 89.03680 [298,] 19417.81 88.59089 [299,] 19417.81 88.59089 [300,] 19400.23 87.92073 [301,] 19382.59 87.25223 [302,] 19382.59 87.25223 [303,] 19370.76 86.80581 [304,] 19370.76 86.80581 [305,] 19370.76 86.80581 [306,] 19334.90 85.46404 [307,] 19346.89 85.07743 [308,] 19340.88 84.85359 [309,] 19322.77 84.18268 [310,] 19316.71 83.95926 [311,] 19310.64 83.73560 [312,] 19310.64 83.73560 [313,] 19310.64 83.73560 [314,] 19273.84 82.39075 [315,] 19249.23 81.50121 [316,] 19218.38 80.39729 [317,] 19212.18 80.17739 [318,] 19212.18 80.17739 [319,] 19205.95 79.95663 [320,] 19199.70 79.73572 [321,] 19199.70 79.73572 [322,] 19199.70 79.73572 [323,] 19187.09 79.29136 [324,] 19180.76 79.06931 [325,] 19180.76 78.63818 [326,] 19187.12 78.42922 [327,] 19174.35 77.98299 [328,] 19161.54 77.53763 [329,] 19161.54 77.53763 [330,] 19148.65 77.09182 [331,] 19148.65 77.09182 [332,] 19155.13 76.87849 [333,] 19161.64 76.66514 [334,] 19161.64 76.66514 [335,] 19148.55 76.21408 [336,] 19148.55 76.21408 [337,] 19155.13 75.99842 [338,] 19155.13 75.55568 [339,] 19148.51 75.32851 [340,] 19148.51 75.32851 [341,] 19148.51 75.32851 [342,] 19148.51 75.32851 [343,] 19121.71 74.41489 [344,] 19121.71 74.41489 [345,] 19121.71 74.41489 [346,] 19121.71 74.41489 [347,] 19121.71 73.96313 [348,] 19101.32 73.27535 [349,] 19101.32 73.27535 [350,] 19108.16 73.05128 [351,] 19108.16 73.05128 [352,] 19094.41 72.59012 [353,] 19094.41 72.59012 [354,] 19094.41 72.59012 [355,] 19080.54 72.12749 [356,] 19080.54 72.12749 [357,] 19066.60 71.66473 [358,] 19066.60 71.66473 [359,] 19066.60 71.66473 [360,] 19066.60 71.66473 [361,] 19066.60 71.66473 [362,] 19066.60 70.73303 [363,] 19066.60 70.73303 [364,] 19066.60 70.73303 [365,] 19066.60 70.73303 [366,] 19052.30 70.26281 [367,] 19052.30 70.26281 [368,] 19052.30 70.26281 [369,] 19066.71 69.79055 [370,] 19052.26 69.31687 [371,] 19052.26 69.31687 [372,] 19052.26 69.31687 [373,] 19052.26 69.31687 [374,] 19052.26 69.31687 [375,] 19052.26 69.31687 [376,] 19052.26 69.31687 [377,] 19052.26 69.31687 [378,] 19059.64 69.07549 [379,] 19059.64 69.07549 [380,] 19067.07 68.83348 [381,] 19052.18 68.34737 [382,] 19052.18 68.34737 [383,] 19052.18 68.34737 [384,] 19052.18 68.34737 [385,] 19044.66 68.10271 [386,] 19037.12 67.85806 [387,] 19037.12 67.85806 [388,] 19037.12 67.85806 [389,] 19037.12 67.85806 [390,] 19037.12 67.85806 [391,] 19021.85 67.36430 [392,] 19029.51 67.11385 [393,] 19037.18 66.86338 [394,] 19037.18 66.86338 [395,] 19037.18 66.86338 [396,] 19044.92 66.61161 [397,] 19029.41 66.11138 [398,] 19029.41 66.11138 [399,] 19037.20 65.85783 [400,] 19037.20 65.85783 [401,] 19037.20 65.85783 [402,] 19037.20 65.85783 [403,] 19021.46 65.35209 [404,] 19021.46 65.35209 [405,] 19021.46 65.35209 [406,] 19029.39 65.09421 [407,] 19037.34 64.83630 [408,] 19029.37 64.58075 [409,] 19021.38 64.32520 [410,] 19021.38 64.32520 [411,] 19037.44 63.80511 [412,] 19053.53 63.28621 [413,] 19053.53 63.28621 [414,] 19045.45 63.02715 [415,] 19053.55 62.24688 [416,] 19069.80 61.72695 [417,] 19069.80 61.72695 [418,] 19069.80 61.72695 [419,] 19053.43 61.20367 [420,] 19053.43 61.20367 [421,] 19053.43 61.20367 [422,] 19053.43 60.15225 [423,] 19053.43 60.15225 [424,] 19053.43 60.15225 [425,] 19070.04 59.62411 [426,] 19070.04 59.62411 > (tri <- trimean(x)) [,1] [,2] [1,] 20446.48 148.60334 [2,] 20440.80 148.36610 [3,] 20435.17 148.13211 [4,] 20429.52 147.89574 [5,] 20423.94 147.67094 [6,] 20418.28 147.44787 [7,] 20412.60 147.22664 [8,] 20406.89 147.00584 [9,] 20401.23 146.78842 [10,] 20395.56 146.56875 [11,] 20389.91 146.34965 [12,] 20384.24 146.12827 [13,] 20378.51 145.90590 [14,] 20372.79 145.68265 [15,] 20367.06 145.45852 [16,] 20361.32 145.23204 [17,] 20355.58 145.00887 [18,] 20349.84 144.78482 [19,] 20344.09 144.55842 [20,] 20338.39 144.33549 [21,] 20332.72 144.11739 [22,] 20327.15 143.90572 [23,] 20321.60 143.69472 [24,] 20316.12 143.48728 [25,] 20310.67 143.28197 [26,] 20305.26 143.07738 [27,] 20299.86 142.87351 [28,] 20294.57 142.67613 [29,] 20289.28 142.48231 [30,] 20284.02 142.28929 [31,] 20278.76 142.09563 [32,] 20273.52 141.90278 [33,] 20268.31 141.71073 [34,] 20263.10 141.51805 [35,] 20257.88 141.32329 [36,] 20252.63 141.12856 [37,] 20247.35 140.93241 [38,] 20242.07 140.73560 [39,] 20236.84 140.54106 [40,] 20231.59 140.34440 [41,] 20226.36 140.15073 [42,] 20221.12 139.95493 [43,] 20215.83 139.75769 [44,] 20210.55 139.55977 [45,] 20205.28 139.36117 [46,] 20200.00 139.16189 [47,] 20194.71 138.96040 [48,] 20189.45 138.75973 [49,] 20184.16 138.55683 [50,] 20178.86 138.35169 [51,] 20173.54 138.14426 [52,] 20168.31 137.94223 [53,] 20163.06 137.74248 [54,] 20157.79 137.54048 [55,] 20152.55 137.33933 [56,] 20147.26 137.13891 [57,] 20142.03 136.94090 [58,] 20136.77 136.74066 [59,] 20131.48 136.53888 [60,] 20126.19 136.33639 [61,] 20120.88 136.13392 [62,] 20115.58 135.93073 [63,] 20110.25 135.72523 [64,] 20104.95 135.52058 [65,] 20099.66 135.31519 [66,] 20094.43 135.11384 [67,] 20089.21 134.91178 [68,] 20083.92 134.70890 [69,] 20078.64 134.50531 [70,] 20073.34 134.29936 [71,] 20068.04 134.09268 [72,] 20062.77 133.88686 [73,] 20057.43 133.68017 [74,] 20052.15 133.47600 [75,] 20046.86 133.27109 [76,] 20041.56 133.06381 [77,] 20036.26 132.85576 [78,] 20030.90 132.64681 [79,] 20025.52 132.43542 [80,] 20020.14 132.22323 [81,] 20014.74 132.00855 [82,] 20009.35 131.79305 [83,] 20003.93 131.57501 [84,] 19998.52 131.35611 [85,] 19993.09 131.13464 [86,] 19987.69 130.91399 [87,] 19982.26 130.69073 [88,] 19976.82 130.46483 [89,] 19971.35 130.23626 [90,] 19965.89 130.00672 [91,] 19960.41 129.77446 [92,] 19954.91 129.53943 [93,] 19949.34 129.30320 [94,] 19943.75 129.06414 [95,] 19938.15 128.82222 [96,] 19932.52 128.57740 [97,] 19926.85 128.33307 [98,] 19921.15 128.08580 [99,] 19915.44 127.83554 [100,] 19909.73 127.58410 [101,] 19904.04 127.33334 [102,] 19898.33 127.07952 [103,] 19892.60 126.82262 [104,] 19886.88 126.56448 [105,] 19881.15 126.30508 [106,] 19875.38 126.04334 [107,] 19869.59 125.77838 [108,] 19863.77 125.51294 [109,] 19857.96 125.24616 [110,] 19852.13 124.97609 [111,] 19846.27 124.70266 [112,] 19840.39 124.42585 [113,] 19834.52 124.14757 [114,] 19828.64 123.86782 [115,] 19822.70 123.58630 [116,] 19816.75 123.30327 [117,] 19810.78 123.01668 [118,] 19804.82 122.72852 [119,] 19798.83 122.43672 [120,] 19792.82 122.14123 [121,] 19786.81 121.84702 [122,] 19780.77 121.55203 [123,] 19774.71 121.25327 [124,] 19768.63 120.95070 [125,] 19762.55 120.64637 [126,] 19756.47 120.34025 [127,] 19750.37 120.03020 [128,] 19744.27 119.71829 [129,] 19738.14 119.40234 [130,] 19732.01 119.08447 [131,] 19725.96 118.77121 [132,] 19719.89 118.45387 [133,] 19713.76 118.13643 [134,] 19707.64 117.81705 [135,] 19701.49 117.49347 [136,] 19695.34 117.16787 [137,] 19689.19 116.84022 [138,] 19683.00 116.51237 [139,] 19676.85 116.18698 [140,] 19670.73 115.86180 [141,] 19664.63 115.53686 [142,] 19658.51 115.20757 [143,] 19652.37 114.87387 [144,] 19646.33 114.54726 [145,] 19640.43 114.23233 [146,] 19634.62 113.92555 [147,] 19628.76 113.61563 [148,] 19623.00 113.31297 [149,] 19617.21 113.00624 [150,] 19611.43 112.69769 [151,] 19605.73 112.39417 [152,] 19599.98 112.08756 [153,] 19594.20 111.78010 [154,] 19588.48 111.47544 [155,] 19582.99 111.18971 [156,] 19577.48 110.90011 [157,] 19572.00 110.61118 [158,] 19566.60 110.32752 [159,] 19561.18 110.04001 [160,] 19555.76 109.75089 [161,] 19550.31 109.45802 [162,] 19544.87 109.16665 [163,] 19539.65 108.89558 [164,] 19534.40 108.62078 [165,] 19529.18 108.34685 [166,] 19524.00 108.07380 [167,] 19518.82 107.79931 [168,] 19513.61 107.52104 [169,] 19508.55 107.25297 [170,] 19503.46 106.98119 [171,] 19498.53 106.72199 [172,] 19493.70 106.46846 [173,] 19488.84 106.21141 [174,] 19484.44 105.99172 [175,] 19480.13 105.77785 [176,] 19475.86 105.56542 [177,] 19471.57 105.35004 [178,] 19467.32 105.13611 [179,] 19463.34 104.94333 [180,] 19459.32 104.74898 [181,] 19455.26 104.55303 [182,] 19451.20 104.35651 [183,] 19447.18 104.16161 [184,] 19443.15 103.96396 [185,] 19439.18 103.77015 [186,] 19435.24 103.57800 [187,] 19431.32 103.38533 [188,] 19427.43 103.19435 [189,] 19423.67 103.01166 [190,] 19420.08 102.84160 [191,] 19416.56 102.67553 [192,] 19413.06 102.50920 [193,] 19409.59 102.34477 [194,] 19406.17 102.18225 [195,] 19402.84 102.02592 [196,] 19399.49 101.86730 [197,] 19396.13 101.70635 [198,] 19392.76 101.54303 [199,] 19389.34 101.37853 [200,] 19385.88 101.21284 [201,] 19382.40 101.04469 [202,] 19378.97 100.87850 [203,] 19375.51 100.70984 [204,] 19372.08 100.54091 [205,] 19368.65 100.37171 [206,] 19365.24 100.20224 [207,] 19361.81 100.03024 [208,] 19358.39 99.85794 [209,] 19354.99 99.68535 [210,] 19351.60 99.51245 [211,] 19348.19 99.33695 [212,] 19344.80 99.16113 [213,] 19341.39 98.98265 [214,] 19338.03 98.80615 [215,] 19334.65 98.62698 [216,] 19331.28 98.44744 [217,] 19327.90 98.26518 [218,] 19324.59 98.08726 [219,] 19321.26 97.90663 [220,] 19317.92 97.72324 [221,] 19314.61 97.54212 [222,] 19311.27 97.35932 [223,] 19307.94 97.17613 [224,] 19304.71 97.00355 [225,] 19301.53 96.83307 [226,] 19298.36 96.66234 [227,] 19295.21 96.49135 [228,] 19292.29 96.33669 [229,] 19289.35 96.17959 [230,] 19286.40 96.02001 [231,] 19283.43 95.85791 [232,] 19280.48 95.69564 [233,] 19277.57 95.53560 [234,] 19274.66 95.37302 [235,] 19271.91 95.22222 [236,] 19269.21 95.07374 [237,] 19266.50 94.92286 [238,] 19263.77 94.76955 [239,] 19261.06 94.61617 [240,] 19258.37 94.46271 [241,] 19255.63 94.31222 [242,] 19252.88 94.15927 [243,] 19250.18 94.00869 [244,] 19247.50 93.85809 [245,] 19244.80 93.70501 [246,] 19242.13 93.55187 [247,] 19239.50 93.40114 [248,] 19236.86 93.24792 [249,] 19234.20 93.09216 [250,] 19231.56 92.93632 [251,] 19228.91 92.77790 [252,] 19226.28 92.61938 [253,] 19223.67 92.46076 [254,] 19221.04 92.29949 [255,] 19218.40 92.13553 [256,] 19215.78 91.97141 [257,] 19213.15 91.80456 [258,] 19210.50 91.63491 [259,] 19207.87 91.46504 [260,] 19205.22 91.29233 [261,] 19202.60 91.11937 [262,] 19199.93 90.94521 [263,] 19197.24 90.76810 [264,] 19194.51 90.58972 [265,] 19191.76 90.40831 [266,] 19189.03 90.22658 [267,] 19186.29 90.04175 [268,] 19183.60 89.85933 [269,] 19180.89 89.67379 [270,] 19178.20 89.48788 [271,] 19175.50 89.29879 [272,] 19172.79 89.10645 [273,] 19170.12 88.91651 [274,] 19167.51 88.72900 [275,] 19164.89 88.53824 [276,] 19162.25 88.34418 [277,] 19159.67 88.15254 [278,] 19157.11 87.96046 [279,] 19154.56 87.76795 [280,] 19152.01 87.57207 [281,] 19149.44 87.37277 [282,] 19146.78 87.17382 [283,] 19144.18 86.97732 [284,] 19141.57 86.77736 [285,] 19138.94 86.57387 [286,] 19136.29 86.36677 [287,] 19133.67 86.15906 [288,] 19131.03 85.94766 [289,] 19128.38 85.73250 [290,] 19125.74 85.51662 [291,] 19123.17 85.30315 [292,] 19120.61 85.08899 [293,] 19118.08 84.87413 [294,] 19115.56 84.66378 [295,] 19113.07 84.45801 [296,] 19110.78 84.26719 [297,] 19108.59 84.08211 [298,] 19106.57 83.90879 [299,] 19104.61 83.73828 [300,] 19102.64 83.56461 [301,] 19100.77 83.39674 [302,] 19098.99 83.23471 [303,] 19097.21 83.06963 [304,] 19095.49 82.90745 [305,] 19093.76 82.74219 [306,] 19092.02 82.57382 [307,] 19090.50 82.42028 [308,] 19088.89 82.26831 [309,] 19087.30 82.11639 [310,] 19085.83 81.97049 [311,] 19084.38 81.82470 [312,] 19082.96 81.67904 [313,] 19081.53 81.53053 [314,] 19080.09 81.37910 [315,] 19078.88 81.24263 [316,] 19077.81 81.11511 [317,] 19076.93 80.99929 [318,] 19076.08 80.88390 [319,] 19075.23 80.76608 [320,] 19074.41 80.64867 [321,] 19073.62 80.53165 [322,] 19072.83 80.41216 [323,] 19072.03 80.29017 [324,] 19071.31 80.17141 [325,] 19070.63 80.05302 [326,] 19069.94 79.93756 [327,] 19069.20 79.82216 [328,] 19068.54 79.71010 [329,] 19067.96 79.60138 [330,] 19067.37 79.49026 [331,] 19066.86 79.38249 [332,] 19066.35 79.27232 [333,] 19065.79 79.16235 [334,] 19065.19 79.05257 [335,] 19064.58 78.94033 [336,] 19064.05 78.83152 [337,] 19063.52 78.72026 [338,] 19062.95 78.60919 [339,] 19062.37 78.50132 [340,] 19061.83 78.39399 [341,] 19061.28 78.28418 [342,] 19060.73 78.17185 [343,] 19060.18 78.05693 [344,] 19059.79 77.95163 [345,] 19059.40 77.84385 [346,] 19059.01 77.73354 [347,] 19058.61 77.62064 [348,] 19058.21 77.51110 [349,] 19057.94 77.40822 [350,] 19057.66 77.30285 [351,] 19057.34 77.19787 [352,] 19057.02 77.09034 [353,] 19056.79 76.98647 [354,] 19056.55 76.88006 [355,] 19056.31 76.77104 [356,] 19056.15 76.66568 [357,] 19056.00 76.55771 [358,] 19055.93 76.45341 [359,] 19055.86 76.34650 [360,] 19055.79 76.23693 [361,] 19055.73 76.12463 [362,] 19055.66 76.00955 [363,] 19055.59 75.90453 [364,] 19055.52 75.79684 [365,] 19055.45 75.68641 [366,] 19055.37 75.57318 [367,] 19055.39 75.46371 [368,] 19055.41 75.35143 [369,] 19055.43 75.23629 [370,] 19055.36 75.12492 [371,] 19055.38 75.01741 [372,] 19055.40 74.90708 [373,] 19055.42 74.79387 [374,] 19055.44 74.67772 [375,] 19055.46 74.55855 [376,] 19055.48 74.43630 [377,] 19055.50 74.31089 [378,] 19055.52 74.18226 [379,] 19055.50 74.05392 [380,] 19055.47 73.92226 [381,] 19055.40 73.79085 [382,] 19055.42 73.66336 [383,] 19055.44 73.53253 [384,] 19055.46 73.39829 [385,] 19055.48 73.26054 [386,] 19055.55 73.12297 [387,] 19055.67 72.98556 [388,] 19055.79 72.84454 [389,] 19055.92 72.69982 [390,] 19056.04 72.55131 [391,] 19056.16 72.39891 [392,] 19056.39 72.25034 [393,] 19056.57 72.10189 [394,] 19056.70 71.95356 [395,] 19056.83 71.80128 [396,] 19056.96 71.64494 [397,] 19057.04 71.48860 [398,] 19057.22 71.33623 [399,] 19057.41 71.17975 [400,] 19057.54 71.02329 [401,] 19057.68 70.86259 [402,] 19057.82 70.69755 [403,] 19057.95 70.52805 [404,] 19058.20 70.36250 [405,] 19058.45 70.19242 [406,] 19058.70 70.01772 [407,] 19058.89 69.84278 [408,] 19059.04 69.66761 [409,] 19059.24 69.49206 [410,] 19059.50 69.31611 [411,] 19059.76 69.13528 [412,] 19059.91 68.95881 [413,] 19059.96 68.78673 [414,] 19060.00 68.60978 [415,] 19060.10 68.43250 [416,] 19060.15 68.26432 [417,] 19060.08 68.10068 [418,] 19060.01 67.93228 [419,] 19059.94 67.75898 [420,] 19059.99 67.59029 [421,] 19060.03 67.41664 [422,] 19060.08 67.23791 [423,] 19060.13 67.07357 [424,] 19060.17 66.90432 [425,] 19060.22 66.73000 [426,] 19060.15 66.56048 > (midr <- midrange(x)) [1] 24050 > midm <- array(NA,dim=8) > for (j in 1:8) midm[j] <- midmean(x,j) > midm [1] 19010.21 19010.21 19010.21 19010.21 19010.21 19010.21 19010.21 19016.05 > postscript(file="/var/www/rcomp/tmp/19coh1324588320.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > 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() null device 1 > postscript(file="/var/www/rcomp/tmp/24rjv1324588320.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > 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() null device 1 > > #Note: the /var/www/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/rcomp/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="/var/www/rcomp/tmp/396eq1324588320.tab") > > try(system("convert tmp/19coh1324588320.ps tmp/19coh1324588320.png",intern=TRUE)) character(0) > try(system("convert tmp/24rjv1324588320.ps tmp/24rjv1324588320.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 13.920 0.140 14.052