R version 3.0.2 (2013-09-25) -- "Frisbee Sailing"
Copyright (C) 2013 The R Foundation for Statistical Computing
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(15086
+ ,20559
+ ,39383
+ ,68032
+ ,42346
+ ,40184
+ ,29183
+ ,18049
+ ,24377
+ ,48033
+ ,61410
+ ,40637
+ ,37914
+ ,22481
+ ,17061
+ ,21921
+ ,34043
+ ,36152
+ ,38929
+ ,36713
+ ,27928
+ ,14178
+ ,19331
+ ,35030
+ ,37834
+ ,45257
+ ,39650
+ ,33001
+ ,22268
+ ,34283
+ ,63733
+ ,54548
+ ,31666
+ ,32654
+ ,28409
+ ,16234
+ ,19544
+ ,33375
+ ,30839
+ ,27581
+ ,28596
+ ,28542
+ ,14285
+ ,20799
+ ,28729
+ ,30278
+ ,31266
+ ,38341
+ ,24938
+ ,15700
+ ,17382
+ ,34229
+ ,29584
+ ,33909
+ ,29397
+ ,24644
+ ,14124
+ ,18770
+ ,32681
+ ,30518
+ ,34310
+ ,24217
+ ,15806
+ ,12362
+ ,13751
+ ,25312
+ ,23549
+ ,30305
+ ,30411
+ ,15593
+ ,17382
+ ,15806
+ ,29183
+ ,33348
+ ,33589
+ ,30198
+ ,23469
+ ,15593
+ ,16474
+ ,27341
+ ,26673
+ ,33268
+ ,26406
+ ,23656
+ ,15086
+ ,18316
+ ,27528
+ ,30305
+ ,34550
+ ,34710
+ ,23843
+ ,12362
+ ,21200
+ ,27261
+ ,31586
+ ,33322
+ ,32814
+ ,21974
+ ,17115
+ ,14445
+ ,26673
+ ,31746
+ ,28249
+ ,29103
+ ,21974
+ ,17622
+ ,21066
+ ,28382
+ ,32147
+ ,35885
+ ,40824
+ ,28622
+ ,24297
+ ,30465
+ ,47927
+ ,51665
+ ,57058
+ ,48407
+ ,30011
+ ,21413
+ ,30518
+ ,41171
+ ,37487
+ ,44456
+ ,41465
+ ,28409
+ ,34870
+ ,29557
+ ,44936
+ ,44028
+ ,45310
+ ,39302
+ ,28516
+ ,19865
+ ,30091
+ ,38101
+ ,44402
+ ,52359
+ ,71182
+ ,60796
+ ,43735
+ ,63947
+ ,89071
+ ,86321
+ ,90967
+ ,54441
+ ,36339
+ ,29290
+ ,65842
+ ,99404
+ ,83865
+ ,84452
+ ,76309
+ ,69554
+ ,44616
+ ,67337
+ ,94144
+ ,84425
+ ,76763
+ ,73559
+ ,60956
+ ,46271
+ ,74066
+ ,91634
+ ,63947
+ ,56017
+ ,53293
+ ,54895
+ ,48781
+ ,60449
+ ,99965
+ ,84746
+ ,81088
+ ,77777
+ ,60315
+ ,46218
+ ,62985
+ ,92596
+ ,78525
+ ,72811
+ ,62318
+ ,48140
+ ,45443
+ ,66910
+ ,90032
+ ,91661
+ ,103516
+ ,100285
+ ,72384
+ ,51985
+ ,72224
+ ,99351
+ ,108616
+ ,102128
+ ,97669
+ ,71716
+ ,51584
+ ,63413
+ ,88137
+ ,86535
+ ,91955
+ ,84853
+ ,54388
+ ,40931
+ ,58900
+ ,76576
+ ,80447
+ ,60315
+ ,52012
+ ,37487
+ ,19972
+ ,25338
+ ,23549
+ ,19785
+ ,29637
+ ,34416
+ ,26914
+ ,25071
+ ,24003
+ ,19358
+ ,28783
+ ,46992
+ ,49796
+ ,40904
+ ,33722
+ ,42640
+ ,56097
+ ,62184
+ ,63039
+ ,61730
+ ,44349
+ ,32307
+ ,45123
+ ,63466
+ ,71903
+ ,73211
+ ,67845
+ ,51237
+ ,31506
+ ,45043
+ ,74386
+ ,76602
+ ,72704
+ ,71770
+ ,52225
+ ,40958
+ ,55856
+ ,79700
+ ,81889
+ ,82049
+ ,66643
+ ,54655
+ ,36739
+ ,53026
+ ,78098
+ ,81248
+ ,84372
+ ,75614
+ ,63012
+ ,38795
+ ,55670
+ ,81168
+ ,85787
+ ,88591
+ ,75401
+ ,63413
+ ,42907
+ ,62291
+ ,86241
+ ,85200
+ ,84559
+ ,75801
+ ,72304
+ ,50944
+ ,70034
+ ,104691
+ ,88404
+ ,93397
+ ,88030
+ ,60582
+ ,50890
+ ,69847
+ ,90166
+ ,84158
+ ,82503
+ ,83811
+ ,66296
+ ,42640
+ ,55563
+ ,82450
+ ,84799
+ ,82370
+ ,83384
+ ,60075
+ ,41091
+ ,55589
+ ,79513
+ ,87656
+ ,74760
+ ,64214
+ ,47019
+ ,34229
+ ,45150
+ ,64187
+ ,73879
+ ,76175
+ ,83785
+ ,52786
+ ,43681
+ ,61864
+ ,74573
+ ,79913
+ ,78765
+ ,79112
+ ,62718
+ ,40557
+ ,57058
+ ,88804
+ ,89205
+ ,82904
+ ,78445
+ ,62825
+ ,57832
+ ,70675
+ ,102795
+ ,107094
+ ,105625
+ ,98390
+ ,74360
+ ,64641
+ ,82423
+ ,116759
+ ,101941
+ ,101140
+ ,98790
+ ,65949
+ ,49635
+ ,77323
+ ,109203
+ ,102902
+ ,93637
+ ,89365
+ ,72250
+ ,53267
+ ,73238
+ ,99725
+ ,69100
+ ,66323
+ ,80447
+ ,80954
+ ,48194
+ ,52546
+ ,91821
+ ,84372
+ ,90994
+ ,72410
+ ,62745
+ ,47366
+ ,63947
+ ,82690
+ ,89712
+ ,79085
+ ,80527
+ ,61303
+ ,40691
+ ,58340
+ ,79166
+ ,82850
+ ,82076
+ ,73292
+ ,60823
+ ,51157
+ ,61410
+ ,64294
+ ,67658
+ ,66323
+ ,66563
+ ,52145
+ ,36980
+ ,50703
+ ,68112
+ ,68379
+ ,65735
+ ,60716
+ ,50463
+ ,33749
+ ,39436
+ ,66990
+ ,70541
+ ,75027
+ ,69233
+ ,53213
+ ,31186
+ ,41145
+ ,67604
+ ,68699
+ ,71690
+ ,60022
+ ,51718
+ ,32814
+ ,42400
+ ,71316
+ ,70141
+ ,64694
+ ,60636
+ ,37620
+ ,35431
+ ,41145
+ ,69393
+ ,72891
+ ,62638
+ ,63439
+ ,54014
+ ,37567
+ ,38875
+ ,64854
+ ,75214
+ ,72944
+ ,71796
+ ,60689
+ ,53000
+ ,50757
+ ,71289
+ ,75080
+ ,77377
+ ,68779
+ ,49155
+ ,30678
+ ,49342
+ ,67604
+ ,63733
+ ,71209
+ ,59461
+ ,49609
+ ,30839
+ ,39302
+ ,67685
+ ,67925
+ ,66056
+ ,60743
+ ,49422
+ ,30972
+ ,34924
+ ,56818
+ ,62024
+ ,59274
+ ,59007
+ ,43494
+ ,26486
+ ,41252
+ ,63146
+ ,70114
+ ,57779
+ ,61944
+ ,48728
+ ,32708
+ ,40237
+ ,67845
+ ,68913
+ ,66697
+ ,68272
+ ,51317
+ ,30812
+ ,37300
+ ,63599
+ ,72918
+ ,74386
+ ,73292
+ ,56577
+ ,41652
+ ,58900
+ ,86188
+ ,88350
+ ,85787
+ ,86535
+ ,67978
+ ,42720
+ ,64187
+ ,98897
+ ,108429
+ ,92249
+ ,111820
+ ,77163
+ ,59648
+ ,83651
+ ,103569
+ ,108749
+ ,96948
+ ,90113
+ ,75614
+ ,52733
+ ,78178
+ ,103302
+ ,96040
+ ,98843
+ ,93183
+ ,75267
+ ,54815
+ ,93904
+ ,110725
+ ,107735
+ ,100793
+ ,102795
+ ,73772
+ ,63092
+ ,82476
+ ,162042
+ ,108215
+ ,100846
+ ,104183
+ ,82156
+ ,59728
+ ,97028
+ ,148746
+ ,101380
+ ,100953
+ ,100205
+ ,73292
+ ,57138
+ ,63039
+ ,135236
+ ,104370
+ ,109844
+ ,104290
+ ,78231
+ ,58153
+ ,88938
+ ,132459
+ ,107654
+ ,113021
+ ,117587
+ ,85387
+ ,64347
+ ,88778
+ ,130777
+ ,159185
+ ,139988
+ ,136731
+ ,90273
+ ,61277
+ ,81275
+ ,82503
+ ,74333
+ ,77350
+ ,102662
+ ,76923
+ ,62291
+ ,97455
+ ,127733
+ ,101861
+ ,87763
+ ,80954
+ ,72918
+ ,67952
+ ,96574
+ ,148105
+ ,151576
+ ,132592
+ ,117213
+ ,102982
+ ,92622
+ ,114196
+ ,145515
+ ,157664
+ ,129522
+ ,121886
+ ,99271
+ ,76923
+ ,103649
+ ,111232
+ ,123621
+ ,104077
+ ,112354
+ ,69260
+ ,47286
+ ,66910
+ ,78738
+ ,62905
+ ,44082
+ ,25392
+ ,31666
+ ,30732
+ ,35751
+ ,47579
+ ,44002
+ ,30064
+ ,28809
+ ,44375
+ ,39249
+ ,48007
+ ,73452
+ ,91154
+ ,71342
+ ,63226
+ ,58820
+ ,50623
+ ,63573
+ ,75347
+ ,80714
+ ,88644
+ ,82370
+ ,63439
+ ,44776
+ ,62024
+ ,86081
+ ,95212
+ ,88244
+ ,80287
+ ,61330
+ ,47206
+ ,64400
+ ,89579
+ ,97722
+ ,109203
+ ,120551
+ ,80234
+ ,52679
+ ,76923
+ ,99164
+ ,100499
+ ,101620
+ ,94385
+ ,75241
+ ,58420
+ ,72677
+ ,76148
+ ,70568
+ ,76469
+ ,71449
+ ,49368
+ ,33135
+ ,62131
+ ,77003
+ ,81836
+ ,75107
+ ,112968
+ ,84399
+ ,65015
+ ,83758
+ ,116225
+ ,104984
+ ,115451
+ ,115718
+ ,82930
+ ,61170
+ ,88377
+ ,118254
+ ,124208
+ ,124502
+ ,109871
+ ,81782
+ ,59034
+ ,79593
+ ,117667
+ ,107788
+ ,108135
+ ,104958
+ ,78685
+ ,57378
+ ,76496
+ ,104370
+ ,102875
+ ,99591
+ ,97802
+ ,79299
+ ,57138
+ ,78738
+ ,123007
+ ,117293
+ ,147731
+ ,156409
+ ,85173
+ ,66430
+ ,85387
+ ,123728
+ ,114703
+ ,113769
+ ,96307
+ ,79433
+ ,58607
+ ,75828
+ ,108429
+ ,104157
+ ,99538
+ ,90353
+ ,69874
+ ,53694
+ ,66910
+ ,102261
+ ,104584
+ ,109550
+ ,99858
+ ,83598
+ ,64748
+ ,86642
+ ,84719
+ ,79566
+ ,93236
+ ,81649
+ ,61891
+ ,44723
+ ,62291
+ ,78365
+ ,127493
+ ,123568
+ ,114356
+ ,88911
+ ,63092
+ ,85574
+ ,115104
+ ,115638
+ ,107975
+ ,117533
+ ,80474
+ ,58259
+ ,78231
+ ,70408
+ ,73986
+ ,78471
+ ,83037
+ ,52813
+ ,34123
+ ,45791
+ ,64454
+ ,61196
+ ,63039
+ ,59781
+ ,50997
+ ,31346
+ ,37861
+ ,53614
+ ,66857
+ ,73478
+ ,82183
+ ,51531
+ ,35564
+ ,42346
+ ,75054
+ ,81088
+ ,76255
+ ,71503
+ ,57992
+ ,46298
+ ,58313
+ ,92943
+ ,97055
+ ,93664
+ ,79993
+ ,67364
+ ,53106
+ ,63813
+ ,87843
+ ,95239
+ ,91661
+ ,85947
+ ,70408
+ ,47152
+ ,56043
+ ,79753
+ ,82530
+ ,76148
+ ,71583
+ ,59060
+ ,38582
+ ,54521
+ ,81702
+ ,81755
+ ,77136
+ ,64748
+ ,51932
+ ,34683
+ ,39730
+ ,68325
+ ,71049
+ ,68779
+ ,72758
+ ,51878
+ ,35484
+ ,47499
+ ,68405
+ ,77350
+ ,59461
+ ,67578
+ ,56524
+ ,35030
+ ,41839
+ ,68112
+ ,61570
+ ,76335
+ ,68432
+ ,46832
+ ,36766
+ ,42320
+ ,66563
+ ,67498
+ ,71903
+ ,63813
+ ,53000
+ ,32467
+ ,41065
+ ,64187
+ ,69740
+ ,68539
+ ,59327
+ ,49876
+ ,35111
+ ,48487
+ ,71503
+ ,68726
+ ,66803
+ ,65735
+ ,54068
+ ,30732
+ ,47606
+ ,72117
+ ,80207
+ ,77136
+ ,67231
+ ,54468
+ ,39169
+ ,53000
+ ,78551
+ ,76148
+ ,78738
+ ,74119
+ ,54602
+ ,41679
+ ,52279
+ ,79379
+ ,90219
+ ,83571
+ ,86081
+ ,58019
+ ,41465
+ ,56871
+ ,83090
+ ,96280
+ ,94091
+ ,97161
+ ,70141
+ ,50303
+ ,66002
+ ,84746
+ ,88404
+ ,119936
+ ,97669
+ ,78712
+ ,59114
+ ,85200
+ ,101380
+ ,103382
+ ,111713
+ ,106373
+ ,75134
+ ,63012
+ ,77056
+ ,100712
+ ,117800
+ ,118254
+ ,102448
+ ,79272
+ ,68592
+ ,93877
+ ,133206
+ ,147571
+ ,144687
+ ,111846
+ ,87069
+ ,65976
+ ,95880
+ ,150828
+ ,128614
+ ,122526
+ ,105946
+ ,82984
+ ,66216
+ ,98897
+ ,128080
+ ,140122
+ ,125223
+ ,122473
+ ,85467
+ ,66350
+ ,86722
+ ,106266
+ ,123381
+ ,128294
+ ,117080
+ ,84292
+ ,70007
+ ,101647
+ ,151336
+ ,149200
+ ,143646
+ ,106319
+ ,78284
+ ,58553
+ ,88564
+ ,114303
+ ,120337
+ ,123354
+ ,117133
+ ,101674
+ ,60876
+ ,87816
+ ,118495
+ ,114169
+ ,117720
+ ,120871
+ ,84826
+ ,64774
+ ,91955
+ ,115184
+ ,103783
+ ,96948
+ ,89979
+ ,84933
+ ,72143
+ ,87896
+ ,122286
+ ,131631
+ ,132806
+ ,139508
+ ,102074
+ ,74466
+ ,104717
+ ,132272
+ ,130323
+ ,138066
+ ,119990
+ ,102261
+ ,78605
+ ,104584
+ ,127599
+ ,111579
+ ,123381
+ ,103890
+ ,61811
+ ,47366
+ ,53507
+ ,75614
+ ,60262
+ ,47980
+ ,28703
+ ,20986
+ ,29290
+ ,37594
+ ,48888
+ ,50757
+ ,56417
+ ,30465
+ ,30411
+ ,38822
+ ,50383
+ ,73051
+ ,68325
+ ,84078
+ ,93477
+ ,72090
+ ,56177
+ ,77243
+ ,102181
+ ,93877
+ ,90433
+ ,88617
+ ,65655
+ ,67177
+ ,62665
+ ,95212
+ ,103543
+ ,109871
+ ,105305
+ ,78738
+ ,56898
+ ,109497
+ ,168344
+ ,123915
+ ,120711
+ ,113128
+ ,83838
+ ,67151
+ ,85574
+ ,116999
+ ,136117
+ ,127065
+ ,111473
+ ,80100
+ ,74093
+ ,102021
+ ,146690
+ ,150508
+ ,135289
+ ,116012
+ ,78498
+ ,59701
+ ,86294
+ ,107468
+ ,115371
+ ,118388
+ ,110485
+ ,100659
+ ,79139
+ ,109310
+ ,139935
+ ,126051
+ ,141270
+ ,118949
+ ,97455
+ ,77777
+ ,107040
+ ,149760
+ ,147731
+ ,159132
+ ,124369
+ ,98016
+ ,73185
+ ,102528
+ ,242516
+ ,293460
+ ,320987
+ ,268736
+ ,239873
+ ,189143
+ ,248550
+ ,269670
+ ,278801
+ ,272634
+ ,181240
+ ,90566
+ ,67017
+ ,92809
+ ,134755
+ ,137291
+ ,138199
+ ,125704
+ ,107414
+ ,75588
+ ,100606
+ ,120791
+ ,124075
+ ,101460
+ ,105198
+ ,70007
+ ,58446
+ ,70808
+ ,101220
+ ,127599
+ ,120497
+ ,110752
+ ,90113
+ ,70541
+ ,84559
+ ,120524
+ ,118922
+ ,123648
+ ,113448
+ ,84132
+ ,65228
+ ,101460
+ ,125303
+ ,141029
+ ,130136
+ ,124502
+ ,97241
+ ,71556
+ ,97188
+ ,141430
+ ,134621
+ ,134354
+ ,113261
+ ,74520
+ ,70835
+ ,87736
+ ,128908
+ ,118575
+ ,118014
+ ,111232
+ ,79593
+ ,56043
+ ,81221
+ ,115024
+ ,91661
+ ,110298
+ ,99591
+ ,70488
+ ,50757
+ ,69473
+ ,117186
+ ,82637
+ ,91207
+ ,76682
+ ,66964
+ ,44536
+ ,64267
+ ,99458
+ ,89338
+ ,104397
+ ,90486
+ ,70702
+ ,58633
+ ,68672
+ ,93477
+ ,110485
+ ,155875
+ ,86828
+ ,74653
+ ,57512
+ ,88484
+ ,78979
+ ,84158
+ ,76389
+ ,77136
+ ,55189
+ ,36686
+ ,50703
+ ,75428
+ ,78979
+ ,78151
+ ,69260
+ ,56844
+ ,39436
+ ,51985
+ ,73585
+ ,67177
+ ,65442
+ ,61757
+ ,50890
+ ,29050
+ ,48647
+ ,70488
+ ,63199
+ ,62825
+ ,66803
+ ,49555
+ ,30892
+ ,40103
+ ,59434
+ ,67604
+ ,70194
+ ,63760
+ ,53106
+ ,34310
+ ,44883
+ ,68058
+ ,74333
+ ,71156
+ ,67337
+ ,52279
+ ,37727
+ ,43254
+ ,74493
+ ,77991
+ ,68699
+ ,66483
+ ,62318
+ ,35885
+ ,49902
+ ,78231
+ ,67044
+ ,71476
+ ,68112
+ ,52492
+ ,41492
+ ,46805
+ ,69847
+ ,73478
+ ,78632
+ ,74813
+ ,60182
+ ,41305
+ ,39116
+ ,72010
+ ,79246
+ ,71690
+ ,69420
+ ,57939
+ ,46298
+ ,60422
+ ,81035
+ ,82637
+ ,79726
+ ,86348
+ ,66296
+ ,33695
+ ,52199
+ ,72517
+ ,72277
+ ,73078
+ ,76763
+ ,65655
+ ,34657
+ ,56417
+ ,68699
+ ,77163
+ ,79700
+ ,78151
+ ,59781
+ ,44829
+ ,46992
+ ,77030
+ ,92596
+ ,84933
+ ,83918
+ ,63306
+ ,47286
+ ,77403
+ ,102768
+ ,118361
+ ,107708
+ ,94705
+ ,76122
+ ,45417
+ ,79806
+ ,106346
+ ,95826
+ ,93797
+ ,99351
+ ,83838
+ ,54174
+ ,84559
+ ,116492
+ ,114543
+ ,123407
+ ,114677
+ ,110378
+ ,96040
+ ,141857
+ ,233465
+ ,146022
+ ,142124
+ ,118067
+ ,83010
+ ,75695
+ ,104157
+ ,150161
+ ,174832
+ ,139935
+ ,108749
+ ,109977
+ ,78872
+ ,126451
+ ,161295)
> 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] 75160.31
> sqrtn <- sqrt(length(x))
> (armse <- sd(x) / sqrtn)
[1] 1007.687
> (armose <- arm / armse)
[1] 74.58696
> (geo <- geomean(x))
[1] 67117.73
> (har <- harmean(x))
[1] 58493.56
> (qua <- quamean(x))
[1] 83121.74
> (win <- winmean(x))
           [,1]      [,2]
  [1,] 75138.14 1003.5654
  [2,] 75116.77  999.4462
  [3,] 75102.78  996.9894
  [4,] 75093.41  995.4671
  [5,] 75090.08  994.8548
  [6,] 74993.33  980.2071
  [7,] 74962.94  975.2475
  [8,] 74945.91  972.9064
  [9,] 74903.15  966.4833
 [10,] 74546.29  924.9545
 [11,] 74477.25  918.1318
 [12,] 74416.36  912.4325
 [13,] 74348.45  906.5840
 [14,] 74282.23  900.5761
 [15,] 74276.11  899.7200
 [16,] 74256.49  897.2342
 [17,] 74256.50  897.1409
 [18,] 74239.10  895.3329
 [19,] 74219.90  893.9009
 [20,] 74215.17  893.0679
 [21,] 74149.70  887.4664
 [22,] 74150.18  886.9275
 [23,] 74149.18  885.8446
 [24,] 74153.83  884.8703
 [25,] 74147.39  884.3584
 [26,] 74142.89  883.5841
 [27,] 74135.96  882.4867
 [28,] 74127.52  881.6975
 [29,] 74115.06  880.5566
 [30,] 74120.20  879.2468
 [31,] 74126.19  878.9534
 [32,] 74126.89  878.4402
 [33,] 74105.60  876.7673
 [34,] 74090.99  875.3732
 [35,] 74082.70  874.1380
 [36,] 74073.43  871.8523
 [37,] 74043.99  869.7642
 [38,] 73997.43  866.7882
 [39,] 73998.27  865.8125
 [40,] 73991.38  864.6163
 [41,] 74018.71  862.7316
 [42,] 74013.27  862.0927
 [43,] 73981.87  860.1363
 [44,] 73980.91  859.6634
 [45,] 73985.77  859.2255
 [46,] 73991.69  858.9472
 [47,] 73983.63  857.5708
 [48,] 73936.13  854.3355
 [49,] 73934.04  853.8700
 [50,] 73913.59  851.4883
 [51,] 73902.67  849.5518
 [52,] 73882.53  847.7670
 [53,] 73857.48  845.2107
 [54,] 73856.31  845.0234
 [55,] 73837.40  843.6698
 [56,] 73877.08  841.2209
 [57,] 73868.50  840.3422
 [58,] 73823.62  836.8582
 [59,] 73804.62  835.7745
 [60,] 73805.92  834.6590
 [61,] 73816.43  833.5192
 [62,] 73811.09  832.8109
 [63,] 73788.06  830.5529
 [64,] 73746.79  827.9742
 [65,] 73741.19  825.8450
 [66,] 73748.31  824.5364
 [67,] 73722.36  822.3994
 [68,] 73690.22  820.5033
 [69,] 73673.89  819.6193
 [70,] 73661.88  818.3789
 [71,] 73651.13  817.6557
 [72,] 73634.15  816.4400
 [73,] 73627.80  815.9519
 [74,] 73632.63  815.7375
 [75,] 73627.80  815.3269
 [76,] 73604.91  813.7899
 [77,] 73568.46  811.7116
 [78,] 73558.47  809.7246
 [79,] 73539.77  808.4250
 [80,] 73519.09  806.8574
 [81,] 73513.88  806.5878
 [82,] 73473.34  803.8278
 [83,] 73473.34  803.8278
 [84,] 73471.58  803.0477
 [85,] 73471.51  802.0018
 [86,] 73464.17  801.4503
 [87,] 73456.68  800.7168
 [88,] 73469.93  798.8795
 [89,] 73467.99  798.4226
 [90,] 73467.99  798.2378
 [91,] 73471.95  797.6983
 [92,] 73465.95  796.8365
 [93,] 73466.02  796.6499
 [94,] 73466.02  796.6499
 [95,] 73472.07  796.1919
 [96,] 73445.24  794.8411
 [97,] 73411.90  792.7756
 [98,] 73407.71  792.5669
 [99,] 73397.03  791.6398
