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(-14.81039
+ ,-3.28274
+ ,-4.85329
+ ,-25.62907
+ ,-2.05413
+ ,-0.96706
+ ,9.96965
+ ,-6.81537
+ ,-3.85900
+ ,7.68954
+ ,-39.33000
+ ,-15.16607
+ ,-0.13154
+ ,1.55439
+ ,-1.21081
+ ,2.59094
+ ,4.25931
+ ,-0.34667
+ ,-21.90876
+ ,-4.51886
+ ,18.64956
+ ,-8.09796
+ ,27.19089
+ ,-4.93466
+ ,5.54314
+ ,14.67181
+ ,4.35485
+ ,-35.36298
+ ,21.09116
+ ,8.54014
+ ,-21.07678
+ ,0.70893
+ ,3.25688
+ ,-13.23521
+ ,-8.35929
+ ,-4.65221
+ ,-5.06926
+ ,-30.82941
+ ,7.17970
+ ,1.29943
+ ,-3.55868
+ ,-1.94040
+ ,-5.68342
+ ,10.24959
+ ,10.51949
+ ,-5.62233
+ ,17.90070
+ ,22.44965
+ ,9.43289
+ ,32.61133
+ ,-6.81711
+ ,-0.65314
+ ,4.75309
+ ,-13.44323
+ ,-4.80908
+ ,-0.04435
+ ,7.05527
+ ,2.01473
+ ,2.62446
+ ,-0.63780
+ ,-10.64846
+ ,-11.45022
+ ,6.12018
+ ,-8.18714
+ ,-4.23073
+ ,39.26613
+ ,4.87331
+ ,-15.42067
+ ,34.03244
+ ,3.53316
+ ,26.54658
+ ,2.94135
+ ,-3.46076
+ ,-0.05704
+ ,7.35246
+ ,-7.07442
+ ,-6.55026
+ ,-22.66062
+ ,7.13510
+ ,0.38424
+ ,-4.83275
+ ,-2.25130
+ ,-31.85808
+ ,10.77965
+ ,-31.82712
+ ,6.39089
+ ,0.01278
+ ,-2.19794
+ ,14.51179
+ ,3.65711
+ ,10.82303
+ ,9.07264
+ ,-11.37350
+ ,10.46728
+ ,-14.97325
+ ,30.67256
+ ,-18.92760
+ ,15.30698
+ ,-34.03504
+ ,-7.24521
+ ,10.61113
+ ,0.77871
+ ,-0.10387
+ ,-16.26651
+ ,-11.33622
+ ,2.84565
+ ,3.50492
+ ,6.10891
+ ,34.72095
+ ,4.36559
+ ,39.31253
+ ,2.45698
+ ,36.57613
+ ,-4.15898
+ ,-1.42988
+ ,-13.17862
+ ,0.68764
+ ,0.64259
+ ,9.81365
+ ,0.58821
+ ,9.06359
+ ,-21.45345
+ ,8.82159
+ ,-14.05202
+ ,-17.35022
+ ,-9.66790
+ ,3.34815
+ ,1.07363
+ ,-1.42596
+ ,-21.61176
+ ,4.76899
+ ,-7.45835
+ ,17.52726
+ ,-12.23760
+ ,-0.31758
+ ,-15.25193
+ ,-19.92141
+ ,-6.64785
+ ,-12.81371
+ ,0.64526
+ ,6.08556
+ ,5.70685
+ ,2.45812
+ ,-4.70240
+ ,-10.01302
+ ,15.30227
+ ,6.66024
+ ,-27.41800
+ ,-6.75046
+ ,-5.71698
+ ,28.08686
+ ,5.69599
+ ,2.03406
+ ,13.22204
+ ,8.88650
+ ,18.84657
+ ,25.04196
+ ,10.71160
+ ,1.13755
+ ,7.45398
+ ,-1.49701
+ ,12.55633
+ ,-14.74332
+ ,-1.99871
+ ,-10.90802
+ ,4.10447
+ ,-1.50770
+ ,-47.28901
+ ,13.63685
+ ,10.28751
+ ,-2.36258
+ ,-23.95407
+ ,60.80184
+ ,4.20974
+ ,-1.68746
+ ,29.84660
+ ,32.39293
+ ,7.24087
+ ,-5.94028
+ ,6.67141
+ ,16.13431
+ ,6.33344
+ ,-18.60496
+ ,-19.82071
+ ,16.01946
+ ,28.64296
+ ,60.62326
+ ,-17.61620
+ ,-21.91786
+ ,3.66145
+ ,-9.01293
+ ,1.02439
+ ,-14.67385
+ ,-5.32970
+ ,-2.78519
+ ,3.96277
+ ,-6.08144
+ ,-40.03321
+ ,31.72985
+ ,-15.47080
+ ,-1.91390
+ ,-35.43998
+ ,-8.19525
+ ,-30.12629
+ ,69.17736
+ ,-16.80747
+ ,-8.38222
+ ,-0.33511
+ ,-12.91487
+ ,-3.86879
+ ,-2.06974
+ ,-5.61551
+ ,11.97373
+ ,-8.03816
+ ,-6.42337
+ ,-8.63532
+ ,-4.51673
+ ,4.60513
+ ,1.38954
+ ,5.13314
+ ,-5.77963
+ ,11.59462
+ ,-3.09008
+ ,9.48020
+ ,-5.30047
+ ,6.13748
+ ,-21.69430
+ ,0.09540
+ ,-3.44514
+ ,10.83315
+ ,2.85724
+ ,15.20208
+ ,-8.97473
+ ,1.11793
+ ,-0.96618
+ ,-6.49897
+ ,-4.79496
+ ,-1.66741
+ ,-10.54649
+ ,-3.67153
+ ,-13.90539)
> 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] -8.298755e-08
> sqrtn <- sqrt(length(x))
> (armse <- sd(x) / sqrtn)
[1] 1.036054
> (armose <- arm / armse)
[1] -8.009964e-08
> (geo <- geomean(x))
[1] NaN
Warning message:
In log(x) : NaNs produced
> (har <- harmean(x))
[1] 8.696741
> (qua <- quamean(x))
[1] 16.05048
> (win <- winmean(x))
               [,1]      [,2]
 [1,] -0.0046462241 1.0215830
 [2,] -0.0002924481 1.0202707
 [3,] -0.2171477593 0.9566651
 [4,] -0.2166398755 0.9563367
 [5,] -0.2448983817 0.9428270
 [6,] -0.2368872614 0.9275439
 [7,] -0.2559862241 0.9243002
 [8,] -0.2700409959 0.9124442
 [9,] -0.2519394191 0.9075879
