R version 2.7.2 (2008-08-25)
Copyright (C) 2008 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.171742608
+ ,-0.166771574
+ ,-0.161414962
+ ,-0.156036741
+ ,-0.150896429
+ ,-0.146148477
+ ,-0.141861106
+ ,-0.138041489
+ ,-0.134659209
+ ,-0.131664423
+ ,-0.129000169
+ ,-0.126609769
+ ,-0.124440755
+ ,-0.122446621
+ ,-0.120587352
+ ,-0.118829328
+ ,-0.117144949
+ ,-0.115512116
+ ,-0.113913662
+ ,-0.112336736
+ ,-0.11077218
+ ,-0.109213906
+ ,-0.107658285
+ ,-0.106103586
+ ,-0.104549458
+ ,-0.102996482
+ ,-0.101445791
+ ,-0.099898764
+ ,-0.098356785
+ ,-0.096821068
+ ,-0.095292541
+ ,-0.093771779
+ ,-0.092258981
+ ,-0.090753982
+ ,-0.089256295
+ ,-0.087765174
+ ,-0.086279694
+ ,-0.084798837
+ ,-0.083321583
+ ,-0.081847005
+ ,-0.080374348
+ ,-0.078903103
+ ,-0.077433068
+ ,-0.075964391
+ ,-0.074497596
+ ,-0.073033588
+ ,-0.071573645
+ ,-0.070119386
+ ,-0.068672724
+ ,-0.067235811
+ ,-0.06581096
+ ,-0.064400572
+ ,-0.063007045
+ ,-0.061632691
+ ,-0.060279647
+ ,-0.058949796
+ ,-0.057644694
+ ,-0.056365505
+ ,-0.055112952
+ ,-0.053887281
+ ,-0.052688237
+ ,-0.051515057
+ ,-0.050366484
+ ,-0.049240785
+ ,-0.048135791
+ ,-0.047048946
+ ,-0.045977366
+ ,-0.044917907
+ ,-0.043867232
+ ,-0.042821894
+ ,-0.041778406
+ ,-0.040733318
+ ,-0.039683286
+ ,-0.038625135
+ ,-0.037555919
+ ,-0.036472968
+ ,-0.035373926
+ ,-0.034256782
+ ,-0.033119887
+ ,-0.031961962
+ ,-0.030782094
+ ,-0.029579729
+ ,-0.028354652
+ ,-0.027106954
+ ,-0.025837007
+ ,-0.024545422
+ ,-0.023233011
+ ,-0.021900741
+ ,-0.020549693
+ ,-0.019181022
+ ,-0.017795916
+ ,-0.016395557
+ ,-0.014981093
+ ,-0.013553607
+ ,-0.012114096
+ ,-0.010663457
+ ,-0.009202472
+ ,-0.007731803
+ ,-0.006251996
+ ,-0.004763485
+ ,-0.0032666
+ ,-0.001761583
+ ,-0.000248607
+ ,0.001272204
+ ,0.00280075
+ ,0.004336931
+ ,0.005880618
+ ,0.007431626
+ ,0.008989688
+ ,0.010554426
+ ,0.012125323
+ ,0.013701706
+ ,0.015282717
+ ,0.016867302
+ ,0.018454193
+ ,0.020041903
+ ,0.021628719
+ ,0.023212709
+ ,0.024791729
+ ,0.02636344
+ ,0.027925328
+ ,0.029474738
+ ,0.031008908
+ ,0.032525005
+ ,0.03402018
+ ,0.035491611
+ ,0.03693656
+ ,0.038352428
+ ,0.039736813
+ ,0.041087566
+ ,0.042402847
+ ,0.043681187
+ ,0.044921534
+ ,0.046123318
+ ,0.047286493
+ ,0.048411595
+ ,0.04949979
+ ,0.050552922
+ ,0.051573562
+ ,0.052565052
+ ,0.053531549
+ ,0.054478062
+ ,0.055410479
+ ,0.05633559
+ ,0.057261084
+ ,0.058195538
+ ,0.05914837
+ ,0.060129767
+ ,0.061150572
+ ,0.062222132
+ ,0.063356098
+ ,0.064564184
+ ,0.065857871
+ ,0.067248073
+ ,0.068744771
+ ,0.070356616
+ ,0.072090529
+ ,0.073951308
+ ,0.075941275
+ ,0.078059967
+ ,0.080303931
+ ,0.082666614
+ ,0.085138392
+ ,0.087706761
+ ,0.090356687
+ ,0.09307113
+ ,0.095831738
+ ,0.098619689
+ ,0.101416639
+ ,0.104205746
+ ,0.106972694
+ ,0.109706651
+ ,0.112401075
+ ,0.11505428
+ ,0.117669692
+ ,0.120255712
+ ,0.122825156
+ ,0.125394269
+ ,0.127981369
+ ,0.13060524
+ ,0.133283479
+ ,0.136031066
+ ,0.138859509
+ ,0.141776921
+ ,0.144789371
+ ,0.147903689
+ ,0.151131651
+ ,0.154494962
+ ,0.158029827
+ ,0.161789098
+ ,0.165839243
+ ,0.170249032
+ ,0.175067431
+ ,0.180290421
+ ,0.185820847
+ ,0.191432001
+ ,0.196753033)
> 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: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.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] 1.528408e-05
> sqrtn <- sqrt(length(x))
> (armse <- sd(x) / sqrtn)
[1] 0.006206152
> (armose <- arm / armse)
[1] 0.00246273
> (geo <- geomean(x))
[1] NaN
Warning message:
In log(x) : NaNs produced
> (har <- harmean(x))
[1] -0.05005614
> (qua <- quamean(x))
[1] 0.08688613
> (win <- winmean(x))
               [,1]        [,2]
 [1,]  1.350744e-05 0.006198326
 [2,]  1.092325e-05 0.006182115
 [3,]  8.605411e-06 0.006158618
 [4,]  6.926670e-06 0.006129670
 [5,]  5.138675e-06 0.006097111
 [6,]  1.410208e-06 0.006062391
 [7,] -6.781142e-06 0.006026458
 [8,] -2.209042e-05 0.005989824
 [9,] -4.676408e-05 0.005952695
