R version 2.7.0 (2008-04-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(0.025749385 + ,0.02374787 + ,0.021734644 + ,0.019846432 + ,0.018140909 + ,0.016621471 + ,0.015270188 + ,0.014070406 + ,0.013016966 + ,0.012114996 + ,0.011375257 + ,0.010806475 + ,0.01041044 + ,0.010179372 + ,0.01009575 + ,0.010135252 + ,0.010270244 + ,0.010473029 + ,0.010718977 + ,0.010986893 + ,0.011260233 + ,0.011526183 + ,0.011775741 + ,0.012002907 + ,0.012204698 + ,0.012380026 + ,0.012529858 + ,0.012656553 + ,0.012763782 + ,0.012855895 + ,0.012937608 + ,0.013013837 + ,0.013089429 + ,0.013168544 + ,0.013254564 + ,0.013349985 + ,0.013456339 + ,0.013573675 + ,0.013700963 + ,0.013836316 + ,0.013976347 + ,0.014117471 + ,0.014255109 + ,0.014384556 + ,0.014500777 + ,0.014598919 + ,0.01467451 + ,0.014723505 + ,0.014742413 + ,0.014728713 + ,0.014680949 + ,0.01459833 + ,0.014481446 + ,0.014331633 + ,0.014151607 + ,0.01394507 + ,0.013716369 + ,0.013470887 + ,0.013214505 + ,0.012953678 + ,0.012695236 + ,0.012445944 + ,0.012212696 + ,0.012001897 + ,0.011819503 + ,0.011670936 + ,0.011560378 + ,0.011491247 + ,0.011465717 + ,0.01148472 + ,0.011547978 + ,0.011654156 + ,0.011800748 + ,0.011984507 + ,0.012201422 + ,0.012446842 + ,0.012715908 + ,0.013003495 + ,0.01330457 + ,0.013613992 + ,0.013926996 + ,0.014239117 + ,0.014546304 + ,0.014844989 + ,0.015132148 + ,0.015405269 + ,0.015662525 + ,0.015902616 + ,0.01612482 + ,0.016328967 + ,0.016515394 + ,0.016684831 + ,0.016838382 + ,0.016977485 + ,0.017103717 + ,0.017218825 + ,0.017324554 + ,0.017422675 + ,0.017514828 + ,0.017602504 + ,0.017687027 + ,0.017769428 + ,0.017850499 + ,0.01793075 + ,0.018010355 + ,0.018089198 + ,0.018166812 + ,0.018242448 + ,0.018314982 + ,0.018383012 + ,0.018444798 + ,0.018498321 + ,0.018541308 + ,0.018571245 + ,0.018585367 + ,0.018580929 + ,0.018555004 + ,0.018504726 + ,0.018427342 + ,0.018320256 + ,0.018181244 + ,0.018008336 + ,0.017800172 + ,0.017555895 + ,0.017275153 + ,0.01695839 + ,0.016606633 + ,0.016221732 + ,0.015806088 + ,0.01536282 + ,0.014895774 + ,0.014409298 + ,0.013908231 + ,0.013398015 + ,0.012884423 + ,0.012373828 + ,0.011872755 + ,0.011388287 + ,0.01092791 + ,0.010499544 + ,0.010111815 + ,0.009773851 + ,0.009495565 + ,0.009287581 + ,0.009160983 + ,0.009127086 + ,0.009196839 + ,0.009380001 + ,0.009684446 + ,0.010115235 + ,0.010674166 + ,0.01135928 + ,0.012164855 + ,0.013081277 + ,0.014095168 + ,0.015189115 + ,0.016342299 + ,0.017530232 + ,0.018725319 + ,0.019897617 + ,0.021015315 + ,0.022046049 + ,0.022957976 + ,0.023721479 + ,0.024310322 + ,0.024703403 + ,0.024886415 + ,0.024852666 + ,0.024604672 + ,0.024154591 + ,0.023524617 + ,0.022746938 + ,0.02186292 + ,0.020922174 + ,0.019980527 + ,0.019097146 + ,0.018329924 + ,0.017729391 + ,0.017330373 + ,0.017144205 + ,0.017153293 + ,0.01731132 + ,0.017549121 + ,0.017785239 + ,0.017938518 + ,0.017941391 + ,0.017753954 + ,0.017378803 + ,0.016876717 + ,0.016383702 + ,0.016120959 + ,0.016377277 + ,0.017432385 + ,0.019434723 + ,0.022319625 + ,0.025817482 + ,0.029505023) > 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] 0.01550756 > sqrtn <- sqrt(length(x)) > (armse <- sd(x) / sqrtn) [1] 0.0002808938 > (armose <- arm / armse) [1] 55.20791 > (geo <- geomean(x)) [1] 0.01503586 > (har <- harmean(x)) [1] 0.01458896 > (qua <- quamean(x)) [1] 0.01599841 > (win <- winmean(x)) [,1] [,2] [1,] 0.01548901 0.0002767074 [2,] 0.01548869 0.0002765338 [3,] 0.01547693 0.0002739782 [4,] 0.01547812 0.0002736435 [5,] 0.01547726 0.0002726545 [6,] 0.01548001 0.0002715018 [7,] 0.01547273 0.0002693853 [8,] 0.01547947 0.0002669536 [9,] 0.01546163 0.0002638485 [10,] 0.01546046 0.0002636162 [11,] 0.01545059 0.0002617567 [12,] 0.01541876 0.0002561685 [13,] 0.01541083 0.0002534715 [14,] 0.01539042 0.0002480386 [15,] 0.01537436 0.0002446168 [16,] 0.01536164 0.0002423437 [17,] 0.01536564 0.0002393096 [18,] 0.01530401 0.0002302702 [19,] 0.01530346 0.0002282851 [20,] 0.01522019 0.0002156359 [21,] 0.01521764 0.0002140088 [22,] 0.01524245 0.0002103566 [23,] 0.01520595 0.0002040143 [24,] 0.01516677 0.0001995845 [25,] 0.01512124 0.0001948245 [26,] 0.01511298 0.0001921066 [27,] 0.01511498 0.0001917986 [28,] 0.01511453 0.0001915822 [29,] 0.01511728 0.0001908674 [30,] 0.01511852 0.0001903577 [31,] 0.01511471 0.0001896445 [32,] 0.01512890 0.0001881034 [33,] 0.01512275 0.0001870226 [34,] 0.01513782 0.0001850594 [35,] 0.01513439 0.0001839361 [36,] 0.01512812 0.0001827486 [37,] 0.01513630 0.0001816660 [38,] 0.01515684 0.0001796157 [39,] 0.01514593 0.0001780242 [40,] 0.01513370 0.0001769081 [41,] 0.01515403 0.0001745545 [42,] 0.01515913 0.0001731276 [43,] 0.01515583 0.0001714359 [44,] 0.01513895 0.0001698434 [45,] 0.01514032 0.0001696426 [46,] 0.01516231 0.0001650181 [47,] 0.01516310 0.0001648317 [48,] 0.01517727 0.0001632919 [49,] 0.01515753 0.0001615653 [50,] 0.01516583 0.0001586865 [51,] 0.01519476 0.0001556036 [52,] 0.01520080 0.0001544017 [53,] 0.01520220 0.0001535887 [54,] 0.01520859 0.0001519377 [55,] 0.01522248 0.0001488386 [56,] 0.01520656 0.0001461613 [57,] 0.01520846 0.0001437920 [58,] 0.01521120 0.0001432451 [59,] 0.01522046 0.0001415782 [60,] 0.01521892 0.0001409365 [61,] 0.01519436 0.0001387476 [62,] 0.01521155 0.0001368821 [63,] 0.01520012 0.0001355214 [64,] 0.01521009 0.0001321986 [65,] 0.01522334 0.0001308500 > (tri <- trimean(x)) [,1] [,2] [1,] 0.01550756 0.0002725236 [2,] 0.01546850 0.0002680564 [3,] 0.01542635 0.0002633771 [4,] 0.01542635 0.0002593577 [5,] 0.01539051 0.0002551586 [6,] 0.01537204 0.0002509073 [7,] 0.01535267 0.0002465953 [8,] 0.01535267 0.0002423720 [9,] 0.01531399 0.0002382492 [10,] 0.01529573 0.0002343127 [11,] 0.01527718 0.0002301159 [12,] 0.01525923 0.0002258584 [13,] 0.01524392 0.0002220087 [14,] 0.01522895 0.0002181959 [15,] 0.01521534 0.0002147273 [16,] 0.01521534 0.0002113752 [17,] 0.01519068 0.0002080087 [18,] 0.01517809 0.0002046885 [19,] 0.01516942 0.0002020227 [20,] 0.01516057 0.0001993340 [21,] 0.01515678 0.0001975994 [22,] 0.01515305 0.0001958627 [23,] 0.01514775 0.0001942853 [24,] 0.01514440 0.0001930960 [25,] 0.01514315 0.0001921488 [26,] 0.01514434 0.0001914692 [27,] 0.01514600 0.0001909115 [28,] 0.01514761 0.0001903036 [29,] 0.01514928 0.0001896337 [30,] 0.01515087 0.0001889322 [31,] 0.01515244 0.0001881803 [32,] 0.01515244 0.0001873869 [33,] 0.01515544 0.0001866076 [34,] 0.01515695 0.0001858069 [35,] 0.01515782 0.0001850451 [36,] 0.01515888 0.0001842620 [37,] 0.01516025 0.0001834583 [38,] 0.01516130 0.0001826237 [39,] 0.01516149 0.0001818234 [40,] 0.01516217 0.0001810239 [41,] 0.01516339 0.0001801884 [42,] 0.01516378 0.0001794033 [43,] 0.01516398 0.0001786022 [44,] 0.01516432 0.0001778008 [45,] 0.01516538 0.0001769893 [46,] 0.01516643 0.0001760603 [47,] 0.01516660 0.0001753381 [48,] 0.01516675 0.0001744985 [49,] 0.01516631 0.0001736329 [50,] 0.01516667 0.0001727516 [51,] 0.01516671 0.0001719433 [52,] 0.01516554 0.0001712279 [53,] 0.01516554 0.0001704589 [54,] 0.01516248 0.0001695956 [55,] 0.01516055 0.0001686964 [56,] 0.01515794 0.0001678798 [57,] 0.01515794 0.0001671202 [58,] 0.01515363 0.0001663909 [59,] 0.01515116 0.0001655232 [60,] 0.01514815 0.0001646020 [61,] 0.01514506 0.0001635175 [62,] 0.01514287 0.0001624084 [63,] 0.01513980 0.0001612242 [64,] 0.01513980 0.0001599025 [65,] 0.01513371 0.0001586368 > (midr <- midrange(x)) [1] 0.01931605 > midm <- array(NA,dim=8) > for (j in 1:8) midm[j] <- midmean(x,j) > midm [1] 0.01513892 0.01516631 0.01516631 0.01516631 0.01516631 0.01513911 0.01516631 [8] 0.01516631 > postscript(file="/var/www/html/rcomp/tmp/1ank01224942933.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/23b111224942933.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/3m2a91224942933.tab") > > system("convert tmp/1ank01224942933.ps tmp/1ank01224942933.png") > system("convert tmp/23b111224942933.ps tmp/23b111224942933.png") > > > proc.time() user system elapsed 3.751 0.947 3.898