R version 2.13.0 (2011-04-13)
Copyright (C) 2011 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(1.439
+ ,1.444
+ ,1.435
+ ,1.430
+ ,1.427
+ ,1.453
+ ,1.448
+ ,1.456
+ ,1.449
+ ,1.437
+ ,1.437
+ ,1.428
+ ,1.413
+ ,1.406
+ ,1.414
+ ,1.415
+ ,1.409
+ ,1.407
+ ,1.400
+ ,1.397
+ ,1.391
+ ,1.394
+ ,1.398
+ ,1.385
+ ,1.369
+ ,1.368
+ ,1.376
+ ,1.374
+ ,1.372
+ ,1.357
+ ,1.361
+ ,1.365
+ ,1.373
+ ,1.357
+ ,1.352
+ ,1.363
+ ,1.358
+ ,1.355
+ ,1.349
+ ,1.357
+ ,1.353
+ ,1.355
+ ,1.364
+ ,1.367
+ ,1.358
+ ,1.366
+ ,1.356
+ ,1.361
+ ,1.366
+ ,1.377
+ ,1.371
+ ,1.372
+ ,1.376
+ ,1.366
+ ,1.355
+ ,1.347
+ ,1.352
+ ,1.334
+ ,1.336
+ ,1.335
+ ,1.347
+ ,1.348
+ ,1.348
+ ,1.347
+ ,1.340
+ ,1.334
+ ,1.330
+ ,1.338
+ ,1.359
+ ,1.358
+ ,1.362
+ ,1.354
+ ,1.354
+ ,1.343
+ ,1.349
+ ,1.337
+ ,1.334
+ ,1.331
+ ,1.332
+ ,1.329
+ ,1.325
+ ,1.326
+ ,1.332
+ ,1.324
+ ,1.309
+ ,1.292
+ ,1.273
+ ,1.275
+ ,1.297
+ ,1.270
+ ,1.267
+ ,1.259
+ ,1.249
+ ,1.235
+ ,1.243
+ ,1.227
+ ,1.233
+ ,1.250
+ ,1.236
+ ,1.222
+ ,1.231
+ ,1.226
+ ,1.238
+ ,1.231
+ ,1.216
+ ,1.222
+ ,1.227
+ ,1.206
+ ,1.196
+ ,1.194
+ ,1.201
+ ,1.205
+ ,1.213
+ ,1.225
+ ,1.226
+ ,1.228
+ ,1.236
+ ,1.237
+ ,1.239
+ ,1.226
+ ,1.227
+ ,1.226
+ ,1.229
+ ,1.234
+ ,1.220
+ ,1.227
+ ,1.233
+ ,1.255
+ ,1.253
+ ,1.258
+ ,1.257
+ ,1.266
+ ,1.264
+ ,1.257
+ ,1.257
+ ,1.270
+ ,1.283
+ ,1.300
+ ,1.296
+ ,1.284
+ ,1.282
+ ,1.285
+ ,1.290
+ ,1.293
+ ,1.303
+ ,1.299
+ ,1.307
+ ,1.303
+ ,1.307
+ ,1.322
+ ,1.321
+ ,1.318
+ ,1.318
+ ,1.325
+ ,1.313
+ ,1.302
+ ,1.279
+ ,1.280
+ ,1.282
+ ,1.286
+ ,1.288
+ ,1.284
+ ,1.271
+ ,1.270
+ ,1.261
+ ,1.261
+ ,1.269
+ ,1.271
+ ,1.270
+ ,1.268
+ ,1.280
+ ,1.282
+ ,1.283
+ ,1.287
+ ,1.274
+ ,1.270
+ ,1.272
+ ,1.273
+ ,1.280
+ ,1.285
+ ,1.299
+ ,1.308
+ ,1.306
+ ,1.307
+ ,1.312
+ ,1.336
+ ,1.332
+ ,1.341
+ ,1.348
+ ,1.346
+ ,1.361
+ ,1.365
+ ,1.373
+ ,1.371
+ ,1.378
+ ,1.386
+ ,1.397
+ ,1.387
+ ,1.394
+ ,1.383
+ ,1.396
+ ,1.410
+ ,1.409
+ ,1.390
+ ,1.386
+ ,1.386
+ ,1.402
+ ,1.393
+ ,1.403
+ ,1.391
+ ,1.380
+ ,1.386
+ ,1.386
+ ,1.393
+ ,1.402
+ ,1.401
+ ,1.424
+ ,1.408
+ ,1.392
+ ,1.395
+ ,1.377
+ ,1.370
+ ,1.371
+ ,1.363
+ ,1.361
+ ,1.348
+ ,1.365
+ ,1.367
+ ,1.365
+ ,1.350
+ ,1.334
+ ,1.332
+ ,1.323
+ ,1.315
+ ,1.300
+ ,1.312
+ ,1.316
+ ,1.325
+ ,1.328
+ ,1.336
+ ,1.320
+ ,1.321
+ ,1.324
+ ,1.327
+ ,1.344
+ ,1.336
+ ,1.324
+ ,1.326
+ ,1.315
+ ,1.316
+ ,1.311
+ ,1.306
+ ,1.310
+ ,1.314
+ ,1.320
+ ,1.314
+ ,1.328
+ ,1.336)
> 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] 1.325771
> sqrtn <- sqrt(length(x))
> (armse <- sd(x) / sqrtn)
[1] 0.00373962
> (armose <- arm / armse)
[1] 354.5203
> (geo <- geomean(x))
[1] 1.324408
> (har <- harmean(x))
[1] 1.323036
> (qua <- quamean(x))
[1] 1.327126
> (win <- winmean(x))
          [,1]        [,2]
 [1,] 1.325767 0.003737007
 [2,] 1.325775 0.003727816
 [3,] 1.325810 0.003720358
 [4,] 1.325764 0.003710598
 [5,] 1.325802 0.003682189
 [6,] 1.325826 0.003668443
 [7,] 1.325934 0.003656003
 [8,] 1.325934 0.003641785
 [9,] 1.325760 0.003621855
