R version 2.11.1 (2010-05-31)
Copyright (C) 2010 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.

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(-4.5741188614556E-14
+ ,0.0069756137364543
+ ,2.9309887850104E-14
+ ,-0.014000228673371
+ ,1.9095836023553E-14
+ ,1.9095836023553E-14
+ ,1.9095836023553E-14
+ ,1.9095836023553E-14
+ ,1.9095836023553E-14
+ ,1.9095836023553E-14
+ ,1.9095836023553E-14
+ ,1.9095836023553E-14
+ ,1.9095836023553E-14
+ ,0.014000228673419
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,-0.014000228673371
+ ,1.9095836023553E-14
+ ,0.014000228673419
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,-0.014000228673371
+ ,0.014000228673419
+ ,0.00099255591272929
+ ,0.0059347355197779
+ ,-3.6859404417555E-14
+ ,-3.6859404417555E-14
+ ,-3.6859404417555E-14
+ ,-3.6859404417555E-14
+ ,-0.0059347355198369
+ ,-2.2204460492503E-14
+ ,0.0059347355197779
+ ,-3.6859404417555E-14
+ ,-3.6859404417555E-14
+ ,-3.6859404417555E-14
+ ,-3.6859404417555E-14
+ ,-0.0059347355198369
+ ,0.0059347355197779
+ ,-0.0059347355198369
+ ,-2.2204460492503E-14
+ ,-2.2204460492503E-14
+ ,-2.2204460492503E-14
+ ,-2.2204460492503E-14
+ ,0.0059347355197779
+ ,-3.6859404417555E-14
+ ,-0.0059347355198369
+ ,-2.2204460492503E-14
+ ,0.0059347355197779
+ ,-3.6859404417555E-14
+ ,-3.6859404417555E-14
+ ,-3.6859404417555E-14
+ ,-3.6859404417555E-14
+ ,-0.0059347355198369
+ ,-2.2204460492503E-14
+ ,-2.2204460492503E-14
+ ,-2.2204460492503E-14
+ ,-2.2204460492503E-14
+ ,0.0059347355197779
+ ,-3.6859404417555E-14
+ ,-3.6859404417555E-14
+ ,-3.6859404417555E-14
+ ,-3.6859404417555E-14
+ ,0.030113980247763
+ ,-1.9984014443253E-14
+ ,0.0057252064780799
+ ,4.0412118096356E-14
+ ,4.0412118096356E-14
+ ,4.0412118096356E-14
+ ,4.0412118096356E-14
+ ,-0.0057252064780595
+ ,-1.9984014443253E-14
+ ,-1.9984014443253E-14
+ ,0.0057252064780799
+ ,-0.0057252064780595
+ ,0.0057252064780799
+ ,0.027218949241341
+ ,2.6201263381154E-14
+ ,2.6201263381154E-14
+ ,2.6201263381154E-14
+ ,0.0092166551049262
+ ,2.2204460492503E-15
+ ,2.2204460492503E-15
+ ,-0.0092166551048978
+ ,2.6201263381154E-14
+ ,0.0092166551049262
+ ,2.2204460492503E-15
+ ,-0.0092166551048978
+ ,2.6201263381154E-14
+ ,2.6201263381154E-14
+ ,0.0092166551049262
+ ,0.027150989065902
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,0.0071174677688512
+ ,-1.2878587085652E-14
+ ,-0.007117467768913
+ ,-4.8849813083507E-14
+ ,0.0071174677688512
+ ,-0.007117467768913
+ ,0.0071174677688512
+ ,-1.2878587085652E-14
+ ,-1.2878587085652E-14
+ ,-1.2878587085652E-14
+ ,-0.007117467768913
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,0.0071174677688512
+ ,-1.2878587085652E-14
+ ,-0.007117467768913
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,0.0071174677688512
+ ,-0.007117467768913
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,0.0071174677688512
+ ,-1.2878587085652E-14
+ ,-1.2878587085652E-14
+ ,-1.2878587085652E-14
+ ,-1.2878587085652E-14
+ ,-0.007117467768913
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,0.0071174677688512
+ ,-1.2878587085652E-14
+ ,-1.2878587085652E-14
+ ,-1.2878587085652E-14
+ ,-1.2878587085652E-14
+ ,-1.2878587085652E-14
+ ,-1.2878587085652E-14
+ ,-1.2878587085652E-14
+ ,-0.007117467768913
+ ,-4.8849813083507E-14
+ ,0.0071174677688512
+ ,-0.007117467768913
+ ,0.0071174677688512
+ ,-1.2878587085652E-14
+ ,-0.007117467768913
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,0.0071174677688512
+ ,-1.2878587085652E-14
+ ,-1.2878587085652E-14
+ ,-1.2878587085652E-14
+ ,-0.007117467768913
+ ,-4.8849813083507E-14
+ ,-4.8849813083507E-14
+ ,0.014184634991951
+ ,-5.3290705182008E-15
+ ,-5.3290705182008E-15
+ ,-5.3290705182008E-15
+ ,-5.3290705182008E-15
+ ,-5.3290705182008E-15
+ ,-5.3290705182008E-15
+ ,-5.3290705182008E-15
+ ,-5.3290705182008E-15
+ ,-5.3290705182008E-15
+ ,-5.3290705182008E-15
+ ,-5.3290705182008E-15
+ ,-5.3290705182008E-15
+ ,-5.3290705182008E-15
+ ,0.00087989446939485
+ ,-4.4408920985006E-14
+ ,0.0052631700442554
+ ,-0.0052631700443193
+ ,-4.4408920985006E-14
+ ,-4.4408920985006E-14
+ ,-4.4408920985006E-14
+ ,-4.4408920985006E-14
+ ,0.0052631700442554
+ ,-1.9539925233403E-14
+ ,-0.0052631700443193
+ ,-4.4408920985006E-14
+ ,-4.4408920985006E-14
+ ,0.0052631700442554
+ ,-1.9539925233403E-14
+ ,-0.0052631700443193
+ ,-4.4408920985006E-14
+ ,0.0052631700442554
+ ,-0.0052631700443193
+ ,-4.4408920985006E-14
+ ,-4.4408920985006E-14
+ ,-4.4408920985006E-14
+ ,-4.4408920985006E-14
+ ,0.0052631700442554
+ ,-1.9539925233403E-14
+ ,-1.9539925233403E-14
+ ,-1.9539925233403E-14
+ ,-1.9539925233403E-14
+ ,-1.9539925233403E-14
+ ,-0.0052631700443193
+ ,-4.4408920985006E-14
+ ,-4.4408920985006E-14
+ ,0.0052631700442554
+ ,-1.9539925233403E-14
+ ,-0.0052631700443193
+ ,0.0052631700442554
+ ,-1.9539925233403E-14
+ ,-1.9539925233403E-14
+ ,-1.9539925233403E-14
+ ,-0.0052631700443193
+ ,0.0052631700442554
+ ,-1.9539925233403E-14
+ ,-0.0052631700443193
+ ,-4.4408920985006E-14
+ ,-4.4408920985006E-14
+ ,-4.4408920985006E-14
+ ,0.0052631700442554
+ ,-1.9539925233403E-14
+ ,0.0069747449270805
+ ,8.4376949871512E-15
+ ,0.013805223564709
+ ,3.5527136788005E-14
+ ,3.5527136788005E-14
+ ,3.5527136788005E-14
+ ,3.5527136788005E-14
+ ,3.5527136788005E-14
+ ,3.5527136788005E-14
+ ,3.5527136788005E-14
+ ,3.5527136788005E-14
+ ,0.0068317942917355
+ ,-0.006831794291668
+ ,3.5527136788005E-14
+ ,3.5527136788005E-14
+ ,3.5527136788005E-14
+ ,3.5527136788005E-14
+ ,0.0068317942917355
+ ,3.1974423109205E-14
+ ,-0.006831794291668
+ ,0.0068317942917355
+ ,3.1974423109205E-14
+ ,3.1974423109205E-14
+ ,-0.006831794291668
+ ,3.5527136788005E-14
+ ,3.5527136788005E-14
+ ,0.0051282163669355
+ ,1.5987211554602E-14
+ ,1.5987211554602E-14
+ ,1.5987211554602E-14
+ ,1.5987211554602E-14
+ ,1.5987211554602E-14
+ ,1.5987211554602E-14
+ ,0.0067969675439161
+ ,-0.0067969675438708
+ ,1.5987211554602E-14
+ ,1.5987211554602E-14
+ ,1.5987211554602E-14
+ ,0.0067969675439161
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,0.011784648168229
+ ,-1.9539925233403E-14
+ ,-1.9539925233403E-14
+ ,-1.9539925233403E-14
+ ,-1.9539925233403E-14
+ ,-1.9539925233403E-14
+ ,-1.9539925233403E-14
+ ,0.02397799870668
+ ,1.9984014443253E-14
+ ,0.0081367413930198
+ ,-4.1744385725906E-14
+ ,-4.1744385725906E-14
+ ,-4.1744385725906E-14
+ ,-0.0081367413930415
+ ,1.9984014443253E-14
+ ,1.9984014443253E-14
+ ,0.0081367413930198
+ ,-4.1744385725906E-14
+ ,-0.0081367413930415
+ ,0.03846628082782
+ ,2.3980817331903E-14
+ ,2.3980817331903E-14
+ ,2.3980817331903E-14
+ ,0.0062696130136239
+ ,2.8421709430404E-14
+ ,2.8421709430404E-14
+ ,2.8421709430404E-14
+ ,2.8421709430404E-14
+ ,-0.0062696130135715
+ ,0.0062696130136239
+ ,2.8421709430404E-14
+ ,2.8421709430404E-14
+ ,-0.0062696130135715
+ ,2.3980817331903E-14
+ ,2.3980817331903E-14
+ ,2.3980817331903E-14
+ ,2.3980817331903E-14
+ ,2.3980817331903E-14
+ ,0.012500162764224
+ ,-7.5495165674511E-15
+ ,-7.5495165674511E-15
+ ,0.0061919702478925
+ ,-2.8865798640254E-14
+ ,-2.8865798640254E-14
+ ,0.0061538655743711
+ ,-6.6613381477509E-15
+ ,-6.6613381477509E-15
+ ,-0.0061538655744067
+ ,-2.8865798640254E-14
+ ,-2.8865798640254E-14
+ ,-2.8865798640254E-14
+ ,-2.8865798640254E-14
+ ,0.010744539282971
+ ,-0.0061255934267059
+ ,-2.3536728122053E-14
+ ,-2.3536728122053E-14
+ ,-2.3536728122053E-14
+ ,-2.3536728122053E-14
+ ,0.0061255934266766
+ ,-5.7731597280508E-15
+ ,-5.7731597280508E-15
+ ,-5.7731597280508E-15
+ ,-5.7731597280508E-15
+ ,-5.7731597280508E-15
+ ,-5.7731597280508E-15
+ ,0.023388467086494
+ ,1.2878587085652E-14
+ ,0.0044642931287129
+ ,2.7089441800854E-14
+ ,0.011074310299127
+ ,0.0073153215763337
+ ,-0.0073153215763395
+ ,3.3750779948605E-14
+ ,3.3750779948605E-14
+ ,3.3750779948605E-14
+ ,3.3750779948605E-14
+ ,3.3750779948605E-14
+ ,3.3750779948605E-14
+ ,3.3750779948605E-14
+ ,3.3750779948605E-14
+ ,3.3750779948605E-14
+ ,3.3750779948605E-14
+ ,3.3750779948605E-14
+ ,3.3750779948605E-14
+ ,3.3750779948605E-14
+ ,3.3750779948605E-14
+ ,3.3750779948605E-14
+ ,3.3750779948605E-14
+ ,3.3750779948605E-14
+ ,3.3750779948605E-14
+ ,3.3750779948605E-14
+ ,3.3750779948605E-14
+ ,3.3750779948605E-14
+ ,3.3750779948605E-14
+ ,0.0073153215763337
+ ,-3.9523939676656E-14
+ ,-3.9523939676656E-14
+ ,-3.9523939676656E-14
+ ,-3.9523939676656E-14
+ ,-3.9523939676656E-14
+ ,-3.9523939676656E-14
+ ,-3.9523939676656E-14
+ ,-3.9523939676656E-14
+ ,-3.9523939676656E-14
+ ,-0.0073153215763395
+ ,3.3750779948605E-14
+ ,0.0073153215763337
+ ,-3.9523939676656E-14
+ ,-3.9523939676656E-14
+ ,-3.9523939676656E-14
+ ,-0.0073153215763395
+ ,0.026803487956034
+ ,1.0214051826551E-14
+ ,1.0214051826551E-14
+ ,0.0057020824489102
+ ,-0.0057020824488978
+ ,0.0057020824489102
+ ,2.2204460492503E-15
+ ,0.0099010709827021
+ ,-9.3258734068513E-15
+ ,-9.3258734068513E-15
+ ,-0.0056457454126093
+ ,-3.6415315207705E-14
+ ,0.0056457454125636
+ ,-9.3258734068513E-15
+ ,-9.3258734068513E-15
+ ,-0.0056457454126093
+ ,-3.6415315207705E-14
+ ,-3.6415315207705E-14
+ ,-3.6415315207705E-14
+ ,-3.6415315207705E-14
+ ,-3.6415315207705E-14
+ ,-3.6415315207705E-14
+ ,-3.6415315207705E-14
+ ,-3.6415315207705E-14
+ ,0.0056457454125636
+ ,-9.3258734068513E-15
+ ,-9.3258734068513E-15
+ ,-9.3258734068513E-15
+ ,-9.3258734068513E-15
+ ,-9.3258734068513E-15
+ ,-9.3258734068513E-15
+ ,-9.3258734068513E-15
+ ,-9.3258734068513E-15
+ ,-9.3258734068513E-15
+ ,-9.3258734068513E-15
+ ,-9.3258734068513E-15
+ ,-9.3258734068513E-15
+ ,0.0056140498327908
+ ,-0.0056140498327761
+ ,-9.3258734068513E-15
+ ,0.0056140498327908
+ ,2.3980817331903E-14
+ ,2.3980817331903E-14
+ ,-0.0056140498327761
+ ,-9.3258734068513E-15
+ ,0.0056140498327908
+ ,-0.0056140498327761
+ ,-9.3258734068513E-15
+ ,-9.3258734068513E-15
+ ,-9.3258734068513E-15
+ ,-9.3258734068513E-15
+ ,-9.3258734068513E-15
+ ,-9.3258734068513E-15
+ ,-9.3258734068513E-15
+ ,-9.3258734068513E-15
+ ,-9.3258734068513E-15
+ ,0.033221048076791
+ ,0.00068050359891414
+ ,9.3258734068513E-15
+ ,9.3258734068513E-15
+ ,9.3258734068513E-15
+ ,0.0061037829380095
+ ,-8.4376949871512E-15
+ ,-8.4376949871512E-15
+ ,-8.4376949871512E-15
+ ,-8.4376949871512E-15
+ ,-8.4376949871512E-15
+ ,-8.4376949871512E-15
+ ,-8.4376949871512E-15
+ ,-8.4376949871512E-15
+ ,-0.0067842865369085
+ ,0.0067842865369143
+ ,-0.0067842865369085
+ ,0.0067842865369143
+ ,-0.0067842865369085
+ ,1.4210854715202E-14
+ ,1.4210854715202E-14
+ ,1.4210854715202E-14
+ ,1.4210854715202E-14
+ ,0.0067842865369143
+ ,-8.4376949871512E-15
+ ,-8.4376949871512E-15
+ ,-8.4376949871512E-15
+ ,-8.4376949871512E-15
+ ,-8.4376949871512E-15
+ ,-8.4376949871512E-15
+ ,-8.4376949871512E-15
+ ,-8.4376949871512E-15
+ ,-0.0067842865369085
+ ,0.0067842865369143
+ ,-8.4376949871512E-15
+ ,-0.0067842865369085
+ ,1.4210854715202E-14
+ ,0.0067842865369143
+ ,-8.4376949871512E-15
+ ,-8.4376949871512E-15
+ ,-8.4376949871512E-15
+ ,-8.4376949871512E-15
+ ,-8.4376949871512E-15
+ ,-8.4376949871512E-15
+ ,-8.4376949871512E-15
+ ,-0.0067842865369085
+ ,1.4210854715202E-14
+ ,0.0067842865369143
+ ,-0.0067842865369085
+ ,1.4210854715202E-14
+ ,1.4210854715202E-14
+ ,1.4210854715202E-14
+ ,1.4210854715202E-14
+ ,1.4210854715202E-14
+ ,1.4210854715202E-14
+ ,0.025539022452614
+ ,-3.9968028886506E-15
+ ,-3.9968028886506E-15
+ ,-3.9968028886506E-15
+ ,-3.9968028886506E-15
+ ,-0.0053227005239038
+ ,-6.2172489379009E-15
+ ,-6.2172489379009E-15
+ ,0.019815707826494
+ ,0.0052185376095291
+ ,-0.0052185376095042
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,2.9309887850104E-14
+ ,0.0052185376095291
+ ,-0.0052185376095042
+ ,0.0039164540922294
+ ,0.005198192341894
+ ,3.0642155479654E-14
+ ,-0.0051981923418691
+ ,0.019355442952994
+ ,3.7747582837255E-14
+ ,3.7747582837255E-14
+ ,0.005098800084038
+ ,1.3766765505352E-14
+ ,1.3766765505352E-14
+ ,1.3766765505352E-14
+ ,1.3766765505352E-14
+ ,1.3766765505352E-14
+ ,1.3766765505352E-14
+ ,1.3766765505352E-14
+ ,1.3766765505352E-14
+ ,-0.0050988000839864
+ ,0.005098800084038
+ ,1.3766765505352E-14
+ ,1.3766765505352E-14
+ ,1.3766765505352E-14
+ ,1.3766765505352E-14
+ ,1.3766765505352E-14
+ ,-0.0050988000839864
+ ,3.7747582837255E-14
+ ,0.005098800084038
+ ,1.3766765505352E-14
+ ,-0.0050988000839864
+ ,3.7747582837255E-14
+ ,3.7747582837255E-14
+ ,3.7747582837255E-14
+ ,3.7747582837255E-14
+ ,3.7747582837255E-14
+ ,0.00063877357653785
+ ,4.3076653355456E-14
+ ,0.0038240964384433
+ ,3.9523939676656E-14
+ ,3.9523939676656E-14
+ ,3.9523939676656E-14
+ ,3.9523939676656E-14
+ ,3.9523939676656E-14
+ ,3.9523939676656E-14
+ ,3.9523939676656E-14
+ ,3.9523939676656E-14
+ ,3.9523939676656E-14
+ ,3.9523939676656E-14
+ ,3.9523939676656E-14
+ ,3.9523939676656E-14
+ ,3.9523939676656E-14
+ ,-0.0038240964383607
+ ,4.3076653355456E-14
+ ,4.3076653355456E-14
+ ,4.3076653355456E-14
+ ,4.3076653355456E-14
+ ,0.0038240964384433
+ ,3.9523939676656E-14
+ ,3.9523939676656E-14
+ ,3.9523939676656E-14
+ ,-0.0038240964383607
+ ,4.3076653355456E-14
+ ,4.3076653355456E-14
+ ,4.3076653355456E-14
+ ,4.3076653355456E-14
+ ,0.0038240964384433
+ ,-0.0038240964383607
+ ,4.3076653355456E-14
+ ,4.3076653355456E-14
+ ,4.3076653355456E-14
+ ,0.0038240964384433
+ ,0.0094967475372392
+ ,0.0062814276885823
+ ,-2.1316282072803E-14
+ ,-0.0062814276886214
+ ,-1.7763568394003E-14
+ ,-1.7763568394003E-14
+ ,-1.7763568394003E-14
+ ,0.0062814276885823
+ ,-2.1316282072803E-14
+ ,-2.1316282072803E-14
+ ,-0.0062814276886214
+ ,-1.7763568394003E-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] 0.000733088
> sqrtn <- sqrt(length(x))
> (armse <- sd(x) / sqrtn)
[1] 0.0001997335
> 
> (armose <- arm / armse)
[1] 3.670330
> (geo <- geomean(x))
[1] NaN
Warning message:
In log(x) : NaNs produced
> (har <- harmean(x))
[1] -8.471833e-14
> (qua <- quamean(x))
[1] 0.005062649
> (win <- winmean(x))
                [,1]         [,2]
  [1,]  7.247622e-04 1.973927e-04
  [2,]  7.148985e-04 1.949186e-04
  [3,]  7.238916e-04 1.894079e-04
  [4,]  7.234601e-04 1.893120e-04
  [5,]  7.292729e-04 1.880228e-04
  [6,]  7.172304e-04 1.854135e-04
  [7,]  7.090125e-04 1.811274e-04
  [8,]  7.015264e-04 1.796112e-04
  [9,]  6.504870e-04 1.698945e-04
 [10,]  6.463218e-04 1.683618e-04
 [11,]  5.560378e-04 1.539613e-04
 [12,]  5.525253e-04 1.534695e-04
 [13,]  5.525253e-04 1.534695e-04
 [14,]  5.525253e-04 1.534695e-04
 [15,]  5.478823e-04 1.528259e-04
 [16,]  5.147379e-04 1.484103e-04
 [17,]  4.954304e-04 1.459835e-04
 [18,]  4.751350e-04 1.435442e-04
 [19,]  4.651895e-04 1.423897e-04
 [20,]  4.474818e-04 1.386380e-04
 [21,]  4.340043e-04 1.371997e-04
 [22,]  4.242233e-04 1.361841e-04
 [23,]  4.254948e-04 1.360766e-04
 [24,]  4.259779e-04 1.360358e-04
 [25,]  3.831241e-04 1.318279e-04
 [26,]  3.831241e-04 1.318279e-04
 [27,]  3.479204e-04 1.286650e-04
 [28,]  3.479204e-04 1.286650e-04
 [29,]  3.479204e-04 1.286650e-04
 [30,]  3.384988e-04 1.278624e-04
 [31,]  3.632427e-04 1.257267e-04
 [32,]  3.632427e-04 1.257267e-04
 [33,]  3.638615e-04 1.256747e-04
 [34,]  3.638615e-04 1.256747e-04
 [35,]  3.702919e-04 1.251384e-04
 [36,]  3.719075e-04 1.250047e-04
 [37,]  3.831166e-04 1.240879e-04
 [38,]  3.831166e-04 1.240879e-04
 [39,]  3.831166e-04 1.240879e-04
 [40,]  3.741100e-04 1.233160e-04
 [41,]  3.740535e-04 1.233112e-04
 [42,]  3.784920e-04 1.213764e-04
 [43,]  3.784920e-04 1.213764e-04
 [44,]  3.801070e-04 1.212474e-04
 [45,]  3.816435e-04 1.207169e-04
 [46,]  3.816435e-04 1.207169e-04
 [47,]  3.830620e-04 1.204495e-04
 [48,]  3.830620e-04 1.204495e-04
 [49,]  3.830620e-04 1.204495e-04
 [50,]  4.061850e-04 1.186466e-04
 [51,]  4.110041e-04 1.182779e-04
 [52,]  4.110041e-04 1.182779e-04
 [53,]  3.687001e-04 1.147317e-04
 [54,]  3.687001e-04 1.147317e-04
 [55,]  3.676687e-04 1.146472e-04
 [56,]  3.676687e-04 1.146472e-04
 [57,]  3.606439e-04 1.140743e-04
 [58,]  3.571358e-04 1.137897e-04
 [59,]  3.586680e-04 1.132475e-04
 [60,]  3.565908e-04 1.130795e-04
 [61,]  3.421926e-04 1.116113e-04
 [62,]  3.519741e-04 1.108429e-04
 [63,]  3.519741e-04 1.108429e-04
 [64,]  3.519741e-04 1.108429e-04
 [65,]  4.834911e-04 1.012215e-04
 [66,]  4.834911e-04 1.012215e-04
 [67,]  4.612079e-04 9.932868e-05
 [68,]  8.739675e-04 8.050212e-05
 [69,]  8.739675e-04 8.050212e-05
 [70,]  8.713982e-04 8.025611e-05
 [71,]  8.713982e-04 8.025611e-05
 [72,]  8.649597e-04 7.964081e-05
 [73,]  8.649597e-04 7.964081e-05
 [74,]  8.612367e-04 7.928575e-05
 [75,]  8.612367e-04 7.928575e-05
 [76,]  8.612367e-04 7.928575e-05
 [77,]  8.183514e-04 7.522722e-05
 [78,]  8.183514e-04 7.522722e-05
 [79,]  8.183514e-04 7.522722e-05
 [80,]  8.183514e-04 7.522722e-05
 [81,]  8.183514e-04 7.522722e-05
 [82,]  8.183514e-04 7.522722e-05
 [83,]  8.183514e-04 7.522722e-05
 [84,]  8.183514e-04 7.522722e-05
 [85,]  8.183514e-04 7.522722e-05
 [86,]  8.122587e-04 7.465521e-05
 [87,]  8.122587e-04 7.465521e-05
 [88,]  8.094168e-04 7.438860e-05
 [89,]  7.995313e-04 7.346195e-05
 [90,]  7.953290e-04 7.306837e-05
 [91,]  7.953290e-04 7.306837e-05
 [92,]  7.953290e-04 7.306837e-05
 [93,]  7.016637e-04 6.432614e-05
 [94,]  6.199226e-04 5.674933e-05
 [95,]  6.059957e-04 5.546389e-05
 [96,]  6.059957e-04 5.546389e-05
 [97,]  6.059957e-04 5.546389e-05
 [98,]  6.059957e-04 5.546389e-05
 [99,]  1.610393e-04 1.454415e-05
