R version 2.5.1 (2007-06-27)
Copyright (C) 2007 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(221
+ ,213
+ ,203
+ ,195
+ ,185
+ ,183
+ ,230
+ ,250
+ ,251
+ ,244
+ ,232
+ ,225
+ ,217
+ ,206
+ ,193
+ ,180
+ ,173
+ ,169
+ ,215
+ ,239
+ ,237
+ ,230
+ ,224
+ ,213
+ ,214
+ ,197
+ ,186
+ ,177
+ ,169
+ ,159
+ ,203
+ ,217
+ ,216
+ ,206
+ ,197
+ ,195
+ ,189
+ ,175
+ ,163
+ ,153
+ ,146
+ ,144
+ ,184
+ ,194
+ ,195
+ ,196
+ ,184
+ ,185
+ ,179
+ ,163
+ ,153
+ ,143
+ ,135
+ ,132
+ ,171
+ ,177
+ ,177
+ ,171
+ ,162
+ ,158
+ ,151
+ ,139
+ ,127
+ ,117
+ ,109
+ ,109
+ ,149
+ ,157
+ ,153
+ ,145
+ ,138
+ ,134
+ ,129
+ ,116
+ ,109
+ ,97
+ ,90
+ ,89
+ ,125
+ ,136
+ ,132
+ ,125
+ ,118
+ ,115
+ ,111
+ ,101
+ ,92
+ ,85
+ ,79
+ ,81
+ ,120
+ ,128
+ ,129
+ ,124
+ ,118
+ ,117
+ ,114
+ ,109
+ ,101
+ ,95
+ ,89
+ ,90
+ ,128
+ ,139
+ ,139
+ ,135
+ ,129
+ ,128
+ ,125
+ ,117
+ ,109
+ ,104
+ ,98
+ ,100
+ ,135
+ ,146
+ ,145
+ ,140
+ ,136
+ ,139
+ ,137
+ ,131
+ ,124
+ ,119
+ ,114
+ ,116
+ ,155
+ ,167
+ ,169
+ ,163
+ ,159
+ ,158
+ ,155
+ ,147
+ ,137
+ ,129
+ ,121
+ ,119
+ ,153
+ ,166
+ ,165
+ ,158
+ ,150
+ ,147
+ ,143
+ ,138
+ ,130
+ ,124
+ ,116
+ ,114
+ ,151
+ ,164
+ ,164
+ ,160
+ ,152
+ ,151
+ ,150
+ ,143
+ ,134
+ ,128
+ ,120
+ ,119
+ ,155
+ ,169
+ ,167
+ ,158
+ ,150
+ ,146
+ ,142
+ ,136
+ ,130
+ ,124
+ ,116
+ ,116
+ ,148
+ ,160
+ ,163
+ ,151
+ ,141
+ ,139
+ ,134
+ ,130
+ ,123
+ ,117
+ ,110
+ ,105
+ ,137
+ ,160
+ ,155
+ ,142
+ ,132
+ ,131
+ ,128
+ ,123
+ ,118
+ ,112
+ ,105
+ ,102
+ ,131
+ ,149
+ ,145
+ ,132
+ ,122
+ ,119
+ ,116
+ ,111
+ ,104
+ ,100
+ ,93
+ ,91
+ ,119
+ ,139
+ ,134
+ ,124
+ ,113
+ ,109
+ ,109
+ ,106
+ ,101
+ ,98
+ ,93
+ ,91
+ ,122
+ ,139
+ ,140
+ ,132
+ ,117
+ ,114
+ ,113
+ ,110
+ ,107
+ ,103
+ ,98
+ ,98
+ ,137
+ ,148
+ ,147
+ ,139
+ ,130
+ ,128
+ ,127
+ ,123
+ ,118
+ ,114
+ ,108
+ ,111
+ ,151
+ ,159
+ ,158
+ ,148
+ ,138
+ ,137
+ ,136
+ ,133
+ ,126
+ ,120
+ ,114
+ ,116
+ ,153
+ ,162
+ ,161
+ ,149
+ ,139
+ ,135
+ ,130
+ ,127
+ ,122
+ ,117
+ ,112
+ ,113
+ ,149
+ ,157
+ ,157
+ ,147
+ ,137
+ ,132
+ ,125
+ ,123
+ ,117
+ ,114
+ ,111
+ ,112
+ ,144
+ ,150
+ ,149
+ ,134
+ ,123
+ ,116
+ ,117
+ ,111
+ ,105
+ ,102
+ ,95
+ ,93
+ ,124
+ ,130
+ ,124)
> 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: Wessa, P., (2007), Central Tendency (v1.0.2) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_centraltendency.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] 140.8687
> sqrtn <- sqrt(length(x))
> (armse <- sd(x) / sqrtn)
[1] 1.967071
> (armose <- arm / armse)
[1] 71.61343
> (geo <- geomean(x))
[1] 137.1557
> (har <- harmean(x))
[1] 133.7376
> (qua <- quamean(x))
[1] 144.8769
> (win <- winmean(x))
           [,1]      [,2]
  [1,] 140.8721 1.9657325
  [2,] 140.8586 1.9556579
  [3,] 140.8485 1.9430775
  [4,] 140.8215 1.9385217
  [5,] 140.7542 1.9232092
  [6,] 140.7138 1.9167918
  [7,] 140.7374 1.9147043
  [8,] 140.6027 1.8939580
  [9,] 140.6027 1.8867552
 [10,] 140.5354 1.8689547
 [11,] 140.3872 1.8478026
 [12,] 140.3872 1.8478026
 [13,] 140.4310 1.8342354
 [14,] 140.3838 1.8277041
 [15,] 140.4343 1.8124319
 [16,] 140.4343 1.8007458
 [17,] 140.4343 1.8007458
 [18,] 140.0101 1.7447351
 [19,] 140.0101 1.7447351
 [20,] 139.9428 1.7084312
 [21,] 139.9428 1.7084312
 [22,] 139.5724 1.6486906
 [23,] 139.5724 1.6486906
 [24,] 139.4916 1.6392303
 [25,] 139.4916 1.6228412
 [26,] 139.4916 1.6228412
 [27,] 139.5825 1.6158167
 [28,] 139.5825 1.5977786
 [29,] 139.4848 1.5865976
 [30,] 139.1818 1.5337758
 [31,] 138.8687 1.4999692
 [32,] 138.7609 1.4885969
 [33,] 138.8721 1.4801717
 [34,] 138.8721 1.4596490
 [35,] 138.9899 1.4510508
 [36,] 138.9899 1.4297284
 [37,] 138.6162 1.3915109
 [38,] 138.4882 1.3787325
 [39,] 138.2256 1.3529806
 [40,] 138.2256 1.3529806
 [41,] 138.2256 1.3529806
 [42,] 137.9428 1.3259334
 [43,] 137.7980 1.2882773
 [44,] 137.5017 1.2613046
 [45,] 137.6532 1.2502674
 [46,] 137.3434 1.2227602
 [47,] 137.3434 1.2227602
 [48,] 137.3434 1.2227602
 [49,] 137.3434 1.2227602
 [50,] 137.1751 1.1814847
 [51,] 137.1751 1.1814847
 [52,] 137.0000 1.1666667
 [53,] 137.0000 1.1389496
 [54,] 136.8182 1.1239710
 [55,] 136.8182 1.1239710
 [56,] 136.8182 1.0952674
 [57,] 136.8182 1.0952674
 [58,] 136.8182 1.0952674
 [59,] 136.8182 1.0952674
 [60,] 136.6162 1.0790816
 [61,] 136.6162 1.0790816
 [62,] 136.4074 1.0626247
 [63,] 136.4074 1.0311361
 [64,] 136.6229 1.0161907
 [65,] 136.6229 1.0161907
 [66,] 136.4007 0.9990629
 [67,] 136.4007 0.9990629
 [68,] 136.4007 0.9990629
 [69,] 136.1684 0.9814550
 [70,] 136.1684 0.9814550
 [71,] 136.1684 0.9814550
 [72,] 136.4108 0.9647998
 [73,] 136.4108 0.9647998
 [74,] 136.1616 0.9461105
 [75,] 136.1616 0.9461105
 [76,] 136.1616 0.9461105
 [77,] 135.6431 0.9081371
 [78,] 135.6431 0.9081371
 [79,] 135.6431 0.9081371
 [80,] 135.9125 0.8896333
 [81,] 135.3670 0.8507529
 [82,] 135.3670 0.8507529
 [83,] 135.3670 0.8507529
 [84,] 135.6498 0.8314310
 [85,] 135.6498 0.8314310
 [86,] 135.3603 0.8111885
 [87,] 135.0673 0.7910704
 [88,] 135.0673 0.7910704
 [89,] 135.3670 0.7706938
 [90,] 135.3670 0.7706938
 [91,] 135.3670 0.7706938
 [92,] 135.3670 0.7288202
 [93,] 135.6801 0.7081726
 [94,] 135.6801 0.7081726
 [95,] 135.6801 0.7081726
 [96,] 135.6801 0.6652371
 [97,] 135.6801 0.6652371
 [98,] 135.6801 0.6652371
 [99,] 135.6801 0.6652371
