R version 2.8.0 (2008-10-20)
Copyright (C) 2008 The R Foundation for Statistical Computing
ISBN 3-900051-07-0

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

  Natural language support but running in an English locale

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(0.00099950033305429
+ ,0.00099850232957088
+ ,-0.0009985023296184
+ ,-0.0030015030037291
+ ,0.001001502337127
+ ,-0.004012041489962
+ ,-0.0040282028286014
+ ,0.0020161297152037
+ ,0.0040201059166187
+ ,0.001002506349653
+ ,-0.0030105391528727
+ ,0.0010045204259987
+ ,-6.6613381477509E-15
+ ,0.0050075217305934
+ ,-0.0009995003331289
+ ,-4.5741188614556E-14
+ ,-0.001000500333546
+ ,0.001000500333538
+ ,-4.5741188614556E-14
+ ,-4.5741188614556E-14
+ ,-0.001000500333546
+ ,-0.0010015023370622
+ ,-0.001002506349673
+ ,0.002004008686753
+ ,0.001000500333538
+ ,-0.002002002670646
+ ,2.7089441800854E-14
+ ,-0.004016069654873
+ ,0.0020100509280176
+ ,-6.6613381477509E-15
+ ,0.0020060187268931
+ ,-0.001002506349673
+ ,0.002004008686753
+ ,-0.0020040086867623
+ ,-4.7073456244107E-14
+ ,-4.7073456244107E-14
+ ,-4.7073456244107E-14
+ ,0.001002506349653
+ ,0.0030015030037274
+ ,-2.8865798640254E-14
+ ,-0.0020000006666292
+ ,-0.0020040086867623
+ ,-4.7073456244107E-14
+ ,0.001002506349653
+ ,-0.001002506349673
+ ,-4.7073456244107E-14
+ ,0.001002506349653
+ ,-0.001002506349673
+ ,-0.0010035123772467
+ ,-6.6613381477509E-15
+ ,-0.0010045204260067
+ ,0.0020080328031984
+ ,0.001002506349653
+ ,-0.0020060187268727
+ ,-6.6613381477509E-15
+ ,-0.0010045204260067
+ ,-1.3322676295502E-15
+ ,-1.3322676295502E-15
+ ,-0.0010055305020016
+ ,1.7319479184152E-14
+ ,1.7319479184152E-14
+ ,0.0040160696549174
+ ,2.7089441800854E-14
+ ,2.7089441800854E-14
+ ,-0.0030105391528727
+ ,0.0010045204259987
+ ,-0.002010050928007
+ ,1.7319479184152E-14
+ ,-0.0010065426113823
+ ,1.9095836023553E-14
+ ,1.9095836023553E-14
+ ,0.002012073113419
+ ,0.0010045204259987
+ ,-0.002010050928007
+ ,1.7319479184152E-14
+ ,0.0020100509280176
+ ,-0.0030165935394066
+ ,0.0010065426114187
+ ,1.7319479184152E-14
+ ,0.0010055305020176
+ ,-0.0020120731134012
+ ,-0.001007556760281
+ ,1.8651746813703E-14
+ ,1.8651746813703E-14
+ ,0.0020140993717184
+ ,-0.0020140993716824
+ ,0.0050276626769183
+ ,0.002004008686753
+ ,0.0020000006666381
+ ,-2.8865798640254E-14
+ ,0.00099850232957088
+ ,-0.0019980026627184
+ ,-0.001000500333546
+ ,3.7747582837255E-14
+ ,0.001000500333538
+ ,-4.5741188614556E-14
+ ,0.00099950033305429
+ ,-0.0009995003331289
+ ,0.00099950033305429
+ ,0.0049825713444709
+ ,-0.00099453016649331
+ ,-0.00099552024148508
+ ,0.001990050408017
+ ,0.00099354205890689
+ ,-0.00099354205887092
+ ,0.00099354205890689
+ ,-0.00099354205887092
+ ,-0.0019900504079935
+ ,0.00099552024151706
+ ,-0.0019920325312848
+ ,-0.0019960086467439
+ ,0.001996008646671
+ ,0.00099651228975572
+ ,1.6875389974302E-14
+ ,0.00099552024151706
+ ,0.0029806281381153
+ ,-0.0029806281381224
+ ,1.5099033134902E-14
+ ,0.0019880722254153
+ ,2.9309887850104E-14
+ ,0.00099255591272929
+ ,-2.2204460492503E-14
+ ,-0.00099255591272218
+ ,-0.00099354205887092
+ ,0.00099354205890689
+ ,0.00099255591272929
+ ,0.0019821612039777
+ ,-1.3766765505352E-14
+ ,-0.00099058948171393
+ ,-0.00099157172231745
+ ,0.0019821612039777
+ ,-0.0019821612040136
+ ,-2.2204460492503E-14
+ ,-0.0029806281381224
+ ,-0.00099552024148508
+ ,0.001990050408017
+ ,0.00099354205890689
+ ,-0.00099354205887092
+ ,-0.00099453016649331
+ ,1.5099033134902E-14
+ ,1.5099033134902E-14
+ ,-0.0029895388483849
+ ,-0.0009985023296184
+ ,0.00099850232957088
+ ,0.00099750631708151
+ ,-4.3964831775156E-14
+ ,-4.3964831775156E-14
+ ,-0.00099750631714413
+ ,0.00099750631708151
+ ,0.00099651228975572
+ ,0.00099552024151706
+ ,1.5099033134902E-14
+ ,-0.0029895388483849
+ ,0.00099750631708151
+ ,-4.3964831775156E-14
+ ,-4.3964831775156E-14
+ ,-0.00099750631714413
+ ,-0.0009985023296184
+ ,-0.0009995003331289
+ ,-4.5741188614556E-14
+ ,0.00099950033305429
+ ,-0.0009995003331289
+ ,-0.002002002670646
+ ,0.001001502337127
+ ,3.7747582837255E-14
+ ,0.0020000006666381
+ ,-2.8865798640254E-14
+ ,-0.0009995003331289
+ ,0.001998002662654
+ ,-0.0009985023296184
+ ,-2.8865798640254E-14
+ ,0.00099850232957088
+ ,-0.0009985023296184
+ ,0.001996008646671
+ ,-0.0029955089798439
+ ,-4.5741188614556E-14
+ ,-4.5741188614556E-14
+ ,-4.5741188614556E-14
+ ,-0.001000500333546
+ ,0.001000500333538
+ ,0.0029955089797542
+ ,-4.3964831775156E-14
+ ,-4.3964831775156E-14
+ ,-4.3964831775156E-14
+ ,-4.3964831775156E-14
+ ,0.0019920325312559
+ ,0.00099453016651507
+ ,-0.00099453016649331
+ ,0.00099453016651507
+ ,0.00099354205890689
+ ,-0.0019880722253709
+ ,1.5099033134902E-14
+ ,-0.00099552024148508
+ ,0.00099552024151706
+ ,0.0019880722254153
+ ,2.9309887850104E-14
+ ,0.00099255591272929
+ ,-2.2204460492503E-14
+ ,-2.2204460492503E-14
+ ,-0.0019860979716224
+ ,-0.0029865626977932
+ ,0.00099651228975572
+ ,0.001990050408017
+ ,0.0019860979716069
+ ,-0.00099255591272218
+ ,2.9309887850104E-14
+ ,-0.0029835924668711
+ ,-0.00099651228978281
+ ,-0.0019960086467439
+ ,-0.0009995003331289
+ ,0.00099950033305429
+ ,0.001996008646671
+ ,-0.0019960086467439
+ ,0.00099850232957088
+ ,-1.8651746813703E-14
+ ,-1.8651746813703E-14
+ ,0.00099750631708151
+ ,-4.3964831775156E-14
+ ,0.00099651228975572
+ ,-0.001994018606883
+ ,-0.0009985023296184
+ ,-2.8865798640254E-14
+ ,-0.0009995003331289
+ ,0.00099950033305429
+ ,0.001996008646671
+ ,-4.3964831775156E-14
+ ,-0.00099750631714413
+ ,-1.8651746813703E-14
+ ,-0.0029985029962187
+ ,0.001000500333538
+ ,0.001998002662654
+ ,-0.0009985023296184
+ ,-2.8865798640254E-14
+ ,-2.8865798640254E-14
+ ,-2.8865798640254E-14
+ ,-0.0009995003331289
+ ,-0.0030045090203461
+ ,-4.7073456244107E-14
+ ,0.001002506349653
+ ,0.001001502337127
+ ,-0.0020040086867623
+ ,0.001002506349653
+ ,2.7089441800854E-14
+ ,-0.001002506349673
+ ,0.001002506349653
+ ,-0.001002506349673
+ ,-0.0010035123772467
+ ,-6.6613381477509E-15
+ ,0.001003512377193
+ ,-0.0010035123772467
+ ,0.001003512377193
+ ,-4.7073456244107E-14
+ ,-4.7073456244107E-14
+ ,-0.0010035123772467
+ ,-6.6613381477509E-15
+ ,0.0020060187268931
+ ,-0.0020060187268727
+ ,-0.0010045204260067
+ ,-0.0010055305020016
+ ,-0.0010065426113823
+ ,-0.001007556760281
+ ,1.8651746813703E-14
+ ,-0.0020181641562815
+ ,0.0010095912013557
+ ,-0.0020202027072966
+ ,0.0030287756621794
+ ,1.8651746813703E-14
+ ,0.0020140993717184
+ ,0.0010055305020176
+ ,-0.0020120731134012
+ ,0.0030165935394191
+ ,-0.0010045204260067
+ ,-1.3322676295502E-15
+ ,-1.3322676295502E-15
+ ,0.0020080328031984
+ ,-4.7073456244107E-14
+ ,0.001002506349653
+ ,2.7089441800854E-14
+ ,-0.001002506349673
+ ,0.0030045090202533
+ ,-0.001000500333546
+ ,-0.0010015023370622
+ ,0.0030015030037274
+ ,-0.0009995003331289
+ ,-0.002002002670646
+ ,0.0020020026706273
+ ,-4.5741188614556E-14
+ ,0.001998002662654
+ ,-0.0019980026627184
+ ,-0.001000500333546
+ ,-0.0010015023370622
+ ,-0.001002506349673
+ ,0.002004008686753
+ ,3.7747582837255E-14
+ ,0.0020000006666381
+ ,-2.8865798640254E-14
+ ,-2.8865798640254E-14
+ ,-0.0020000006666292
+ ,3.7747582837255E-14
+ ,-0.0020040086867623
+ ,0.001002506349653
+ ,-0.0020060187268727
+ ,-6.6613381477509E-15
+ ,-6.6613381477509E-15
+ ,0.001003512377193
+ ,-4.7073456244107E-14
+ ,-4.7073456244107E-14
+ ,0.001002506349653
+ ,2.7089441800854E-14
+ ,0.0020020026706273
+ ,-0.001000500333546
+ ,3.7747582837255E-14
+ ,-0.0010015023370622
+ ,0.0020020026706273
+ ,-0.001000500333546
+ ,-0.0010015023370622
+ ,0.0040000053333271
+ ,-0.0009985023296184
+ ,0.00099850232957088
+ ,0.00099750631708151
+ ,-4.3964831775156E-14
+ ,-4.3964831775156E-14
+ ,-4.3964831775156E-14
+ ,0.00099651228975572
+ ,0.00099552024151706
+ ,-0.0019920325312848
+ ,-4.3964831775156E-14
+ ,-4.3964831775156E-14
+ ,0.0019920325312559
+ ,1.5099033134902E-14
+ ,1.5099033134902E-14
+ ,0.00099453016651507
+ ,0.00099354205890689
+ ,2.9309887850104E-14
+ ,0.00099255591272929
+ ,-2.2204460492503E-14
+ ,-0.00099255591272218
+ ,0.00099255591272929
+ ,0.0019821612039777
+ ,0.00098960918518642
+ ,-0.0029717703891796
+ ,-2.2204460492503E-14
+ ,-0.00099255591272218
+ ,-0.00099354205887092
+ ,0.00099354205890689
+ ,2.9309887850104E-14
+ ,-0.00099354205887092
+ ,0.0019860979716069
+ ,0.00099157172227748
+ ,-1.7763568394003E-14
+ ,-1.7763568394003E-14
+ ,-1.7763568394003E-14
+ ,-0.00099157172231745
+ ,0.00099157172227748
+ ,-1.7763568394003E-14
+ ,0.0009905894816824
+ ,-1.3766765505352E-14
+ ,0.00098960918518642
+ ,-0.00098960918517976
+ ,-1.3766765505352E-14
+ ,0.001978240012086
+ ,-0.00098863082691913
+ ,0.00098863082692002
+ ,0.00098765440128057
+ ,-0.0019762852281917
+ ,0.0029629651306204
+ ,0.001970443987263
+ ,-3.5971225997855E-14
+ ,0.00098376791016408
+ ,-0.00098376791016852
+ ,-3.5971225997855E-14
+ ,-0.0019704439873358
+ ,0.002954211897463
+ ,-0.00098376791016852
+ ,0.001966568972064
+ ,2.3536728122053E-14
+ ,0.0039215736532237
+ ,4.1744385725906E-14
+ ,0.00097799518794162
+ ,-3.5083047578155E-14
+ ,-0.0039177327289353
+ ,-3.3306690738755E-14
+ ,-3.3306690738755E-14
+ ,-3.3306690738755E-14
+ ,-0.00098183611223313
+ ,0.0029426210542236
+ ,-0.0019607849419745
+ ,-0.00098183611223313
+ ,2.3536728122053E-14
+ ,2.3536728122053E-14
+ ,0.00098183611222336
+ ,-0.00098183611223313
+ ,2.3536728122053E-14
+ ,2.3536728122053E-14
+ ,2.3536728122053E-14
+ ,2.3536728122053E-14
+ ,2.3536728122053E-14
+ ,2.3536728122053E-14
+ ,2.3536728122053E-14
+ ,2.3536728122053E-14
+ ,0.0019627091678234
+ ,-0.00098087305562533
+ ,0.0019607849419669
+ ,-0.0039254221160743
+ ,-0.0019685045726687
+ ,3.5527136788005E-15
+ ,3.5527136788005E-15
+ ,0.00098473666250376
+ ,0.0029484050842639
+ ,-0.00098183611223313
+ ,0.00098183611222336
+ ,-0.0029484050843331
+ ,-3.5971225997855E-14
+ ,-0.00098473666253618
+ ,-0.00098570732479608
+ ,-0.00098667990243717
+ ,0.00098667990240831
+ ,-0.00098667990243717
+ ,0.00098667990240831
+ ,-3.6859404417555E-14
+ ,0.002954211897463
+ ,3.1530333899354E-14
+ ,0.00098280106193149
+ ,2.3536728122053E-14
+ ,-0.00098280106187643
+ ,3.1530333899354E-14
+ ,-0.00098376791016852
+ ,0.0029484050842639
+ ,0.00098087305556671
+ ,0.00097991188637492
+ ,2.5757174171304E-14
+ ,-0.00097991188637447
+ ,-2.5313084961454E-14
+ ,-0.0019627091678251
+ ,0.0019627091678234
+ ,-2.5313084961454E-14
+ ,-0.00098087305562533
+ ,0.00098087305556671
+ ,-0.00098087305562533
+ ,-0.00098183611223313
+ ,2.3536728122053E-14
+ ,0.0029426210542236
+ ,2.5757174171304E-14
+ ,2.5757174171304E-14
+ ,2.5757174171304E-14
+ ,2.5757174171304E-14
+ ,-0.00097991188637447
+ ,0.00097991188637492
+ ,-0.0019607849419745
+ ,0.00098087305556671
+ ,-2.5313084961454E-14
+ ,-0.00098087305562533
+ ,0.00098087305556671
+ ,-0.0019627091678251
+ ,2.3536728122053E-14
+ ,-0.00098280106187643
+ ,3.1530333899354E-14
+ ,3.1530333899354E-14
+ ,-0.00098376791016852
+ ,-0.0019704439873358
+ ,0.003937012959363
+ ,0.0019627091678234
+ ,-0.00098087305562533
+ ,-0.0029484050843331
+ ,0.0029484050842639
+ ,-3.3306690738755E-14
+ ,-0.0019646371741331
+ ,3.1530333899354E-14
+ ,3.1530333899354E-14
+ ,0.0019646371741313
+ ,-0.00098183611223313
+ ,2.3536728122053E-14
+ ,2.3536728122053E-14
+ ,2.3536728122053E-14
+ ,2.3536728122053E-14
+ ,2.3536728122053E-14
+ ,0.0048995688411235
+ ,-3.5083047578155E-14
+ ,0.00097703964786477
+ ,3.8191672047105E-14
+ ,0.00097608597303811
+ ,0.00097513415818273
+ ,0.0019474202843766
+ ,-0.00097323608651889
+ ,-0.001949318356067
+ ,-1.7319479184152E-14
+ ,-1.7319479184152E-14
+ ,0.00097513415818273
+ ,-2.3536728122053E-14
+ ,0.00097418419787676
+ ,-0.001949318356067
+ ,0.001949318356083
+ ,-0.001949318356067
+ ,-0.0019531256209171
+ ,-3.5083047578155E-14
+ ,-3.5083047578155E-14
+ ,-3.5083047578155E-14
+ ,-0.00097799518793495
+ ,0.0019550348358415
+ ,0.00097608597303811
+ ,-1.7319479184152E-14
+ ,-0.0019531256209171
+ ,-0.00097799518793495
+ ,0.00097799518794162
+ ,-0.00097799518793495
+ ,-0.0009789525989583
+ ,0.0019569477869257
+ ,-0.00097799518793495
+ ,4.1744385725906E-14
+ ,4.1744385725906E-14
+ ,4.1744385725906E-14
+ ,4.1744385725906E-14
+ ,-0.0009789525989583
+ ,0.00097895259902581
+ ,4.1744385725906E-14
+ ,-0.0009789525989583
+ ,0.00097895259902581
+ ,4.1744385725906E-14
+ ,4.1744385725906E-14
+ ,-0.0009789525989583
+ ,2.5757174171304E-14
+ ,-0.00097991188637447
+ ,0.00097991188637492
+ ,2.5757174171304E-14
+ ,0.00097895259902581
+ ,-0.0009789525989583
+ ,-0.0019607849419745
+ ,-3.3306690738755E-14
+ ,0.00098087305556671
+ ,0.001958864485375
+ ,4.1744385725906E-14
+ ,0.0068259650704414
+ ,-0.00097228981895814
+ ,-0.004875680063519
+ ,-0.00097799518793495
+ ,4.1744385725906E-14
+ ,4.1744385725906E-14
+ ,0.0019550348358415
+ ,3.8191672047105E-14
+ ,0.0019512201312382
+ ,-2.3536728122053E-14
+ ,-0.00097513415822359
+ ,0.0029225544425828
+ ,-0.00097323608651889
+ ,0.0009732360865331
+ ,0.00097228981898079
+ ,4.1744385725906E-14
+ ,0.00097134538964205
+ ,-0.00097134538959009
+ ,4.1744385725906E-14
+ ,0.0019417481829418
+ ,-0.001941748182869
+ ,4.1744385725906E-14
+ ,0.0019417481829418
+ ,-0.00097040279326865
+ ,-0.00097134538959009
+ ,0.00097134538964205
+ ,0.00097040279330995
+ ,3.1086244689504E-14
+ ,-0.00097040279326865
+ ,0.0019398648178099
+ ,0.00096852307808337
+ ,-4.7073456244107E-14
+ ,-4.7073456244107E-14
+ ,-0.00096852307814732
+ ,-1.6875389974302E-14
+ ,-1.6875389974302E-14
+ ,-1.6875389974302E-14
+ ,0.00096852307808337
+ ,0.0019342365798529
+ ,-0.0019342365798778
+ ,-0.00096852307814732
+ ,-0.0019398648178166
+ ,0.0038759738447101
+ ,-0.00096758594878299
+ ,-4.7073456244107E-14
+ ,-0.0019379851026473
+ ,0.00096946202453108
+ ,0.0019361090268832
+ ,-0.0019361090268832
+ ,0.0029027596579834
+ ,2.2204460492503E-14
+ ,2.2204460492503E-14
+ ,-0.00096665063107793
+ ,0.00096665063111701
+ ,2.2204460492503E-14
+ ,-0.00096665063107793
+ ,1.6875389974302E-14
+ ,0.00096665063111701
+ ,-0.00096665063107793
+ ,0.0019323677510168
+ ,-0.0019323677510368
+ ,1.6875389974302E-14
+ ,0.00096665063111701
+ ,0.00096571711992199
+ ,-3.6859404417555E-14
+ ,0.0009647854100634
+ ,-0.00096478541013578
+ ,-0.0019323677510368
+ ,0.0019323677510168
+ ,-0.00096571711993665
+ ,2.2204460492503E-14
+ ,-0.00096665063107793
+ ,0.0019323677510168
+ ,0.0009647854100634
+ ,-0.00096478541013578
+ ,0.0009647854100634
+ ,-3.5527136788005E-14
+ ,-0.00096478541013578
+ ,-3.6859404417555E-14
+ ,0.0009647854100634
+ ,0.0009638554962641
+ ,0.00096292737335757
+ ,-0.00096292737343617
+ ,-4.2632564145606E-14
+ ,-0.0019286409064425
+ ,-3.6859404417555E-14
+ ,0.0009647854100634
+ ,0.0019267828696643
+ ,-3.5971225997855E-14
+ ,0.00096200103616395
+ ,-2.7089441800854E-14
+ ,0.001921230177873
+ ,-0.00096015369832081
+ ,0.00096015369832259
+ ,0.00095923268747899
+ ,-0.00095923268748077
+ ,-2.0872192862953E-14
+ ,-0.00096015369832081
+ ,2.2648549702353E-14)
> ylimmax = ''
> ylimmin = ''
> main = 'Robustness of Central Tendency'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> geomean <- function(x) {
+ 
+ return(exp(mean(log(x))))
+ }
> harmean <- function(x) {
+ return(1/mean(1/x))
+ }
> quamean <- function(x) {
+ 
+ return(sqrt(mean(x*x)))
+ }
> winmean <- function(x) {
+ x <-sort(x[!is.na(x)])
+ n<-length(x)
+ denom <- 3
+ 
+ nodenom <- n/denom
+ if (nodenom>40) denom <- n/40
+ sqrtn = sqrt(n)
+ roundnodenom = floor(nodenom)
+ win <- array(NA,dim=c(roundnodenom,2))
+ for (j in 1:roundnodenom) {
+ 
+ win[j,1] <- (j*x[j+1]+sum(x[(j+1):(n-j)])+j*x[n-j])/n
+ win[j,2] <- sd(c(rep(x[j+1],j),x[(j+1):(n-j)],rep(x[n-j],j)))/sqrtn
+ }
+ return(win)
+ }
> trimean <- function(x) {
+ 
+ x <-sort(x[!is.na(x)])
+ n<-length(x)
+ denom <- 3
+ nodenom <- n/denom
+ if (nodenom>40) denom <- n/40
+ sqrtn = sqrt(n)
+ 
+ roundnodenom = floor(nodenom)
+ tri <- array(NA,dim=c(roundnodenom,2))
+ for (j in 1:roundnodenom) {
+ tri[j,1] <- mean(x,trim=j/n)
+ tri[j,2] <- sd(x[(j+1):(n-j)]) / sqrt(n-j*2)
+ }
+ 
+ return(tri)
+ }
> midrange <- function(x) {
+ return((max(x)+min(x))/2)
+ }
> q1 <- function(data,n,p,i,f) {
+ 
+ np <- n*p;
+ i <<- floor(np)
+ f <<- np - i
+ qvalue <- (1-f)*data[i] + f*data[i+1]
+ }
> q2 <- function(data,n,p,i,f) {
+ 
+ np <- (n+1)*p
+ i <<- floor(np)
+ f <<- np - i
+ qvalue <- (1-f)*data[i] + f*data[i+1]
+ }
> q3 <- function(data,n,p,i,f) {
+ 
+ np <- n*p
+ i <<- floor(np)
+ f <<- np - i
+ if (f==0) {
+ qvalue <- data[i]
+ } else {
+ 
+ qvalue <- data[i+1]
+ }
+ }
> q4 <- function(data,n,p,i,f) {
+ np <- n*p
+ i <<- floor(np)
+ 
+ f <<- np - i
+ if (f==0) {
+ qvalue <- (data[i]+data[i+1])/2
+ } else {
+ qvalue <- data[i+1]
+ }
+ 
+ }
> q5 <- function(data,n,p,i,f) {
+ np <- (n-1)*p
+ i <<- floor(np)
+ f <<- np - i
+ if (f==0) {
+ 
+ qvalue <- data[i+1]
+ } else {
+ qvalue <- data[i+1] + f*(data[i+2]-data[i+1])
+ }
+ }
> q6 <- function(data,n,p,i,f) {
+ 
+ np <- n*p+0.5
+ i <<- floor(np)
+ f <<- np - i
+ qvalue <- data[i]
+ }
> q7 <- function(data,n,p,i,f) {
+ 
+ np <- (n+1)*p
+ i <<- floor(np)
+ f <<- np - i
+ if (f==0) {
+ qvalue <- data[i]
+ } else {
+ 
+ qvalue <- f*data[i] + (1-f)*data[i+1]
+ }
+ }
> q8 <- function(data,n,p,i,f) {
+ np <- (n+1)*p
+ i <<- floor(np)
+ 
+ f <<- np - i
+ if (f==0) {
+ qvalue <- data[i]
+ } else {
+ if (f == 0.5) {
+ qvalue <- (data[i]+data[i+1])/2
+ 
+ } else {
+ if (f < 0.5) {
+ qvalue <- data[i]
+ } else {
+ qvalue <- data[i+1]
+ }
+ 
+ }
+ }
+ }
> midmean <- function(x,def) {
+ x <-sort(x[!is.na(x)])
+ n<-length(x)
+ 
+ if (def==1) {
+ qvalue1 <- q1(x,n,0.25,i,f)
+ qvalue3 <- q1(x,n,0.75,i,f)
+ }
+ if (def==2) {
+ qvalue1 <- q2(x,n,0.25,i,f)
+ 
+ qvalue3 <- q2(x,n,0.75,i,f)
+ }
+ if (def==3) {
+ qvalue1 <- q3(x,n,0.25,i,f)
+ qvalue3 <- q3(x,n,0.75,i,f)
+ }
+ 
+ if (def==4) {
+ qvalue1 <- q4(x,n,0.25,i,f)
+ qvalue3 <- q4(x,n,0.75,i,f)
+ }
+ if (def==5) {
+ qvalue1 <- q5(x,n,0.25,i,f)
+ 
+ qvalue3 <- q5(x,n,0.75,i,f)
+ }
+ if (def==6) {
+ qvalue1 <- q6(x,n,0.25,i,f)
+ qvalue3 <- q6(x,n,0.75,i,f)
+ }
+ 
+ if (def==7) {
+ qvalue1 <- q7(x,n,0.25,i,f)
+ qvalue3 <- q7(x,n,0.75,i,f)
+ }
+ if (def==8) {
+ qvalue1 <- q8(x,n,0.25,i,f)
+ 
+ qvalue3 <- q8(x,n,0.75,i,f)
+ }
+ midm <- 0
+ myn <- 0
+ roundno4 <- round(n/4)
+ round3no4 <- round(3*n/4)
+ 
+ for (i in 1:n) {
+ if ((x[i]>=qvalue1) & (x[i]<=qvalue3)){
+ midm = midm + x[i]
+ myn = myn + 1
+ }
+ }
+ 
+ midm = midm / myn
+ return(midm)
+ }
> (arm <- mean(x))
[1] 6.378062e-05
> sqrtn <- sqrt(length(x))
> (armse <- sd(x) / sqrtn)
[1] 5.836973e-05
> 
> (armose <- arm / armse)
[1] 1.092700
> (geo <- geomean(x))
[1] NaN
Warning message:
In log(x) : NaNs produced
> (har <- harmean(x))
[1] -1.961217e-13
> (qua <- quamean(x))
[1] 0.001465294
> (win <- winmean(x))
               [,1]         [,2]
  [1,] 6.227137e-05 5.774507e-05
  [2,] 6.224595e-05 5.773201e-05
  [3,] 6.214632e-05 5.771372e-05
  [4,] 6.216929e-05 5.758167e-05
  [5,] 5.525045e-05 5.672037e-05
  [6,] 6.379429e-05 5.586123e-05
  [7,] 6.368306e-05 5.583530e-05
  [8,] 6.288316e-05 5.574628e-05
  [9,] 6.274874e-05 5.571441e-05
 [10,] 6.207265e-05 5.563094e-05
 [11,] 4.733268e-05 5.417088e-05
 [12,] 4.715767e-05 5.414552e-05
 [13,] 4.703150e-05 5.411282e-05
 [14,] 4.696470e-05 5.410702e-05
 [15,] 4.703556e-05 5.410070e-05
 [16,] 4.695877e-05 5.408077e-05
 [17,] 4.663721e-05 5.403891e-05
 [18,] 4.613255e-05 5.399546e-05
 [19,] 4.613571e-05 5.394903e-05
 [20,] 4.687746e-05 5.388327e-05
 [21,] 4.668390e-05 5.386668e-05
 [22,] 7.909732e-05 5.137214e-05
 [23,] 7.917174e-05 5.136731e-05
 [24,] 7.910625e-05 5.133771e-05
 [25,] 7.918665e-05 5.133250e-05
 [26,] 7.835851e-05 5.125925e-05
 [27,] 7.759682e-05 5.117899e-05
 [28,] 3.819105e-05 4.814958e-05
 [29,] 3.828320e-05 4.813094e-05
 [30,] 3.828320e-05 4.813094e-05
 [31,] 3.818349e-05 4.812443e-05
 [32,] 3.818288e-05 4.811085e-05
 [33,] 3.818288e-05 4.811085e-05
 [34,] 3.807396e-05 4.810375e-05
 [35,] 3.807396e-05 4.810375e-05
 [36,] 3.807350e-05 4.808853e-05
 [37,] 3.807350e-05 4.808853e-05
 [38,] 3.795226e-05 4.808064e-05
 [39,] 3.807619e-05 4.807229e-05
 [40,] 3.807619e-05 4.807229e-05
 [41,] 3.820622e-05 4.806353e-05
 [42,] 3.807249e-05 4.805483e-05
 [43,] 3.820859e-05 4.804567e-05
 [44,] 3.820859e-05 4.804567e-05
 [45,] 3.806559e-05 4.803638e-05
 [46,] 3.821089e-05 4.802660e-05
 [47,] 3.835906e-05 4.801664e-05
 [48,] 3.820683e-05 4.800675e-05
 [49,] 3.836100e-05 4.799639e-05
 [50,] 3.851799e-05 4.798585e-05
 [51,] 3.851639e-05 4.796464e-05
 [52,] 3.884133e-05 4.794285e-05
 [53,] 3.933566e-05 4.790976e-05
 [54,] 3.983634e-05 4.787630e-05
 [55,] 3.948922e-05 4.785377e-05
 [56,] 3.966161e-05 4.784226e-05
 [57,] 3.983218e-05 4.780729e-05
 [58,] 4.000968e-05 4.779546e-05
 [59,] 4.000968e-05 4.779546e-05
 [60,] 4.000454e-05 4.777104e-05
 [61,] 4.000454e-05 4.777104e-05
 [62,] 3.981025e-05 4.775845e-05
 [63,] 4.057618e-05 4.770750e-05
 [64,] 4.017626e-05 4.768159e-05
 [65,] 4.056907e-05 4.765549e-05
 [66,] 4.056907e-05 4.765549e-05
 [67,] 4.015206e-05 4.762850e-05
 [68,] 4.012768e-05 4.751989e-05
 [69,] 3.990955e-05 4.747882e-05
 [70,] 3.990376e-05 4.745113e-05
 [71,] 3.989791e-05 4.742311e-05
 [72,] 4.011190e-05 4.740894e-05
 [73,] 4.032845e-05 4.739461e-05
 [74,] 4.010243e-05 4.738004e-05
 [75,] 4.031747e-05 4.733595e-05
 [76,] 1.512012e-04 4.107069e-05
 [77,] 1.509674e-04 4.105436e-05
 [78,] 1.510929e-04 4.104872e-05
 [79,] 1.506146e-04 4.101532e-05
 [80,] 1.505016e-04 4.099271e-05
 [81,] 1.502576e-04 4.097568e-05
 [82,] 1.496508e-04 4.091835e-05
 [83,] 1.496508e-04 4.091835e-05
 [84,] 1.493996e-04 4.090087e-05
 [85,] 1.488929e-04 4.086562e-05
 [86,] 1.487749e-04 4.084167e-05
 [87,] 1.485168e-04 4.082374e-05
 [88,] 1.485168e-04 4.082374e-05
 [89,] 1.485168e-04 4.082374e-05
 [90,] 1.478627e-04 4.076190e-05
 [91,] 1.470607e-04 4.070629e-05
 [92,] 1.393242e-05 3.272256e-05
 [93,] 1.348239e-05 3.270082e-05
 [94,] 1.333138e-05 3.269353e-05
 [95,] 1.333138e-05 3.269353e-05
 [96,] 1.317746e-05 3.268610e-05
 [97,] 1.317746e-05 3.268610e-05
 [98,] 1.333364e-05 3.267839e-05
 [99,] 1.317524e-05 3.267075e-05
