R version 2.12.0 (2010-10-15)
Copyright (C) 2010 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-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(684.28
+ ,722.57
+ ,695.96
+ ,688.13
+ ,720.76
+ ,737.26
+ ,736.9
+ ,727.38
+ ,728.83
+ ,715.58
+ ,735.93
+ ,758.69
+ ,758.69
+ ,740.99
+ ,719.44
+ ,721.24
+ ,737.38
+ ,756.52
+ ,745.32
+ ,733.76
+ ,738.46
+ ,736.65
+ ,737.5
+ ,725.22
+ ,722.45
+ ,720.44
+ ,733.6
+ ,662.93
+ ,763.87
+ ,746.67
+ ,712.59
+ ,655.02
+ ,660.2
+ ,633.59
+ ,627.08
+ ,635.63
+ ,668.87
+ ,656.7
+ ,631.9
+ ,610.83
+ ,632.98
+ ,631.06
+ ,644.3
+ ,647.55
+ ,655.14
+ ,674.4
+ ,676.21
+ ,670.43
+ ,683.43
+ ,703.42
+ ,708.48
+ ,714.5
+ ,706.19
+ ,694.15
+ ,658.15
+ ,648.4
+ ,630.82
+ ,634.19
+ ,658.15
+ ,636.11
+ ,640.33
+ ,601.8
+ ,611.07
+ ,634.31
+ ,625.76
+ ,596.86
+ ,604.09
+ ,587.11
+ ,583.13
+ ,579.16
+ ,551.47
+ ,541.83
+ ,569.53
+ ,564.71
+ ,572.06
+ ,578.56
+ ,595.78
+ ,567.84
+ ,534.25
+ ,519.08
+ ,531.96
+ ,551.83
+ ,560.62
+ ,579.52
+ ,593.13
+ ,626.72
+ ,715.1
+ ,832.38
+ ,837.2
+ ,879.46
+ ,881.99
+ ,1044.9
+ ,924.73
+ ,987.95
+ ,966.15
+ ,1016.72
+ ,1107.27
+ ,1296.67
+ ,1379.75
+ ,1543.03
+ ,1570.12
+ ,1538.81
+ ,1484.63
+ ,1451.27
+ ,1414.91
+ ,1456.93
+ ,1319.19
+ ,1267.89
+ ,1349.77
+ ,1240.2
+ ,1189.51
+ ,1117.87
+ ,1080.06
+ ,1054.77
+ ,1019.47
+ ,1049.96
+ ,1060.79
+ ,1070.42
+ ,1075.24
+ ,1080.06
+ ,1074.04
+ ,1062
+ ,1064.4
+ ,1071.63
+ ,1057.18
+ ,898.24
+ ,895.83
+ ,951.22
+ ,936.77
+ ,901.85
+ ,888.61
+ ,870.55
+ ,887.41
+ ,596.02
+ ,586.39
+ ,596.02
+ ,721
+ ,777.51
+ ,723.48
+ ,680.64
+ ,613.45
+ ,558.18
+ ,641.49
+ ,652.19
+ ,619.37
+ ,655.83
+ ,667.93
+ ,667.7
+ ,663.07
+ ,633.89
+ ,595.28
+ ,568.94
+ ,572.72
+ ,535.26
+ ,508.04
+ ,512.94
+ ,495.22
+ ,469.37
+ ,469.37
+ ,429.69
+ ,468.13
+ ,470.06
+ ,464.93
+ ,450.74
+ ,423.51
+ ,454.84
+ ,497.77
+ ,465.45
+ ,542.31
+ ,606.03
+ ,609.58
+ ,645.79
+ ,719.63
+ ,779.41
+ ,773.5
+ ,806.82
+ ,876.1
+ ,824.64
+ ,881.7
+ ,878.41
+ ,904.18
+ ,892.34
+ ,887.13
+ ,867.85
+ ,839.28
+ ,826.06
+ ,751.11
+ ,789.25
+ ,732.98
+ ,622.07
+ ,600.95
+ ,590.53
+ ,584.39
+ ,525.31
+ ,573.83
+ ,597.67
+ ,743.54
+ ,701.36
+ ,671.43
+ ,751.65
+ ,738.33
+ ,681.61
+ ,616.97
+ ,632.94
+ ,677.73
+ ,730.96
+ ,719.66
+ ,764.21
+ ,805
+ ,829.35
+ ,826.26
+ ,765.93
+ ,801.91
+ ,769.16
+ ,739.89
+ ,688.07
+ ,636.11
+ ,631.72
+ ,625.92
+ ,627.82
+ ,606.13
+ ,595.3
+ ,583.14
+ ,500.19
+ ,462.89
+ ,417.47
+ ,472.27
+ ,474.81
+ ,489.07
+ ,493.14
+ ,626.64
+ ,680.43
+ ,620.3
+ ,676.74
+ ,690.03
+ ,631.04
+ ,623.26
+ ,619.83
+ ,631.74
+ ,648.77
+ ,724.21
+ ,727.09
+ ,767.31
+ ,801.42
+ ,817.72
+ ,764.33
+ ,746.93
+ ,717.29
+ ,695.9
+ ,688.38
+ ,663.53
+ ,688.39
+ ,716.14
+ ,733.28
+ ,688.23
+ ,760.63
+ ,716.77
+ ,683.81
+ ,630.79
+ ,617.91
+ ,524.19
+ ,441.98
+ ,466.09
+ ,501.44
+ ,599.55
+ ,621.99
+ ,607.57
+ ,614.56
+ ,619.09
+ ,603.14
+ ,569.12
+ ,575.72
+ ,642.41
+ ,748.18
+ ,768.39
+ ,763.16
+ ,800.63
+ ,778.49
+ ,733.7
+ ,740.17
+ ,678.33
+ ,697.09
+ ,678.54
+ ,670.87
+ ,674.52
+ ,664.39
+ ,646.52
+ ,661.25
+ ,729.36
+ ,721.52
+ ,708.75
+ ,706.12
+ ,676.93
+ ,708.15
+ ,717.64
+ ,714.35
+ ,703.6
+ ,718.96
+ ,736.03
+ ,717.96
+ ,730.77
+ ,734.66
+ ,728.58
+ ,729.18
+ ,717.13
+ ,684.6
+ ,692.06
+ ,668.84
+ ,647.4
+ ,622.57
+ ,593.69
+ ,583.24
+ ,639.17
+ ,709.73
+ ,698.75
+ ,697.7
+ ,732.25
+ ,714.49
+ ,704.85
+ ,717.56
+ ,704.91
+ ,690.74
+ ,790.34
+ ,790.69
+ ,893.66
+ ,875.62
+ ,914.02
+ ,940.37
+ ,989.08
+ ,987.73
+ ,936.34
+ ,914.11
+ ,943.08
+ ,1080.64
+ ,1102.71
+ ,1097.26
+ ,1157.37
+ ,1154.25
+ ,1136
+ ,1133.42
+ ,1062.44
+ ,1041.1
+ ,1047.97
+ ,941.71
+ ,896.79
+ ,734.91
+ ,646.17
+ ,666.39
+ ,625.45
+ ,586.26
+ ,617.15
+ ,686.86
+ ,761.34
+ ,741.73
+ ,763.62
+ ,822.57
+ ,867.58
+ ,944.85
+ ,953.95
+ ,970.08
+ ,1003.22
+ ,972.81
+ ,1003.86
+ ,991.83
+ ,919.13
+ ,919.52
+ ,916.17
+ ,936.18
+ ,960.65)
> par20 = ''
> par19 = ''
> par18 = ''
> par17 = ''
> par16 = ''
> par15 = ''
> par14 = ''
> par13 = ''
> par12 = ''
> par11 = ''
> par10 = ''
> par9 = ''
> par8 = ''
> par7 = ''
> par6 = ''
> par5 = ''
> par4 = ''
> par3 = ''
> par2 = ''
> par1 = ''
> 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] 750.0261
> sqrtn <- sqrt(length(x))
> (armse <- sd(x) / sqrtn)
[1] 10.47143
> (armose <- arm / armse)
[1] 71.62597
> (geo <- geomean(x))
[1] 727.0037
> (har <- harmean(x))
[1] 706.8421
> (qua <- quamean(x))
[1] 776.4603
> (win <- winmean(x))
           [,1]      [,2]
  [1,] 749.9691 10.454647
  [2,] 749.9797 10.447127
  [3,] 749.6391 10.351231
  [4,] 749.4338 10.286618
  [5,] 749.4127 10.267971
  [6,] 748.9523 10.150232
  [7,] 748.3240 10.030745
  [8,] 747.6853  9.920678
  [9,] 746.9551  9.798769
 [10,] 746.4001  9.698984
 [11,] 745.5791  9.566297
 [12,] 744.6786  9.435222
 [13,] 742.9171  9.187781
 [14,] 741.7815  9.024125
 [15,] 741.7580  8.999940
 [16,] 741.5850  8.853868
 [17,] 741.6536  8.825064
 [18,] 740.9965  8.727034
 [19,] 740.5820  8.653553
 [20,] 740.4660  8.615232
 [21,] 740.2270  8.574578
 [22,] 739.6296  8.434503
 [23,] 739.8989  8.407934
 [24,] 740.2982  8.378962
 [25,] 740.3179  8.318464
 [26,] 740.3122  8.303660
 [27,] 740.6225  8.250578
 [28,] 740.7044  8.228680
 [29,] 740.3107  8.172017
 [30,] 740.6855  8.118863
 [31,] 740.6889  8.112201
 [32,] 741.3783  8.049034
 [33,] 741.0876  8.012257
 [34,] 741.4507  7.951362
 [35,] 741.2259  7.888227
 [36,] 741.4308  7.842963
 [37,] 741.4368  7.791215
 [38,] 741.1587  7.743057
 [39,] 738.8917  7.505393
 [40,] 738.6380  7.472446
 [41,] 737.4902  7.312173
 [42,] 737.4925  7.300352
 [43,] 736.2946  7.162405
 [44,] 736.1920  7.116828
 [45,] 736.4005  7.082427
 [46,] 736.4479  7.075256
 [47,] 734.5934  6.891481
 [48,] 734.7079  6.829655
 [49,] 734.1873  6.780882
 [50,] 733.4556  6.711140
 [51,] 732.6885  6.616797
 [52,] 732.5673  6.565907
 [53,] 731.6711  6.482492
 [54,] 731.5174  6.452968
 [55,] 731.8230  6.403998
 [56,] 732.0142  6.362362
 [57,] 731.5446  6.307901
 [58,] 731.7269  6.287135
 [59,] 731.7046  6.284684
 [60,] 729.9208  6.117173
 [61,] 729.0992  6.040646
 [62,] 729.0337  6.035035
 [63,] 728.6717  5.983301
 [64,] 728.4549  5.944523
 [65,] 728.7702  5.923461
 [66,] 727.2607  5.760330
 [67,] 726.9919  5.715900
 [68,] 726.5736  5.646082
 [69,] 726.4801  5.613164
 [70,] 726.6660  5.576431
 [71,] 726.2677  5.540950
 [72,] 726.2911  5.503307
 [73,] 725.9509  5.419708
 [74,] 725.9609  5.385498
 [75,] 725.9528  5.378030
 [76,] 725.3843  5.263734
 [77,] 725.5554  5.245488
 [78,] 725.5914  5.178054
 [79,] 725.4051  5.157720
 [80,] 725.0691  5.108019
 [81,] 725.2227  5.084823
 [82,] 724.1583  4.991103
 [83,] 723.6544  4.936891
 [84,] 723.6999  4.925904
 [85,] 717.5703  4.398775
 [86,] 717.1041  4.361311
 [87,] 716.0856  4.269735
 [88,] 715.5276  4.206626
 [89,] 715.3105  4.120567
 [90,] 715.3373  4.112520
 [91,] 715.0266  4.084604
 [92,] 714.6900  4.036403
 [93,] 713.4878  3.946950
 [94,] 710.8028  3.744993
 [95,] 710.5248  3.700846
 [96,] 710.4936  3.598496
 [97,] 710.3726  3.589136
 [98,] 710.2212  3.571168
 [99,] 707.5598  3.390186
