R version 2.6.0 (2007-10-03)
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(-183.923544514060
+ ,-177.072609114252
+ ,-228.635109114254
+ ,-237.447609114242
+ ,-127.760109114251
+ ,-193.010109114250
+ ,-220.635109114248
+ ,-164.510109114256
+ ,-268.322609114246
+ ,-333.697609114249
+ ,-34.2601091142503
+ ,-154.885109114249
+ ,-97.7452805255622
+ ,101.105654874272
+ ,2.54315487427317
+ ,-43.2693451257275
+ ,-163.581845125727
+ ,-162.831845125727
+ ,46.5431548742728
+ ,26.6681548742733
+ ,-107.144345125727
+ ,42.4806548742729
+ ,76.918154874273
+ ,196.293154874273
+ ,201.43298346296
+ ,12.2839188627955
+ ,-0.278581137204542
+ ,42.9089188627948
+ ,87.5964188627952
+ ,84.3464188627952
+ ,57.721418862795
+ ,173.846418862796
+ ,-185.966081137205
+ ,47.6589188627952
+ ,89.0964188627953
+ ,-68.5285811372049
+ ,272.611247451482
+ ,146.462182851318
+ ,162.899682851318
+ ,10.0871828513171
+ ,279.774682851318
+ ,212.524682851317
+ ,248.899682851317
+ ,-41.9753171486821
+ ,-5.78781714868267
+ ,52.8371828513174
+ ,274.274682851318
+ ,414.649682851317
+ ,310.789511440004
+ ,362.64044683984
+ ,26.0779468398400
+ ,403.265446839839
+ ,327.95294683984
+ ,193.702946839840
+ ,317.07794683984
+ ,202.20294683984
+ ,321.39044683984
+ ,178.015446839840
+ ,16.4529468398398
+ ,-68.1720531601603
+ ,-157.032224571473
+ ,-76.1812891716376
+ ,-81.7437891716377
+ ,-134.556289171638
+ ,77.1312108283621
+ ,199.881210828362
+ ,105.256210828362
+ ,198.381210828362
+ ,262.568710828362
+ ,196.193710828362
+ ,11.6312108283621
+ ,-145.993789171638
+ ,-166.853960582951
+ ,-202.003025183115
+ ,43.4344748168846
+ ,-113.378025183116
+ ,-113.690525183116
+ ,-155.940525183116
+ ,-210.565525183116
+ ,-124.440525183115
+ ,-64.2530251831158
+ ,-298.628025183116
+ ,-154.190525183116
+ ,23.1844748168843
+ ,-249.675696594429
+ ,118.175238805407
+ ,-180.387261194593
+ ,-79.1997611945937
+ ,-81.5122611945933
+ ,-246.762261194593
+ ,-105.387261194593
+ ,-319.262261194593
+ ,-72.0747611945935
+ ,-90.4497611945934
+ ,-80.0122611945933
+ ,119.362738805407
+ ,-53.4974326059064
+ ,-114.646497206071
+ ,-155.208997206071
+ ,-50.0214972060714
+ ,-196.333997206071
+ ,-14.5839972060711
+ ,-82.2089972060712
+ ,17.9160027939294
+ ,-162.896497206071
+ ,-132.271497206071
+ ,-16.8339972060710
+ ,81.5410027939288
+ ,275.680831382616
+ ,-32.4682332175485
+ ,17.9692667824515
+ ,27.1567667824508
+ ,-123.155733217549
+ ,108.594266782451
+ ,67.9692667824511
+ ,34.0942667824517
+ ,-13.7182332175489
+ ,-113.093233217549
+ ,54.3442667824512
+ ,149.719266782451
+ ,153.859095371138
+ ,-28.2899692290262
+ ,238.147530770974
+ ,50.3350307709731
+ ,8.02253077097358
+ ,-61.2274692290265
+ ,-140.852469229027
+ ,-28.7274692290261
+ ,9.46003077097336
+ ,-121.914969229026
+ ,41.5225307709736
+ ,115.897530770973
+ ,27.0373593596605
+ ,-91.111705240504
+ ,3.32579475949605
+ ,-29.4867052405046
+ ,-73.7992052405041
+ ,50.9507947594957
+ ,-86.6742052405044
+ ,-9.54920524050376
+ ,-66.3617052405043
+ ,73.2632947594957
+ ,-216.299205240504
+ ,-128.924205240504
+ ,-142.784376651817
+ ,27.0665587480184
+ ,60.5040587480183
+ ,35.6915587480177
+ ,16.3790587480182
+ ,-64.870941251982
+ ,115.504058748018
+ ,-30.3709412519815
+ ,87.816558748018
+ ,205.441558748018
+ ,-64.1209412519819
+ ,-322.745941251982
+ ,-139.606112663295
+ ,35.2448227365406
+ ,-4.31767726345942
+ ,17.8698227365400
+ ,2.55732273654045
+ ,129.307322736540
+ ,-16.3176772634598
+ ,164.807322736541
+ ,21.9948227365402
+ ,138.619822736540
+ ,87.0573227365404
+ ,51.4323227365403
+ ,-80.4278486747727
+ ,-105.191879672279
+ ,5.24562032772056
+ ,68.4331203277199
+ ,-0.879379672279542
+ ,-105.129379672280
+ ,-82.7543796722797
+ ,-132.629379672279
+ ,102.558120327720
+ ,23.1831203277203
+ ,-180.379379672280
+ ,-267.00437967228
+ ,30.1354489164073
+ ,23.986384316243
+ ,90.4238843162428
+ ,31.6113843162422
+ ,81.2988843162427
+ ,25.0488843162426
+ ,-13.5761156837575
+ ,33.548884316243
+ ,140.736384316242
+ ,132.361384316243
+ ,94.7988843162427
+ ,4.17388431624253)
> 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] -1.323944e-15
> sqrtn <- sqrt(length(x))
> (armse <- sd(x) / sqrtn)
[1] 10.61887
> (armose <- arm / armse)
[1] -1.246784e-16
> (geo <- geomean(x))
[1] NaN
Warning message:
In log(x) : NaNs produced
> (har <- harmean(x))
[1] -66.34536
> (qua <- quamean(x))
[1] 146.7557
> (win <- winmean(x))
               [,1]      [,2]
 [1,]  -0.002252959 10.597662
 [2,]  -0.389141709 10.511457
 [3,]  -0.608723958 10.367590
 [4,]  -0.114079874 10.254262
 [5,]  -0.192055669 10.231234
 [6,]   0.152952071 10.127553
 [7,]  -0.871577889  9.939375
 [8,]  -0.654044531  9.864732
 [9,]  -0.306871806  9.803922