[100,] 73364.83  790.0416
[101,] 73295.30  785.4159
[102,] 73293.16  784.9015
[103,] 73286.53  784.5782
[104,] 73279.83  783.6336
[105,] 73279.83  783.4230
[106,] 73277.52  783.3111
[107,] 73268.31  782.4429
[108,] 73245.09  780.6801
[109,] 73259.13  779.6333
[110,] 73178.80  775.1366
[111,] 73183.54  774.7104
[112,] 73166.67  772.5977
[113,] 73166.67  771.9371
[114,] 73164.20  771.1541
[115,] 73161.70  771.0360
[116,] 73159.17  770.2403
[117,] 73196.95  768.6092
[118,] 73194.38  767.1173
[119,] 73204.64  766.2204
[120,] 73202.03  764.4741
[121,] 73196.86  763.9954
[122,] 73194.31  763.6377
[123,] 73196.89  762.8179
[124,] 73191.49  762.5652
[125,] 73186.16  760.6329
[126,] 73191.64  759.6756
[127,] 73202.48  758.9695
[128,] 73202.58  758.4786
[129,] 73199.78  758.1074
[130,] 73194.12  757.5925
[131,] 73191.38  755.4598
[132,] 73174.27  753.6676
[133,] 73148.57  752.2238
[134,] 73128.50  751.0418
[135,] 73113.94  748.8399
[136,] 73119.85  747.8200
[137,] 73108.05  746.5046
[138,] 73110.94  745.6032
[139,] 73090.01  744.6487
[140,] 73087.08  743.9822
[141,] 73081.06  743.4403
[142,] 73044.36  741.7756
[143,] 73053.57  741.1270
[144,] 73053.57  739.7723
[145,] 73044.23  738.2631
[146,] 73041.06  737.8542
[147,] 73031.59  737.1501
[148,] 73047.44  736.2076
[149,] 73005.93  733.7828
[150,] 72979.96  731.5157
[151,] 72957.23  730.5036
[152,] 72950.86  729.3660
[153,] 72977.10  727.1369
[154,] 72977.10  726.5720
[155,] 72910.46  722.7724
[156,] 72870.14  718.9984
[157,] 72883.79  718.1509
[158,] 72870.18  717.5566
[159,] 72887.21  715.3917
[160,] 72897.64  713.8013
[161,] 72911.38  710.5901
[162,] 72914.90  710.4455
[163,] 72855.32  707.5770
[164,] 72855.32  707.2766
[165,] 72851.87  704.7352
[166,] 72894.64  702.9856
[167,] 72905.39  701.3414
[168,] 72894.57  700.8757
[169,] 72861.78  698.5584
[170,] 72850.97  697.7872
[171,] 72854.55  697.6415
[172,] 72865.63  696.8791
[173,] 72839.58  694.5260
[174,] 72835.94  694.0562
[175,] 72817.05  692.6329
[176,] 72828.39  690.9148
[177,] 72862.59  689.5322
[178,] 72832.07  685.3847
[179,] 72862.77  684.1507
[180,] 72847.40  683.1797
[181,] 72827.88  681.7215
[182,] 72835.79  681.4041
[183,] 72831.81  678.9748
[184,] 72827.81  678.1632
[185,] 72815.89  676.6880
[186,] 72795.83  675.1997
[187,] 72787.85  674.8666
[188,] 72796.02  674.2061
[189,] 72795.87  673.5414
[190,] 72767.41  672.3564
[191,] 72792.02  670.6985
[192,] 72754.92  668.1528
[193,] 72804.49  665.5033
[194,] 72712.96  660.7038
[195,] 72717.04  660.1998
[196,] 72763.28  658.0280
[197,] 72767.56  657.1825
[198,] 72725.15  654.7527
[199,] 72695.03  651.8145
[200,] 72656.38  649.2139
[201,] 72643.44  648.3405
[202,] 72613.02  646.7614
[203,] 72626.26  645.9008
[204,] 72590.95  644.1388
[205,] 72595.57  643.6152
[206,] 72577.82  642.9027
[207,] 72582.16  642.7320
[208,] 72564.40  640.9430
[209,] 72568.78  640.4102
[210,] 72595.83  639.3474
[211,] 72582.24  638.8034
[212,] 72568.59  637.8926
[213,] 72591.57  636.6386
[214,] 72596.22  636.4565
[215,] 72610.07  634.8202
[216,] 72661.20  630.2652
[217,] 72647.05  629.3432
[218,] 72623.52  628.4065
[219,] 72618.94  627.4759
[220,] 72656.85  625.6462
[221,] 72652.04  625.4551
[222,] 72642.39  623.9507
[223,] 72661.60  622.0848
[224,] 72666.47  617.3904
[225,] 72695.46  615.1509
[226,] 72724.57  613.6568
[227,] 72719.82  612.7018
[228,] 72768.83  610.8376
[229,] 72768.65  610.4585
[230,] 72759.02  609.3024
[231,] 72783.76  606.4373
[232,] 72773.67  605.6635
[233,] 72743.65  604.0889
[234,] 72753.83  603.7045
[235,] 72753.83  602.5371
[236,] 72758.96  601.5675
[237,] 72769.07  600.3930
[238,] 72769.07  599.6115
[239,] 72763.88  598.6230
[240,] 72758.85  597.6238
[241,] 72732.85  595.8152
[242,] 72748.44  594.8261
[243,] 72758.81  594.0322
[244,] 72732.68  592.6041
[245,] 72753.79  591.8131
[246,] 72748.44  590.8011
[247,] 72769.72  590.0055
[248,] 72742.96  588.5645
[249,] 72796.69  585.3382
[250,] 72845.00  580.6841
[251,] 72834.09  579.4455
[252,] 72828.81  578.8244
[253,] 72812.31  578.1824
[254,] 72905.16  573.9450
[255,] 72899.82  573.3180
[256,] 72910.74  571.2524
[257,] 72866.46  569.5319
[258,] 72855.45  568.6822
[259,] 72833.14  565.7304
[260,] 72822.04  564.4519
[261,] 72810.91  562.7594
[262,] 72782.64  561.6711
[263,] 72799.58  561.0519
[264,] 72788.31  560.6185
[265,] 72799.62  558.9264
[266,] 72805.19  557.8573
[267,] 72799.60  557.2082
[268,] 72868.87  554.6915
[269,] 72863.02  553.6094
[270,] 72845.63  552.5033
[271,] 72793.04  550.0657
[272,] 72816.48  549.2165
[273,] 72816.48  548.3328
[274,] 72851.77  546.1880
[275,] 72780.92  542.1687
[276,] 72733.36  537.7290
[277,] 72703.70  535.2711
[278,] 72697.66  534.1668
[279,] 72709.79  532.8508
[280,] 72769.98  530.6874
[281,] 72763.87  527.3432
[282,] 72769.78  527.1319
[283,] 72733.32  524.8543
[284,] 72751.61  523.3109
[285,] 72757.81  522.6354
[286,] 72739.39  521.5039
[287,] 72770.35  519.9430
[288,] 72770.35  518.5857
[289,] 72776.40  518.3705
[290,] 72782.71  511.7724
[291,] 72738.89  508.7696
[292,] 72751.35  507.8670
[293,] 72732.48  504.4222
[294,] 72751.42  503.7538
[295,] 72713.41  502.3398
[296,] 72713.41  501.4075
[297,] 72573.04  496.2171
[298,] 72566.57  495.9787
[299,] 72598.58  494.8455
[300,] 72476.12  490.3539
[301,] 72411.65  487.0633
[302,] 72365.94  484.4790
[303,] 72392.05  482.6250
[304,] 72365.86  480.2746
[305,] 72378.87  479.3390
[306,] 72431.60  477.4757
[307,] 72424.92  476.2984
[308,] 72412.03  474.3988
[309,] 72418.49  473.6918
[310,] 72418.49  470.8463
[311,] 72432.02  470.3712
[312,] 72425.23  469.1796
[313,] 72384.66  467.7201
[314,] 72378.09  466.9990
[315,] 72350.95  463.6322
[316,] 72330.59  461.9306
[317,] 72289.50  459.9950
[318,] 72296.67  459.2721
[319,] 72296.67  459.2721
[320,] 72227.88  455.3720
[321,] 72186.27  451.4770
[322,] 72200.27  450.9851
[323,] 72200.01  448.5421
[324,] 72172.36  446.5785
[325,] 72186.22  446.0922
[326,] 72193.31  445.8439
[327,] 72235.44  443.8718
[328,] 72122.67  439.9151
[329,] 72108.63  439.4246
[330,] 72073.03  437.6997
[331,] 72087.15  436.7018
[332,] 71979.96  432.9808
[333,] 71987.20  431.2320
[334,] 71987.47  430.2302
[335,] 71972.90  429.2393
[336,] 72009.15  426.4654
[337,] 72023.53  424.9455
[338,] 72030.88  423.6883
[339,] 72103.76  420.1508
[340,] 72118.54  419.6382
[341,] 72125.41  417.8619
[342,] 72169.74  415.3213
[343,] 72110.64  412.2738
[344,] 72081.28  410.7496
[345,] 72036.56  408.2037
[346,] 72036.56  407.6864
[347,] 72014.49  405.8955
[348,] 71991.79  404.1021
[349,] 71999.38  401.2772
[350,] 72014.31  400.2409
[351,] 72067.16  396.3547
[352,] 72165.79  392.4785
[353,] 72135.09  390.9290
[354,] 72150.48  389.8803
[355,] 72196.22  387.8193
[356,] 72234.63  385.9910
[357,] 72219.10  384.9545
[358,] 72196.05  384.1700
[359,] 72211.65  383.6414
[360,] 72227.02  382.5909
[361,] 72288.93  379.9707
[362,] 72258.03  378.9193
[363,] 72258.03  376.8039
[364,] 72242.21  375.2173
[365,] 72273.65  371.5039
[366,] 72265.69  370.7181
[367,] 72265.69  370.1806
[368,] 72257.99  369.3811
[369,] 72297.50  367.5191
[370,] 72273.67  366.7130
[371,] 72122.23  360.0161
[372,] 72050.04  357.6068
[373,] 72090.29  354.1294
[374,] 72106.55  351.9780
[375,] 72154.86  348.2152
[376,] 72170.90  347.6817
[377,] 72146.62  344.7157
[378,] 72154.53  343.9074
[379,] 72154.53  343.3606
[380,] 72179.31  341.4640
[381,] 72195.57  339.8491
[382,] 72179.27  338.2094
[383,] 72187.60  337.9344
[384,] 72170.90  335.7468
[385,] 72129.36  333.8432
[386,] 72162.62  332.7461
[387,] 72113.07  330.0057
[388,] 72121.20  329.7379
[389,] 72146.25  326.7087
[390,] 72146.25  325.0585
[391,] 72146.25  325.0585
[392,] 72121.00  322.0154
[393,] 72129.55  321.7351
[394,] 72129.23  321.1843
[395,] 72070.07  318.1285
[396,] 72121.09  316.4571
[397,] 72112.46  315.0692
[398,] 72138.09  313.6671
[399,] 72138.09  313.1014
[400,] 72121.02  312.5455
[401,] 72181.40  310.5785
[402,] 72189.82  309.7361
[403,] 72146.66  307.7616
[404,] 72163.90  307.2013
[405,] 72163.90  307.2013
[406,] 72207.70  303.5103
[407,] 72216.22  302.6611
[408,] 72242.83  301.2483
[409,] 72260.29  300.1099
[410,] 72277.78  299.5465
[411,] 72251.31  298.6859
[412,] 72242.35  296.1008
[413,] 72251.33  295.8121
[414,] 72287.00  294.1074
> (tri <- trimean(x))
           [,1]     [,2]
  [1,] 75012.70 988.3290
  [2,] 74886.86 972.7373
  [3,] 74771.34 958.9496
  [4,] 74660.14 945.7329
  [5,] 74550.95 932.6421
  [6,] 74442.07 919.4002
  [7,] 74349.15 908.5917
  [8,] 74260.32 898.3655
  [9,] 74173.36 888.2756
 [10,] 74090.95 878.8107
 [11,] 74044.59 873.8650
 [12,] 74004.48 869.5407
 [13,] 73969.43 865.6887
 [14,] 73939.60 862.2893
 [15,] 73914.52 859.3265
 [16,] 73889.77 856.3942
 [17,] 73866.21 853.6026
 [18,] 73842.56 850.7864
 [19,] 73819.84 848.0526
 [20,] 73798.08 845.3743
 [21,] 73776.50 842.7141
 [22,] 73758.08 840.3299
 [23,] 73739.57 837.9482
 [24,] 73721.04 835.5949
 [25,] 73702.25 833.2632
 [26,] 73683.67 830.9303
 [27,] 73665.20 828.6071
 [28,] 73646.95 826.3067
 [29,] 73628.94 824.0148
 [30,] 73611.33 821.7451
 [31,] 73593.47 819.5028
 [32,] 73575.36 817.2474
 [33,] 73557.15 814.9870
 [34,] 73539.57 812.7637
 [35,] 73522.38 810.5664
 [36,] 73505.39 808.3885
 [37,] 73488.61 806.2649
 [38,] 73472.62 804.1883
 [39,] 73457.89 802.1865
 [40,] 73443.08 800.1942
 [41,] 73428.40 798.2178
 [42,] 73412.96 796.2770
 [43,] 73397.60 794.3345
 [44,] 73382.98 792.4288
 [45,] 73368.33 790.5162
 [46,] 73353.51 788.5952
 [47,] 73338.50 786.6610
 [48,] 73323.62 784.7436
 [49,] 73309.77 782.8930
 [50,] 73295.91 781.0347
 [51,] 73282.45 779.2184
 [52,] 73269.18 777.4320
 [53,] 73256.29 775.6715
 [54,] 73243.86 773.9555
 [55,] 73231.42 772.2255
 [56,] 73219.31 770.5095
 [57,] 73206.38 768.8324
 [58,] 73193.56 767.1577
 [59,] 73181.56 765.5457
 [60,] 73169.87 763.9411
 [61,] 73158.11 762.3440
 [62,] 73146.12 760.7546
 [63,] 73134.19 759.1634
 [64,] 73122.62 757.6040
 [65,] 73111.72 756.0831
 [66,] 73100.89 754.5906
 [67,] 73089.89 753.1089
 [68,] 73079.29 751.6553
 [69,] 73069.18 750.2247
 [70,] 73059.31 748.7964
 [71,] 73049.59 747.3773
 [72,] 73040.00 745.9568
 [73,] 73030.65 744.5446
 [74,] 73021.37 743.1263
 [75,] 73011.97 741.6963
 [76,] 73002.62 740.2583
 [77,] 72993.57 738.8335
 [78,] 72985.03 737.4321
 [79,] 72976.61 736.0520
 [80,] 72968.42 734.6805
 [81,] 72960.51 733.3224
 [82,] 72952.64 731.9540
 [83,] 72945.31 730.6199
 [84,] 72937.95 729.2709
 [85,] 72930.59 727.9204
 [86,] 72923.20 726.5729
 [87,] 72915.89 725.2198
 [88,] 72908.65 723.8640
 [89,] 72901.20 722.5243
 [90,] 72893.75 721.1770
 [91,] 72886.28 719.8175
 [92,] 72878.72 718.4514
 [93,] 72871.21 717.0842
 [94,] 72863.68 715.7043
 [95,] 72856.11 714.3086
 [96,] 72848.44 712.9042
 [97,] 72841.08 711.5060
 [98,] 72834.09 710.1254
 [99,] 72834.09 708.7323