[100,] 707.6438  3.372912
[101,] 707.3509  3.352737
[102,] 704.6751  3.172583
[103,] 704.7086  3.138091
[104,] 704.4437  3.119777
[105,] 703.4762  3.036883
[106,] 702.3157  2.954079
[107,] 702.1794  2.934895
[108,] 701.8984  2.913350
[109,] 701.8807  2.863982
[110,] 701.5468  2.826093
[111,] 701.5107  2.823914
[112,] 702.3363  2.759559
[113,] 702.6150  2.732870
[114,] 702.8312  2.702097
[115,] 702.5508  2.650109
[116,] 702.9217  2.600130
[117,] 702.7790  2.534233
[118,] 702.9005  2.526813
[119,] 702.3136  2.478124
[120,] 701.0160  2.367450
[121,] 700.8882  2.354144
[122,] 700.2005  2.281051
[123,] 699.9071  2.249696
> (tri <- trimean(x))
           [,1]      [,2]
  [1,] 748.6977 10.247711
  [2,] 747.4123 10.031033
  [3,] 746.1074  9.807949
  [4,] 744.9041  9.610081
  [5,] 743.7402  9.421529
  [6,] 742.5675  9.228463
  [7,] 741.4614  9.050151
  [8,] 740.4366  8.884823
  [9,] 739.4841  8.729441
 [10,] 738.6064  8.584617
 [11,] 737.7776  8.446513
 [12,] 737.0190  8.318520
 [13,] 737.0190  8.199869
 [14,] 735.7843  8.102239
 [15,] 735.3180  8.016638
 [16,] 734.8479  7.929864
 [17,] 734.3841  7.852274
 [18,] 733.9102  7.773985
 [19,] 733.4713  7.700431
 [20,] 733.0516  7.629594
 [21,] 732.6332  7.558769
 [22,] 732.2227  7.488038
 [23,] 731.8380  7.424113
 [24,] 731.4352  7.359350
 [25,] 731.4352  7.293774
 [26,] 731.0080  7.229321
 [27,] 730.1358  7.163145
 [28,] 729.6779  7.097359
 [29,] 729.2106  7.030106
 [30,] 728.7536  6.963392
 [31,] 728.2755  6.896793
 [32,] 727.7910  6.827717
 [33,] 727.2740  6.758975
 [34,] 726.7608  6.689307
 [35,] 726.2276  6.619748
 [36,] 725.6952  6.550510
 [37,] 725.1484  6.480479
 [38,] 724.5940  6.409960
 [39,] 724.0413  6.338775
 [40,] 723.5551  6.277046
 [41,] 723.0703  6.214210
 [42,] 722.6149  6.156679
 [43,] 722.1530  6.097079
 [44,] 721.7212  6.041671
 [45,] 721.2862  5.985886
 [46,] 720.8388  5.929060
 [47,] 720.3835  5.869872
 [48,] 719.9748  5.816965
 [49,] 719.5569  5.764358
 [50,] 719.5569  5.711606
 [51,] 718.7518  5.659688
 [52,] 718.7518  5.609724
 [53,] 717.9883  5.559646
 [54,] 717.6233  5.511077
 [55,] 717.2567  5.461469
 [56,] 716.8764  5.411488
 [57,] 716.4853  5.360790
 [58,] 716.0999  5.310078
 [59,] 715.7038  5.257635
 [60,] 715.3019  5.202584
 [61,] 714.9379  5.152708
 [62,] 714.5883  5.103856
 [63,] 714.2345  5.052700
 [64,] 713.8836  5.001244
 [65,] 713.5321  4.948833
 [66,] 713.1671  4.894428
 [67,] 712.8318  4.844875
 [68,] 712.4971  4.794674
 [69,] 712.1664  4.744960
 [70,] 711.8321  4.693982
 [71,] 711.4877  4.641729
 [72,] 711.1463  4.588248
 [73,] 710.7982  4.533405
 [74,] 710.4516  4.479366
 [75,] 710.0985  4.423728
 [76,] 709.7390  4.365106
 [77,] 709.3857  4.308615
 [78,] 709.0219  4.249465
 [79,] 708.6504  4.189898
 [80,] 708.2760  4.127648
 [81,] 707.9018  4.064040
 [82,] 707.5169  3.997387
 [83,] 707.1480  3.931660
 [84,] 706.7829  3.864600
 [85,] 706.4094  3.793532
 [86,] 706.1635  3.747427
 [87,] 705.9228  3.700322
 [88,] 705.6994  3.655007
 [89,] 705.4837  3.610017
 [90,] 705.2681  3.566398
 [91,] 705.0473  3.520221
 [92,] 704.8286  3.472422
 [93,] 704.6125  3.423945
 [94,] 704.4179  3.377207
 [95,] 704.2779  3.338504
 [96,] 704.1408  3.299426
 [97,] 704.0013  3.262891
 [98,] 703.8612  3.224241
 [99,] 703.7211  3.183770