[10,] 1.325798 0.003600418
[11,] 1.325798 0.003591099
[12,] 1.325659 0.003575983
[13,] 1.325205 0.003529260
[14,] 1.325151 0.003523912
[15,] 1.325151 0.003511887
[16,] 1.324965 0.003494023
[17,] 1.324899 0.003487812
[18,] 1.324899 0.003487812
[19,] 1.324899 0.003472910
[20,] 1.324899 0.003457328
[21,] 1.324981 0.003432372
[22,] 1.324725 0.003409196
[23,] 1.324814 0.003382303
[24,] 1.324814 0.003382303
[25,] 1.324814 0.003363524
[26,] 1.324814 0.003344115
[27,] 1.324709 0.003315040
[28,] 1.324601 0.003305748
[29,] 1.324713 0.003294064
[30,] 1.324713 0.003272152
[31,] 1.324713 0.003249642
[32,] 1.325085 0.003188897
[33,] 1.325853 0.003113592
[34,] 1.325853 0.003089785
[35,] 1.326260 0.003051335
[36,] 1.326399 0.003013588
[37,] 1.326543 0.002975204
[38,] 1.326543 0.002975204
[39,] 1.326391 0.002962518
[40,] 1.326081 0.002910019
[41,] 1.326081 0.002882692
[42,] 1.326407 0.002853432
[43,] 1.326407 0.002853432
[44,] 1.326919 0.002808310
[45,] 1.327267 0.002778137
[46,] 1.327267 0.002748259
[47,] 1.327085 0.002703186
[48,] 1.326713 0.002643010
[49,] 1.326523 0.002597353
[50,] 1.326329 0.002582482
[51,] 1.326329 0.002582482
[52,] 1.326128 0.002567171
[53,] 1.326128 0.002567171
[54,] 1.325919 0.002518011
[55,] 1.325705 0.002502253
[56,] 1.325922 0.002483853
[57,] 1.325922 0.002448999
[58,] 1.325922 0.002448999
[59,] 1.325922 0.002413129
[60,] 1.326155 0.002393724
[61,] 1.327101 0.002316174
[62,] 1.327101 0.002279141
[63,] 1.326857 0.002261345
[64,] 1.326609 0.002243446
[65,] 1.326860 0.002184901
[66,] 1.326860 0.002184901
[67,] 1.326601 0.002166430
[68,] 1.326864 0.002145389
[69,] 1.326864 0.002145389
[70,] 1.326864 0.002104651
[71,] 1.326864 0.002104651
[72,] 1.327143 0.002082607
[73,] 1.327143 0.002082607
[74,] 1.327143 0.002039832
[75,] 1.327143 0.001996678
[76,] 1.327438 0.001973704
[77,] 1.327736 0.001906614
[78,] 1.328039 0.001839164
[79,] 1.328345 0.001815895
[80,] 1.329275 0.001746167
[81,] 1.329589 0.001722973
[82,] 1.329589 0.001631545
[83,] 1.329267 0.001609085
[84,] 1.329593 0.001585346
[85,] 1.329593 0.001585346
[86,] 1.329926 0.001513986
> (tri <- trimean(x))
          [,1]        [,2]
 [1,] 1.325777 0.003698483
 [2,] 1.325787 0.003658074
 [3,] 1.325794 0.003620633
 [4,] 1.325788 0.003584087
 [5,] 1.325794 0.003548486
 [6,] 1.325793 0.003517673
 [7,] 1.325787 0.003488014
 [8,] 1.325764 0.003458994
 [9,] 1.325742 0.003430679
