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(12.182
+ ,12.242
+ ,12.368
+ ,12.552
+ ,12.377
+ ,12.350
+ ,12.333
+ ,12.427
+ ,12.281
+ ,12.381
+ ,12.569
+ ,12.684
+ ,12.689
+ ,12.580
+ ,12.264
+ ,12.259
+ ,12.205
+ ,12.255
+ ,12.039
+ ,11.893
+ ,11.741
+ ,12.149
+ ,12.096
+ ,12.146
+ ,11.946
+ ,11.976
+ ,12.392
+ ,12.102
+ ,12.362
+ ,12.547
+ ,12.532
+ ,12.422
+ ,12.304
+ ,12.216
+ ,12.267
+ ,12.652
+ ,12.605
+ ,12.626
+ ,12.613
+ ,12.603
+ ,12.575
+ ,12.527
+ ,12.580
+ ,12.325
+ ,12.304
+ ,12.372
+ ,12.617
+ ,12.627
+ ,12.851
+ ,12.825
+ ,12.721
+ ,12.765
+ ,12.848
+ ,12.891
+ ,12.870
+ ,12.831
+ ,12.818
+ ,13.013
+ ,13.057
+ ,12.969
+ ,13.011
+ ,12.815
+ ,12.861
+ ,12.746
+ ,12.872
+ ,12.825
+ ,12.891
+ ,12.993
+ ,12.985
+ ,13.026
+ ,12.825
+ ,12.598
+ ,12.621
+ ,12.480
+ ,12.543
+ ,12.594
+ ,12.647
+ ,12.638
+ ,12.503
+ ,12.392
+ ,12.389
+ ,12.603
+ ,12.210
+ ,12.278
+ ,12.286
+ ,12.090
+ ,12.145
+ ,12.307
+ ,12.270
+ ,12.159
+ ,12.023
+ ,12.062
+ ,11.844
+ ,11.842
+ ,11.805
+ ,11.809
+ ,11.453
+ ,11.346
+ ,11.345
+ ,11.382
+ ,11.216
+ ,11.289
+ ,11.225
+ ,11.382
+ ,11.148
+ ,11.226
+ ,11.104
+ ,11.051
+ ,10.962
+ ,11.238
+ ,11.437
+ ,11.495
+ ,11.458
+ ,11.603
+ ,11.630
+ ,11.341
+ ,11.369
+ ,11.133
+ ,11.398
+ ,11.578
+ ,11.380
+ ,11.326
+ ,11.286
+ ,11.604
+ ,11.655
+ ,11.432
+ ,11.730
+ ,11.782
+ ,11.633
+ ,11.532
+ ,11.611
+ ,11.660
+ ,11.478
+ ,11.346
+ ,11.415
+ ,11.427
+ ,11.626
+ ,11.384
+ ,11.412
+ ,11.500
+ ,11.713
+ ,11.546
+ ,11.506
+ ,11.532
+ ,11.186
+ ,11.225
+ ,11.515
+ ,11.234
+ ,11.264
+ ,11.429
+ ,11.416
+ ,10.906
+ ,11.057
+ ,10.609
+ ,11.028
+ ,11.394
+ ,11.016
+ ,10.850
+ ,10.827
+ ,11.019
+ ,11.140
+ ,10.372
+ ,10.847
+ ,10.826
+ ,10.484
+ ,10.323
+ ,9.955
+ ,9.437
+ ,9.262
+ ,8.569
+ ,8.462
+ ,9.389
+ ,9.302
+ ,8.577
+ ,8.975
+ ,8.852
+ ,9.263
+ ,9.028
+ ,8.520
+ ,8.683
+ ,8.376
+ ,8.179
+ ,9.062
+ ,9.005
+ ,9.179
+ ,9.326
+ ,9.324
+ ,9.617
+ ,9.134
+ ,8.696
+ ,8.947
+ ,8.864
+ ,8.685
+ ,8.822
+ ,8.495
+ ,8.273
+ ,8.421
+ ,7.996
+ ,7.552
+ ,8.048
+ ,8.445
+ ,8.464
+ ,8.725
+ ,8.827
+ ,8.154
+ ,8.409
+ ,8.587
+ ,8.376
+ ,8.638
+ ,8.934
+ ,8.693
+ ,8.750
+ ,8.563
+ ,8.629
+ ,8.566
+ ,8.922
+ ,8.824
+ ,8.607
+ ,8.573
+ ,8.519
+ ,8.417
+ ,8.469
+ ,8.516
+ ,8.488
+ ,8.666
+ ,8.772
+ ,9.027
+ ,8.955
+ ,8.997
+ ,8.770
+ ,8.739
+ ,8.599
+ ,8.475
+ ,8.446
+ ,8.196
+ ,8.216
+ ,8.280
+ ,7.949
+ ,8.224
+ ,8.109
+ ,8.078
+ ,8.117
+ ,8.176
+ ,8.373
+ ,8.149
+ ,8.001
+ ,7.937
+ ,8.070
+ ,7.955
+ ,8.056)
> 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] 10.84238
> sqrtn <- sqrt(length(x))
> (armse <- sd(x) / sqrtn)
[1] 0.1072105
> (armose <- arm / armse)
[1] 101.1317
> (geo <- geomean(x))
[1] 10.70196
> (har <- harmean(x))
[1] 10.55468
> (qua <- quamean(x))
[1] 10.97357
> (win <- winmean(x))
          [,1]       [,2]
 [1,] 10.84380 0.10702136
 [2,] 10.84379 0.10700242
 [3,] 10.84384 0.10699265
 [4,] 10.84420 0.10689864
 [5,] 10.84414 0.10687506
 [6,] 10.84489 0.10672465
 [7,] 10.84293 0.10652949
 [8,] 10.84338 0.10648253
 [9,] 10.84298 0.10639983