[100,] 1.317524e-05 3.267075e-05
[101,] 1.317524e-05 3.267075e-05
[102,] 1.301236e-05 3.266291e-05
[103,] 1.317617e-05 3.265482e-05
[104,] 1.317617e-05 3.265482e-05
[105,] 1.317617e-05 3.265482e-05
[106,] 1.317617e-05 3.265482e-05
[107,] 1.317617e-05 3.265482e-05
[108,] 1.317617e-05 3.265482e-05
[109,] 1.317617e-05 3.265482e-05
[110,] 1.317617e-05 3.265482e-05
[111,] 1.335237e-05 3.264613e-05
[112,] 1.335237e-05 3.264613e-05
[113,] 1.317228e-05 3.263745e-05
[114,] 1.317228e-05 3.263745e-05
[115,] 1.317228e-05 3.263745e-05
[116,] 1.298778e-05 3.262858e-05
[117,] 1.298778e-05 3.262858e-05
[118,] 1.298778e-05 3.262858e-05
[119,] 1.298778e-05 3.262858e-05
[120,] 1.317788e-05 3.261920e-05
[121,] 1.298582e-05 3.260996e-05
[122,] 1.298582e-05 3.260996e-05
[123,] 1.298582e-05 3.260996e-05
[124,] 1.298582e-05 3.260996e-05
[125,] 1.298582e-05 3.260996e-05
[126,] 1.298582e-05 3.260996e-05
[127,] 1.278463e-05 3.260029e-05
[128,] 1.298700e-05 3.259031e-05
[129,] 1.298700e-05 3.259031e-05
[130,] 1.298700e-05 3.259031e-05
[131,] 1.319369e-05 3.258012e-05
[132,] 1.340155e-05 3.256988e-05
[133,] 1.319128e-05 3.255977e-05
[134,] 1.319128e-05 3.255977e-05
[135,] 1.340344e-05 3.254933e-05
[136,] 1.340344e-05 3.254933e-05
[137,] 1.340344e-05 3.254933e-05
[138,] 1.340214e-05 3.252821e-05
[139,] 1.340214e-05 3.252821e-05
[140,] 1.340214e-05 3.252821e-05
[141,] 1.340214e-05 3.252821e-05
[142,] 1.340214e-05 3.252821e-05
[143,] 1.317696e-05 3.251739e-05
[144,] 1.340237e-05 3.250630e-05
[145,] 1.340237e-05 3.250630e-05
[146,] 1.340237e-05 3.250630e-05
[147,] 1.340237e-05 3.250630e-05
[148,] 1.340099e-05 3.248376e-05
[149,] 1.340099e-05 3.248376e-05
[150,] 1.363485e-05 3.247225e-05
[151,] 1.363298e-05 3.244933e-05
[152,] 1.386903e-05 3.243773e-05
[153,] 1.434282e-05 3.241446e-05
[154,] 1.434282e-05 3.241446e-05
[155,] 1.458211e-05 3.240272e-05
[156,] 1.482246e-05 3.239092e-05
[157,] 1.506388e-05 3.237908e-05
[158,] 1.481656e-05 3.236720e-05
[159,] 1.481656e-05 3.236720e-05
[160,] 1.481656e-05 3.236720e-05
[161,] 1.506365e-05 3.235508e-05
[162,] 1.506365e-05 3.235508e-05
[163,] 1.505867e-05 3.233061e-05
[164,] 1.505867e-05 3.233061e-05
[165,] 1.480142e-05 3.231826e-05
[166,] 1.454311e-05 3.230587e-05
[167,] 1.454311e-05 3.230587e-05
[168,] 1.428222e-05 3.229335e-05
[169,] 1.427863e-05 3.226812e-05
[170,] 1.401567e-05 3.225551e-05
[171,] 1.401567e-05 3.225551e-05
[172,] 1.348514e-05 3.223008e-05
[173,] 1.321912e-05 3.221734e-05
[174,] 1.321755e-05 3.219153e-05
[175,] 1.294951e-05 3.217870e-05
[176,] 1.294951e-05 3.217870e-05
[177,] 1.294845e-05 3.215253e-05
[178,] 1.294845e-05 3.215253e-05
[179,] 1.294845e-05 3.215253e-05
[180,] 1.294845e-05 3.215253e-05
[181,] 1.294845e-05 3.215253e-05
[182,] 1.294736e-05 3.212567e-05
[183,] 1.294736e-05 3.212567e-05
[184,] 1.294736e-05 3.212567e-05
[185,] 1.266567e-05 3.211220e-05
[186,] 1.266567e-05 3.211220e-05
[187,] 1.351489e-05 3.207056e-05
[188,] 1.379560e-05 3.202913e-05
[189,] 1.379560e-05 3.202913e-05
[190,] 1.379280e-05 3.200136e-05
[191,] 1.379000e-05 3.197350e-05
[192,] 1.379000e-05 3.197350e-05
[193,] 1.378717e-05 3.194540e-05
[194,] 1.378717e-05 3.194540e-05
[195,] 1.407495e-05 3.190284e-05
[196,] 1.377999e-05 3.188873e-05
[197,] 1.377713e-05 3.186023e-05
[198,] 1.377427e-05 3.183164e-05
[199,] 1.377427e-05 3.183164e-05
[200,] 1.347503e-05 3.181733e-05
[201,] 1.317488e-05 3.180298e-05
[202,] 1.317314e-05 3.177392e-05
[203,] 1.347335e-05 3.175922e-05
[204,] 1.286704e-05 3.173024e-05
[205,] 1.286704e-05 3.173024e-05
[206,] 1.347459e-05 3.170048e-05
[207,] 1.407921e-05 3.164120e-05
[208,] 1.377160e-05 3.162649e-05
[209,] 1.407714e-05 3.161154e-05
[210,] 3.338214e-04 1.829824e-05
> (tri <- trimean(x))
               [,1]         [,2]
  [1,] 6.087819e-05 5.701466e-05
  [2,] 5.947611e-05 5.626456e-05
  [3,] 5.807787e-05 5.550076e-05
  [4,] 5.670428e-05 5.472218e-05
  [5,] 5.531599e-05 5.395786e-05
  [6,] 5.532936e-05 5.336401e-05
  [7,] 5.388647e-05 5.291438e-05
  [8,] 5.245049e-05 5.245848e-05
  [9,] 5.110805e-05 5.200456e-05
 [10,] 4.977223e-05 5.154407e-05
 [11,] 4.849768e-05 5.108243e-05
 [12,] 4.860779e-05 5.076546e-05
 [13,] 4.873383e-05 5.044415e-05
 [14,] 4.887087e-05 5.011886e-05
 [15,] 4.901383e-05 4.978699e-05
 [16,] 4.915278e-05 4.944833e-05
 [17,] 4.929772e-05 4.910369e-05
 [18,] 4.946371e-05 4.875442e-05
 [19,] 4.966065e-05 4.840032e-05
 [20,] 4.985875e-05 4.804130e-05
 [21,] 5.001847e-05 4.767833e-05
 [22,] 5.018918e-05 4.730802e-05
 [23,] 4.877167e-05 4.707603e-05
 [24,] 4.734092e-05 4.683931e-05
 [25,] 4.590326e-05 4.659904e-05
 [26,] 4.445215e-05 4.635380e-05
 [27,] 4.302580e-05 4.610696e-05
 [28,] 4.162048e-05 4.585870e-05
 [29,] 4.175538e-05 4.574867e-05
 [30,] 4.188771e-05 4.563708e-05
 [31,] 4.202098e-05 4.552306e-05
 [32,] 4.215876e-05 4.540684e-05
 [33,] 4.229755e-05 4.528864e-05
 [34,] 4.243732e-05 4.516786e-05
 [35,] 4.258170e-05 4.504472e-05
 [36,] 4.272711e-05 4.491888e-05
 [37,] 4.287358e-05 4.479085e-05
 [38,] 4.302111e-05 4.466001e-05
 [39,] 4.317335e-05 4.452657e-05
 [40,] 4.332306e-05 4.439048e-05
 [41,] 4.347386e-05 4.425138e-05
 [42,] 4.362210e-05 4.410949e-05
 [43,] 4.362210e-05 4.396474e-05
 [44,] 4.392560e-05 4.381707e-05
 [45,] 4.407718e-05 4.366608e-05
 [46,] 4.423362e-05 4.351201e-05
 [47,] 4.438751e-05 4.335479e-05
 [48,] 4.453883e-05 4.319432e-05
 [49,] 4.469505e-05 4.303053e-05
 [50,] 4.484871e-05 4.286333e-05
 [51,] 4.499978e-05 4.269264e-05
 [52,] 4.515204e-05 4.251869e-05
 [53,] 4.529795e-05 4.234140e-05
 [54,] 4.543372e-05 4.216103e-05
 [55,] 4.555930e-05 4.197750e-05
 [56,] 4.569353e-05 4.179035e-05
 [57,] 4.582504e-05 4.159914e-05
 [58,] 4.595391e-05 4.140447e-05
 [59,] 4.608001e-05 4.120553e-05
 [60,] 4.620711e-05 4.100182e-05
 [61,] 4.633531e-05 4.079394e-05
 [62,] 4.646453e-05 4.058102e-05
 [63,] 4.659869e-05 4.036329e-05
 [64,] 4.671866e-05 4.014173e-05
 [65,] 4.684746e-05 3.991547e-05
 [66,] 4.696965e-05 3.968436e-05
 [67,] 4.709283e-05 3.944747e-05
 [68,] 4.722495e-05 3.920539e-05
 [69,] 4.735859e-05 3.896037e-05
 [70,] 4.749739e-05 3.871028e-05
 [71,] 4.763744e-05 3.845456e-05
 [72,] 4.777875e-05 3.819302e-05
 [73,] 4.791735e-05 3.792505e-05
 [74,] 4.805323e-05 3.765042e-05
 [75,] 4.819425e-05 3.736888e-05
 [76,] 4.833267e-05 3.708106e-05
 [77,] 4.654123e-05 3.695330e-05
 [78,] 4.473871e-05 3.682252e-05
 [79,] 4.291876e-05 3.668831e-05
 [80,] 4.109144e-05 3.655134e-05
 [81,] 4.109144e-05 3.641120e-05
 [82,] 3.739764e-05 3.626767e-05
 [83,] 3.553895e-05 3.612172e-05
 [84,] 3.366416e-05 3.597172e-05
 [85,] 3.177717e-05 3.581803e-05
 [86,] 3.177717e-05 3.566101e-05
 [87,] 2.797190e-05 3.550022e-05
 [88,] 2.604919e-05 3.533540e-05
 [89,] 2.410946e-05 3.516591e-05
 [90,] 2.215249e-05 3.499161e-05
 [91,] 2.215249e-05 3.481401e-05
 [92,] 1.821888e-05 3.463286e-05
 [93,] 1.828499e-05 3.464077e-05
 [94,] 1.835860e-05 3.464857e-05
 [95,] 1.843517e-05 3.465594e-05
 [96,] 1.851244e-05 3.466272e-05
 [97,] 1.859275e-05 3.466904e-05
 [98,] 1.867379e-05 3.467474e-05
 [99,] 1.875325e-05 3.467996e-05
