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(553.12
+ ,568.15
+ ,552.75
+ ,510.65
+ ,524.93
+ ,532.95
+ ,540.84
+ ,540.22
+ ,553.75
+ ,529.69
+ ,525.93
+ ,527.31
+ ,527.31
+ ,512.03
+ ,502.76
+ ,496.62
+ ,492.23
+ ,495.12
+ ,469.93
+ ,492.36
+ ,497.87
+ ,480.21
+ ,462.29
+ ,456.03
+ ,456.22
+ ,460.41
+ ,466.59
+ ,441.37
+ ,455.31
+ ,426.96
+ ,419.7
+ ,419.7
+ ,416.44
+ ,404.04
+ ,388.63
+ ,397.65
+ ,390.38
+ ,378.1
+ ,384.87
+ ,419.19
+ ,427.96
+ ,413.81
+ ,408.04
+ ,410.3
+ ,405.66
+ ,400.9
+ ,387
+ ,388.25
+ ,390
+ ,416.44
+ ,436.11
+ ,428.79
+ ,424.16
+ ,409.12
+ ,392.01
+ ,388.37
+ ,373.97
+ ,358.93
+ ,371.96
+ ,353.92
+ ,364.95
+ ,340.52
+ ,353.67
+ ,361.19
+ ,364.7
+ ,359.43
+ ,371.21
+ ,385.24
+ ,389.63
+ ,433.23
+ ,407.79
+ ,400.9
+ ,385.87
+ ,406.67
+ ,406.04
+ ,418.19
+ ,429.22
+ ,420.95
+ ,402.78
+ ,391.01
+ ,416.94
+ ,397.14
+ ,406.67
+ ,419.44
+ ,422.2
+ ,435.98
+ ,470.81
+ ,504.51
+ ,497.25
+ ,508.52
+ ,522.43
+ ,551.62
+ ,537.59
+ ,559.76
+ ,558.26
+ ,563.27
+ ,558.01
+ ,563.27
+ ,564.15
+ ,582.81
+ ,592.96
+ ,602.73
+ ,581.06
+ ,595.47
+ ,605.36
+ ,615.39
+ ,602.48
+ ,565.65
+ ,565.65
+ ,566.9
+ ,558.63
+ ,547.48
+ ,543.72
+ ,517.42
+ ,526.31
+ ,512.4
+ ,496.12
+ ,503.63
+ ,501.13
+ ,499.88
+ ,501.13
+ ,494.86
+ ,488.6
+ ,482.34
+ ,447.26
+ ,440.99
+ ,418.44
+ ,418.44
+ ,412.18
+ ,394.64
+ ,334.5
+ ,328.24
+ ,319.47
+ ,323.23
+ ,328.24
+ ,365.82
+ ,371.88
+ ,340.6
+ ,337.48
+ ,337.22
+ ,313.96
+ ,308.39
+ ,307.57
+ ,295.4
+ ,297.04
+ ,341.89
+ ,341.18
+ ,352.12
+ ,367.36
+ ,364.88
+ ,363.09
+ ,351.25
+ ,349.29
+ ,335.47
+ ,320.1
+ ,310.7
+ ,312.39
+ ,309.68
+ ,309.67
+ ,328.92
+ ,337.01
+ ,327.79
+ ,324.38
+ ,313.93
+ ,310.83
+ ,316.62
+ ,325.5
+ ,320.03
+ ,320.1
+ ,338.13
+ ,379.25
+ ,376.82
+ ,398.37
+ ,394.46
+ ,411.95
+ ,425.91
+ ,444.48
+ ,433.45
+ ,446.07
+ ,426.04
+ ,447.06
+ ,467.94
+ ,516.73
+ ,520.27
+ ,516.49
+ ,519.41
+ ,537.66
+ ,547.44
+ ,507.04
+ ,495.99
+ ,436.58
+ ,453.11
+ ,456.77
+ ,450.38
+ ,439.18
+ ,416.56
+ ,440.06
+ ,447.61
+ ,420.32
+ ,417.67
+ ,404.38
+ ,416.41
+ ,419.48
+ ,417.72
+ ,408.62
+ ,442.94
+ ,425.82
+ ,451.19
+ ,467.49
+ ,478.76
+ ,478.56
+ ,427.63
+ ,448.81
+ ,435.41
+ ,434.67
+ ,413.62
+ ,399.02
+ ,406.64
+ ,384.83
+ ,379.81
+ ,355.7
+ ,348.24
+ ,308.83
+ ,296.93
+ ,280.11
+ ,286.85
+ ,294.93
+ ,294.77
+ ,299.37
+ ,287.1
+ ,297.46
+ ,298.88
+ ,288.74
+ ,288.32
+ ,286.32
+ ,254.34
+ ,247.09
+ ,247.29
+ ,255.49
+ ,267.26
+ ,276.44
+ ,260.07
+ ,267.1
+ ,273.81
+ ,290.37
+ ,293.98
+ ,302.36
+ ,289.92
+ ,283.27
+ ,279.87
+ ,267.66
+ ,286.88
+ ,309.69
+ ,323.95
+ ,315.36
+ ,327.52
+ ,325.69
+ ,326.92
+ ,328.26
+ ,348.94
+ ,340.53
+ ,330.29
+ ,335.91
+ ,376.13
+ ,444.11
+ ,516.35
+ ,529.16
+ ,525.07
+ ,519.78
+ ,548.84
+ ,539.68
+ ,534.99
+ ,584.34
+ ,664.34
+ ,691.81
+ ,689.34
+ ,725.81
+ ,734.89
+ ,681.58
+ ,685.72
+ ,633.01
+ ,680.28
+ ,684.95
+ ,653.47
+ ,647.28
+ ,602.73
+ ,589.76
+ ,588.41
+ ,613.51
+ ,611.93
+ ,587.69
+ ,554.63
+ ,533.09
+ ,560.59
+ ,553.05
+ ,528.97
+ ,500.93
+ ,508.86
+ ,537.86
+ ,547.3
+ ,556.94
+ ,549.33
+ ,545.18
+ ,543.69
+ ,543.49
+ ,553.32
+ ,563.57
+ ,531.87
+ ,517.85
+ ,500.83
+ ,481.51
+ ,479.73
+ ,496.6
+ ,520.76
+ ,528.71
+ ,515.46
+ ,522
+ ,515.63
+ ,522.87
+ ,534.29
+ ,521.83
+ ,524.66
+ ,640.02
+ ,644.21
+ ,715.51
+ ,706.96
+ ,724.75
+ ,742.55
+ ,788.25
+ ,787.58
+ ,761.64
+ ,732.83
+ ,765.24
+ ,807.45
+ ,828.19
+ ,817.45
+ ,840.86
+ ,843.77
+ ,835.3
+ ,823.97
+ ,781.49
+ ,778.46
+ ,793.63
+ ,717.57
+ ,656.35
+ ,510.28
+ ,450.98
+ ,477.83
+ ,464.59
+ ,449.72
+ ,460.47
+ ,516.6
+ ,599.31
+ ,609.15
+ ,609.3
+ ,655.49
+ ,690.91
+ ,743.59
+ ,756.8
+ ,780.4
+ ,842.74
+ ,818.72
+ ,847.84
+ ,817.97
+ ,763.94
+ ,762.81
+ ,777.29
+ ,786.94
+ ,766.28)
> 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] 475.2487
> sqrtn <- sqrt(length(x))
> (armse <- sd(x) / sqrtn)
[1] 7.251475
> (armose <- arm / armse)
[1] 65.5382
> (geo <- geomean(x))
[1] 456.1261
> (har <- harmean(x))
[1] 438.2187
> (qua <- quamean(x))
[1] 495.189
> (win <- winmean(x))
           [,1]     [,2]
  [1,] 475.2382 7.249897
  [2,] 475.2708 7.245911
  [3,] 475.2649 7.243037
  [4,] 475.2543 7.230772
  [5,] 475.2532 7.210258
  [6,] 475.1872 7.200975
  [7,] 475.0951 7.187356
  [8,] 475.2122 7.174930
  [9,] 475.2637 7.168421
 [10,] 475.0856 7.126708
 [11,] 474.6808 7.074832
 [12,] 474.6086 7.045901
 [13,] 474.6925 7.035173
 [14,] 474.6883 7.030779
 [15,] 474.4680 7.004127
 [16,] 474.4302 6.997811
 [17,] 474.3971 6.983134
 [18,] 474.3605 6.974907
 [19,] 473.8543 6.904430
 [20,] 473.8224 6.896186
 [21,] 473.9538 6.872961
 [22,] 473.9336 6.861899
 [23,] 473.8706 6.852859
 [24,] 473.5864 6.814953
 [25,] 472.7950 6.707942
 [26,] 472.7295 6.699360
 [27,] 472.1998 6.636348
 [28,] 472.1512 6.611876
 [29,] 471.6380 6.550487
 [30,] 471.7949 6.524140
 [31,] 471.6294 6.430396
 [32,] 471.5219 6.406944
 [33,] 470.7966 6.325860
 [34,] 469.4780 6.181393
 [35,] 469.3936 6.172993
 [36,] 469.2414 6.158018
 [37,] 468.9797 6.115785
 [38,] 468.9138 6.107215
 [39,] 468.7225 6.061525
 [40,] 468.7485 6.036424
 [41,] 466.9808 5.870922
 [42,] 466.2307 5.777418
 [43,] 466.2773 5.758068
 [44,] 466.3763 5.712721
 [45,] 465.6897 5.641330
 [46,] 465.3157 5.607389
 [47,] 464.7820 5.561406
 [48,] 464.2773 5.455344
 [49,] 462.0331 5.257412
 [50,] 461.8366 5.233125
 [51,] 461.7731 5.204935
 [52,] 461.4292 5.174109
 [53,] 461.5843 5.159874
 [54,] 461.1175 5.110821
 [55,] 460.7657 5.078037
 [56,] 460.8340 5.073182
 [57,] 460.7954 5.070248
 [58,] 460.3003 5.032405
 [59,] 459.7918 4.979097
 [60,] 459.6065 4.933106
 [61,] 459.7734 4.845153
 [62,] 459.7096 4.817214
 [63,] 459.6618 4.803038
 [64,] 459.2715 4.747843
 [65,] 459.0390 4.726017
 [66,] 458.7725 4.700556
 [67,] 456.5464 4.530840
 [68,] 456.7565 4.484291
 [69,] 456.5247 4.468603
 [70,] 456.5379 4.467667
 [71,] 456.3609 4.440703
 [72,] 456.3863 4.423491
 [73,] 457.5832 4.332300
 [74,] 457.7236 4.322686
 [75,] 457.2500 4.281883
 [76,] 457.4827 4.243117
 [77,] 457.4285 4.215365
 [78,] 457.6779 4.188085
 [79,] 457.6779 4.180982
 [80,] 457.8318 4.139926
 [81,] 458.0338 4.059754
 [82,] 457.9493 4.039764
 [83,] 458.2485 4.007384
 [84,] 458.6355 3.976092
 [85,] 458.9902 3.950915
 [86,] 458.9622 3.943674
 [87,] 458.7123 3.925427
 [88,] 458.3737 3.876992
 [89,] 458.6269 3.845432
 [90,] 459.2343 3.764172
 [91,] 459.3896 3.753059
 [92,] 459.3747 3.749571
 [93,] 459.3469 3.683631
 [94,] 459.5253 3.625595
 [95,] 459.6952 3.614021
 [96,] 459.9762 3.590039
 [97,] 459.5818 3.527495
 [98,] 459.5659 3.508001
 [99,] 460.7679 3.416641