[10,] -0.1670759336 0.8885766
[11,] -0.1336815768 0.8712148
[12,] -0.0914057261 0.8552506
[13,] -0.0865612033 0.8379037
[14,] -0.0757178838 0.8285164
[15,] -0.1309172614 0.8206365
[16,] -0.1594550207 0.8131699
[17,] -0.2597677178 0.7982515
[18,] -0.4415602490 0.7721654
[19,] -0.5189651452 0.7558022
[20,] -0.6093568465 0.7234887
[21,] -0.6177490041 0.7205955
[22,] -0.6045809544 0.7040960
[23,] -0.6094290871 0.6968862
[24,] -0.6496803734 0.6716566
[25,] -0.6340029876 0.6675256
[26,] -0.6523141079 0.6535486
[27,] -0.5922363071 0.6473161
[28,] -0.5114290041 0.6369341
[29,] -0.5692051867 0.6298459
[30,] -0.5681197095 0.6257830
[31,] -0.6696196266 0.6135807
[32,] -0.6990954772 0.6057508
[33,] -0.7679504564 0.5949410
[34,] -0.8406808299 0.5864310
[35,] -0.8856492946 0.5804606
[36,] -0.9065083817 0.5612766
[37,] -0.8855504149 0.5589502
[38,] -0.8195187137 0.5513392
[39,] -0.7968679668 0.5471740
[40,] -0.8041509544 0.5448261
[41,] -0.7748708299 0.5392375
[42,] -0.7663401245 0.5367926
[43,] -0.6956238589 0.5244748
[44,] -0.5587929876 0.5108896
[45,] -0.5967386307 0.5050342
[46,] -0.6193988797 0.5018136
[47,] -0.6009206639 0.4888182
[48,] -0.5586468050 0.4834890
[49,] -0.6111601660 0.4754229
[50,] -0.5023593361 0.4657345
[51,] -0.4668011203 0.4563555
[52,] -0.3394852697 0.4434576
[53,] -0.3929800830 0.4374748
[54,] -0.5075205809 0.4153141
[55,] -0.5035176763 0.4062203
[56,] -0.5217792531 0.4038687
[57,] -0.5093740664 0.3986003
[58,] -0.5221436929 0.3972525
[59,] -0.5112299170 0.3946187
[60,] -0.5162166390 0.3918245
[61,] -0.4666193361 0.3724200
[62,] -0.4146602490 0.3679246
[63,] -0.4404248963 0.3588648
[64,] -0.3873500415 0.3524245
[65,] -0.4397330290 0.3482090
[66,] -0.4266946058 0.3464821
[67,] -0.4013013278 0.3440716
[68,] -0.3803539419 0.3414740
[69,] -0.4740965975 0.3318806
[70,] -0.4552924481 0.3299778
[71,] -0.3995883817 0.3189458
[72,] -0.4799057261 0.3063833
[73,] -0.5099478008 0.2967981
[74,] -0.5227427386 0.2929602
[75,] -0.5172468880 0.2918168
[76,] -0.5446415768 0.2868953
[77,] -0.6189960996 0.2810440
[78,] -0.5299692946 0.2739836
[79,] -0.5517057676 0.2709624
[80,] -0.4914103320 0.2641701
> (tri <- trimean(x))
             [,1]      [,2]
 [1,] -0.09158314 0.9838732
 [2,] -0.17998734 0.9431846
 [3,] -0.27212877 0.8998312
 [4,] -0.29108502 0.8785317
 [5,] -0.31050199 0.8559585
 [6,] -0.32431026 0.8352863
 [7,] -0.33977938 0.8165508
 [8,] -0.35260107 0.7972291
 [9,] -0.36375408 0.7786782
