R version 2.15.2 (2012-10-26) -- "Trick or Treat"
Copyright (C) 2012 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i686-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(210907
+ ,120982
+ ,176508
+ ,179321
+ ,123185
+ ,52746
+ ,385534
+ ,33170
+ ,101645
+ ,149061
+ ,165446
+ ,237213
+ ,173326
+ ,133131
+ ,258873
+ ,180083
+ ,324799
+ ,230964
+ ,236785
+ ,135473
+ ,202925
+ ,215147
+ ,344297
+ ,153935
+ ,132943
+ ,174724
+ ,174415
+ ,225548
+ ,223632
+ ,124817
+ ,221698
+ ,210767
+ ,170266
+ ,260561
+ ,84853
+ ,294424
+ ,101011
+ ,215641
+ ,325107
+ ,7176
+ ,167542
+ ,106408
+ ,96560
+ ,265769
+ ,269651
+ ,149112
+ ,175824
+ ,152871
+ ,111665
+ ,116408
+ ,362301
+ ,78800
+ ,183167
+ ,277965
+ ,150629
+ ,168809
+ ,24188
+ ,329267
+ ,65029
+ ,101097
+ ,218946
+ ,244052
+ ,341570
+ ,103597
+ ,233328
+ ,256462
+ ,206161
+ ,311473
+ ,235800
+ ,177939
+ ,207176
+ ,196553
+ ,174184
+ ,143246
+ ,187559
+ ,187681
+ ,119016
+ ,182192
+ ,73566
+ ,194979
+ ,167488
+ ,143756
+ ,275541
+ ,243199
+ ,182999
+ ,135649
+ ,152299
+ ,120221
+ ,346485
+ ,145790
+ ,193339
+ ,80953
+ ,122774
+ ,130585
+ ,112611
+ ,286468
+ ,241066
+ ,148446
+ ,204713
+ ,182079
+ ,140344
+ ,220516
+ ,243060
+ ,162765
+ ,182613
+ ,232138
+ ,265318
+ ,85574
+ ,310839
+ ,225060
+ ,232317
+ ,144966
+ ,43287
+ ,155754
+ ,164709
+ ,201940
+ ,235454
+ ,220801
+ ,99466
+ ,92661
+ ,133328
+ ,61361
+ ,125930
+ ,100750
+ ,224549
+ ,82316
+ ,102010
+ ,101523
+ ,243511
+ ,22938
+ ,41566
+ ,152474
+ ,61857
+ ,99923
+ ,132487
+ ,317394
+ ,21054
+ ,209641
+ ,22648
+ ,31414
+ ,46698
+ ,131698
+ ,91735
+ ,244749
+ ,184510
+ ,79863
+ ,128423
+ ,97839
+ ,38214
+ ,151101
+ ,272458
+ ,172494
+ ,108043
+ ,328107
+ ,250579
+ ,351067
+ ,158015
+ ,98866
+ ,85439
+ ,229242
+ ,351619
+ ,84207
+ ,120445
+ ,324598
+ ,131069
+ ,204271
+ ,165543
+ ,141722
+ ,116048
+ ,250047
+ ,299775
+ ,195838
+ ,173260
+ ,254488
+ ,104389
+ ,136084
+ ,199476
+ ,92499
+ ,224330
+ ,135781
+ ,74408
+ ,81240
+ ,14688
+ ,181633
+ ,271856
+ ,7199
+ ,46660
+ ,17547
+ ,133368
+ ,95227
+ ,152601
+ ,98146
+ ,79619
+ ,59194
+ ,139942
+ ,118612
+ ,72880
+ ,65475
+ ,99643
+ ,71965
+ ,77272
+ ,49289
+ ,135131
+ ,108446
+ ,89746
+ ,44296
+ ,77648
+ ,181528
+ ,134019
+ ,124064
+ ,92630
+ ,121848
+ ,52915
+ ,81872
+ ,58981
+ ,53515
+ ,60812
+ ,56375
+ ,65490
+ ,80949
+ ,76302
+ ,104011
+ ,98104
+ ,67989
+ ,30989
+ ,135458
+ ,73504
+ ,63123
+ ,61254
+ ,74914
+ ,31774
+ ,81437
+ ,87186
+ ,50090
+ ,65745
+ ,56653
+ ,158399
+ ,46455
+ ,73624
+ ,38395
+ ,91899
+ ,139526
+ ,52164
+ ,51567
+ ,70551
+ ,84856
+ ,102538
+ ,86678
+ ,85709
+ ,34662
+ ,150580
+ ,99611
+ ,19349
+ ,99373
+ ,86230
+ ,30837
+ ,31706
+ ,89806
+ ,62088
+ ,40151
+ ,27634
+ ,76990
+ ,37460
+ ,54157
+ ,49862
+ ,84337
+ ,64175
+ ,59382
+ ,119308
+ ,76702
+ ,103425
+ ,70344
+ ,43410
+ ,104838
+ ,62215
+ ,69304
+ ,53117
+ ,19764
+ ,86680
+ ,84105
+ ,77945
+ ,89113
+ ,91005
+ ,40248
+ ,64187
+ ,50857
+ ,56613
+ ,62792
+ ,72535)
> 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] 138088.1
> sqrtn <- sqrt(length(x))
> (armse <- sd(x) / sqrtn)
[1] 4843.494
> (armose <- arm / armse)
[1] 28.51001
> (geo <- geomean(x))
[1] 112398.6
> (har <- harmean(x))
[1] 83412.91
> (qua <- quamean(x))
[1] 160700.4
> (win <- winmean(x))
          [,1]     [,2]
 [1,] 138007.7 4829.874
 [2,] 137985.6 4813.476
 [3,] 138009.6 4809.984
 [4,] 137971.1 4798.167
 [5,] 137940.4 4791.865
 [6,] 137910.6 4781.162
 [7,] 137651.2 4734.917
 [8,] 137627.1 4729.739
 [9,] 137572.6 4713.489