[10,] 10.84414 0.10626491
[11,] 10.84410 0.10619831
[12,] 10.84515 0.10600411
[13,] 10.84526 0.10596572
[14,] 10.84554 0.10576844
[15,] 10.84536 0.10572309
[16,] 10.84644 0.10561317
[17,] 10.84780 0.10547649
[18,] 10.84788 0.10538094
[19,] 10.85137 0.10499398
[20,] 10.84793 0.10464140
[21,] 10.85415 0.10376416
[22,] 10.85221 0.10357856
[23,] 10.84927 0.10336669
[24,] 10.85196 0.10302934
[25,] 10.84956 0.10272620
[26,] 10.84945 0.10265005
[27,] 10.85107 0.10233628
[28,] 10.84995 0.10223951
[29,] 10.85169 0.10205652
[30,] 10.85133 0.10199210
[31,] 10.85146 0.10189931
[32,] 10.85171 0.10179164
[33,] 10.85237 0.10155768
[34,] 10.85305 0.10144990
[35,] 10.85599 0.10117614
[36,] 10.85570 0.10108614
[37,] 10.85526 0.10103145
[38,] 10.85967 0.10028065
[39,] 10.86014 0.10023762
[40,] 10.85982 0.10013843
[41,] 10.85949 0.10001056
[42,] 10.85730 0.09975341
[43,] 10.85816 0.09953701
[44,] 10.85957 0.09929574
[45,] 10.85903 0.09902958
[46,] 10.86216 0.09859822
[47,] 10.85934 0.09813964
[48,] 10.86030 0.09735593
[49,] 10.85324 0.09636673
[50,] 10.85264 0.09626502
[51,] 10.84816 0.09571888
[52,] 10.84878 0.09566247
[53,] 10.85429 0.09506692
[54,] 10.85559 0.09468326
[55,] 10.85713 0.09440925
[56,] 10.86049 0.09393587
[57,] 10.86003 0.09383622
[58,] 10.87024 0.09271339
[59,] 10.86788 0.09248888
[60,] 10.86452 0.09216314
[61,] 10.86867 0.09149752
[62,] 10.86718 0.09094960
[63,] 10.88104 0.08961340
[64,] 10.88411 0.08934390
[65,] 10.88281 0.08874978
[66,] 10.88360 0.08848107
[67,] 10.88816 0.08796136
[68,] 10.89197 0.08730118
[69,] 10.89335 0.08705629
[70,] 10.89867 0.08646698
[71,] 10.89753 0.08635230
[72,] 10.90617 0.08542914
[73,] 10.92340 0.08337073
[74,] 10.92902 0.08173621
[75,] 10.95212 0.07948834
[76,] 10.95091 0.07936585
[77,] 10.95584 0.07789202
[78,] 10.95552 0.07685428
[79,] 10.95300 0.07660195
[80,] 10.97220 0.07482522
[81,] 10.98742 0.07348510
[82,] 11.03236 0.06760690
[83,] 11.14258 0.05815022
> (tri <- trimean(x))
          [,1]       [,2]
 [1,] 10.84672 0.10688179
 [2,] 10.84969 0.10673001
 [3,] 10.85271 0.10657539
 [4,] 10.85576 0.10641120
 [5,] 10.85877 0.10625912
 [6,] 10.86185 0.10609897
 [7,] 10.86484 0.10595349
 [8,] 10.86819 0.10582641
 [9,] 10.87153 0.10569330