[10,] -8.224921e-05 0.005915106
[11,] -1.290166e-04 0.005877020
[12,] -1.865986e-04 0.005838386
[13,] -2.537971e-04 0.005799155
[14,] -3.289951e-04 0.005759277
[15,] -4.104991e-04 0.005718689
[16,] -4.968510e-04 0.005677295
[17,] -5.870637e-04 0.005634958
[18,] -6.807566e-04 0.005591501
[19,] -7.781846e-04 0.005546714
[20,] -8.801701e-04 0.005500369
[21,] -9.879598e-04 0.005452238
[22,] -1.103030e-03 0.005402118
[23,] -1.226869e-03 0.005349848
[24,] -1.360767e-03 0.005295329
[25,] -1.505621e-03 0.005238533
[26,] -1.661788e-03 0.005179511
[27,] -1.828985e-03 0.005118393
[28,] -2.006242e-03 0.005055384
[29,] -2.191906e-03 0.004990751
[30,] -2.383696e-03 0.004924811
[31,] -2.578799e-03 0.004857911
[32,] -2.773990e-03 0.004790417
[33,] -2.965780e-03 0.004722690
[34,] -3.150567e-03 0.004655076
[35,] -3.324796e-03 0.004587890
[36,] -3.485097e-03 0.004521406
[37,] -3.628422e-03 0.004455852
[38,] -3.752151e-03 0.004391404
[39,] -3.854182e-03 0.004328185
[40,] -3.932989e-03 0.004266269
[41,] -3.987655e-03 0.004205681
[42,] -4.017889e-03 0.004146402
[43,] -4.024005e-03 0.004088377
[44,] -4.006898e-03 0.004031522
[45,] -3.967992e-03 0.003975726
[46,] -3.909183e-03 0.003920859
[47,] -3.832768e-03 0.003866781
[48,] -3.741373e-03 0.003813343
[49,] -3.637874e-03 0.003760391
[50,] -3.525322e-03 0.003707774
[51,] -3.406868e-03 0.003655342
[52,] -3.285692e-03 0.003602951
[53,] -3.164933e-03 0.003550464
[54,] -3.047631e-03 0.003497748
[55,] -2.936672e-03 0.003444681
[56,] -2.834738e-03 0.003391145
[57,] -2.744264e-03 0.003337031
[58,] -2.667403e-03 0.003282233
[59,] -2.605998e-03 0.003226653
[60,] -2.561557e-03 0.003170198
[61,] -2.535242e-03 0.003112777
[62,] -2.527855e-03 0.003054306
[63,] -2.539840e-03 0.002994703
[64,] -2.571285e-03 0.002933891
[65,] -2.621933e-03 0.002871798
> (tri <- trimean(x))
               [,1]        [,2]
 [1,]  1.528408e-05 0.006124496
 [2,] -1.128178e-04 0.006045923
 [3,] -3.720721e-04 0.005971170
 [4,] -3.720721e-04 0.005900412
 [5,] -6.389864e-04 0.005833434
 [6,] -7.761676e-04 0.005769811
 [7,] -9.156784e-04 0.005709039
 [8,] -9.156784e-04 0.005650629
 [9,] -1.199371e-03 0.005594148