[100,]  1.431565e-04 1.292076e-05
[101,]  1.111907e-04 1.002858e-05
[102,]  1.044344e-04 9.419004e-06
[103,]  5.286776e-17 1.394668e-15
[104,]  5.286776e-17 1.394668e-15
[105,]  5.286776e-17 1.394668e-15
[106,]  5.286776e-17 1.394668e-15
[107,]  5.286776e-17 1.394668e-15
[108,]  5.286776e-17 1.394668e-15
[109,]  5.286776e-17 1.394668e-15
[110,]  5.286776e-17 1.394668e-15
[111,]  5.286776e-17 1.394668e-15
[112,]  5.286776e-17 1.394668e-15
[113,]  5.286776e-17 1.394668e-15
[114,]  5.286776e-17 1.394668e-15
[115,] -4.335157e-16 1.371220e-15
[116,]  1.388660e-16 1.339566e-15
[117,]  3.862871e-16 1.326186e-15
[118,]  3.862871e-16 1.326186e-15
[119,]  2.185201e-16 1.318184e-15
[120,]  2.185201e-16 1.318184e-15
[121,]  2.185201e-16 1.318184e-15
[122,]  2.185201e-16 1.318184e-15
[123,]  2.185201e-16 1.318184e-15
[124,]  2.185201e-16 1.318184e-15
[125,]  2.185201e-16 1.318184e-15
[126,]  2.185201e-16 1.318184e-15
[127,]  2.185201e-16 1.318184e-15
[128,]  2.185201e-16 1.318184e-15
[129,]  2.185201e-16 1.318184e-15
[130,]  2.185201e-16 1.318184e-15
[131,]  2.185201e-16 1.318184e-15
[132,]  2.185201e-16 1.318184e-15
[133,]  2.185201e-16 1.318184e-15
[134,]  7.852625e-16 1.288053e-15
[135,]  4.046146e-16 1.270054e-15
[136,]  4.046146e-16 1.270054e-15
[137,]  4.046146e-16 1.270054e-15
[138,]  8.909980e-16 1.244666e-15
[139,]  8.909980e-16 1.244666e-15
[140,]  8.909980e-16 1.244666e-15
[141,]  8.909980e-16 1.244666e-15
[142,]  8.909980e-16 1.244666e-15
[143,]  3.869920e-16 1.221272e-15
[144,]  3.869920e-16 1.221272e-15
[145,]  3.869920e-16 1.221272e-15
[146,]  3.869920e-16 1.221272e-15
[147,]  3.869920e-16 1.221272e-15
[148,]  3.869920e-16 1.221272e-15
[149,]  3.869920e-16 1.221272e-15
[150,]  1.021405e-15 1.188715e-15
[151,]  1.021405e-15 1.188715e-15
[152,]  1.021405e-15 1.188715e-15
[153,]  1.021405e-15 1.188715e-15
[154,]  1.021405e-15 1.188715e-15
[155,]  1.021405e-15 1.188715e-15
[156,]  1.021405e-15 1.188715e-15
[157,]  5.787258e-16 1.168509e-15
[158,]  5.787258e-16 1.168509e-15
[159,]  5.787258e-16 1.168509e-15
[160,]  5.787258e-16 1.168509e-15
[161,]  5.787258e-16 1.168509e-15
[162,]  5.787258e-16 1.168509e-15
[163,]  5.787258e-16 1.168509e-15
[164,]  5.787258e-16 1.168509e-15
[165,]  5.787258e-16 1.168509e-15
[166,]  5.787258e-16 1.168509e-15
[167,]  6.964447e-16 1.162523e-15
[168,]  6.964447e-16 1.162523e-15
[169,]  6.964447e-16 1.162523e-15
[170,]  6.964447e-16 1.162523e-15
[171,]  6.964447e-16 1.162523e-15
[172,]  6.964447e-16 1.162523e-15
[173,]  6.964447e-16 1.162523e-15
[174,]  6.964447e-16 1.162523e-15
[175,]  6.964447e-16 1.162523e-15
[176,]  2.805516e-15 1.058691e-15
[177,]  2.805516e-15 1.058691e-15
[178,]  2.805516e-15 1.058691e-15
[179,]  2.805516e-15 1.058691e-15
[180,]  2.297985e-15 1.035332e-15
[181,]  2.297985e-15 1.035332e-15
[182,]  3.837495e-15 9.636584e-16
[183,]  3.450503e-15 9.458312e-16
[184,]  3.061396e-15 9.281907e-16
[185,]  3.061396e-15 9.281907e-16
[186,]  3.454732e-15 9.104172e-16
[187,]  3.454732e-15 9.104172e-16
[188,]  3.454732e-15 9.104172e-16
[189,]  3.454732e-15 9.104172e-16
[190,]  3.454732e-15 9.104172e-16
[191,]  3.454732e-15 9.104172e-16
[192,]  3.454732e-15 9.104172e-16
[193,]  3.454732e-15 9.104172e-16
[194,]  3.454732e-15 9.104172e-16
[195,]  3.454732e-15 9.104172e-16
[196,]  3.731054e-15 8.981002e-16
[197,]  3.731054e-15 8.981002e-16
[198,]  3.731054e-15 8.981002e-16
[199,]  4.151882e-15 8.795898e-16
[200,]  4.151882e-15 8.795898e-16
[201,]  4.151882e-15 8.795898e-16
[202,]  4.294272e-15 8.733951e-16
[203,]  4.294272e-15 8.733951e-16
[204,]  4.294272e-15 8.733951e-16
[205,]  4.294272e-15 8.733951e-16
[206,]  4.294272e-15 8.733951e-16
[207,]  4.294272e-15 8.733951e-16
[208,]  4.294272e-15 8.733951e-16
[209,]  4.294272e-15 8.733951e-16
[210,]  4.294272e-15 8.733951e-16
> (tri <- trimean(x))
               [,1]         [,2]
  [1,] 6.964640e-04 1.896809e-04
  [2,] 6.679849e-04 1.815273e-04
  [3,] 6.443026e-04 1.742962e-04
  [4,] 6.174317e-04 1.688010e-04
  [5,] 5.904971e-04 1.630657e-04
  [6,] 5.622030e-04 1.573527e-04
  [7,] 5.357779e-04 1.518995e-04
  [8,] 5.103852e-04 1.469669e-04
  [9,] 4.857898e-04 1.420432e-04
 [10,] 4.668901e-04 1.383427e-04
 [11,] 4.482977e-04 1.346913e-04
 [12,] 4.381152e-04 1.326305e-04
 [13,] 4.281706e-04 1.305593e-04
 [14,] 4.181599e-04 1.284236e-04
 [15,] 4.080825e-04 1.262196e-04
 [16,] 3.982638e-04 1.240027e-04
 [17,] 3.905689e-04 1.220973e-04
 [18,] 3.840268e-04 1.203317e-04
 [19,] 3.786403e-04 1.187025e-04
 [20,] 3.737762e-04 1.171115e-04
 [21,] 3.698277e-04 1.157377e-04
 [22,] 3.665422e-04 1.144228e-04
 [23,] 3.637139e-04 1.131381e-04
 [24,] 3.608062e-04 1.118258e-04
 [25,] 3.578566e-04 1.104803e-04
 [26,] 3.567550e-04 1.093444e-04
 [27,] 3.556457e-04 1.081791e-04
 [28,] 3.559597e-04 1.071551e-04
 [29,] 3.562760e-04 1.061053e-04
 [30,] 3.565944e-04 1.050287e-04
 [31,] 3.572635e-04 1.039646e-04
 [32,] 3.570488e-04 1.029782e-04
 [33,] 3.568326e-04 1.019667e-04
 [34,] 3.565938e-04 1.009315e-04
 [35,] 3.563533e-04 9.986913e-05
 [36,] 3.559037e-04 9.880329e-05
 [37,] 3.554000e-04 9.771496e-05
 [38,] 3.545481e-04 9.663818e-05
 [39,] 3.536901e-04 9.553203e-05
 [40,] 3.528258e-04 9.439517e-05
 [41,] 3.522140e-04 9.326016e-05
 [42,] 3.515994e-04 9.209307e-05
 [43,] 3.515994e-04 9.097524e-05
 [44,] 3.501109e-04 8.982535e-05
 [45,] 3.493156e-04 8.864726e-05
 [46,] 3.484743e-04 8.745642e-05
 [47,] 3.476268e-04 8.622967e-05
 [48,] 3.467373e-04 8.497630e-05
 [49,] 3.458411e-04 8.368355e-05
 [50,] 3.449382e-04 8.234922e-05
 [51,] 3.434766e-04 8.104598e-05
 [52,] 3.418908e-04 7.971503e-05
 [53,] 3.402928e-04 7.833928e-05
 [54,] 3.396459e-04 7.706467e-05
 [55,] 3.389941e-04 7.574725e-05
 [56,] 3.383600e-04 7.438787e-05
 [57,] 3.377210e-04 7.298025e-05
 [58,] 3.372281e-04 7.154538e-05
 [59,] 3.368057e-04 7.006891e-05
 [60,] 3.363480e-04 6.855871e-05
 [61,] 3.359296e-04 6.699597e-05
 [62,] 3.358018e-04 6.543332e-05
 [63,] 3.354757e-04 6.383920e-05
 [64,] 3.351470e-04 6.217665e-05
 [65,] 3.348158e-04 6.043952e-05
 [66,] 3.319222e-04 5.904036e-05
 [67,] 3.290052e-04 5.758242e-05
 [68,] 3.264889e-04 5.614504e-05
 [69,] 3.161794e-04 5.528709e-05
 [70,] 3.057859e-04 5.439621e-05
 [71,] 2.953545e-04 5.348017e-05
 [72,] 2.848373e-04 5.252737e-05
 [73,] 2.743495e-04 5.156050e-05
 [74,] 2.637747e-04 5.055307e-05
 [75,] 2.531779e-04 4.951704e-05
 [76,] 2.424923e-04 4.843506e-05
 [77,] 2.317170e-04 4.730341e-05
 [78,] 2.215909e-04 4.629454e-05
 [79,] 2.113791e-04 4.523970e-05
 [80,] 2.010804e-04 4.413501e-05
 [81,] 2.010804e-04 4.297603e-05
 [82,] 1.802176e-04 4.175762e-05
 [83,] 1.696514e-04 4.047374e-05
 [84,] 1.589937e-04 3.911724e-05
 [85,] 1.482433e-04 3.767952e-05
 [86,] 1.482433e-04 3.615007e-05
 [87,] 1.265575e-04 3.454770e-05
 [88,] 1.156204e-04 3.282804e-05
 [89,] 1.046316e-04 3.098699e-05
 [90,] 9.370058e-05 2.904367e-05
 [91,] 9.370058e-05 2.693853e-05
 [92,] 7.167637e-05 2.459941e-05
 [93,] 6.051545e-05 2.195014e-05
 [94,] 5.068907e-05 1.965842e-05
 [95,] 4.201844e-05 1.767562e-05
 [96,] 3.347949e-05 1.550428e-05
 [97,] 2.486220e-05 1.290962e-05
 [98,] 1.616549e-05 9.550513e-06
 [99,] 7.388258e-06 3.741822e-06
