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(0.0269573944933363
+ ,0.0252543244856407
+ ,0.0234663563270153
+ ,0.0216705100164706
+ ,0.019898101790598
+ ,0.0181658374794494
+ ,0.0164973985000073
+ ,0.0149320643170635
+ ,0.0135200236899052
+ ,0.0123129179368339
+ ,0.0113530379065303
+ ,0.0106667931378226
+ ,0.0102555318123724
+ ,0.0100957841561224
+ ,0.0101439448988436
+ ,0.0103452471356505
+ ,0.0106454297830635
+ ,0.0109956496156571
+ ,0.0113567511168418
+ ,0.0116995743651474
+ ,0.0120038016021861
+ ,0.0122577893985301
+ ,0.0124566565619037
+ ,0.0126007410198417
+ ,0.0126938786083839
+ ,0.0127440867974884
+ ,0.0127603271073115
+ ,0.0127518866517497
+ ,0.012730091554938
+ ,0.0127033824211204
+ ,0.0126820474023146
+ ,0.0126728978301890
+ ,0.0126835863914231
+ ,0.0127200720059746
+ ,0.0127867341866925
+ ,0.0128848132514507
+ ,0.0130168589419847
+ ,0.0131824919279555
+ ,0.0133792638239547
+ ,0.0136034381117238
+ ,0.0138502575294115
+ ,0.0141130119758550
+ ,0.0143843141006615
+ ,0.0146561929705374
+ ,0.0149204497755614
+ ,0.0151679099347271
+ ,0.0153906854905040
+ ,0.0155801772476093
+ ,0.0157301860536652
+ ,0.0158344564820207
+ ,0.0158878669513172
+ ,0.0158880790329646
+ ,0.0158340406664293
+ ,0.0157246203179212
+ ,0.0155622721024531
+ ,0.0153502961637555
+ ,0.0150927825940047
+ ,0.0147967443636773
+ ,0.0144685093678446
+ ,0.0141153426366991
+ ,0.0137473184287137
+ ,0.0133719165690060
+ ,0.0129991295856794
+ ,0.0126365920656438
+ ,0.0122936356537062
+ ,0.0119770702098744
+ ,0.0116946193511517
+ ,0.0114515525892263
+ ,0.0112529429085974
+ ,0.0111014197365865
+ ,0.010999659627745
+ ,0.0109477734358669
+ ,0.0109464198725327
+ ,0.0109931007694523
+ ,0.0110853666437776
+ ,0.0112205638661813
+ ,0.0113934494815145
+ ,0.0116009018973178
+ ,0.0118377710108104
+ ,0.0120995035156501
+ ,0.0123815303981410
+ ,0.0126789822256679
+ ,0.0129872303156284
+ ,0.0133021864947266
+ ,0.0136196681678444
+ ,0.0139359989666353
+ ,0.0142476384269921
+ ,0.0145523172936869
+ ,0.0148462726424358
+ ,0.0151286665705089
+ ,0.0153968574410418
+ ,0.0156514568224672
+ ,0.0158888458140516
+ ,0.0161107155798414
+ ,0.0163170168398825
+ ,0.0165072626562962
+ ,0.0166835964786100
+ ,0.0168461130125906
+ ,0.0169958250296389
+ ,0.0171360070011546
+ ,0.0172662956343239
+ ,0.0173893712601468
+ ,0.0175061026742448
+ ,0.0176186506159529
+ ,0.0177271432733880
+ ,0.0178325809743100
+ ,0.0179354792646741
+ ,0.0180353593201978
+ ,0.0181320155244646
+ ,0.0182230128596788
+ ,0.0183077662975503
+ ,0.0183828934562379
+ ,0.0184463808146833
+ ,0.0184952850532389
+ ,0.0185246727014246
+ ,0.0185319229683955
+ ,0.0185138103656993
+ ,0.0184666748473341
+ ,0.0183855781746181
+ ,0.0182702374435566
+ ,0.0181164266371800
+ ,0.0179232845951886
+ ,0.0176897336530753
+ ,0.0174158042740473
+ ,0.0171020082981199
+ ,0.0167501453461214
+ ,0.0163625277823769
+ ,0.0159423552840085
+ ,0.0154943700489033
+ ,0.0150223367314834
+ ,0.0145319815828029
+ ,0.014029792621004
+ ,0.0135203292570985
+ ,0.0130112917810563
+ ,0.0125093891760375
+ ,0.0120211829058130
+ ,0.0115538848908445
+ ,0.0111151347888026
+ ,0.0107123284371675
+ ,0.0103529021597579
+ ,0.0100453762165042
+ ,0.0097974605179079
+ ,0.00961840694887808
+ ,0.00951648107667327
+ ,0.00950056056842361
+ ,0.00957855807009369
+ ,0.00975759221647046
+ ,0.0100447763823385
+ ,0.0104431663572113
+ ,0.01095514391496
+ ,0.0115801587363620
+ ,0.0123143847965651
+ ,0.0131521161314924
+ ,0.0140846710486466
+ ,0.0150986330666529
+ ,0.0161783375672453
+ ,0.0173052540392484
+ ,0.0184564955443016
+ ,0.0196074476016893
+ ,0.0207299446723496
+ ,0.0217947446358812
+ ,0.0227713498747957
+ ,0.023629841040952
+ ,0.0243411998027985
+ ,0.0248804178650263
+ ,0.0252257777054356
+ ,0.0253629722171458
+ ,0.0252836819176938
+ ,0.0249884910510709
+ ,0.0244879883611272
+ ,0.023802452480233
+ ,0.0229617545827468
+ ,0.0220061083343501
+ ,0.0209839523582308
+ ,0.0199516220897085
+ ,0.0189704331793296
+ ,0.0181027036227214
+ ,0.0174061270858631
+ ,0.0169250032337177
+ ,0.0166821398953368
+ ,0.0166699999186468
+ ,0.0168486913540258
+ ,0.0171506672653043
+ ,0.0174905503729930
+ ,0.0177773514719896
+ ,0.0179306304679235
+ ,0.0178932968299711
+ ,0.0176501245954888
+ ,0.0172482830632603
+ ,0.0168193889721676
+ ,0.0166009255666754
+ ,0.0169256553778928
+ ,0.0181401809021045
+ ,0.0204502954920460
+ ,0.023808545265327
+ ,0.0279169589287818
+ ,0.0322865204651344)
> 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] 0.01561351
> sqrtn <- sqrt(length(x))
> (armse <- sd(x) / sqrtn)
[1] 0.0002960484
> (armose <- arm / armse)
[1] 52.73971
> (geo <- geomean(x))
[1] 0.01511069
> (har <- harmean(x))
[1] 0.01464902
> (qua <- quamean(x))
[1] 0.01615425
> (win <- winmean(x))
            [,1]         [,2]
 [1,] 0.01559141 0.0002904444
 [2,] 0.01558230 0.0002883421
 [3,] 0.01555862 0.0002836908
 [4,] 0.01555984 0.0002831090
 [5,] 0.01556010 0.0002828731
 [6,] 0.01556677 0.0002819533
 [7,] 0.01555836 0.0002804912
 [8,] 0.01555602 0.0002795367
 [9,] 0.01554029 0.0002763105