> (tri <- trimean(x))
           [,1]      [,2]
  [1,] 140.7051 1.9333986
  [2,] 140.5358 1.8994577
  [3,] 140.3711 1.8693045
  [4,] 140.2076 1.8423469
  [5,] 140.0488 1.8153850
  [6,] 139.9018 1.7907087
  [7,] 139.7597 1.7661285
  [8,] 139.6121 1.7407053
  [9,] 139.4803 1.7173227
 [10,] 139.3466 1.6938198
 [11,] 139.2182 1.6715064
 [12,] 139.1026 1.6506984
 [13,] 138.9852 1.6288701
 [14,] 138.8625 1.6073970
 [15,] 138.7416 1.5855237
 [16,] 138.6151 1.5640085
 [17,] 138.4867 1.5424775
 [18,] 138.3563 1.5198261
 [19,] 138.2510 1.5009534
 [20,] 138.1440 1.4811132
 [21,] 138.0392 1.4631684
 [22,] 137.9328 1.4442928
 [23,] 137.8446 1.4290048
 [24,] 137.7550 1.4129308
 [25,] 137.6680 1.3967120
 [26,] 137.5796 1.3807796
 [27,] 137.4897 1.3640025
 [28,] 137.3942 1.3467698
 [29,] 137.2971 1.3297973
 [30,] 137.2025 1.3126544
 [31,] 137.1191 1.2980371
 [32,] 137.0472 1.2847974
 [33,] 136.9784 1.2715516
 [34,] 136.9039 1.2580537
 [35,] 136.8282 1.2450228
 [36,] 136.7467 1.2317280
 [37,] 136.6637 1.2189154
 [38,] 136.5928 1.2076589
 [39,] 136.5251 1.1965051
 [40,] 136.4654 1.1861981
 [41,] 136.4047 1.1752889
 [42,] 136.3427 1.1637364
 [43,] 136.2891 1.1530514
 [44,] 136.2392 1.1437855
 [45,] 136.1981 1.1354378
 [46,] 136.1512 1.1271222
 [47,] 136.1133 1.1197656
 [48,] 136.0746 1.1119420
 [49,] 136.0352 1.1036202
 [50,] 135.9949 1.0947663
 [51,] 135.9590 1.0874974
 [52,] 135.9223 1.0797465
 [53,] 135.8901 1.0722673
 [54,] 135.8571 1.0656850
 [55,] 135.8289 1.0594231
 [56,] 135.8000 1.0527174
 [57,] 135.8000 1.0469717
 [58,] 135.7403 1.0408005
 [59,] 135.7095 1.0341737
 [60,] 135.6780 1.0270582
 [61,] 135.6514 1.0202826
 [62,] 135.6243 1.0129982
 [63,] 135.6023 1.0060525
 [64,] 135.5799 1.0001967
 [65,] 135.5509 0.9945983
 [66,] 135.5212 0.9885370
 [67,] 135.4969 0.9829022
 [68,] 135.4720 0.9767905
 [69,] 135.4465 0.9701643
 [70,] 135.4268 0.9639596
 [71,] 135.4065 0.9572206
 [72,] 135.3856 0.9499029
 [73,] 135.3576 0.9428408
 [74,] 135.3289 0.9351570
 [75,] 135.3061 0.9279120
 [76,] 135.2828 0.9200150
 [77,] 135.2587 0.9114067
 [78,] 135.2482 0.9043493
 [79,] 135.2374 0.8966267
 [80,] 135.2263 0.8881763
 [81,] 135.2074 0.8800407
 [82,] 135.2030 0.8735734
 [83,] 135.1985 0.8664436
 [84,] 135.1938 0.8585862
 [85,] 135.1938 0.8511531
 [86,] 135.1680 0.8429393
 [87,] 135.1626 0.8352470
 [88,] 135.1653 0.8280838
 [89,] 135.1681 0.8201215
 [90,] 135.1624 0.8126938
 [91,] 135.1565 0.8044111
 [92,] 135.1504 0.7951756
 [93,] 135.1441 0.7879983
 [94,] 135.1284 0.7814209
 [95,] 135.1121 0.7739963
 [96,] 135.0952 0.7656217
 [97,] 135.0777 0.7595686
 [98,] 135.0594 0.7526615
 [99,] 135.0404 0.7447915
> (midr <- midrange(x))
[1] 165
> midm <- array(NA,dim=8)
> for (j in 1:8) midm[j] <- midmean(x,j)
> midm
[1] 135.0861 135.0861 135.0861 135.0861 135.0861 135.0861 135.0861 135.0861
> postscript(file="/var/www/html/rcomp/tmp/1ca1z1192983530.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/2x0n61192983530.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 
> 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/3xp0o1192983530.tab") 
> 
> system("convert tmp/1ca1z1192983530.ps tmp/1ca1z1192983530.png")
> system("convert tmp/2x0n61192983530.ps tmp/2x0n61192983530.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.157   0.363   1.383