[100,] 703.7211  3.150564
[101,] 703.5479  3.115792
[102,] 703.4637  3.079552
[103,] 703.4368  3.050597
[104,] 703.4368  3.021244
[105,] 703.3854  2.990603
[106,] 703.3834  2.962118
[107,] 703.4074  2.935847
[108,] 703.4350  2.908465
[109,] 703.4698  2.879998
[110,] 703.5059  2.851963
[111,] 703.5506  2.823696
[112,] 703.5974  2.793105
[113,] 703.6264  2.763809
[114,] 703.6499  2.733563
[115,] 703.6499  2.702532
[116,] 703.6951  2.671957
[117,] 703.7133  2.641789
[118,] 703.7355  2.613056
[119,] 703.7355  2.582060
[120,] 703.7901  2.551308
[121,] 703.8572  2.524823
[122,] 703.9297  2.496507
[123,] 704.0214  2.470158
> (midr <- midrange(x))
[1] 993.795
> midm <- array(NA,dim=8)
> for (j in 1:8) midm[j] <- midmean(x,j)
> midm
[1] 704.1887 704.8286 704.8286 704.8286 704.8286 704.4044 704.8286 704.8286
> postscript(file="/var/www/rcomp/tmp/1c5xf1321397463.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))
Error: object 'ylimmin' not found
Execution halted