[10,] 0.01553850 0.0002745252
[11,] 0.01551377 0.0002692629
[12,] 0.01551386 0.0002691589
[13,] 0.01550843 0.0002667997
[14,] 0.01551119 0.0002636300
[15,] 0.01547439 0.0002576691
[16,] 0.01546262 0.0002550404
[17,] 0.01541679 0.0002437473
[18,] 0.01539760 0.0002410965
[19,] 0.01538633 0.0002394161
[20,] 0.01532048 0.0002299980
[21,] 0.01529368 0.0002266136
[22,] 0.01526289 0.0002228035
[23,] 0.01521468 0.0002150929
[24,] 0.01521012 0.0002141711
[25,] 0.01517497 0.0002099532
[26,] 0.01510481 0.0001998923
[27,] 0.01504915 0.0001937260
[28,] 0.01506235 0.0001922199
[29,] 0.01506129 0.0001920192
[30,] 0.01506406 0.0001912120
[31,] 0.01506870 0.0001899084
[32,] 0.01508367 0.0001881540
[33,] 0.01508638 0.0001875786
[34,] 0.01507946 0.0001862827
[35,] 0.01509564 0.0001846640
[36,] 0.01508281 0.0001833393
[37,] 0.01510172 0.0001802862
[38,] 0.01511948 0.0001770123
[39,] 0.01511345 0.0001755250
[40,] 0.01511177 0.0001747449
[41,] 0.01512638 0.0001731321
[42,] 0.01515680 0.0001698625
[43,] 0.01516163 0.0001689168
[44,] 0.01515089 0.0001672137
[45,] 0.01512841 0.0001651916
[46,] 0.01514296 0.0001637374
[47,] 0.01515913 0.0001620513
[48,] 0.01516467 0.0001603267
[49,] 0.01517229 0.0001571098
[50,] 0.01516737 0.0001551441
[51,] 0.01516377 0.0001532508
[52,] 0.01515551 0.0001522770
[53,] 0.01514567 0.0001513062
[54,] 0.01513747 0.0001505441
[55,] 0.01510892 0.0001476823
[56,] 0.01510720 0.0001470911
[57,] 0.01509040 0.0001449090
[58,] 0.01509050 0.0001444297
[59,] 0.01508968 0.0001436702
[60,] 0.01506643 0.0001413910
[61,] 0.01505698 0.0001402004
[62,] 0.01505962 0.0001390511
[63,] 0.01505977 0.0001339470
[64,] 0.01508828 0.0001308257
[65,] 0.01508099 0.0001296098
> (tri <- trimean(x))
            [,1]         [,2]
 [1,] 0.01561351 0.0002847896
 [2,] 0.01555935 0.0002787240
 [3,] 0.01549793 0.0002733992
 [4,] 0.01549793 0.0002694950
 [5,] 0.01545498 0.0002654911
 [6,] 0.01543259 0.0002612622
 [7,] 0.01540852 0.0002569201
 [8,] 0.01540852 0.0002525391
 [9,] 0.01536173 0.0002479966