[10,] 137681.2 4702.192
[11,] 137795.5 4691.355
[12,] 137502.7 4650.083
[13,] 137255.4 4613.209
[14,] 137238.9 4608.061
[15,] 136668.1 4534.343
[16,] 136449.2 4491.562
[17,] 136068.9 4428.391
[18,] 135713.6 4353.352
[19,] 135603.8 4331.496
[20,] 135403.0 4306.762
[21,] 135486.9 4292.022
[22,] 135326.4 4273.032
[23,] 135122.3 4231.531
[24,] 135227.8 4216.621
[25,] 134826.9 4172.240
[26,] 134754.8 4150.362
[27,] 134731.2 4111.939
[28,] 134559.8 4090.894
[29,] 134171.4 4051.074
[30,] 134385.1 4025.611
[31,] 133878.3 3965.154
[32,] 133826.4 3955.817
[33,] 133852.2 3943.428
[34,] 133899.0 3933.798
[35,] 133954.5 3926.932
[36,] 133778.6 3897.843
[37,] 133306.9 3849.587
[38,] 133277.2 3842.393
[39,] 133198.0 3826.100
[40,] 133239.0 3815.233
[41,] 133252.0 3764.804
[42,] 133139.7 3748.823
[43,] 133119.0 3745.956
[44,] 133294.7 3704.667
[45,] 133059.7 3677.909
[46,] 132501.7 3622.902
[47,] 132654.9 3599.619
[48,] 132643.4 3586.980
[49,] 132624.4 3582.425
[50,] 132589.5 3565.786
[51,] 132288.9 3532.852
[52,] 132150.4 3517.120
[53,] 132204.0 3505.256
[54,] 131972.4 3475.426
[55,] 131543.7 3407.386
[56,] 131450.3 3399.032
[57,] 130780.1 3316.780
[58,] 130841.5 3308.272
[59,] 130614.7 3288.830
[60,] 130155.9 3242.731
[61,] 130415.3 3192.658
[62,] 130386.7 3148.049
[63,] 130517.1 3124.935
[64,] 130264.9 3097.515
[65,] 130361.3 3058.304
[66,] 129928.8 3004.188
[67,] 129331.1 2944.778
[68,] 129309.7 2921.704
[69,] 129119.4 2904.541
[70,] 128736.3 2872.458
[71,] 127538.8 2752.552
[72,] 127634.5 2741.813
[73,] 127214.9 2660.435
[74,] 126973.5 2628.007
[75,] 127004.6 2619.811
[76,] 126977.3 2607.372
[77,] 126965.3 2592.444
[78,] 127014.9 2584.899
[79,] 127126.7 2560.523
[80,] 127324.4 2543.563
[81,] 126987.8 2509.244
[82,] 127079.7 2473.348
[83,] 126684.0 2444.269
[84,] 126351.4 2408.667
[85,] 126208.2 2390.377
[86,] 126010.3 2358.427
[87,] 126051.0 2343.084
[88,] 126525.4 2302.961
[89,] 126292.6 2282.003
[90,] 126312.5 2277.940
[91,] 126233.8 2250.317
[92,] 125525.5 2199.897
[93,] 125244.2 2154.846
[94,] 124876.0 2123.185
[95,] 124902.6 2119.095
[96,] 124429.6 2063.077
> (tri <- trimean(x))
          [,1]     [,2]
 [1,] 137682.0 4778.367
 [2,] 137351.7 4724.496
 [3,] 137028.0 4676.948
 [4,] 136691.5 4628.474
 [5,] 136360.2 4581.068
 [6,] 136030.4 4532.893
 [7,] 135701.1 4484.541
 [8,] 135406.2 4441.890
 [9,] 135110.1 4398.015
