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(358
+ ,361
+ ,360
+ ,361
+ ,360
+ ,362
+ ,375
+ ,377
+ ,380
+ ,382
+ ,382
+ ,382
+ ,384
+ ,383
+ ,381
+ ,378
+ ,379
+ ,377
+ ,391
+ ,397
+ ,393
+ ,393
+ ,394
+ ,392
+ ,405
+ ,399
+ ,384
+ ,369
+ ,360
+ ,347
+ ,352
+ ,351
+ ,348
+ ,346
+ ,344
+ ,347
+ ,351
+ ,346
+ ,342
+ ,338
+ ,336
+ ,334
+ ,338
+ ,338
+ ,337
+ ,344
+ ,341
+ ,348
+ ,352
+ ,345
+ ,343
+ ,339
+ ,336
+ ,334
+ ,344
+ ,340
+ ,338
+ ,340
+ ,338
+ ,341
+ ,342
+ ,337
+ ,332
+ ,326
+ ,321
+ ,316
+ ,327
+ ,324
+ ,316
+ ,315
+ ,313
+ ,316
+ ,316
+ ,313
+ ,312
+ ,303
+ ,299
+ ,295
+ ,306
+ ,309
+ ,300
+ ,299
+ ,298
+ ,301
+ ,302
+ ,298
+ ,294
+ ,289
+ ,285
+ ,284
+ ,298
+ ,300
+ ,295
+ ,297
+ ,299
+ ,306
+ ,309
+ ,310
+ ,305
+ ,303
+ ,301
+ ,301
+ ,316
+ ,322
+ ,317
+ ,322
+ ,323
+ ,331
+ ,336
+ ,334
+ ,334
+ ,335
+ ,332
+ ,336
+ ,353
+ ,360
+ ,357
+ ,360
+ ,365
+ ,377
+ ,384
+ ,388
+ ,388
+ ,390
+ ,391
+ ,395
+ ,415
+ ,425
+ ,425
+ ,423
+ ,427
+ ,434
+ ,439
+ ,438
+ ,435
+ ,434
+ ,434
+ ,435
+ ,449
+ ,456
+ ,452
+ ,448
+ ,445
+ ,452
+ ,457
+ ,454
+ ,446
+ ,444
+ ,440
+ ,441
+ ,458
+ ,466
+ ,464
+ ,464
+ ,457
+ ,465
+ ,472
+ ,461
+ ,451
+ ,445
+ ,435
+ ,426
+ ,444
+ ,451
+ ,440
+ ,432
+ ,429
+ ,434
+ ,438
+ ,435
+ ,431
+ ,427
+ ,420
+ ,424
+ ,440
+ ,447
+ ,436
+ ,427
+ ,421
+ ,427
+ ,427
+ ,424
+ ,417
+ ,410
+ ,402
+ ,400
+ ,417
+ ,424
+ ,414
+ ,398
+ ,390
+ ,395
+ ,398
+ ,393
+ ,386
+ ,378
+ ,374
+ ,373
+ ,392
+ ,402
+ ,387
+ ,379
+ ,370
+ ,373
+ ,377
+ ,370
+ ,358
+ ,357
+ ,349
+ ,348
+ ,369
+ ,381
+ ,368
+ ,361
+ ,351
+ ,351
+ ,358
+ ,354
+ ,347
+ ,345
+ ,343
+ ,340
+ ,362
+ ,370
+ ,373
+ ,371
+ ,354
+ ,357
+ ,363
+ ,364
+ ,363
+ ,358
+ ,357
+ ,357
+ ,380
+ ,378
+ ,376
+ ,380
+ ,379
+ ,384
+ ,392
+ ,394
+ ,392
+ ,396
+ ,392
+ ,396
+ ,419
+ ,421
+ ,420
+ ,418
+ ,410
+ ,418
+ ,426
+ ,428
+ ,430
+ ,424
+ ,423
+ ,427
+ ,441
+ ,449
+ ,452
+ ,462
+ ,455
+ ,461
+ ,461
+ ,463
+ ,462
+ ,456
+ ,455
+ ,456
+ ,472
+ ,472
+ ,471
+ ,465
+ ,459
+ ,465
+ ,468
+ ,467
+ ,463
+ ,460
+ ,462
+ ,461
+ ,476
+ ,476
+ ,471
+ ,453
+ ,443
+ ,442
+ ,444
+ ,438
+ ,427
+ ,424
+ ,416
+ ,406
+ ,431
+ ,434
+ ,418)
> 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] 385.3502
> sqrtn <- sqrt(length(x))
> (armse <- sd(x) / sqrtn)
[1] 2.991095
> (armose <- arm / armse)
[1] 128.8325
> (geo <- geomean(x))
[1] 381.86
> (har <- harmean(x))
[1] 378.3334
> (qua <- quamean(x))
[1] 388.7711
> (win <- winmean(x))
           [,1]     [,2]
  [1,] 385.3535 2.990711
  [2,] 385.3535 2.985016
  [3,] 385.4040 2.979645
  [4,] 385.4175 2.978256
  [5,] 385.4007 2.976612
  [6,] 385.4411 2.972508
  [7,] 385.3939 2.963382
  [8,] 385.3670 2.960859
  [9,] 385.3367 2.958052
 [10,] 385.3367 2.951623
 [11,] 385.3367 2.951623
 [12,] 385.3367 2.951623
 [13,] 385.3367 2.943345
 [14,] 385.3367 2.943345
 [15,] 385.3367 2.933881
 [16,] 385.3367 2.933881
 [17,] 385.2795 2.928789
 [18,] 385.3401 2.922924
 [19,] 385.4040 2.916790
 [20,] 385.3367 2.910846
 [21,] 385.4781 2.897454
 [22,] 385.5522 2.890535
 [23,] 385.5522 2.890535
 [24,] 385.7138 2.861176
 [25,] 385.6296 2.853829
 [26,] 385.6296 2.838346
 [27,] 385.7205 2.814312
 [28,] 385.8148 2.806008
 [29,] 385.7172 2.797680
 [30,] 385.9192 2.780104
 [31,] 386.0236 2.771151
 [32,] 385.9158 2.762003
 [33,] 385.9158 2.762003
 [34,] 385.8013 2.752375
 [35,] 385.6835 2.742556
 [36,] 385.6835 2.722174
 [37,] 386.1818 2.680461
 [38,] 386.3098 2.669999
 [39,] 386.1785 2.659135
 [40,] 386.3131 2.648205
 [41,] 386.1751 2.614484
 [42,] 386.4579 2.591978
 [43,] 386.4579 2.568853
 [44,] 386.9024 2.511409
 [45,] 386.9024 2.487859
 [46,] 386.7475 2.475487
 [47,] 387.0640 2.452093
 [48,] 386.9024 2.439252
 [49,] 386.9024 2.439252
 [50,] 386.9024 2.439252
 [51,] 386.9024 2.413190
 [52,] 386.9024 2.386819
 [53,] 386.7239 2.372965
 [54,] 386.7239 2.372965
 [55,] 386.5387 2.358720
 [56,] 386.7273 2.345142
 [57,] 386.7273 2.345142
 [58,] 386.7273 2.316269
 [59,] 386.5286 2.301190
 [60,] 386.5286 2.301190
 [61,] 386.5286 2.301190
 [62,] 386.1111 2.269914
 [63,] 386.1111 2.239081
 [64,] 386.3266 2.223839
 [65,] 386.3266 2.223839
 [66,] 386.3266 2.223839
 [67,] 386.3266 2.191388
 [68,] 386.3266 2.191388
 [69,] 386.5589 2.175231
 [70,] 386.5589 2.175231
 [71,] 386.7980 2.158766
 [72,] 386.3131 2.123233
 [73,] 386.3131 2.088534
 [74,] 386.3131 2.088534
 [75,] 386.0606 2.070354
 [76,] 386.0606 2.034551
 [77,] 385.8013 2.016142
 [78,] 385.8013 1.979727
 [79,] 385.8013 1.979727
 [80,] 386.0707 1.961516
 [81,] 386.0707 1.961516
 [82,] 386.0707 1.961516
 [83,] 386.3502 1.942795
 [84,] 386.3502 1.942795
 [85,] 386.0640 1.922638
 [86,] 386.3535 1.903354
 [87,] 386.6465 1.844202
 [88,] 386.6465 1.844202
 [89,] 386.3468 1.823221
 [90,] 386.3468 1.823221
 [91,] 386.6532 1.803240
 [92,] 386.6532 1.803240
 [93,] 386.9663 1.782999
 [94,] 386.9663 1.740514
 [95,] 386.9663 1.740514
 [96,] 387.2896 1.633965
 [97,] 387.2896 1.633965
 [98,] 386.9596 1.611034
 [99,] 386.9596 1.611034