[10,]   0.023157601  9.741110
[11,]  -0.303784901  9.628517
[12,]  -0.799744145  9.466618
[13,]  -0.948001838  9.300903
[14,]  -2.402967876  9.011702
[15,]  -2.696658189  8.914183
[16,]  -2.379540184  8.816482
[17,]  -2.266864427  8.788561
[18,]  -2.080816551  8.735274
[19,]  -2.228474108  8.717265
[20,]  -2.101524670  8.654449
[21,]  -0.994736680  8.537128
[22,]  -1.011570406  8.475485
[23,]  -2.779603887  8.243205
[24,]  -3.215063894  8.174980
[25,]  -4.383611296  8.042551
[26,]  -3.856572247  7.933909
[27,]  -4.974384635  7.780441
[28,]  -5.471428474  7.705499
[29,]  -5.914463262  7.648922
[30,]  -6.700590544  7.545410
[31,]  -5.718895714  7.376806
[32,]  -6.227066697  7.217817
[33,]  -6.419936943  7.132453
[34,]  -7.960248039  6.938235
[35,]  -7.256178652  6.826758
[36,]  -7.321953378  6.750616
[37,]  -7.328811782  6.736674
[38,]  -8.033889907  6.543360
[39,]  -8.475475497  6.457589
[40,]  -8.345997699  6.339329
[41,]  -8.381802975  6.284710
[42,]  -9.489991912  6.134379
[43,]  -8.841974740  5.892917
[44,]  -8.927108651  5.845573
[45,]  -9.153833678  5.812430
[46,]  -9.138344117  5.801421
[47,]  -7.814072126  5.654573
[48,]  -8.052527111  5.554944
[49,]  -8.718629678  5.488362
[50,]  -8.765405323  5.481437
[51,]  -7.911042258  5.207604
[52,]  -6.172151605  5.042260
[53,]  -6.998321145  4.939680
[54,]  -7.294932592  4.731352
[55,]  -6.304940606  4.622255
[56,]  -8.323223064  4.428173
[57,]  -9.011210644  4.349298
[58,]  -9.961451259  4.261883
[59,] -10.091334661  4.196141
[60,] -10.400482360  4.150382
[61,] -10.295329790  4.115809
[62,]  -9.519455321  4.015219
[63,]  -9.615933251  3.879181
[64,]  -9.413039898  3.800766
> (tri <- trimean(x))
            [,1]      [,2]
 [1,] -0.4260635 10.356801
 [2,] -0.8588914 10.098127
 [3,] -1.1013428  9.868962
 [4,] -1.2726885  9.679170
 [5,] -1.5782556  9.509815
 [6,] -1.8739783  9.334054
 [7,] -2.2383703  9.167682
 [8,] -2.4513769  9.024919
 [9,] -2.6992848  8.884428