[100,] 5.114353e-06 2.976283e-06
[101,] 3.082423e-06 2.178143e-06
[102,] 1.499468e-06 1.499468e-06
[103,] 1.480996e-15 1.362479e-15
[104,] 1.501695e-15 1.360151e-15
[105,] 1.522592e-15 1.357742e-15
[106,] 1.543688e-15 1.355250e-15
[107,] 1.564987e-15 1.352673e-15
[108,] 1.586493e-15 1.350008e-15
[109,] 1.608207e-15 1.347253e-15
[110,] 1.630132e-15 1.344405e-15
[111,] 1.652273e-15 1.341462e-15
[112,] 1.674632e-15 1.338419e-15
[113,] 1.697212e-15 1.335275e-15
[114,] 1.720017e-15 1.332026e-15
[115,] 1.743050e-15 1.328670e-15
[116,] 1.773009e-15 1.325710e-15
[117,] 1.795421e-15 1.323385e-15
[118,] 1.814679e-15 1.321278e-15
[119,] 1.834134e-15 1.319090e-15
[120,] 1.856065e-15 1.316986e-15
[121,] 1.878223e-15 1.314800e-15
[122,] 1.900610e-15 1.312527e-15
[123,] 1.923230e-15 1.310166e-15
[124,] 1.946087e-15 1.307714e-15
[125,] 1.969185e-15 1.305167e-15
[126,] 1.992527e-15 1.302523e-15
[127,] 2.016118e-15 1.299779e-15
[128,] 2.039961e-15 1.296931e-15
[129,] 2.064060e-15 1.293977e-15
[130,] 2.088420e-15 1.290911e-15
[131,] 2.113044e-15 1.287732e-15
[132,] 2.137938e-15 1.284434e-15
[133,] 2.163105e-15 1.281015e-15
[134,] 2.188550e-15 1.277469e-15
[135,] 2.206877e-15 1.274525e-15
[136,] 2.230370e-15 1.271865e-15
[137,] 2.254127e-15 1.269098e-15
[138,] 2.278153e-15 1.266219e-15
[139,] 2.296143e-15 1.263836e-15
[140,] 2.314339e-15 1.261351e-15
[141,] 2.332744e-15 1.258761e-15
[142,] 2.351363e-15 1.256062e-15
[143,] 2.370197e-15 1.253250e-15
[144,] 2.395744e-15 1.250843e-15
[145,] 2.421592e-15 1.248329e-15
[146,] 2.447746e-15 1.245704e-15
[147,] 2.474211e-15 1.242965e-15
[148,] 2.500993e-15 1.240106e-15
[149,] 2.528098e-15 1.237125e-15
[150,] 2.555532e-15 1.234015e-15
[151,] 2.575176e-15 1.231596e-15
[152,] 2.595061e-15 1.229064e-15
[153,] 2.615192e-15 1.226415e-15
[154,] 2.635573e-15 1.223646e-15
[155,] 2.656209e-15 1.220752e-15
[156,] 2.677104e-15 1.217727e-15
[157,] 2.698264e-15 1.214567e-15
[158,] 2.725350e-15 1.211742e-15
[159,] 2.752784e-15 1.208786e-15
[160,] 2.780571e-15 1.205692e-15
[161,] 2.808720e-15 1.202456e-15
[162,] 2.808720e-15 1.199072e-15
[163,] 2.866128e-15 1.195533e-15
[164,] 2.895403e-15 1.191834e-15
[165,] 2.925068e-15 1.187967e-15
[166,] 2.955131e-15 1.183927e-15
[167,] 2.985600e-15 1.179705e-15
[168,] 3.014973e-15 1.175467e-15
[169,] 3.044749e-15 1.171038e-15
[170,] 3.074935e-15 1.166409e-15
[171,] 3.105541e-15 1.161572e-15
[172,] 3.105541e-15 1.156518e-15
[173,] 3.168045e-15 1.151235e-15
[174,] 3.199962e-15 1.145715e-15
[175,] 3.232335e-15 1.139945e-15
[176,] 3.265174e-15 1.133915e-15
[177,] 3.271135e-15 1.130775e-15
[178,] 3.277184e-15 1.127472e-15
[179,] 3.283321e-15 1.124000e-15
[180,] 3.289550e-15 1.120351e-15
[181,] 3.302499e-15 1.117173e-15
[182,] 3.302499e-15 1.113825e-15
[183,] 3.308801e-15 1.112339e-15
[184,] 3.306939e-15 1.111230e-15
[185,] 3.310173e-15 1.110494e-15
[186,] 3.313456e-15 1.109668e-15
[187,] 3.311587e-15 1.109229e-15
[188,] 3.309688e-15 1.108706e-15
[189,] 3.307760e-15 1.108097e-15
[190,] 3.305800e-15 1.107396e-15
[191,] 3.303809e-15 1.106600e-15
[192,] 3.301785e-15 1.105703e-15
[193,] 3.299728e-15 1.104701e-15
[194,] 3.297638e-15 1.103589e-15
[195,] 3.295512e-15 1.102361e-15
[196,] 3.293351e-15 1.101012e-15
[197,] 3.287389e-15 1.099900e-15
[198,] 3.281326e-15 1.098667e-15
[199,] 3.275158e-15 1.097308e-15
[200,] 3.263090e-15 1.096367e-15
[201,] 3.263090e-15 1.095307e-15
[202,] 3.238314e-15 1.094121e-15
[203,] 3.223612e-15 1.092990e-15
[204,] 3.208645e-15 1.091725e-15
[205,] 3.193405e-15 1.090321e-15
[206,] 3.177886e-15 1.088770e-15
[207,] 3.162080e-15 1.087063e-15
[208,] 3.145978e-15 1.085193e-15
[209,] 3.129572e-15 1.083149e-15
[210,] 3.112854e-15 1.080921e-15
> 
> (midr <- midrange(x))
[1] 0.01223303
> midm <- array(NA,dim=8)
> for (j in 1:8) midm[j] <- midmean(x,j)
> midm
[1] 1.840977e-15 1.840977e-15 1.840977e-15 1.840977e-15 1.840977e-15
[6] 1.840977e-15 1.840977e-15 1.840977e-15
> postscript(file="/var/www/rcomp/tmp/15qkn1288122667.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/rcomp/tmp/2gi1q1288122667.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/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/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/rcomp/tmp/3edb81288122668.tab") 
> try(system("convert tmp/15qkn1288122667.ps tmp/15qkn1288122667.png",intern=TRUE))
character(0)
> try(system("convert tmp/2gi1q1288122667.ps tmp/2gi1q1288122667.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
   5.80    0.44    6.23