[10,] 10.87498 0.10555770
[11,] 10.87836 0.10542499
[12,] 10.88181 0.10528602
[13,] 10.88521 0.10515340
[14,] 10.88868 0.10501047
[15,] 10.89218 0.10487124
[16,] 10.89576 0.10472119
[17,] 10.89932 0.10456517
[18,] 10.90286 0.10440456
[19,] 10.90647 0.10423510
[20,] 10.90992 0.10407849
[21,] 10.91364 0.10393037
[22,] 10.91708 0.10382862
[23,] 10.92070 0.10372443
[24,] 10.92454 0.10361862
[25,] 10.92832 0.10351969
[26,] 10.93230 0.10342403
[27,] 10.93636 0.10331740
[28,] 10.94043 0.10321382
[29,] 10.94464 0.10309920
[30,] 10.94886 0.10297836
[31,] 10.95318 0.10284317
[32,] 10.95759 0.10269431
[33,] 10.96209 0.10253161
[34,] 10.96665 0.10236126
[35,] 10.97129 0.10217529
[36,] 10.97592 0.10198239
[37,] 10.98066 0.10177064
[38,] 10.98553 0.10153657
[39,] 10.99035 0.10131847
[40,] 10.99526 0.10107641
[41,] 11.00030 0.10081146
[42,] 11.00547 0.10052340
[43,] 11.01085 0.10021687
[44,] 11.01633 0.09988895
[45,] 11.02189 0.09953941
[46,] 11.02762 0.09916701
[47,] 11.03338 0.09878038
[48,] 11.03940 0.09837718
[49,] 11.04553 0.09797546
[50,] 11.05207 0.09758245
[51,] 11.05881 0.09715046
[52,] 11.06588 0.09669942
[53,] 11.07313 0.09620153
[54,] 11.08040 0.09568681
[55,] 11.08784 0.09513859
[56,] 11.09543 0.09454717
[57,] 11.10315 0.09392162
[58,] 11.11110 0.09323304
[59,] 11.11897 0.09254637
[60,] 11.12715 0.09179584
[61,] 11.13570 0.09098038
[62,] 11.14439 0.09012114
[63,] 11.15340 0.08920076
[64,] 11.16226 0.08827945
[65,] 11.17132 0.08727335
[66,] 11.18072 0.08619176
[67,] 11.19042 0.08500361
[68,] 11.20032 0.08371869
[69,] 11.21044 0.08233506
[70,] 11.22088 0.08080342
[71,] 11.23154 0.07913710
[72,] 11.24263 0.07726644
[73,] 11.25387 0.07524795
[74,] 11.26496 0.07318056
[75,] 11.27631 0.07099771
[76,] 11.28734 0.06877360
[77,] 11.29886 0.06623386
[78,] 11.31071 0.06348038
[79,] 11.32309 0.06040500
[80,] 11.33610 0.05682147
[81,] 11.34902 0.05286945
[82,] 11.36200 0.04836024
[83,] 11.37396 0.04398949
> (midr <- midrange(x))
[1] 10.3045
> midm <- array(NA,dim=8)
> for (j in 1:8) midm[j] <- midmean(x,j)
> midm
[1] 11.13509 11.14439 11.14439 11.14439 11.15340 11.14439 11.14439 11.14439
> postscript(file="/var/www/html/rcomp/tmp/1uced1240419008.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/2u8op1240419008.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/3g7x91240419008.tab") 
> 
> system("convert tmp/1uced1240419008.ps tmp/1uced1240419008.png")
> system("convert tmp/2u8op1240419008.ps tmp/2u8op1240419008.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.292   0.336   1.410