[10,] 1.325739 0.003403691
[11,] 1.325733 0.003378044
[12,] 1.325726 0.003352160
[13,] 1.325733 0.003326589
[14,] 1.325778 0.003304506
[15,] 1.325829 0.003281745
[16,] 1.325881 0.003258867
[17,] 1.325946 0.003236313
[18,] 1.326018 0.003213028
[19,] 1.326091 0.003188431
[20,] 1.326165 0.003163648
[21,] 1.326241 0.003138675
[22,] 1.326313 0.003114182
[23,] 1.326401 0.003089954
[24,] 1.326486 0.003066263
[25,] 1.326572 0.003041156
[26,] 1.326660 0.003015856
[27,] 1.326750 0.002990359
[28,] 1.326847 0.002965253
[29,] 1.326950 0.002939189
[30,] 1.327051 0.002912284
[31,] 1.327153 0.002885117
[32,] 1.327258 0.002857676
[33,] 1.327349 0.002832433
[34,] 1.327411 0.002810382
[35,] 1.327473 0.002788333
[36,] 1.327522 0.002767197
[37,] 1.327565 0.002746914
[38,] 1.327604 0.002727531
[39,] 1.327644 0.002706782
[40,] 1.327691 0.002685308
[41,] 1.327750 0.002665381
[42,] 1.327810 0.002645605
[43,] 1.327860 0.002626124
[44,] 1.327912 0.002605182
[45,] 1.327946 0.002585372
[46,] 1.327970 0.002565839
[47,] 1.327994 0.002546537
[48,] 1.328025 0.002528330
[49,] 1.328069 0.002512095
[50,] 1.328120 0.002497089
[51,] 1.328179 0.002481603
[52,] 1.328240 0.002464760
[53,] 1.328309 0.002447305
[54,] 1.328380 0.002428313
[55,] 1.328459 0.002410531
[56,] 1.328548 0.002392025
[57,] 1.328632 0.002372993
[58,] 1.328718 0.002354324
[59,] 1.328807 0.002333894
[60,] 1.328899 0.002313778
[61,] 1.328985 0.002293022
[62,] 1.329045 0.002275391
[63,] 1.329106 0.002258309
[64,] 1.329177 0.002240529
[65,] 1.329258 0.002221972
[66,] 1.329333 0.002205333
[67,] 1.329411 0.002186853
[68,] 1.329500 0.002167497
[69,] 1.329583 0.002147521
[70,] 1.329669 0.002125296
[71,] 1.329759 0.002103493
[72,] 1.329851 0.002079175
[73,] 1.329937 0.002053824
[74,] 1.330027 0.002025503
[75,] 1.330120 0.001997260
[76,] 1.330217 0.001969124
[77,] 1.330308 0.001939578
[78,] 1.330392 0.001912174
[79,] 1.330470 0.001887142
[80,] 1.330541 0.001860944
[81,] 1.330583 0.001837694
[82,] 1.330617 0.001813364
[83,] 1.330652 0.001793871
[84,] 1.330700 0.001773395
[85,] 1.330739 0.001752107
[86,] 1.330779 0.001727496
> (midr <- midrange(x))
[1] 1.325
> midm <- array(NA,dim=8)
> for (j in 1:8) midm[j] <- midmean(x,j)
> midm
[1] 1.328130 1.328432 1.328432 1.328432 1.329258 1.328432 1.328432 1.328432
> postscript(file="/var/wessaorg/rcomp/tmp/1hyqs1324629240.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/wessaorg/rcomp/tmp/2dwpr1324629240.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/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/wessaorg/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/wessaorg/rcomp/tmp/39gpr1324629240.tab") 
> 
> try(system("convert tmp/1hyqs1324629240.ps tmp/1hyqs1324629240.png",intern=TRUE))
character(0)
> try(system("convert tmp/2dwpr1324629240.ps tmp/2dwpr1324629240.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.470   0.154   1.999