[10,] -2.9960182  8.743022
[11,] -3.3370074  8.600100
[12,] -3.6521474  8.461753
[13,] -3.9270779  8.333669
[14,] -4.1953624  8.215794
[15,] -4.3470995  8.121891
[16,] -4.4791348  8.031224
[17,] -4.6385977  7.943492
[18,] -4.8103069  7.851554
[19,] -4.9993625  7.757384
[20,] -5.1835767  7.657273
[21,] -5.3808280  7.554779
[22,] -5.6517834  7.454287
[23,] -5.9291561  7.350825
[24,] -6.1117388  7.259543
[25,] -6.2749318  7.166304
[26,] -6.3786842  7.076398
[27,] -6.5136467  6.987594
[28,] -6.5941310  6.903728
[29,] -6.6515827  6.818353
[30,] -6.6885542  6.729664
[31,] -6.6879617  6.641136
[32,] -6.7348520  6.557940
[33,] -6.7590322  6.479553
[34,] -6.7749428  6.400208
[35,] -6.7200782  6.328449
[36,] -6.6955707  6.258139
[37,] -6.6672597  6.186405
[38,] -6.6376655  6.107946
[39,] -6.5757830  6.036700
[40,] -6.4922800  5.964352
[41,] -6.4113905  5.893447
[42,] -6.3259526  5.818705
[43,] -6.1894981  5.747504
[44,] -6.0756171  5.687990
[45,] -5.9536282  5.624706
[46,] -5.8170861  5.555895
[47,] -5.6756306  5.478922
[48,] -5.5846331  5.405291
[49,] -5.4796163  5.330653
[50,] -5.3416636  5.251607
[51,] -5.1955840  5.161904
[52,] -5.0794147  5.086557
[53,] -5.0324993  5.016509
[54,] -4.9477199  4.945658
[55,] -4.8459437  4.884499
[56,] -4.7822784  4.824367
[57,] -4.6266325  4.773174
[58,] -4.4323022  4.720218
[59,] -4.1849591  4.665765
[60,] -3.9180042  4.607617
[61,] -3.6216624  4.542605
[62,] -3.3127559  4.468493
[63,] -3.0215325  4.392676
[64,] -2.7075134  4.318479
> (midr <- midrange(x))
[1] 40.47604
> midm <- array(NA,dim=8)
> for (j in 1:8) midm[j] <- midmean(x,j)
> midm
[1] -6.631640 -5.584633 -6.631640 -5.584633 -5.584633 -6.631640 -5.584633
[8] -5.675631
> postscript(file="/var/www/html/rcomp/tmp/1pl5m1197492050.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/27ddf1197492050.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/3razp1197492051.tab") 
> 
> system("convert tmp/1pl5m1197492050.ps tmp/1pl5m1197492050.png")
> system("convert tmp/27ddf1197492050.ps tmp/27ddf1197492050.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.487   0.368   1.626