> (tri <- trimean(x))
           [,1]     [,2]
  [1,] 385.3864 2.975801
  [2,] 385.4198 2.960260
  [3,] 385.4536 2.947071
  [4,] 385.4706 2.935231
  [5,] 385.4843 2.923244
  [6,] 385.5018 2.911085
  [7,] 385.5124 2.899154
  [8,] 385.5302 2.888176
  [9,] 385.5520 2.877054
 [10,] 385.5776 2.865781
 [11,] 385.6036 2.854754
 [12,] 385.6300 2.843200
 [13,] 385.6568 2.831096
 [14,] 385.6840 2.819227
 [15,] 385.7116 2.806786
 [16,] 385.7396 2.794580
 [17,] 385.7681 2.781779
 [18,] 385.8008 2.768762
 [19,] 385.8301 2.755565
 [20,] 385.8560 2.742178
 [21,] 385.8863 2.728556
 [22,] 385.9091 2.715200
 [23,] 385.9283 2.701644
 [24,] 385.9478 2.687400
 [25,] 385.9595 2.674340
 [26,] 385.9755 2.661064
 [27,] 385.9918 2.648051
 [28,] 386.0041 2.635803
 [29,] 386.0126 2.623387
 [30,] 386.0253 2.610779
 [31,] 386.0298 2.598490
 [32,] 386.0300 2.586025
 [33,] 386.0346 2.573372
 [34,] 386.0393 2.560008
 [35,] 386.0485 2.546417
 [36,] 386.0622 2.532586
 [37,] 386.0762 2.519057
 [38,] 386.0724 2.506950
 [39,] 386.0639 2.494664
 [40,] 386.0599 2.482204
 [41,] 386.0512 2.469551
 [42,] 386.0469 2.457845
 [43,] 386.0332 2.446534
 [44,] 386.0191 2.435655
 [45,] 385.9903 2.426867
 [46,] 385.9610 2.418602
 [47,] 385.9360 2.410340
 [48,] 385.9005 2.402569
 [49,] 385.8693 2.394829
 [50,] 385.8376 2.386493
 [51,] 385.8051 2.377524
 [52,] 385.7720 2.369124
 [53,] 385.7382 2.361320
 [54,] 385.7090 2.353551
 [55,] 385.6791 2.345145
 [56,] 385.6541 2.336749
 [57,] 385.6541 2.328299
 [58,] 385.5912 2.319148
 [59,] 385.5587 2.310634
 [60,] 385.5311 2.302134
 [61,] 385.5029 2.292902
 [62,] 385.4740 2.282884
 [63,] 385.4561 2.273586
 [64,] 385.4379 2.264986
 [65,] 385.4132 2.256339
 [66,] 385.3879 2.246893
 [67,] 385.3620 2.236590
 [68,] 385.3354 2.227004
 [69,] 385.3082 2.216524
 [70,] 385.2739 2.205893
 [71,] 385.2387 2.194274
 [72,] 385.1961 2.182428
 [73,] 385.1656 2.171427
 [74,] 385.1342 2.161212
 [75,] 385.1020 2.149972
 [76,] 385.0759 2.138636
 [77,] 385.0490 2.128131
 [78,] 385.0284 2.117566
 [79,] 385.0072 2.107909
 [80,] 384.9854 2.097172
 [81,] 384.9556 2.086278
 [82,] 384.9248 2.074172
 [83,] 384.8931 2.060741
 [84,] 384.8527 2.046963
 [85,] 384.8527 2.031674
 [86,] 384.7760 2.016060
 [87,] 384.7317 1.999936
 [88,] 384.6777 1.985839
 [89,] 384.6218 1.970079
 [90,] 384.5726 1.953984
 [91,] 384.5217 1.935983
 [92,] 384.4602 1.917220
 [93,] 384.3964 1.896213
 [94,] 384.3211 1.874135
 [95,] 384.2430 1.852707
 [96,] 384.1619 1.828605
 [97,] 384.0680 1.809995
 [98,] 383.9703 1.788880
 [99,] 383.8788 1.767061
> (midr <- midrange(x))
[1] 380
> midm <- array(NA,dim=8)
> for (j in 1:8) midm[j] <- midmean(x,j)
> midm
[1] 384.5503 385.1656 385.1656 385.1656 385.1656 385.1656 385.1656 385.1656
> postscript(file="/var/www/html/rcomp/tmp/1r9c81192986188.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/277k81192986188.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/3euqc1192986189.tab") 
> 
> system("convert tmp/1r9c81192986188.ps tmp/1r9c81192986188.png")
> system("convert tmp/277k81192986188.ps tmp/277k81192986188.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.160   0.367   1.296