[10,] -0.37730226 0.7596710
[11,] -0.40043676 0.7421608
[12,] -0.42736931 0.7258404
[13,] -0.45875195 0.7104648
[14,] -0.49114559 0.6961538
[15,] -0.52503796 0.6819713
[16,] -0.55533560 0.6677250
[17,] -0.58414213 0.6533197
[18,] -0.60657376 0.6394825
[19,] -0.61745724 0.6272902
[20,] -0.62367264 0.6158301
[21,] -0.62453950 0.6065106
[22,] -0.62493508 0.5968390
[23,] -0.62607851 0.5879227
[24,] -0.62698244 0.5790135
[25,] -0.62578911 0.5715226
[26,] -0.62537016 0.5638587
[27,] -0.62403460 0.5567456
[28,] -0.62556881 0.5496204
[29,] -0.63093721 0.5427718
[30,] -0.63377155 0.5359706
[31,] -0.63377155 0.5289912
[32,] -0.63527282 0.5223916
[33,] -0.63252617 0.5158747
[34,] -0.62680936 0.5096394
[35,] -0.61794404 0.5035317
[36,] -0.60703669 0.4973796
[37,] -0.59503192 0.4920597
[38,] -0.58356345 0.4864979
[39,] -0.57438276 0.4810291
[40,] -0.56584335 0.4754198
[41,] -0.55681314 0.4695300
[42,] -0.54864911 0.4635559
[43,] -0.54059019 0.4572714
[44,] -0.53491105 0.4513232
[45,] -0.53404477 0.4458108
[46,] -0.53179134 0.4402337
[47,] -0.52866898 0.4344047
[48,] -0.52611393 0.4289639
[49,] -0.52497168 0.4234162
[50,] -0.52196525 0.4179326
[51,] -0.52264511 0.4126288
[52,] -0.52457131 0.4074900
[53,] -0.53092541 0.4027544
[54,] -0.53564165 0.3979910
[55,] -0.53659969 0.3942895
[56,] -0.53772341 0.3908282
[57,] -0.53826370 0.3871709
[58,] -0.53924088 0.3834963
[59,] -0.53981846 0.3795295
[60,] -0.54078355 0.3753303
[61,] -0.54161277 0.3708856
[62,] -0.54161277 0.3673707
[63,] -0.54852183 0.3637758
[64,] -0.55218124 0.3604301
[65,] -0.55777306 0.3571517
[66,] -0.56178826 0.3537923
[67,] -0.56639850 0.3501440
[68,] -0.57205429 0.3462239
[69,] -0.57865049 0.3420140
[70,] -0.58226614 0.3380733
[71,] -0.58668182 0.3337703
[72,] -0.59322887 0.3298340
[73,] -0.59722168 0.3264511
[74,] -0.60031978 0.3234067
[75,] -0.60309615 0.3202330
[76,] -0.60619573 0.3166673
[77,] -0.60843931 0.3130116
[78,] -0.60805059 0.3093426
[79,] -0.61095723 0.3057683
[80,] -0.61318877 0.3018999
> (midr <- midrange(x))
[1] 10.94417
> midm <- array(NA,dim=8)
> for (j in 1:8) midm[j] <- midmean(x,j)
> midm
[1] -0.6040840 -0.5407836 -0.5407836 -0.5407836 -0.5407836 -0.6027276 -0.5407836
[8] -0.5407836
> postscript(file="/var/fisher/rcomp/tmp/18fm11355078018.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/2dnsk1355078018.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/39oqw1355078018.tab") 
> 
> try(system("convert tmp/18fm11355078018.ps tmp/18fm11355078018.png",intern=TRUE))
character(0)
> try(system("convert tmp/2dnsk1355078018.ps tmp/2dnsk1355078018.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  2.218   0.352   2.556