R version 2.8.1 (2008-12-22) Copyright (C) 2008 The R Foundation for Statistical Computing ISBN 3-900051-07-0 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(3.594 + ,3.617 + ,3.647 + ,3.761 + ,3.767 + ,3.692 + ,3.709 + ,3.716 + ,3.755 + ,3.702 + ,3.698 + ,3.765 + ,3.809 + ,3.777 + ,3.747 + ,3.695 + ,3.719 + ,3.714 + ,3.740 + ,3.632 + ,3.643 + ,3.615 + ,3.747 + ,3.687 + ,3.700 + ,3.578 + ,3.608 + ,3.711 + ,3.599 + ,3.721 + ,3.683 + ,3.668 + ,3.715 + ,3.673 + ,3.701 + ,3.846 + ,3.853 + ,3.807 + ,3.860 + ,3.854 + ,3.842 + ,3.826 + ,3.863 + ,3.763 + ,3.789 + ,3.822 + ,3.877 + ,3.867 + ,3.940 + ,3.864 + ,3.859 + ,3.868 + ,3.887 + ,3.908 + ,3.936 + ,3.891 + ,3.962 + ,3.947 + ,3.909 + ,3.897 + ,3.878 + ,3.841 + ,3.841 + ,3.842 + ,3.836 + ,3.821 + ,3.815 + ,3.805 + ,3.815 + ,3.803 + ,3.722 + ,3.745 + ,3.681 + ,3.684 + ,3.676 + ,3.725 + ,3.734 + ,3.747 + ,3.683 + ,3.715 + ,3.717 + ,3.772 + ,3.655 + ,3.615 + ,3.634 + ,3.540 + ,3.580 + ,3.598 + ,3.577 + ,3.568 + ,3.493 + ,3.446 + ,3.376 + ,3.345 + ,3.336 + ,3.341 + ,3.196 + ,3.191 + ,3.157 + ,3.111 + ,3.054 + ,3.109 + ,3.041 + ,3.017 + ,3.045 + ,3.037 + ,3.043 + ,2.992 + ,2.946 + ,2.883 + ,2.988 + ,2.950 + ,3.054 + ,3.083 + ,3.125 + ,3.175 + ,3.066 + ,3.054 + ,2.969 + ,3.007 + ,3.021 + ,2.982 + ,2.987 + ,2.964 + ,3.084 + ,3.060 + ,3.026 + ,3.083 + ,3.089 + ,3.090 + ,3.065 + ,3.072 + ,3.070 + ,3.057 + ,3.009 + ,2.990 + ,2.984 + ,3.014 + ,2.987 + ,3.025 + ,3.030 + ,3.097 + ,3.118 + ,3.153 + ,3.189 + ,3.173 + ,3.044 + ,3.148 + ,3.123 + ,3.092 + ,3.051 + ,3.054 + ,2.995 + ,2.942 + ,2.957 + ,2.772 + ,2.984 + ,3.060 + ,2.943 + ,2.914 + ,2.869 + ,2.879 + ,2.783 + ,2.473 + ,2.728 + ,2.735 + ,2.743 + ,2.638 + ,2.616 + ,2.380 + ,2.354 + ,2.095 + ,2.248 + ,2.443 + ,2.183 + ,1.978 + ,2.040 + ,2.079 + ,2.155 + ,2.052 + ,1.983 + ,1.908 + ,1.861 + ,1.895 + ,2.002 + ,2.091 + ,2.033 + ,2.132 + ,2.170 + ,2.239 + ,2.162 + ,2.130 + ,2.202 + ,2.171 + ,2.168 + ,2.032 + ,2.125 + ,2.070 + ,2.033 + ,2.042 + ,1.899 + ,1.873 + ,1.805 + ,1.819 + ,1.843 + ,1.839 + ,1.875 + ,1.879 + ,1.764 + ,1.834 + ,1.847 + ,1.835 + ,1.874 + ,1.882 + ,1.885 + ,1.904 + ,1.819 + ,1.874 + ,1.863 + ,1.887 + ,1.861 + ,1.884 + ,1.865 + ,1.853 + ,1.861 + ,1.870 + ,1.899 + ,1.916 + ,1.941 + ,2.003 + ,2.008 + ,1.991 + ,1.975 + ,1.980 + ,1.974 + ,1.973 + ,1.982 + ,1.886 + ,1.884 + ,1.884 + ,1.848 + ,1.782 + ,1.838 + ,1.814 + ,1.846 + ,1.893 + ,1.894 + ,1.929 + ,1.887 + ,1.896 + ,1.869 + ,1.897 + ,1.879 + ,1.919) > ylimmax = '' > ylimmin = '' > main = 'Robustness of Central Tendency' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > 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] 2.931244 > sqrtn <- sqrt(length(x)) > (armse <- sd(x) / sqrtn) [1] 0.04754539 > (armose <- arm / armse) [1] 61.65149 > (geo <- geomean(x)) [1] 2.825468 > (har <- harmean(x)) [1] 2.715319 > (qua <- quamean(x)) [1] 3.027229 > (win <- winmean(x)) [,1] [,2] [1,] 2.931256 0.04753354 [2,] 2.931382 0.04751176 [3,] 2.931441 0.04749788 [4,] 2.931094 0.04745549 [5,] 2.931075 0.04745388 [6,] 2.931169 0.04740023 [7,] 2.931031 0.04738443 [8,] 2.931000 0.04736573 [9,] 2.930717 0.04733717 [10,] 2.930835 0.04731972 [11,] 2.930575 0.04727726 [12,] 2.930575 0.04726927 [13,] 2.930472 0.04725263 [14,] 2.930693 0.04722342 [15,] 2.930988 0.04716712 [16,] 2.930925 0.04716222 [17,] 2.930591 0.04713625 [18,] 2.930661 0.04711806 [19,] 2.930287 0.04706428 [20,] 2.930287 0.04701196 [21,] 2.930370 0.04700458 [22,] 2.930543 0.04697480 [23,] 2.930634 0.04696675 [24,] 2.930161 0.04693063 [25,] 2.929276 0.04684709 [26,] 2.929276 0.04677978 [27,] 2.929169 0.04677177 [28,] 2.928839 0.04669271 [29,] 2.929067 0.04667249 [30,] 2.928358 0.04661945 [31,] 2.928114 0.04660124 [32,] 2.927988 0.04657131 [33,] 2.927858 0.04654049 [34,] 2.926118 0.04639011 [35,] 2.924465 0.04626911 [36,] 2.924606 0.04614233 [37,] 2.924024 0.04607669 [38,] 2.923874 0.04604185 [39,] 2.923720 0.04600613 [40,] 2.923563 0.04596953 [41,] 2.922917 0.04587146 [42,] 2.921594 0.04577690 [43,] 2.922441 0.04570225 [44,] 2.923134 0.04564128 [45,] 2.924197 0.04549162 [46,] 2.923835 0.04537958 [47,] 2.924575 0.04513914 [48,] 2.925142 0.04482163 [49,] 2.930736 0.04424881 [50,] 2.930736 0.04421806 [51,] 2.930535 0.04417254 [52,] 2.930740 0.04409118 [53,] 2.930949 0.04404085 [54,] 2.931161 0.04398960 [55,] 2.931378 0.04397113 [56,] 2.932921 0.04380547 [57,] 2.934717 0.04354864 [58,] 2.934488 0.04349715 [59,] 2.934024 0.04328466 [60,] 2.939457 0.04279075 [61,] 2.939457 0.04275374 [62,] 2.938969 0.04271943 [63,] 2.939961 0.04252180 [64,] 2.939709 0.04242663 [65,] 2.940988 0.04212327 [66,] 2.944886 0.04167956 [67,] 2.946996 0.04146436 [68,] 2.950209 0.04119904 [69,] 2.950752 0.04107132 [70,] 2.957642 0.04029627 [71,] 2.958201 0.04012311 [72,] 2.957350 0.03997716 [73,] 2.960224 0.03917672 [74,] 2.959933 0.03884795 [75,] 2.960524 0.03862166 [76,] 2.958429 0.03838539 [77,] 2.958126 0.03831861 [78,] 2.957205 0.03770004 [79,] 2.962492 0.03717822 [80,] 2.974146 0.03623988 [81,] 2.974783 0.03585391 [82,] 3.006098 0.03293806 [83,] 3.014268 0.03225211 [84,] 3.033780 0.03055051 > (tri <- trimean(x)) [,1] [,2] [1,] 2.931786 0.04752181 [2,] 2.932324 0.04750661 [3,] 2.932806 0.04749922 [4,] 2.933276 0.04749335 [5,] 2.933844 0.04749553 [6,] 2.934426 0.04749485 [7,] 2.935000 0.04750081 [8,] 2.935605 0.04750613 [9,] 2.936225 0.04751087 [10,] 2.936889 0.04751596 [11,] 2.937552 0.04751974 [12,] 2.938252 0.04752466 [13,] 2.938965 0.04752686 [14,] 2.939699 0.04752696 [15,] 2.940429 0.04752588 [16,] 2.941149 0.04752575 [17,] 2.941886 0.04752201 [18,] 2.942661 0.04751609 [19,] 2.943444 0.04750719 [20,] 2.944266 0.04749764 [21,] 2.945104 0.04748713 [22,] 2.945952 0.04747233 [23,] 2.946808 0.04745448 [24,] 2.947675 0.04743188 [25,] 2.948583 0.04740597 [26,] 2.949554 0.04737939 [27,] 2.950545 0.04735098 [28,] 2.951561 0.04731689 [29,] 2.952612 0.04728096 [30,] 2.953675 0.04723952 [31,] 2.954792 0.04719398 [32,] 2.955942 0.04714204 [33,] 2.957122 0.04708396 [34,] 2.958333 0.04701933 [35,] 2.959641 0.04695407 [36,] 2.961044 0.04688628 [37,] 2.962472 0.04681621 [38,] 2.963955 0.04674009 [39,] 2.965477 0.04665582 [40,] 2.967040 0.04656286 [41,] 2.968645 0.04646059 [42,] 2.970312 0.04635153 [43,] 2.972065 0.04623472 [44,] 2.973831 0.04610894 [45,] 2.975616 0.04597270 [46,] 2.977407 0.04583008 [47,] 2.979256 0.04567807 [48,] 2.981127 0.04552319 [49,] 2.983026 0.04536918 [50,] 2.984786 0.04523116 [51,] 2.986592 0.04507860 [52,] 2.988453 0.04491107 [53,] 2.990358 0.04472943 [54,] 2.992308 0.04453071 [55,] 2.994306 0.04431348 [56,] 2.996352 0.04407421 [57,] 2.998407 0.04382010 [58,] 3.000464 0.04355531 [59,] 3.002588 0.04326562 [60,] 3.004791 0.04295850 [61,] 3.006886 0.04265214 [62,] 3.009046 0.04231539 [63,] 3.011289 0.04194480 [64,] 3.013571 0.04154840 [65,] 3.015935 0.04111595 [66,] 3.018336 0.04065881 [67,] 3.020692 0.04018560 [68,] 3.023059 0.03967712 [69,] 3.025405 0.03913373 [70,] 3.027816 0.03853936 [71,] 3.030089 0.03794259 [72,] 3.032427 0.03729099 [73,] 3.034880 0.03657420 [74,] 3.037330 0.03584278 [75,] 3.039885 0.03505062 [76,] 3.042520 0.03418001 [77,] 3.045330 0.03321740 [78,] 3.048265 0.03213122 [79,] 3.051354 0.03095887 [80,] 3.054394 0.02968051 [81,] 3.057163 0.02833419 [82,] 3.060033 0.02682154 [83,] 3.061932 0.02548656 [84,] 3.063628 0.02402808 > (midr <- midrange(x)) [1] 2.863 > midm <- array(NA,dim=8) > for (j in 1:8) midm[j] <- midmean(x,j) > midm [1] 3.005906 3.011289 3.011289 3.011289 3.013571 3.011289 3.011289 3.011289 > postscript(file="/var/www/html/rcomp/tmp/1ehjx1240419618.ps",horizontal=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/html/rcomp/tmp/2fcp51240419618.ps",horizontal=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/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/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/html/rcomp/tmp/3l3ey1240419618.tab") > > system("convert tmp/1ehjx1240419618.ps tmp/1ehjx1240419618.png") > system("convert tmp/2fcp51240419618.ps tmp/2fcp51240419618.png") > > > proc.time() user system elapsed 1.306 0.359 1.413