[100,] 1.883580e-05 3.468468e-05
[101,] 1.891912e-05 3.468872e-05
[102,] 1.900323e-05 3.469208e-05
[103,] 1.909050e-05 3.469488e-05
[104,] 1.917622e-05 3.469713e-05
[105,] 1.926276e-05 3.469862e-05
[106,] 1.935013e-05 3.469935e-05
[107,] 1.943833e-05 3.469928e-05
[108,] 1.952739e-05 3.469841e-05
[109,] 1.961732e-05 3.469669e-05
[110,] 1.970812e-05 3.469412e-05
[111,] 1.979981e-05 3.469065e-05
[112,] 1.988994e-05 3.468646e-05
[113,] 1.998097e-05 3.468132e-05
[114,] 2.007539e-05 3.467540e-05
[115,] 2.017077e-05 3.466848e-05
[116,] 2.026710e-05 3.466055e-05
[117,] 2.036693e-05 3.465175e-05
[118,] 2.046778e-05 3.464187e-05
[119,] 2.056965e-05 3.463087e-05
[120,] 2.067258e-05 3.461874e-05
[121,] 2.077399e-05 3.460562e-05
[122,] 2.087904e-05 3.459147e-05
[123,] 2.098518e-05 3.457609e-05
[124,] 2.109244e-05 3.455942e-05
[125,] 2.120083e-05 3.454143e-05
[126,] 2.131036e-05 3.452207e-05
[127,] 2.142106e-05 3.450131e-05
[128,] 2.153561e-05 3.447930e-05
[129,] 2.164872e-05 3.445601e-05
[130,] 2.176304e-05 3.443118e-05
[131,] 2.187861e-05 3.440478e-05
[132,] 2.199273e-05 3.437696e-05
[133,] 2.210538e-05 3.434768e-05
[134,] 2.222202e-05 3.431688e-05
[135,] 2.233996e-05 3.428431e-05
[136,] 2.245645e-05 3.425012e-05
[137,] 2.257425e-05 3.421404e-05
[138,] 2.269338e-05 3.417600e-05
[139,] 2.281388e-05 3.413641e-05
[140,] 2.293576e-05 3.409473e-05
[141,] 2.305904e-05 3.405091e-05
[142,] 2.318375e-05 3.400486e-05
[143,] 2.330990e-05 3.395652e-05
[144,] 2.344043e-05 3.390604e-05
[145,] 2.356960e-05 3.385336e-05
[146,] 2.370029e-05 3.379815e-05
[147,] 2.383255e-05 3.374031e-05
[148,] 2.396638e-05 3.367976e-05
[149,] 2.410184e-05 3.361693e-05
[150,] 2.423895e-05 3.355120e-05
[151,] 2.437474e-05 3.348274e-05
[152,] 2.451221e-05 3.341173e-05
[153,] 2.464836e-05 3.333780e-05
[154,] 2.478015e-05 3.326112e-05
[155,] 2.491358e-05 3.318102e-05
[156,] 2.504563e-05 3.309766e-05
[157,] 2.517628e-05 3.301091e-05
[158,] 2.530551e-05 3.292066e-05
[159,] 2.543956e-05 3.282675e-05
[160,] 2.557534e-05 3.272874e-05
[161,] 2.571288e-05 3.262647e-05
[162,] 2.571288e-05 3.252009e-05
[163,] 2.598703e-05 3.240910e-05
[164,] 2.612689e-05 3.229397e-05
[165,] 2.626862e-05 3.217387e-05
[166,] 2.641554e-05 3.204892e-05
[167,] 2.656777e-05 3.191892e-05
[168,] 2.672206e-05 3.178331e-05
[169,] 2.688182e-05 3.164219e-05
[170,] 2.704383e-05 3.149569e-05
[171,] 2.721147e-05 3.134319e-05
[172,] 2.721147e-05 3.118407e-05
[173,] 2.756068e-05 3.101875e-05
[174,] 2.774588e-05 3.084658e-05
[175,] 2.793375e-05 3.066763e-05
[176,] 2.812779e-05 3.048118e-05
[177,] 2.832464e-05 3.028646e-05
[178,] 2.852438e-05 3.008388e-05
[179,] 2.872706e-05 2.987220e-05
[180,] 2.893274e-05 2.965096e-05
[181,] 2.914149e-05 2.941963e-05
[182,] 2.914149e-05 2.917767e-05
[183,] 2.956849e-05 2.892542e-05
[184,] 2.978689e-05 2.866136e-05
[185,] 3.000865e-05 2.838480e-05
[186,] 3.023756e-05 2.809549e-05
[187,] 3.047005e-05 2.779213e-05
[188,] 3.069494e-05 2.747550e-05
[189,] 3.091967e-05 2.714471e-05
[190,] 3.114799e-05 2.679718e-05
[191,] 3.138003e-05 2.643288e-05
[192,] 3.161588e-05 2.605067e-05
[193,] 3.185560e-05 2.564801e-05
[194,] 3.209931e-05 2.522461e-05
[195,] 3.234709e-05 2.477755e-05
[196,] 3.259513e-05 2.430688e-05
[197,] 3.285139e-05 2.380917e-05
[198,] 3.311207e-05 2.328271e-05
[199,] 3.337728e-05 2.272477e-05
[200,] 3.364711e-05 2.213061e-05
[201,] 3.364711e-05 2.149709e-05
[202,] 3.421359e-05 2.081965e-05
[203,] 3.450654e-05 2.009373e-05
[204,] 3.480058e-05 1.931184e-05
[205,] 3.510847e-05 1.846653e-05
[206,] 3.542201e-05 1.754521e-05
[207,] 3.573275e-05 1.653599e-05
[208,] 3.604071e-05 1.542249e-05
[209,] 3.635886e-05 1.417405e-05
[210,] 3.667870e-05 1.274889e-05
> 
> (midr <- midrange(x))
[1] 0.0009751425
> midm <- array(NA,dim=8)
> for (j in 1:8) midm[j] <- midmean(x,j)
> midm
[1] 2.210210e-05 4.322036e-05 4.322036e-05 4.322036e-05 2.210210e-05
[6] 4.322036e-05 4.322036e-05 4.322036e-05
> postscript(file="/var/www/html/freestat/rcomp/tmp/1oxxk1290187063.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> lb <- win[,1] - 2*win[,2]
> 
> ub <- win[,1] + 2*win[,2]
> if ((ylimmin == '') | (ylimmax == '')) plot(win[,1],type='b',main=main, xlab='j', pch=19, ylab='Winsorized Mean(j/n)', ylim=c(min(lb),max(ub))) else plot(win[,1],type='l',main=main, xlab='j', pch=19, ylab='Winsorized Mean(j/n)', ylim=c(ylimmin,ylimmax))
> lines(ub,lty=3)
> lines(lb,lty=3)
> grid()
> dev.off()
null device 
          1 
> 
> postscript(file="/var/www/html/freestat/rcomp/tmp/2y7e51290187063.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> lb <- tri[,1] - 2*tri[,2]
> ub <- tri[,1] + 2*tri[,2]
> if ((ylimmin == '') | (ylimmax == '')) plot(tri[,1],type='b',main=main, xlab='j', pch=19, ylab='Trimmed Mean(j/n)', ylim=c(min(lb),max(ub))) else plot(tri[,1],type='l',main=main, xlab='j', pch=19, ylab='Trimmed Mean(j/n)', ylim=c(ylimmin,ylimmax))
> lines(ub,lty=3)
> lines(lb,lty=3)
> 
> grid()
> dev.off()
null device 
          1 
> 
> #Note: the /var/www/html/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/freestat/rcomp/createtable")
> 
> a<-table.start()
> a<-table.row.start(a)
> a<-table.element(a,'Central Tendency - Ungrouped Data',4,TRUE)
> 
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,'Measure',header=TRUE)
> a<-table.element(a,'Value',header=TRUE)
> a<-table.element(a,'S.E.',header=TRUE)
> a<-table.element(a,'Value/S.E.',header=TRUE)
> 
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/arithmetic_mean.htm', 'Arithmetic Mean', 'click to view the definition of the Arithmetic Mean'),header=TRUE)
> a<-table.element(a,arm)
> a<-table.element(a,hyperlink('http://www.xycoon.com/arithmetic_mean_standard_error.htm', armse, 'click to view the definition of the Standard Error of the Arithmetic Mean'))
> a<-table.element(a,armose)
> 
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/geometric_mean.htm', 'Geometric Mean', 'click to view the definition of the Geometric Mean'),header=TRUE)
> a<-table.element(a,geo)
> a<-table.element(a,'')
> a<-table.element(a,'')
> 
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/harmonic_mean.htm', 'Harmonic Mean', 'click to view the definition of the Harmonic Mean'),header=TRUE)
> a<-table.element(a,har)
> a<-table.element(a,'')
> a<-table.element(a,'')
> 
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/quadratic_mean.htm', 'Quadratic Mean', 'click to view the definition of the Quadratic Mean'),header=TRUE)
> a<-table.element(a,qua)
> a<-table.element(a,'')
> a<-table.element(a,'')
> 
> a<-table.row.end(a)
> for (j in 1:length(win[,1])) {
+ a<-table.row.start(a)
+ mylabel <- paste('Winsorized Mean (',j)
+ mylabel <- paste(mylabel,'/')
+ mylabel <- paste(mylabel,length(win[,1]))
+ 
+ mylabel <- paste(mylabel,')')
+ a<-table.element(a,hyperlink('http://www.xycoon.com/winsorized_mean.htm', mylabel, 'click to view the definition of the Winsorized Mean'),header=TRUE)
+ a<-table.element(a,win[j,1])
+ a<-table.element(a,win[j,2])
+ a<-table.element(a,win[j,1]/win[j,2])
+ a<-table.row.end(a)
+ 
+ }
> for (j in 1:length(tri[,1])) {
+ a<-table.row.start(a)
+ mylabel <- paste('Trimmed Mean (',j)
+ mylabel <- paste(mylabel,'/')
+ mylabel <- paste(mylabel,length(tri[,1]))
+ 
+ mylabel <- paste(mylabel,')')
+ a<-table.element(a,hyperlink('http://www.xycoon.com/arithmetic_mean.htm', mylabel, 'click to view the definition of the Trimmed Mean'),header=TRUE)
+ a<-table.element(a,tri[j,1])
+ a<-table.element(a,tri[j,2])
+ a<-table.element(a,tri[j,1]/tri[j,2])
+ a<-table.row.end(a)
+ 
+ }
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/median_1.htm', 'Median', 'click to view the definition of the Median'),header=TRUE)
> a<-table.element(a,median(x))
> a<-table.element(a,'')
> a<-table.element(a,'')
> 
> a<-table.row.end(a)
> a<-table.row.start(a)
> a<-table.element(a,hyperlink('http://www.xycoon.com/midrange.htm', 'Midrange', 'click to view the definition of the Midrange'),header=TRUE)
> a<-table.element(a,midr)
> a<-table.element(a,'')
> a<-table.element(a,'')
> 
> a<-table.row.end(a)
> a<-table.row.start(a)
> mymid <- hyperlink('http://www.xycoon.com/midmean.htm', 'Midmean', 'click to view the definition of the Midmean')
> mylabel <- paste(mymid,hyperlink('http://www.xycoon.com/method_1.htm','Weighted Average at Xnp',''),sep=' - ')
> a<-table.element(a,mylabel,header=TRUE)
> a<-table.element(a,midm[1])
> 
> a<-table.element(a,'')
> a<-table.element(a,'')
> a<-table.row.end(a)
> a<-table.row.start(a)
> mymid <- hyperlink('http://www.xycoon.com/midmean.htm', 'Midmean', 'click to view the definition of the Midmean')
> mylabel <- paste(mymid,hyperlink('http://www.xycoon.com/method_2.htm','Weighted Average at X(n+1)p',''),sep=' - ')
> 
> a<-table.element(a,mylabel,header=TRUE)
> a<-table.element(a,midm[2])
> a<-table.element(a,'')
> a<-table.element(a,'')
> a<-table.row.end(a)
> a<-table.row.start(a)
> 
> mymid <- hyperlink('http://www.xycoon.com/midmean.htm', 'Midmean', 'click to view the definition of the Midmean')
> mylabel <- paste(mymid,hyperlink('http://www.xycoon.com/method_3.htm','Empirical Distribution Function',''),sep=' - ')
> a<-table.element(a,mylabel,header=TRUE)
> a<-table.element(a,midm[3])
> a<-table.element(a,'')
> a<-table.element(a,'')
> 
> a<-table.row.end(a)
> a<-table.row.start(a)
> mymid <- hyperlink('http://www.xycoon.com/midmean.htm', 'Midmean', 'click to view the definition of the Midmean')
> mylabel <- paste(mymid,hyperlink('http://www.xycoon.com/method_4.htm','Empirical Distribution Function - Averaging',''),sep=' - ')
> a<-table.element(a,mylabel,header=TRUE)
> a<-table.element(a,midm[4])
> 
> a<-table.element(a,'')
> a<-table.element(a,'')
> a<-table.row.end(a)
> a<-table.row.start(a)
> mymid <- hyperlink('http://www.xycoon.com/midmean.htm', 'Midmean', 'click to view the definition of the Midmean')
> mylabel <- paste(mymid,hyperlink('http://www.xycoon.com/method_5.htm','Empirical Distribution Function - Interpolation',''),sep=' - ')
> 
> a<-table.element(a,mylabel,header=TRUE)
> a<-table.element(a,midm[5])
> a<-table.element(a,'')
> a<-table.element(a,'')
> a<-table.row.end(a)
> a<-table.row.start(a)
> 
> mymid <- hyperlink('http://www.xycoon.com/midmean.htm', 'Midmean', 'click to view the definition of the Midmean')
> mylabel <- paste(mymid,hyperlink('http://www.xycoon.com/method_6.htm','Closest Observation',''),sep=' - ')
> a<-table.element(a,mylabel,header=TRUE)
> a<-table.element(a,midm[6])
> a<-table.element(a,'')
> a<-table.element(a,'')
> 
> a<-table.row.end(a)
> a<-table.row.start(a)
> mymid <- hyperlink('http://www.xycoon.com/midmean.htm', 'Midmean', 'click to view the definition of the Midmean')
> mylabel <- paste(mymid,hyperlink('http://www.xycoon.com/method_7.htm','True Basic - Statistics Graphics Toolkit',''),sep=' - ')
> a<-table.element(a,mylabel,header=TRUE)
> a<-table.element(a,midm[7])
> 
> a<-table.element(a,'')
> a<-table.element(a,'')
> a<-table.row.end(a)
> a<-table.row.start(a)
> mymid <- hyperlink('http://www.xycoon.com/midmean.htm', 'Midmean', 'click to view the definition of the Midmean')
> mylabel <- paste(mymid,hyperlink('http://www.xycoon.com/method_8.htm','MS Excel (old versions)',''),sep=' - ')
> 
> a<-table.element(a,mylabel,header=TRUE)
> a<-table.element(a,midm[8])
> a<-table.element(a,'')
> a<-table.element(a,'')
> a<-table.row.end(a)
> a<-table.row.start(a)
> 
> a<-table.element(a,'Number of observations',header=TRUE)
> a<-table.element(a,length(x))
> a<-table.element(a,'')
> a<-table.element(a,'')
> a<-table.row.end(a)
> a<-table.end(a)
> 
> table.save(a,file="/var/www/html/freestat/rcomp/tmp/3w26n1290187064.tab") 
> try(system("convert tmp/1oxxk1290187063.ps tmp/1oxxk1290187063.png",intern=TRUE))
character(0)
> try(system("convert tmp/2y7e51290187063.ps tmp/2y7e51290187063.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  5.308   0.631   8.523