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(115.01
+ ,124.56
+ ,128.60
+ ,131.54
+ ,124.93
+ ,126.77
+ ,128.97
+ ,149.91
+ ,149.55
+ ,149.91
+ ,159.47
+ ,162.04
+ ,163.88
+ ,166.82
+ ,173.80
+ ,172.06
+ ,178.21
+ ,174.53
+ ,176.37
+ ,168.65
+ ,166.82
+ ,164.24
+ ,163.51
+ ,164.61
+ ,163.51
+ ,164.61
+ ,167.92
+ ,160.94
+ ,156.53
+ ,159.47
+ ,149.18
+ ,135.58
+ ,124.93
+ ,116.84
+ ,114.27
+ ,114.64
+ ,112.80
+ ,113.17
+ ,111.33
+ ,113.91
+ ,119.42
+ ,120.89
+ ,122.04
+ ,120.15
+ ,118.21
+ ,125.33
+ ,135.37
+ ,146.89
+ ,149.54
+ ,160.84
+ ,169.95
+ ,177.13
+ ,169.35
+ ,159.93
+ ,149.69
+ ,148.67
+ ,136.32
+ ,128.17
+ ,138.74
+ ,140.58
+ ,147.57
+ ,147.83
+ ,155.65
+ ,148.88
+ ,147.90
+ ,141.99
+ ,136.58
+ ,121.82
+ ,127.52
+ ,129.80
+ ,131.29
+ ,135.96
+ ,146.50
+ ,158.65
+ ,153.21
+ ,147.04
+ ,141.04
+ ,140.45
+ ,140.15
+ ,139.30
+ ,137.60
+ ,146.02
+ ,158.79
+ ,167.19
+ ,161.99
+ ,164.62
+ ,156.21
+ ,154.42
+ ,150.39
+ ,148.98
+ ,158.61
+ ,169.98
+ ,190.09
+ ,184.39
+ ,193.67
+ ,203.79
+ ,204.07
+ ,208.92
+ ,206.88
+ ,218.89
+ ,215.52
+ ,251.66
+ ,262.11
+ ,227.27
+ ,202.60
+ ,191.63
+ ,178.71
+ ,178.32
+ ,176.41
+ ,175.70
+ ,175.73
+ ,172.35
+ ,176.61
+ ,183.49
+ ,172.59
+ ,148.39
+ ,138.31
+ ,150.61
+ ,151.74
+ ,151.66
+ ,149.88
+ ,144.62
+ ,137.10
+ ,140.05
+ ,138.92
+ ,130.15
+ ,128.92
+ ,120.64
+ ,118.54
+ ,107.95
+ ,107.93
+ ,126.54
+ ,130.21
+ ,126.21
+ ,125.29
+ ,117.03
+ ,117.34
+ ,113.87
+ ,113.00
+ ,111.41
+ ,103.02
+ ,111.41
+ ,113.19
+ ,108.10
+ ,108.80
+ ,102.16
+ ,105.83
+ ,108.41
+ ,105.70
+ ,105.11
+ ,110.78
+ ,113.51
+ ,108.98
+ ,108.28
+ ,117.49
+ ,128.22
+ ,127.73
+ ,128.01
+ ,132.84
+ ,128.12
+ ,130.28
+ ,129.30
+ ,135.00
+ ,127.23
+ ,123.79
+ ,121.92
+ ,122.03
+ ,123.34
+ ,125.27
+ ,122.53
+ ,125.31
+ ,123.28
+ ,122.56
+ ,123.72
+ ,121.46
+ ,132.03
+ ,149.30
+ ,161.26
+ ,187.84
+ ,190.32
+ ,176.26
+ ,168.98
+ ,149.60
+ ,150.84
+ ,141.81
+ ,138.62
+ ,141.96
+ ,131.35
+ ,131.62
+ ,148.72
+ ,145.62
+ ,147.46
+ ,160.55
+ ,165.57
+ ,166.33
+ ,161.39
+ ,166.28
+ ,166.58
+ ,163.73
+ ,154.74
+ ,150.60
+ ,141.29
+ ,151.03
+ ,150.15
+ ,156.57
+ ,153.87
+ ,153.59
+ ,151.30
+ ,150.99
+ ,140.88
+ ,144.25
+ ,141.93
+ ,143.87
+ ,149.36
+ ,159.71
+ ,167.83
+ ,161.12
+ ,164.44
+ ,167.16
+ ,179.84
+ ,174.44
+ ,180.35
+ ,193.17
+ ,195.16
+ ,202.43
+ ,189.91
+ ,195.98
+ ,212.09
+ ,205.81
+ ,204.31
+ ,196.07
+ ,199.98
+ ,199.10
+ ,198.31
+ ,195.72
+ ,223.04
+ ,238.41
+ ,259.73
+ ,326.54
+ ,335.15
+ ,321.81
+ ,368.62
+ ,369.59
+ ,425.00
+ ,439.72
+ ,362.23
+ ,328.76
+ ,348.55
+ ,328.18
+ ,329.34
+ ,295.55
+ ,237.38
+ ,226.85
+ ,220.14
+ ,239.36
+ ,224.69
+ ,230.98
+ ,233.47
+ ,256.70
+ ,253.41
+ ,224.95
+ ,210.37
+ ,191.09
+ ,198.85
+ ,211.04
+ ,206.25
+ ,201.51
+ ,194.54
+ ,191.07
+ ,192.82
+ ,181.88
+ ,157.67
+ ,195.82
+ ,246.25
+ ,271.69
+ ,270.23
+ ,274.08
+ ,306.53
+ ,326.55
+ ,348.15
+ ,316.75
+ ,336.12
+ ,354.47
+ ,326.43
+ ,303.88
+ ,327.07
+ ,315.92
+ ,289.01
+ ,281.01
+ ,269.03
+ ,274.89
+ ,277.77
+ ,283.88
+ ,266.32
+ ,264.36
+ ,276.19
+ ,345.69
+ ,349.40
+ ,353.42
+ ,358.20
+ ,361.00)
> ylimmax = ''
> ylimmin = ''
> main = 'Robustness of Central Tendency'
> #'GNU S' R Code compiled by R2WASP v. 1.2.291 ()
> #Author: root
> #To cite this work: Wessa, P., (2012), Central Tendency (v1.0.4) 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
> #
> 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] 178.8793
> sqrtn <- sqrt(length(x))
> (armse <- sd(x) / sqrtn)
[1] 3.847429
> (armose <- arm / armse)
[1] 46.4932
> (geo <- geomean(x))
[1] 169.0039
> (har <- harmean(x))
[1] 161.197
> (qua <- quamean(x))
[1] 190.889
> (win <- winmean(x))
           [,1]      [,2]
  [1,] 178.8333 3.8364841
  [2,] 178.4790 3.7649779
  [3,] 178.4752 3.7629653
  [4,] 178.3920 3.7487615
  [5,] 178.4064 3.7432104
  [6,] 178.3510 3.7341673
  [7,] 178.2678 3.7201371
  [8,] 178.2447 3.7154406
  [9,] 178.1283 3.6964686
 [10,] 178.1131 3.6913019
 [11,] 178.1050 3.6886426
 [12,] 178.0787 3.6692053
 [13,] 177.6891 3.6060065
 [14,] 177.6477 3.5991841
 [15,] 177.3582 3.5574482
 [16,] 177.4013 3.5485365
 [17,] 177.3798 3.5431993
 [18,] 177.3236 3.5331864
 [19,] 177.2920 3.5284772
 [20,] 177.3126 3.5270987
 [21,] 177.3301 3.5245055
 [22,] 176.9953 3.4770793
 [23,] 176.6362 3.4221861
 [24,] 176.5995 3.4113764
 [25,] 175.8503 3.3049960
 [26,] 175.7795 3.2653381
 [27,] 175.0493 3.1680591
 [28,] 174.4698 3.0901276
 [29,] 173.9900 3.0288058
 [30,] 173.7757 2.9899552
 [31,] 173.4760 2.9482253
 [32,] 173.4016 2.9226903
 [33,] 173.3391 2.9011118
 [34,] 173.3030 2.8870771
 [35,] 173.0541 2.8531262
 [36,] 172.9477 2.8289112
 [37,] 172.8444 2.8093576
 [38,] 172.5149 2.7697553
 [39,] 172.2752 2.7403277
 [40,] 171.9775 2.7069060
 [41,] 171.7201 2.6669799
 [42,] 171.3015 2.6205020
 [43,] 170.9344 2.5633919
 [44,] 170.6873 2.5355129
 [45,] 169.9353 2.4468468
 [46,] 168.8931 2.3385082
 [47,] 168.8650 2.3159355
 [48,] 168.7597 2.2957736
 [49,] 168.1232 2.2330109
 [50,] 167.7660 2.1892816
 [51,] 167.1408 2.1292145
 [52,] 167.0717 2.1221665
 [53,] 166.7407 2.0907027
 [54,] 166.8519 2.0760093
 [55,] 166.6107 2.0442003
 [56,] 166.1140 1.9921796
 [57,] 165.9644 1.9651258
 [58,] 165.3709 1.9037706
 [59,] 164.7397 1.8427063
 [60,] 164.5862 1.8211072
 [61,] 164.4728 1.8080967
 [62,] 164.1844 1.7822733
 [63,] 163.7679 1.7461305
 [64,] 163.7147 1.7296526
 [65,] 163.6888 1.7172036
 [66,] 163.3709 1.6897330
 [67,] 163.3909 1.6804444
 [68,] 163.4406 1.6677248
 [69,] 163.2480 1.6404221
 [70,] 163.2225 1.6363242
 [71,] 163.0220 1.6179417
 [72,] 162.8976 1.5728073
 [73,] 162.6987 1.5551181
 [74,] 162.6840 1.5471894
 [75,] 162.5693 1.5353888
 [76,] 162.1073 1.4849075
 [77,] 162.2915 1.4692318
 [78,] 162.8097 1.4291911
 [79,] 162.8806 1.4208867
 [80,] 162.7876 1.4058399
 [81,] 162.7230 1.3864669
 [82,] 162.5841 1.3620793
 [83,] 162.5179 1.3470008
 [84,] 162.5653 1.3303068
 [85,] 162.3705 1.2959195
 [86,] 162.4190 1.2714412
 [87,] 162.5029 1.2653561
 [88,] 162.3187 1.2466876
 [89,] 162.3039 1.2382455
 [90,] 162.3637 1.2270835
 [91,] 161.9646 1.1663183
 [92,] 160.9407 1.0874620
 [93,] 160.7553 1.0616345
 [94,] 160.2931 1.0234113
 [95,] 159.9049 0.9836099
 [96,] 159.7933 0.9690758
 [97,] 159.5097 0.9389032
 [98,] 159.5520 0.9193336
 [99,] 159.5553 0.9143348