[10,] 0.01533965 0.0002436261
[11,] 0.01531726 0.0002391846
[12,] 0.01529692 0.0002351325
[13,] 0.01527609 0.0002307812
[14,] 0.01525526 0.0002263816
[15,] 0.01523369 0.0002220040
[16,] 0.01523369 0.0002179713
[17,] 0.01519579 0.0002139061
[18,] 0.01517989 0.0002107136
[19,] 0.01516490 0.0002075372
[20,] 0.01515028 0.0002042709
[21,] 0.01513946 0.0002016518
[22,] 0.01513001 0.0001991336
[23,] 0.01512213 0.0001967532
[24,] 0.01511681 0.0001948455
[25,] 0.01511160 0.0001928519
[26,] 0.01510815 0.0001910379
[27,] 0.01510833 0.0001898636
[28,] 0.01511139 0.0001890370
[29,] 0.01511387 0.0001882305
[30,] 0.01511648 0.0001873460
[31,] 0.01511903 0.0001864209
[32,] 0.01511903 0.0001854866
[33,] 0.01512321 0.0001845737
[34,] 0.01512492 0.0001835950
[35,] 0.01512699 0.0001825948
[36,] 0.01512840 0.0001815948
[37,] 0.01513043 0.0001805689
[38,] 0.01513169 0.0001796408
[39,] 0.01513222 0.0001788298
[40,] 0.01513303 0.0001780123
[41,] 0.01513394 0.0001771346
[42,] 0.01513427 0.0001762511
[43,] 0.01513331 0.0001754797
[44,] 0.01513212 0.0001746563
[45,] 0.01513134 0.0001738316
[46,] 0.01513146 0.0001730273
[47,] 0.01513098 0.0001721993
[48,] 0.01512982 0.0001713606
[49,] 0.01512837 0.0001705115
[50,] 0.01512655 0.0001697657
[51,] 0.01512486 0.0001690327
[52,] 0.01512324 0.0001683058
[53,] 0.01512324 0.0001675081
[54,] 0.01512090 0.0001666293
[55,] 0.01512021 0.0001656392
[56,] 0.01512068 0.0001647080
[57,] 0.01512068 0.0001636346
[58,] 0.01512257 0.0001625457
[59,] 0.01512395 0.0001612743
[60,] 0.01512544 0.0001598198
[61,] 0.01512802 0.0001583064
[62,] 0.01513116 0.0001566081
[63,] 0.01513437 0.0001546918
[64,] 0.01513437 0.0001529778
[65,] 0.01514002 0.0001512686
> (midr <- midrange(x))
[1] 0.02089354
> midm <- array(NA,dim=8)
> for (j in 1:8) midm[j] <- midmean(x,j)
> midm
[1] 0.01510078 0.01512837 0.01512837 0.01512837 0.01512837 0.01510218 0.01512837
[8] 0.01512837
> postscript(file="/var/www/html/rcomp/tmp/15hda1193075209.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/2h8141193075209.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/31s9m1193075210.tab") 
> 
> system("convert tmp/15hda1193075209.ps tmp/15hda1193075209.png")
> system("convert tmp/2h8141193075209.ps tmp/2h8141193075209.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.053   0.345   1.130