[100,] 72820.26 707.3380
[101,] 72813.76 705.9534
[102,] 72808.06 704.6264
[103,] 72808.06 703.2922
[104,] 72796.71 701.9474
[105,] 72791.12 700.6017
[106,] 72785.50 699.2434
[107,] 72785.50 697.8709
[108,] 72774.37 696.4957
[109,] 72769.08 695.1315
[110,] 72763.62 693.7669
[111,] 72759.03 692.4548
[112,] 72754.36 691.1336
[113,] 72749.86 689.8283
[114,] 72745.34 688.5174
[115,] 72740.83 687.2024
[116,] 72736.33 685.8735
[117,] 72731.84 684.5404
[118,] 72726.93 683.2150
[119,] 72722.03 681.8954
[120,] 72717.00 680.5729
[121,] 72711.99 679.2598
[122,] 72707.00 677.9376
[123,] 72702.02 676.6046
[124,] 72696.99 675.2672
[125,] 72692.00 673.9171
[126,] 72687.04 672.5781
[127,] 72682.00 671.2364
[128,] 72676.84 669.8880
[129,] 72671.66 668.5301
[130,] 72666.48 667.1607
[131,] 72661.33 665.7817
[132,] 72656.20 664.4156
[133,] 72651.20 663.0578
[134,] 72646.43 661.7033
[135,] 72641.84 660.3486
[136,] 72637.36 659.0076
[137,] 72632.81 657.6638
[138,] 72628.35 656.3214
[139,] 72623.84 654.9745
[140,] 72619.51 653.6239
[141,] 72615.19 652.2654
[142,] 72610.91 650.8973
[143,] 72606.94 649.5347
[144,] 72602.88 648.1636
[145,] 72598.79 646.7934
[146,] 72594.78 645.4262
[147,] 72590.77 644.0472
[148,] 72586.83 642.6601
[149,] 72582.74 641.2678
[150,] 72578.99 639.8905
[151,] 72575.46 638.5259
[152,] 72572.12 637.1573
[153,] 72568.81 635.7863
[154,] 72565.26 634.4270
[155,] 72561.70 633.0576
[156,] 72558.69 631.7209
[157,] 72556.02 630.4166
[158,] 72553.22 629.1063
[159,] 72550.52 627.7868
[160,] 72547.67 626.4781
[161,] 72544.72 625.1730
[162,] 72541.64 623.8923
[163,] 72538.51 622.5966
[164,] 72535.87 621.3211
[165,] 72533.22 620.0327
[166,] 72530.58 618.7599
[167,] 72530.58 617.4925
[168,] 72524.48 616.2291
[169,] 72521.45 614.9550
[170,] 72518.68 613.6936
[171,] 72515.98 612.4251
[172,] 72513.25 611.1415
[173,] 72510.41 609.8502
[174,] 72507.76 608.5717
[175,] 72507.76 607.2818
[176,] 72502.65 605.9927
[177,] 72500.06 604.7078
[178,] 72497.19 603.4228
[179,] 72494.55 602.1729
[180,] 72491.65 600.9214
[181,] 72488.86 599.6650
[182,] 72486.21 598.4098
[183,] 72486.21 597.1412
[184,] 72480.78 595.8857
[185,] 72478.10 594.6231
[186,] 72475.49 593.3616
[187,] 72473.03 592.1013
[188,] 72470.61 590.8276
[189,] 72468.12 589.5444
[190,] 72465.62 588.2516
[191,] 72463.33 586.9557
[192,] 72460.84 585.6624
[193,] 72458.62 584.3829
[194,] 72456.01 583.1183
[195,] 72454.08 581.8962
[196,] 72452.11 580.6629
[197,] 72449.78 579.4388
[198,] 72449.78 578.2078
[199,] 72445.35 576.9896
[200,] 72443.50 575.7905
[201,] 72441.93 574.6064
[202,] 72440.44 573.4160
[203,] 72439.17 572.2281
[204,] 72437.80 571.0335
[205,] 72436.68 569.8434
[206,] 72436.68 568.6424
[207,] 72434.48 567.4328
[208,] 72433.41 566.2075
[209,] 72432.46 564.9865
[210,] 72431.48 563.7542
[211,] 72430.29 562.5169
[212,] 72429.20 561.2680
[213,] 72428.20 560.0121
[214,] 72428.20 558.7531
[215,] 72425.82 557.4776
[216,] 72424.50 556.2035
[217,] 72422.82 554.9670
[218,] 72421.23 553.7236
[219,] 72419.79 552.4731
[220,] 72418.38 551.2155
[221,] 72416.70 549.9615
[222,] 72415.04 548.6910
[223,] 72413.45 547.4201
[224,] 72411.71 546.1531
[225,] 72409.92 544.9254
[226,] 72407.93 543.7067
[227,] 72405.72 542.4878
[228,] 72403.53 541.2620
[229,] 72403.53 540.0402
[230,] 72398.44 538.8043
[231,] 72395.95 537.5637
[232,] 72393.27 536.3396
[233,] 72390.64 535.1060
[234,] 72388.21 533.8731
[235,] 72385.70 532.6257
[236,] 72383.17 531.3733
[237,] 72383.17 530.1133
[238,] 72377.94 528.8482
[239,] 72375.27 527.5730
[240,] 72372.62 526.2901
[241,] 72369.99 524.9994
[242,] 72367.52 523.7114
[243,] 72364.93 522.4152
[244,] 72362.26 521.1082
[245,] 72362.26 519.7987
[246,] 72357.09 518.4780
[247,] 72354.45 517.1491
[248,] 72351.65 515.8086
[249,] 72349.02 514.4652
[250,] 72346.01 513.1410
[251,] 72342.66 511.8545
[252,] 72339.36 510.5629
[253,] 72336.09 509.2581
[254,] 72332.90 507.9403
[255,] 72329.08 506.6544
[256,] 72325.27 505.3553
[257,] 72321.37 504.0612
[258,] 72317.74 502.7683
[259,] 72314.17 501.4648
[260,] 72310.72 500.1780
[261,] 72307.33 498.8862
[262,] 72303.99 497.5948
[263,] 72300.82 496.2959
[264,] 72297.52 494.9830
[265,] 72294.28 493.6536
[266,] 72290.94 492.3236
[267,] 72287.55 490.9850
[268,] 72284.18 489.6320
[269,] 72280.33 488.2883
[270,] 72276.50 486.9359
[271,] 72272.76 485.5750
[272,] 72269.34 484.2235
[273,] 72265.75 482.8596
[274,] 72262.14 481.4838
[275,] 72258.28 480.1122
[276,] 72254.86 478.7713
[277,] 72251.73 477.4666
[278,] 72248.77 476.1719
[279,] 72245.84 474.8692
[280,] 72242.81 473.5608
[281,] 72239.37 472.2576
[282,] 72235.95 470.9761
[283,] 72232.48 469.6744
[284,] 72229.22 468.3804
[285,] 72225.82 467.0838
[286,] 72222.36 465.7728
[287,] 72218.99 464.4537
[288,] 72215.41 463.1315
[289,] 72211.81 461.8037
[290,] 72208.14 460.4546
[291,] 72204.41 459.1710
[292,] 72200.95 457.9054
[293,] 72197.38 456.6286
[294,] 72193.91 455.3755
[295,] 72190.30 454.1081
[296,] 72186.91 452.8369
[297,] 72183.50 451.5548
[298,] 72180.98 450.3217
[299,] 72178.48 449.0688
[300,] 72175.76 447.8077
[301,] 72173.82 446.5859
[302,] 72172.28 445.3868
[303,] 72171.03 444.2008
[304,] 72169.60 443.0176
[305,] 72168.33 441.8443
[306,] 72166.97 440.6614
[307,] 72165.26 439.4812
[308,] 72163.58 438.2947
[309,] 72161.98 437.1117
[310,] 72160.32 435.9155
[311,] 72158.65 434.7359
[312,] 72156.89 433.5398
[313,] 72155.15 432.3370
[314,] 72153.67 431.1314
[315,] 72152.22 429.9123
[316,] 72150.93 428.7172
[317,] 72149.77 427.5226
[318,] 72148.87 426.3321
[319,] 72148.87 425.1280
[320,] 72146.95 423.8996
[321,] 72146.43 422.7032
[322,] 72146.17 421.5389
[323,] 72145.82 420.3580
[324,] 72145.47 419.1884
[325,] 72145.30 418.0234
[326,] 72145.03 416.8415
[327,] 72144.72 415.6388
[328,] 72144.13 414.4402
[329,] 72144.27 413.2752
[330,] 72144.50 412.0931
[331,] 72144.96 410.9117
[332,] 72145.34 409.7205
[333,] 72146.41 408.5592
[334,] 72146.41 407.3995
[335,] 72147.45 406.2302
[336,] 72149.63 405.0507
[337,] 72150.54 403.8879
[338,] 72151.37 402.7232
[339,] 72152.16 401.5525
[340,] 72152.47 400.4103
[341,] 72152.69 399.2512
[342,] 72152.87 398.0939
[343,] 72152.76 396.9502
[344,] 72153.04 395.8280
[345,] 72153.51 394.7044
[346,] 72154.27 393.5949
[347,] 72155.04 392.4687
[348,] 72155.96 391.3450
[349,] 72157.04 390.2238
[350,] 72157.04 389.1215
[351,] 72159.02 388.0104
[352,] 72159.62 386.9351
[353,] 72159.58 385.8958
[354,] 72159.75 384.8568
[355,] 72159.81 383.8101
[356,] 72159.57 382.7716
[357,] 72159.07 381.7376
[358,] 72158.67 380.6959
[359,] 72158.43 379.6422
[360,] 72158.07 378.5721
[361,] 72157.62 377.4936
[362,] 72156.74 376.4316
[363,] 72156.07 375.3615
[364,] 72155.39 374.3002
[365,] 72155.39 373.2394
[366,] 72154.81 372.2134
[367,] 72153.27 371.1754
[368,] 72152.52 370.1209
[369,] 72151.82 369.0539
[370,] 72150.84 367.9915
[371,] 72150.01 366.9166
[372,] 72150.20 365.9271
[373,] 72150.88 364.9525
[374,] 72151.28 364.0113
[375,] 72151.59 363.0818
[376,] 72151.56 362.1910
[377,] 72151.43 361.2860
[378,] 72151.47 360.4070
[379,] 72151.44 359.5185
[380,] 72151.42 358.6159
[381,] 72151.42 357.7216
[382,] 72150.93 356.8310
[383,] 72150.74 355.9446
[384,] 72150.49 355.0392
[385,] 72150.35 354.1473
[386,] 72150.49 353.2640
[387,] 72150.41 352.3758
[388,] 72150.66 351.5109
[389,] 72150.87 350.6270
[390,] 72150.90 349.7716
[391,] 72150.93 348.9214
[392,] 72150.96 348.0474
[393,] 72151.17 347.2027
[394,] 72151.32 346.3392
[395,] 72151.48 345.4610
[396,] 72151.48 344.6127
[397,] 72152.04 343.7702
[398,] 72152.54 342.9283
[399,] 72152.64 342.0875
[400,] 72152.74 341.2326
[401,] 72152.97 340.3629
[402,] 72152.77 339.5040
[403,] 72152.50 338.6355
[404,] 72152.55 337.7781
[405,] 72152.47 336.9057
[406,] 72152.38 336.0073
[407,] 72151.99 335.1520
[408,] 72151.53 334.2869
[409,] 72150.87 333.4224
[410,] 72150.09 332.5534
[411,] 72149.16 331.6685
[412,] 72149.16 330.7735
[413,] 72147.75 329.9015
[414,] 72146.99 329.0080
> (midr <- midrange(x))
[1] 166674.5
> midm <- array(NA,dim=8)
> for (j in 1:8) midm[j] <- midmean(x,j)
> midm
[1] 72160.32 72160.32 72160.32 72160.32 72192.98 72160.32 72160.32 72160.32
> postscript(file="/var/fisher/rcomp/tmp/11r7e1385468517.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/fisher/rcomp/tmp/2okcq1385468517.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/fisher/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/fisher/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/fisher/rcomp/tmp/37bkb1385468518.tab") 
> 
> try(system("convert tmp/11r7e1385468517.ps tmp/11r7e1385468517.png",intern=TRUE))
character(0)
> try(system("convert tmp/2okcq1385468517.ps tmp/2okcq1385468517.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
 21.604   1.143  22.722