R version 2.12.0 (2010-10-15)
Copyright (C) 2010 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-bit)

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(24.90
+ ,25.06
+ ,25.10
+ ,24.92
+ ,25.46
+ ,25.89
+ ,25.39
+ ,25.38
+ ,25.25
+ ,24.88
+ ,25.00
+ ,25.00
+ ,24.07
+ ,23.60
+ ,23.18
+ ,23.25
+ ,23.04
+ ,22.77
+ ,22.25
+ ,22.41
+ ,22.50
+ ,22.91
+ ,22.88
+ ,21.69
+ ,21.19
+ ,21.56
+ ,22.00
+ ,22.13
+ ,22.27
+ ,22.30
+ ,21.94
+ ,22.40
+ ,22.77
+ ,22.90
+ ,23.03
+ ,23.05
+ ,22.41
+ ,22.26
+ ,21.90
+ ,22.01
+ ,22.62
+ ,22.76
+ ,23.40
+ ,23.63
+ ,24.05
+ ,23.82
+ ,23.71
+ ,23.95
+ ,23.61
+ ,23.98
+ ,23.56
+ ,23.99
+ ,24.33
+ ,24.48
+ ,24.31
+ ,24.38
+ ,24.63
+ ,25.54
+ ,25.75
+ ,25.73
+ ,25.85
+ ,25.78
+ ,25.86
+ ,26.86
+ ,27.36
+ ,27.38
+ ,26.58
+ ,27.65
+ ,27.73
+ ,27.18
+ ,27.32
+ ,27.30
+ ,26.90
+ ,26.70
+ ,26.75
+ ,26.41
+ ,26.29
+ ,27.51
+ ,27.91
+ ,27.70
+ ,27.28
+ ,28.25
+ ,27.62
+ ,27.30
+ ,25.94
+ ,24.99
+ ,25.50
+ ,24.42
+ ,26.58
+ ,25.84
+ ,26.76
+ ,26.74
+ ,26.68
+ ,25.55
+ ,26.40
+ ,25.19
+ ,23.94
+ ,24.20
+ ,24.20
+ ,23.07
+ ,24.07
+ ,25.02
+ ,24.65
+ ,24.68
+ ,24.63
+ ,24.49
+ ,25.05
+ ,24.31
+ ,23.90
+ ,23.68
+ ,24.50
+ ,25.22
+ ,25.48
+ ,26.00
+ ,26.07
+ ,26.06
+ ,26.22
+ ,26.70
+ ,27.20
+ ,26.77
+ ,26.11
+ ,25.43
+ ,24.99
+ ,25.51
+ ,24.00
+ ,23.86
+ ,22.96
+ ,23.41
+ ,23.17
+ ,24.12
+ ,23.87
+ ,24.27
+ ,24.40
+ ,24.16
+ ,25.15
+ ,25.09
+ ,24.60
+ ,24.33
+ ,24.14
+ ,24.36
+ ,25.40
+ ,26.15
+ ,26.77
+ ,26.94
+ ,26.33
+ ,26.24
+ ,26.23
+ ,25.88
+ ,27.00
+ ,26.91
+ ,27.15
+ ,27.78
+ ,28.73
+ ,28.83
+ ,28.68
+ ,27.56
+ ,27.15
+ ,27.41
+ ,27.47
+ ,28.76
+ ,28.47
+ ,27.94
+ ,27.23
+ ,27.01
+ ,26.15
+ ,26.11
+ ,27.20
+ ,27.36
+ ,27.33
+ ,27.43
+ ,28.92
+ ,29.45
+ ,29.01
+ ,29.25
+ ,29.14
+ ,29.64
+ ,30.40
+ ,30.62
+ ,31.25
+ ,31.75
+ ,31.30
+ ,30.70
+ ,31.03
+ ,31.46
+ ,31.28
+ ,31.03
+ ,30.95
+ ,31.17
+ ,31.29
+ ,31.91
+ ,32.10
+ ,31.71
+ ,31.90
+ ,32.02
+ ,32.65
+ ,33.77
+ ,33.51
+ ,34.26
+ ,34.21
+ ,34.13
+ ,34.73
+ ,34.73
+ ,34.57
+ ,34.80
+ ,33.98
+ ,34.40
+ ,34.21
+ ,34.61
+ ,35.25
+ ,35.23
+ ,35.00
+ ,34.52
+ ,33.82
+ ,34.35
+ ,34.81
+ ,34.96
+ ,36.69
+ ,36.42
+ ,36.44
+ ,37.41
+ ,36.40
+ ,36.15
+ ,35.78
+ ,36.95
+ ,36.14
+ ,36.36
+ ,37.31
+ ,37.58
+ ,38.00
+ ,37.23
+ ,37.00
+ ,37.87
+ ,37.70
+ ,36.17
+ ,36.56
+ ,37.70
+ ,38.77
+ ,39.02
+ ,39.88
+ ,39.56
+ ,38.52
+ ,37.20
+ ,38.58
+ ,39.41
+ ,39.08
+ ,38.81
+ ,38.73
+ ,38.70
+ ,39.23
+ ,39.82
+ ,39.97
+ ,40.37
+ ,39.54
+ ,39.21
+ ,39.07
+ ,39.78
+ ,39.40
+ ,38.92)
> ylimmax = ''
> ylimmin = ''
> main = 'Robustness of Central Tendency'
> 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] 28.58155
> sqrtn <- sqrt(length(x))
> (armse <- sd(x) / sqrtn)
[1] 0.3314371
> (armose <- arm / armse)
[1] 86.2352
> (geo <- geomean(x))
[1] 28.12291
> (har <- harmean(x))
[1] 27.70186
> (qua <- quamean(x))
[1] 29.07123
> (win <- winmean(x))
          [,1]      [,2]
 [1,] 28.58143 0.3311047
 [2,] 28.58174 0.3309293
 [3,] 28.58349 0.3306418
 [4,] 28.58349 0.3305113
 [5,] 28.58039 0.3298636
 [6,] 28.58016 0.3297854
 [7,] 28.57988 0.3290812
 [8,] 28.58329 0.3287601
 [9,] 28.57771 0.3279796