[10,] 134816.1 4354.313
[11,] 134506.0 4309.916
[12,] 134179.9 4264.601
[13,] 133875.6 4221.726
[14,] 133587.8 4180.655
[15,] 133296.8 4138.060
[16,] 133044.0 4100.421
[17,] 132802.8 4064.779
[18,] 132583.3 4032.813
[19,] 132383.1 4005.388
[20,] 132186.4 3978.288
[21,] 131998.2 3951.696
[22,] 131802.2 3924.786
[23,] 131611.7 3897.852
[24,] 131428.7 3872.403
[25,] 131428.7 3846.579
[26,] 131062.2 3822.359
[27,] 130887.5 3798.238
[28,] 130710.9 3775.202
[29,] 130539.0 3752.211
[30,] 130380.9 3730.403
[31,] 130211.0 3708.815
[32,] 130059.0 3689.579
[33,] 129906.4 3669.749
[34,] 129750.1 3649.436
[35,] 129589.1 3628.440
[36,] 129422.9 3606.553
[37,] 129260.3 3584.984
[38,] 129111.9 3564.816
[39,] 128961.8 3543.790
[40,] 128811.6 3522.346
[41,] 128657.1 3500.128
[42,] 128499.1 3479.175
[43,] 128341.8 3457.739
[44,] 128182.0 3435.070
[45,] 128013.3 3413.047
[46,] 127848.8 3391.023
[47,] 127698.8 3370.518
[48,] 127540.9 3349.778
[49,] 127380.1 3328.262
[50,] 127380.1 3305.496
[51,] 127050.4 3282.023
[52,] 126889.9 3258.711
[53,] 126730.2 3234.609
[54,] 126565.2 3209.399
[55,] 126403.6 3184.028
[56,] 126251.0 3160.542
[57,] 126097.7 3135.786
[58,] 125960.4 3113.733
[59,] 125818.2 3090.428
[60,] 125679.2 3066.438
[61,] 125679.2 3043.196
[62,] 125410.4 3020.746
[63,] 125268.1 2998.901
[64,] 125118.5 2976.474
[65,] 124972.3 2953.769
[66,] 124819.7 2931.285
[67,] 124675.4 2910.013
[68,] 124544.1 2890.350
[69,] 124410.0 2870.209
[70,] 124277.6 2849.267
[71,] 124152.4 2828.358
[72,] 124057.3 2812.601
[73,] 123956.9 2795.846
[74,] 123865.5 2782.090
[75,] 123778.1 2768.724
[76,] 123687.4 2754.331
[77,] 123687.4 2739.076
[78,] 123499.6 2723.020
[79,] 123400.2 2705.629
[80,] 123294.5 2687.778
[81,] 123179.9 2668.925
[82,] 123071.2 2650.247
[83,] 122956.3 2631.698
[84,] 122849.0 2613.015
[85,] 122747.8 2594.579
[86,] 122647.2 2575.183
[87,] 122549.0 2555.679
[88,] 122446.0 2534.808
[89,] 122325.3 2513.998
[90,] 122207.1 2492.154
[91,] 122083.9 2467.774
[92,] 121958.4 2442.428
[93,] 121849.6 2418.111
[94,] 121745.2 2394.286
[95,] 121647.9 2369.946
[96,] 121545.9 2342.521
> (midr <- midrange(x))
[1] 196355
> midm <- array(NA,dim=8)
> for (j in 1:8) midm[j] <- midmean(x,j)
> midm
[1] 123616.4 124057.3 124057.3 124057.3 124057.3 123717.3 124057.3 124057.3
> postscript(file="/var/fisher/rcomp/tmp/1jg3k1353335556.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/fisher/rcomp/tmp/2hovr1353335556.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/fisher/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/fisher/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/fisher/rcomp/tmp/3w76l1353335556.tab") 
> 
> try(system("convert tmp/1jg3k1353335556.ps tmp/1jg3k1353335556.png",intern=TRUE))
character(0)
> try(system("convert tmp/2hovr1353335556.ps tmp/2hovr1353335556.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  2.567   0.327   2.874