[100,] 159.2065 0.8894389
> (tri <- trimean(x))
           [,1]      [,2]
  [1,] 178.2635 3.7644140
  [2,] 177.6861 3.6884424
  [3,] 177.2816 3.6474801
  [4,] 176.8729 3.6053151
  [5,] 176.4800 3.5651688
  [6,] 176.0788 3.5243697
  [7,] 175.6816 3.4833924
  [8,] 175.2914 3.4428671
  [9,] 174.8987 3.4010825
 [10,] 174.5143 3.3599361
 [11,] 174.1261 3.3174043
 [12,] 173.7330 3.2730354
 [13,] 173.3367 3.2285975
 [14,] 172.9675 3.1887746
 [15,] 172.5962 3.1476060
 [16,] 172.2410 3.1084359
 [17,] 171.8774 3.0680224
 [18,] 171.5097 3.0259509
 [19,] 171.1401 2.9825045
 [20,] 170.7667 2.9370921
 [21,] 170.3863 2.8892485
 [22,] 169.9990 2.8388748
 [23,] 169.6236 2.7895210
 [24,] 169.2609 2.7417312
 [25,] 168.8942 2.6919787
 [26,] 168.5579 2.6475692
 [27,] 168.2194 2.6034538
 [28,] 167.9086 2.5639561
 [29,] 167.6184 2.5277523
 [30,] 167.3439 2.4937509
 [31,] 167.0739 2.4604158
 [32,] 166.8116 2.4279580
 [33,] 166.5479 2.3952584
 [34,] 166.2820 2.3620158
 [35,] 166.0129 2.3276847
 [36,] 165.7485 2.2935093
 [37,] 165.4833 2.2587959
 [38,] 165.2172 2.2231941
 [39,] 164.9580 2.1879675
 [40,] 164.7024 2.1524441
 [41,] 164.4525 2.1168537
 [42,] 164.2066 2.0815368
 [43,] 163.9701 2.0469670
 [44,] 163.7412 2.0137656
 [45,] 163.5160 1.9801901
 [46,] 163.3106 1.9500387
 [47,] 163.1341 1.9246451
 [48,] 162.9551 1.8989035
 [49,] 162.7758 1.8726755
 [50,] 162.6123 1.8485516
 [51,] 162.4564 1.8254341
 [52,] 162.3161 1.8043603
 [53,] 162.3161 1.7823091
 [54,] 162.0406 1.7606541
 [55,] 161.9002 1.7383035
 [56,] 161.7638 1.7162952
 [57,] 161.6387 1.6958538
 [58,] 161.5152 1.6755461
 [59,] 161.4059 1.6574058
 [60,] 161.3119 1.6414935
 [61,] 161.2202 1.6256550
 [62,] 161.1295 1.6094172
 [63,] 161.0447 1.5934918
 [64,] 160.9695 1.5784778
 [65,] 160.8940 1.5632600
 [66,] 160.8175 1.5475915
 [67,] 160.7477 1.5323424
 [68,] 160.6758 1.5164369
 [69,] 160.6007 1.5000011
 [70,] 160.5289 1.4838882
 [71,] 160.4561 1.4667423
 [72,] 160.3868 1.4493605
 [73,] 160.3191 1.4332525
 [74,] 160.2550 1.4169211
 [75,] 160.1895 1.3997081
 [76,] 160.1254 1.3818006
 [77,] 160.0720 1.3656002
 [78,] 160.0122 1.3489073
 [79,] 159.9367 1.3330386
 [80,] 159.8572 1.3162291
 [81,] 159.7778 1.2988890
 [82,] 159.6980 1.2812354
 [83,] 159.6980 1.2636029
 [84,] 159.5405 1.2453348
 [85,] 159.4577 1.2264189
 [86,] 159.3778 1.2081471
 [87,] 159.2939 1.1897478
 [88,] 159.2051 1.1698556
 [89,] 159.1185 1.1494029
 [90,] 159.0295 1.1274396
 [91,] 158.9358 1.1039402
 [92,] 158.8502 1.0830451
 [93,] 158.7907 1.0667068
 [94,] 158.7344 1.0506528
 [95,] 158.6895 1.0360049
 [96,] 158.6541 1.0229517
 [97,] 158.6207 1.0095366
 [98,] 158.5945 0.9970356
 [99,] 158.5659 0.9845782
[100,] 158.5361 0.9709504
> (midr <- midrange(x))
[1] 270.94
> midm <- array(NA,dim=8)
> for (j in 1:8) midm[j] <- midmean(x,j)
> midm
[1] 159.9354 160.1895 160.1895 160.1895 160.1895 160.0011 160.1895 160.1895
> postscript(file="/var/wessaorg/rcomp/tmp/1lwxh1354979293.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/wessaorg/rcomp/tmp/2oc6o1354979293.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/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/wessaorg/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/wessaorg/rcomp/tmp/3roij1354979293.tab") 
> 
> try(system("convert tmp/1lwxh1354979293.ps tmp/1lwxh1354979293.png",intern=TRUE))
character(0)
> try(system("convert tmp/2oc6o1354979293.ps tmp/2oc6o1354979293.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  2.678   0.317   3.167