[10,] -1.341909e-03 0.005539232
[11,] -1.483711e-03 0.005485587
[12,] -1.623950e-03 0.005432976
[13,] -1.761941e-03 0.005381205
[14,] -1.897173e-03 0.005330112
[15,] -2.029308e-03 0.005279561
[16,] -2.029308e-03 0.005229440
[17,] -2.283648e-03 0.005179655
[18,] -2.405763e-03 0.005130136
[19,] -2.524500e-03 0.005080838
[20,] -2.639828e-03 0.005031739
[21,] -2.751652e-03 0.004982846
[22,] -2.859790e-03 0.004934188
[23,] -2.963968e-03 0.004885821
[24,] -3.063825e-03 0.004837820
[25,] -3.158922e-03 0.004790278
[26,] -3.248770e-03 0.004743295
[27,] -3.332857e-03 0.004696982
[28,] -3.410678e-03 0.004651444
[29,] -3.481766e-03 0.004606783
[30,] -3.545723e-03 0.004563089
[31,] -3.602246e-03 0.004520439
[32,] -3.602246e-03 0.004478888
[33,] -3.692369e-03 0.004438475
[34,] -3.725993e-03 0.004399216
[35,] -3.752246e-03 0.004361106
[36,] -3.771493e-03 0.004324122
[37,] -3.784235e-03 0.004288224
[38,] -3.791091e-03 0.004253352
[39,] -3.792787e-03 0.004219439
[40,] -3.790137e-03 0.004186403
[41,] -3.784019e-03 0.004154156
[42,] -3.775360e-03 0.004122605
[43,] -3.765112e-03 0.004091651
[44,] -3.754230e-03 0.004061198
[45,] -3.743657e-03 0.004031145
[46,] -3.734304e-03 0.004001395
[47,] -3.727033e-03 0.003971853
[48,] -3.722645e-03 0.003942424
[49,] -3.721868e-03 0.003913019
[50,] -3.725350e-03 0.003883551
[51,] -3.733646e-03 0.003853933
[52,] -3.747218e-03 0.003824086
[53,] -3.747218e-03 0.003793929
[54,] -3.791553e-03 0.003763386
[55,] -3.822748e-03 0.003732383
[56,] -3.860086e-03 0.003700845
[57,] -3.860086e-03 0.003668702
[58,] -3.953009e-03 0.003635884
[59,] -4.008283e-03 0.003602322
[60,] -4.069091e-03 0.003567948
[61,] -4.135087e-03 0.003532696
[62,] -4.205864e-03 0.003496500
[63,] -4.280959e-03 0.003459296
[64,] -4.280959e-03 0.003421019
[65,] -4.442035e-03 0.003381608
> (midr <- midrange(x))
[1] 0.01250521
> midm <- array(NA,dim=8)
> for (j in 1:8) midm[j] <- midmean(x,j)
> midm
[1] -0.004373416 -0.003721868 -0.003721868 -0.003721868 -0.003721868
[6] -0.004371377 -0.003721868 -0.003721868
> postscript(file="/var/www/html/rcomp/tmp/1pvyu1225038815.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/2q5ij1225038815.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 
> 
> #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> 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/327ud1225038815.tab") 
> 
> system("convert tmp/1pvyu1225038815.ps tmp/1pvyu1225038815.png")
> system("convert tmp/2q5ij1225038815.ps tmp/2q5ij1225038815.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.302   0.405   1.469