[10,] 28.57733 0.3278527
[11,] 28.57306 0.3270610
[12,] 28.57725 0.3266582
[13,] 28.57523 0.3263068
[14,] 28.56981 0.3256333
[15,] 28.56864 0.3244621
[16,] 28.57360 0.3236205
[17,] 28.58019 0.3226437
[18,] 28.57880 0.3223388
[19,] 28.56996 0.3212633
[20,] 28.57384 0.3201056
[21,] 28.53314 0.3149582
[22,] 28.52291 0.3136078
[23,] 28.51221 0.3115497
[24,] 28.51872 0.3111004
[25,] 28.50806 0.3097036
[26,] 28.49194 0.3076919
[27,] 28.48357 0.3063712
[28,] 28.48574 0.3046583
[29,] 28.48349 0.3042057
[30,] 28.46837 0.3010766
[31,] 28.48039 0.2992113
[32,] 28.44938 0.2956041
[33,] 28.45194 0.2925451
[34,] 28.44140 0.2905028
[35,] 28.44004 0.2901244
[36,] 28.44004 0.2896453
[37,] 28.44147 0.2885705
[38,] 28.41791 0.2853134
[39,] 28.43151 0.2839345
[40,] 28.43616 0.2833797
[41,] 28.38054 0.2772829
[42,] 28.29915 0.2681396
[43,] 28.30248 0.2673791
[44,] 28.26496 0.2633415
[45,] 28.26322 0.2623148
[46,] 28.23826 0.2595556
[47,] 28.23826 0.2592602
[48,] 28.23453 0.2573894
[49,] 28.23833 0.2571495
[50,] 28.21508 0.2548720
[51,] 28.21705 0.2534770
[52,] 28.21101 0.2522419
[53,] 28.19047 0.2495910
[54,] 28.18837 0.2480538
[55,] 28.16919 0.2462032
[56,] 28.17353 0.2442099
[57,] 28.18236 0.2436618
[58,] 28.16438 0.2419335
[59,] 28.13465 0.2383675
[60,] 28.09744 0.2348281
[61,] 28.09271 0.2332656
[62,] 28.03504 0.2270758
[63,] 27.82992 0.2074042
[64,] 27.69845 0.1949140
[65,] 27.69341 0.1921592
[66,] 27.66783 0.1895345
[67,] 27.66783 0.1891390
[68,] 27.65465 0.1839956
[69,] 27.65198 0.1825521
[70,] 27.58415 0.1767238
[71,] 27.54562 0.1726232
[72,] 27.55120 0.1718416
[73,] 27.60496 0.1679753
[74,] 27.60209 0.1668813
[75,] 27.58465 0.1645401
[76,] 27.56403 0.1597548
[77,] 27.56403 0.1597548
[78,] 27.54287 0.1575265
[79,] 27.46632 0.1511176
[80,] 27.44771 0.1486611
[81,] 27.38806 0.1423497
[82,] 27.14969 0.1226263
[83,] 27.09822 0.1171823
[84,] 27.03636 0.1119150
[85,] 27.01659 0.1080261
[86,] 26.98659 0.1038286
> (tri <- trimean(x))
          [,1]      [,2]
 [1,] 28.56437 0.3295606
 [2,] 28.54705 0.3279356
 [3,] 28.52929 0.3263186
 [4,] 28.51064 0.3247192
 [5,] 28.49169 0.3230698
 [6,] 28.47309 0.3214784
 [7,] 28.45422 0.3198146
 [8,] 28.43508 0.3181779
 [9,] 28.41517 0.3164971