[100,] 460.2855 3.385150
[101,] 460.3320 3.375615
[102,] 460.4868 3.363898
[103,] 460.0765 3.299913
[104,] 460.2315 3.266602
[105,] 459.9242 3.243572
[106,] 459.9586 3.236632
[107,] 459.9354 3.200133
[108,] 459.4057 3.154862
[109,] 459.3614 3.138686
[110,] 459.4925 3.124061
[111,] 459.7151 3.101452
[112,] 460.0338 3.031886
[113,] 460.0890 3.028648
[114,] 460.5524 2.964897
[115,] 460.5929 2.948534
[116,] 460.5489 2.919160
[117,] 460.7106 2.904581
[118,] 461.2254 2.864840
[119,] 460.6482 2.830229
[120,] 461.1165 2.786751
[121,] 461.3886 2.754882
[122,] 461.4449 2.745168
[123,] 461.5149 2.699830
> (tri <- trimean(x))
           [,1]     [,2]
  [1,] 474.8551 7.192904
  [2,] 474.4679 7.133907
  [3,] 474.0598 7.074917
  [4,] 473.6491 7.014845
  [5,] 473.2367 6.955974
  [6,] 472.8198 6.899552
  [7,] 472.4097 6.842848
  [8,] 472.0087 6.786338
  [9,] 471.5877 6.729556
 [10,] 471.1559 6.671521
 [11,] 470.7380 6.616439
 [12,] 470.3546 6.565206
 [13,] 470.3546 6.515097
 [14,] 469.5804 6.464145
 [15,] 469.1833 6.411672
 [16,] 468.7975 6.359586
 [17,] 468.4097 6.306039
 [18,] 468.0195 6.251609
 [19,] 467.6267 6.195699
 [20,] 467.2591 6.142957
 [21,] 466.8888 6.088742
 [22,] 466.5068 6.033971
 [23,] 466.1212 5.977760
 [24,] 465.7339 5.919884
 [25,] 465.7339 5.862160
 [26,] 465.3554 5.808978
 [27,] 464.6611 5.754160
 [28,] 464.3320 5.700977
 [29,] 464.0006 5.647047
 [30,] 463.6861 5.594519
 [31,] 463.3613 5.541259
 [32,] 463.0386 5.491008
 [33,] 462.7157 5.439935
 [34,] 462.4155 5.391191
 [35,] 462.1592 5.348184
 [36,] 461.9024 5.303804
 [37,] 461.6474 5.258321
 [38,] 461.3978 5.213057
 [39,] 461.1470 5.166275
 [40,] 460.8990 5.119752
 [41,] 460.6467 5.072414
 [42,] 460.4467 5.031283
 [43,] 460.2671 4.992868
 [44,] 460.0836 4.953657
 [45,] 459.8944 4.914794
 [46,] 459.7229 4.877566
 [47,] 459.5597 4.840255
 [48,] 459.4095 4.803402
 [49,] 459.2714 4.769666
 [50,] 459.2714 4.743183
 [51,] 459.1211 4.716515
 [52,] 459.1211 4.689797
 [53,] 458.9844 4.663141
 [54,] 458.9151 4.635810
 [55,] 458.8570 4.609256
 [56,] 458.8072 4.582806
 [57,] 458.7548 4.555269
 [58,] 458.7026 4.526512
 [59,] 458.6621 4.497951
 [60,] 458.6337 4.470218
 [61,] 458.6095 4.443015
 [62,] 458.5807 4.418061
 [63,] 458.5531 4.392970
 [64,] 458.5261 4.367164
 [65,] 458.5082 4.342287
 [66,] 458.4954 4.316996
 [67,] 458.4889 4.291411
 [68,] 458.5348 4.271294
 [69,] 458.5765 4.251912
 [70,] 458.6245 4.232037
 [71,] 458.6729 4.211059
 [72,] 458.7263 4.189955
 [73,] 458.7801 4.168329
 [74,] 458.8075 4.149181
 [75,] 458.8321 4.129274
 [76,] 458.8680 4.109801
 [77,] 458.8993 4.090714
 [78,] 458.9324 4.071566
 [79,] 458.9605 4.052338
 [80,] 458.9892 4.032183
 [81,] 459.0150 4.012452
 [82,] 459.0368 3.994757
 [83,] 459.0609 3.976701
 [84,] 459.0789 3.958778
 [85,] 459.0886 3.940943
 [86,] 459.0908 3.922934
 [87,] 459.0936 3.903976
 [88,] 459.1020 3.884477
 [89,] 459.1180 3.865688
 [90,] 459.1288 3.846923
 [91,] 459.1265 3.830309
 [92,] 459.1207 3.812916
 [93,] 459.1151 3.794361
 [94,] 459.1101 3.777317
 [95,] 459.1009 3.761496
 [96,] 459.0879 3.744905
 [97,] 459.0684 3.728039
 [98,] 459.0571 3.712628
 [99,] 459.0459 3.696787
[100,] 459.0459 3.683690
[101,] 458.9797 3.670827
[102,] 458.9498 3.657171
[103,] 458.9156 3.642760
[104,] 458.9156 3.630018
[105,] 458.8599 3.617560
[106,] 458.8361 3.604929
[107,] 458.8108 3.591314
[108,] 458.7855 3.578091
[109,] 458.7715 3.565728
[110,] 458.7581 3.552804
[111,] 458.7413 3.539189
[112,] 458.7190 3.525232
[113,] 458.6887 3.513332
[114,] 458.6562 3.500151
[115,] 458.6562 3.488772
[116,] 458.5657 3.476798
[117,] 458.5190 3.464894
[118,] 458.4670 3.452213
[119,] 458.4670 3.440053
[120,] 458.3471 3.428299
[121,] 458.2801 3.417338
[122,] 458.2042 3.406555
[123,] 458.1246 3.394666
> (midr <- midrange(x))
[1] 547.465
> midm <- array(NA,dim=8)
> for (j in 1:8) midm[j] <- midmean(x,j)
> midm
[1] 458.6415 459.1207 459.1207 459.1207 459.1207 458.6517 459.1207 459.1207
> postscript(file="/var/www/rcomp/tmp/1vmcz1321397048.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