[10,] 28.39559 0.3148298
[11,] 28.37572 0.3130838
[12,] 28.35594 0.3113319
[13,] 28.33543 0.3095220
[14,] 28.31474 0.3076443
[15,] 28.29412 0.3057250
[16,] 28.27323 0.3038029
[17,] 28.25161 0.3018407
[18,] 28.22914 0.2998427
[19,] 28.20636 0.2977498
[20,] 28.18372 0.2956196
[21,] 28.16042 0.2934481
[22,] 28.13902 0.2915466
[23,] 28.11778 0.2896279
[24,] 28.09671 0.2877383
[25,] 28.07490 0.2857562
[26,] 28.05320 0.2837459
[27,] 28.03186 0.2817468
[28,] 28.01050 0.2797071
[29,] 27.98860 0.2776433
[30,] 27.96636 0.2754711
[31,] 27.94434 0.2733592
[32,] 27.92134 0.2712150
[33,] 27.89917 0.2691614
[34,] 27.87642 0.2671503
[35,] 27.85362 0.2651245
[36,] 27.83038 0.2629753
[37,] 27.80663 0.2606985
[38,] 27.78231 0.2583190
[39,] 27.75833 0.2559737
[40,] 27.73331 0.2535280
[41,] 27.70756 0.2509270
[42,] 27.68322 0.2485295
[43,] 27.66122 0.2465356
[44,] 27.63859 0.2444251
[45,] 27.61673 0.2424037
[46,] 27.59440 0.2402784
[47,] 27.57238 0.2381558
[48,] 27.54981 0.2358736
[49,] 27.52681 0.2335146
[50,] 27.50310 0.2309670
[51,] 27.47955 0.2283571
[52,] 27.45532 0.2256062
[53,] 27.43066 0.2226959
[54,] 27.40600 0.2197119
[55,] 27.38074 0.2165599
[56,] 27.35541 0.2132571
[57,] 27.32924 0.2097768
[58,] 27.30204 0.2059950
[59,] 27.27464 0.2019858
[60,] 27.24739 0.1978675
[61,] 27.22051 0.1936337
[62,] 27.19299 0.1890877
[63,] 27.16644 0.1846059
[64,] 27.14554 0.1813486
[65,] 27.12813 0.1787675
[66,] 27.11032 0.1761175
[67,] 27.09274 0.1733984
[68,] 27.07459 0.1704155
[69,] 27.05625 0.1675090
[70,] 27.03737 0.1643913
[71,] 27.02000 0.1614473
[72,] 27.00325 0.1585289
[73,] 26.98571 0.1553255
[74,] 26.96582 0.1519998
[75,] 26.94528 0.1483581
[76,] 26.92453 0.1444863
[77,] 26.90365 0.1405840
[78,] 26.88196 0.1361527
[79,] 26.86010 0.1313624
[80,] 26.83990 0.1266966
[81,] 26.81948 0.1216432
[82,] 26.80021 0.1166543
[83,] 26.78826 0.1135037
[84,] 26.77756 0.1105418
[85,] 26.76852 0.1077846
[86,] 26.75977 0.1050657
> (midr <- midrange(x))
[1] 30.78
> midm <- array(NA,dim=8)
> for (j in 1:8) midm[j] <- midmean(x,j)
> midm
[1] 27.10713 27.14554 27.14554 27.14554 27.12813 27.14554 27.14554 27.14554
> postscript(file="/var/www/rcomp/tmp/1t8yt1324628841.ps",horizontal=F,onefile=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/rcomp/tmp/2161b1324628841.ps",horizontal=F,onefile=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/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/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/rcomp/tmp/3lsxv1324628841.tab") 
> 
> try(system("convert tmp/1t8yt1324628841.ps tmp/1t8yt1324628841.png",intern=TRUE))
character(0)
> try(system("convert tmp/2161b1324628841.ps tmp/2161b1324628841.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.610   0.050   1.665