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(-4.5741188614556E-14 + ,-4.5741188614556E-14 + ,-4.5741188614556E-14 + ,0.0039920212695539 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,-0.015052968628983 + ,-2.0872192862953E-14 + ,0.015052968628979 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,-0.015052968628983 + ,-2.0872192862953E-14 + ,-2.0872192862953E-14 + ,0.0030287756621794 + ,1.8651746813703E-14 + ,-0.010131798930381 + ,2.5757174171304E-14 + ,2.5757174171304E-14 + ,2.5757174171304E-14 + ,-0.0020387366898742 + ,-0.0061412680221262 + ,0.0061412680220565 + ,-0.0061412680221262 + ,-4.3520742565306E-14 + ,0.0061412680220565 + ,-0.0010209291341265 + ,-1.9095836023553E-14 + ,-1.9095836023553E-14 + ,-1.9095836023553E-14 + ,-1.9095836023553E-14 + ,-1.9095836023553E-14 + ,-0.0082051742392188 + ,-3.3306690738755E-14 + ,0.0082051742391664 + ,-1.9095836023553E-14 + ,-1.9095836023553E-14 + ,-1.9095836023553E-14 + ,-1.9095836023553E-14 + ,-1.9095836023553E-14 + ,-0.0082051742392188 + ,-3.3306690738755E-14 + ,-3.3306690738755E-14 + ,-3.3306690738755E-14 + ,-3.3306690738755E-14 + ,-0.0020618564005335 + ,2.5313084961454E-14 + ,-0.0062112000926748 + ,-3.4638958368305E-14 + ,-3.4638958368305E-14 + ,0.0062112000926655 + ,2.5313084961454E-14 + ,-0.0062112000926748 + ,0.0062112000926655 + ,2.5313084961454E-14 + ,-0.0062112000926748 + ,-3.4638958368305E-14 + ,-3.4638958368305E-14 + ,0.0062112000926655 + ,2.5313084961454E-14 + ,-0.0062112000926748 + ,-3.4638958368305E-14 + ,0.0062112000926655 + ,-0.0062112000926748 + ,-3.4638958368305E-14 + ,-3.4638958368305E-14 + ,-3.4638958368305E-14 + ,-3.4638958368305E-14 + ,0.0062112000926655 + ,-0.0062112000926748 + ,-3.4638958368305E-14 + ,0.0062112000926655 + ,-0.0062112000926748 + ,-3.4638958368305E-14 + ,-3.4638958368305E-14 + ,-3.4638958368305E-14 + ,-3.4638958368305E-14 + ,-3.4638958368305E-14 + ,-3.4638958368305E-14 + ,-3.4638958368305E-14 + ,-3.4638958368305E-14 + ,-3.4638958368305E-14 + ,-3.4638958368305E-14 + ,0.0062112000926655 + ,-0.0062112000926748 + ,-3.4638958368305E-14 + ,-3.4638958368305E-14 + ,-3.4638958368305E-14 + ,-3.4638958368305E-14 + ,-3.4638958368305E-14 + ,-3.4638958368305E-14 + ,-3.4638958368305E-14 + ,0.0062112000926655 + ,2.5313084961454E-14 + ,2.5313084961454E-14 + ,0.010267030639725 + ,-0.010267030639719 + ,2.5313084961454E-14 + ,0.010267030639725 + ,-1.9095836023553E-14 + ,-1.9095836023553E-14 + ,-1.9095836023553E-14 + ,-1.9095836023553E-14 + ,-0.010267030639719 + ,0.010267030639725 + ,-1.9095836023553E-14 + ,-1.9095836023553E-14 + ,-1.9095836023553E-14 + ,-1.9095836023553E-14 + ,-1.9095836023553E-14 + ,-1.9095836023553E-14 + ,-1.9095836023553E-14 + ,-0.010267030639719 + ,2.5313084961454E-14 + ,2.5313084961454E-14 + ,2.5313084961454E-14 + ,2.5313084961454E-14 + ,0.010267030639725 + ,-1.9095836023553E-14 + ,-1.9095836023553E-14 + ,-1.9095836023553E-14 + ,-1.9095836023553E-14 + ,-0.010267030639719 + ,0.010267030639725 + ,-1.9095836023553E-14 + ,-1.9095836023553E-14 + ,-0.0020449904877191 + ,-0.0061601837514913 + ,-3.3306690738755E-14 + ,-3.3306690738755E-14 + ,-3.3306690738755E-14 + ,-3.3306690738755E-14 + ,-3.3306690738755E-14 + ,0.0061601837514664 + ,0.0010230179920088 + ,-2.5757174171304E-14 + ,-2.5757174171304E-14 + ,-2.5757174171304E-14 + ,-2.5757174171304E-14 + ,-2.5757174171304E-14 + ,-2.5757174171304E-14 + ,-2.5757174171304E-14 + ,0.0091603693986739 + ,9.7699626167014E-15 + ,-0.0091603693986899 + ,0.0091603693986739 + ,-0.0091603693986899 + ,-2.5757174171304E-14 + ,0.0091603693986739 + ,9.7699626167014E-15 + ,9.7699626167014E-15 + ,9.7699626167014E-15 + ,-0.0091603693986899 + ,0.0091603693986739 + ,9.7699626167014E-15 + ,9.7699626167014E-15 + ,-0.0091603693986899 + ,-2.5757174171304E-14 + ,-2.5757174171304E-14 + ,-2.5757174171304E-14 + ,-2.5757174171304E-14 + ,-2.5757174171304E-14 + ,0.0091603693986739 + ,9.7699626167014E-15 + ,9.7699626167014E-15 + ,9.7699626167014E-15 + ,9.7699626167014E-15 + ,9.7699626167014E-15 + ,9.7699626167014E-15 + ,9.7699626167014E-15 + ,9.7699626167014E-15 + ,-0.0091603693986899 + ,-2.5757174171304E-14 + ,-2.5757174171304E-14 + ,-2.5757174171304E-14 + ,-0.0051255878488261 + ,-1.3766765505352E-14 + ,-1.3766765505352E-14 + ,-1.3766765505352E-14 + ,-0.010330670387914 + ,-3.4638958368305E-14 + ,-3.4638958368305E-14 + ,-3.4638958368305E-14 + ,0.010330670387865 + ,0.021353124470586 + ,1.7319479184152E-14 + ,0.010010093595117 + ,1.6875389974302E-14 + ,-0.010010093595083 + ,1.7319479184152E-14 + ,1.7319479184152E-14 + ,1.7319479184152E-14 + ,1.7319479184152E-14 + ,0.010010093595117 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,-0.010010093595083 + ,1.7319479184152E-14 + ,0.010010093595117 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,-0.010010093595083 + ,0.010010093595117 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,-0.010010093595083 + ,1.7319479184152E-14 + ,1.7319479184152E-14 + ,1.7319479184152E-14 + ,1.7319479184152E-14 + ,1.7319479184152E-14 + ,1.7319479184152E-14 + ,0.010010093595117 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,-0.010010093595083 + ,1.7319479184152E-14 + ,1.7319479184152E-14 + ,1.7319479184152E-14 + ,1.7319479184152E-14 + ,0.010010093595117 + ,1.6875389974302E-14 + ,-0.010010093595083 + ,0.010010093595117 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,-0.010010093595083 + ,1.7319479184152E-14 + ,1.7319479184152E-14 + ,-0.0010065426113823 + ,0.0060241146034188 + ,3.7747582837255E-14 + ,-0.0060241146033619 + ,0.0060241146034188 + ,3.7747582837255E-14 + ,3.7747582837255E-14 + ,3.7747582837255E-14 + ,3.7747582837255E-14 + ,3.7747582837255E-14 + ,-0.0060241146033619 + ,1.9095836023553E-14 + ,1.9095836023553E-14 + ,1.9095836023553E-14 + ,1.9095836023553E-14 + ,-0.0080890228730812 + ,0.0080890228731025 + ,1.9095836023553E-14 + ,-0.0080890228730812 + ,2.2204460492503E-15 + ,2.2204460492503E-15 + ,2.2204460492503E-15 + ,0.0080890228731025 + ,-0.0080890228730812 + ,2.2204460492503E-15 + ,2.2204460492503E-15 + ,0.0080890228731025 + ,1.9095836023553E-14 + ,1.9095836023553E-14 + ,1.9095836023553E-14 + ,-0.0060606246116812 + ,9.7699626167014E-15 + ,9.7699626167014E-15 + ,9.7699626167014E-15 + ,9.7699626167014E-15 + ,-0.0081383969029902 + ,0.0081383969029809 + ,-0.0081383969029902 + ,-1.9095836023553E-14 + ,0.0081383969029809 + ,-0.0081383969029902 + ,0.0081383969029809 + ,9.7699626167014E-15 + ,9.7699626167014E-15 + ,9.7699626167014E-15 + ,9.7699626167014E-15 + ,-0.0081383969029902 + ,-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 + ,-1.9095836023553E-14 + ,0.0081383969029809 + ,9.7699626167014E-15 + ,9.7699626167014E-15 + ,9.7699626167014E-15 + ,9.7699626167014E-15 + ,9.7699626167014E-15 + ,-0.0081383969029902 + ,-1.9095836023553E-14 + ,0.0081383969029809 + ,9.7699626167014E-15 + ,-0.0081383969029902 + ,-1.9095836023553E-14 + ,0.0020408170348811 + ,2.8421709430404E-14 + ,2.8421709430404E-14 + ,2.8421709430404E-14 + ,2.8421709430404E-14 + ,2.8421709430404E-14 + ,2.8421709430404E-14 + ,2.8421709430404E-14 + ,2.8421709430404E-14 + ,-0.0081883773793718 + ,0.0081883773793865 + ,2.8421709430404E-14 + ,2.8421709430404E-14 + ,-0.0081883773793718 + ,-1.3766765505352E-14 + ,0.0081883773793865 + ,-0.0081883773793718 + ,-1.3766765505352E-14 + ,-1.3766765505352E-14 + ,-1.3766765505352E-14 + ,-1.3766765505352E-14 + ,0.0081883773793865 + ,2.8421709430404E-14 + ,2.8421709430404E-14 + ,-0.0081883773793718 + ,-1.3766765505352E-14 + ,-1.3766765505352E-14 + ,-1.3766765505352E-14 + ,-1.3766765505352E-14 + ,-1.3766765505352E-14 + ,-1.3766765505352E-14 + ,0.0081883773793865 + ,2.8421709430404E-14 + ,-0.0081883773793718 + ,-1.3766765505352E-14 + ,-1.3766765505352E-14 + ,-1.3766765505352E-14 + ,-1.3766765505352E-14 + ,-1.3766765505352E-14 + ,0.0081883773793865 + ,-0.0010198879007715 + ,-0.0081967672042262 + ,-4.7517545453957E-14 + ,0.0081967672041525 + ,-2.6201263381154E-14 + ,-2.6201263381154E-14 + ,-2.6201263381154E-14 + ,-2.6201263381154E-14 + ,-2.6201263381154E-14 + ,-2.6201263381154E-14 + ,-2.6201263381154E-14 + ,-2.6201263381154E-14 + ,-2.6201263381154E-14 + ,-2.6201263381154E-14 + ,-2.6201263381154E-14 + ,0.0010198879007737 + ,2.8421709430404E-14 + ,2.8421709430404E-14 + ,-0.0061349885674717 + ,0.0051151006667443 + ,0.010152371463974 + ,-0.010152371464044 + ,-2.6201263381154E-14 + ,-2.6201263381154E-14 + ,-2.6201263381154E-14 + ,0.010152371463974 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,-0.010152371464044 + ,0.010152371463974 + ,-4.4408920985006E-14 + ,-0.010152371464044 + ,-2.6201263381154E-14 + ,-2.6201263381154E-14 + ,-2.6201263381154E-14 + ,-2.6201263381154E-14 + ,-2.6201263381154E-14 + ,-2.6201263381154E-14 + ,-2.6201263381154E-14 + ,-2.6201263381154E-14 + ,-2.6201263381154E-14 + ,0.0071174677688735 + ,9.7699626167014E-15 + ,9.7699626167014E-15 + ,0.0060606246117101 + ,-0.0060606246116812 + ,0.0060606246117101 + ,1.9095836023553E-14 + ,1.9095836023553E-14 + ,1.9095836023553E-14 + ,1.9095836023553E-14 + ,1.9095836023553E-14 + ,-0.0060606246116812 + ,0.00505306785141 + ,0.010030174359918 + ,-0.010030174359918 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,0.010030174359918 + ,-0.010030174359918 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,0.010030174359918 + ,-0.010030174359918 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,0.010030174359918 + ,-1.8651746813703E-14 + ,-1.8651746813703E-14 + ,-1.8651746813703E-14 + ,-1.8651746813703E-14 + ,-1.8651746813703E-14 + ,-1.8651746813703E-14 + ,-0.010030174359918 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,0.010030174359918 + ,-1.8651746813703E-14 + ,-1.8651746813703E-14 + ,-1.8651746813703E-14 + ,-1.8651746813703E-14 + ,-1.8651746813703E-14 + ,-1.8651746813703E-14 + ,-0.010030174359918 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,0.010030174359918 + ,-1.8651746813703E-14 + ,-0.010030174359918 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,1.8651746813703E-14 + ,0.010030174359918 + ,-1.8651746813703E-14 + ,-1.8651746813703E-14 + ,-0.010030174359918 + ,0.010030174359918 + ,-1.8651746813703E-14 + ,-1.8651746813703E-14 + ,-0.010030174359918 + ,0.010030174359918 + ,-1.8651746813703E-14 + ,-1.8651746813703E-14 + ,-1.8651746813703E-14 + ,-1.8651746813703E-14 + ,-1.8651746813703E-14 + ,-1.8651746813703E-14 + ,-1.8651746813703E-14 + ,-0.010030174359918 + ,0.010030174359918 + ,-1.8651746813703E-14 + ,-1.8651746813703E-14 + ,-1.8651746813703E-14 + ,-1.8651746813703E-14 + ,-1.8651746813703E-14 + ,0.00099750631708151 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,0.0079444310585561 + ,2.0428103653103E-14 + ,2.0428103653103E-14 + ,2.0428103653103E-14 + ,2.0428103653103E-14 + ,-0.0079444310585797 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,0.0079444310585561 + ,-0.0079444310585797 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,-4.3964831775156E-14 + ,0.0029865626977559 + ,-0.0089865806978935 + ,-4.7073456244107E-14 + ,-4.7073456244107E-14 + ,-4.7073456244107E-14 + ,0.0089865806978531 + ,6.6613381477509E-15 + ,6.6613381477509E-15 + ,6.6613381477509E-15 + ,6.6613381477509E-15 + ,-0.0089865806978935 + ,0.0089865806978531 + ,6.6613381477509E-15 + ,-0.0089865806978935 + ,-4.7073456244107E-14 + ,0.0089865806978531 + ,-0.0089865806978935 + ,-4.7073456244107E-14 + ,0.0030045090202533 + ,-4.5741188614556E-14 + ,-4.5741188614556E-14 + ,-4.5741188614556E-14 + ,-4.5741188614556E-14 + ,-4.5741188614556E-14 + ,-4.5741188614556E-14 + ,-4.5741188614556E-14 + ,-4.5741188614556E-14 + ,-4.5741188614556E-14 + ,-4.5741188614556E-14 + ,-0.010050335853546 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,0.010050335853456 + ,-0.010050335853546 + ,-4.4408920985006E-14 + ,0.010050335853456 + ,-4.5741188614556E-14 + ,-4.5741188614556E-14 + ,-4.5741188614556E-14 + ,-4.5741188614556E-14 + ,-4.5741188614556E-14 + ,-4.5741188614556E-14 + ,-4.5741188614556E-14 + ,-4.5741188614556E-14 + ,-4.5741188614556E-14 + ,-0.010050335853546 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,0.010050335853456 + ,-0.010050335853546 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,0.010050335853456 + ,-4.5741188614556E-14 + ,-4.5741188614556E-14 + ,-4.5741188614556E-14 + ,-0.010050335853546 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,0.010050335853456 + ,-4.5741188614556E-14 + ,-4.5741188614556E-14 + ,-0.010050335853546 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,0.010050335853456 + ,-4.5741188614556E-14 + ,-4.5741188614556E-14 + ,-4.5741188614556E-14 + ,-0.010050335853546 + ,0.021978906718755 + ,-1.9539925233403E-14 + ,-0.0099305682026194 + ,-1.8651746813703E-14 + ,-1.8651746813703E-14 + ,-1.8651746813703E-14 + ,0.0099305682025812 + ,-1.9539925233403E-14 + ,0.0078740564308806 + ,-0.005899722127225 + ,0.0058997221271628 + ,-2.5313084961454E-14 + ,-2.5313084961454E-14 + ,-2.5313084961454E-14 + ,-0.005899722127225 + ,-3.6859404417555E-14 + ,0.0058997221271628 + ,-2.5313084961454E-14 + ,-2.5313084961454E-14 + ,-0.005899722127225 + ,-0.0019743343037373 + ,-1.9539925233403E-14 + ,-1.9539925233403E-14 + ,-1.9539925233403E-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] 1.893424e-05 > sqrtn <- sqrt(length(x)) > (armse <- sd(x) / sqrtn) [1] 0.0001863188 > (armose <- arm / armse) [1] 0.1016228 > (geo <- geomean(x)) [1] NaN Warning message: In log(x) : NaNs produced > (har <- harmean(x)) [1] 2.547918e-13 > (qua <- quamean(x)) [1] 0.00467289 > (win <- winmean(x)) [,1] [,2] [1,] 1.794093e-05 1.861352e-04 [2,] 1.293186e-05 1.813394e-04 [3,] -9.252226e-06 1.787923e-04 [4,] -9.656287e-06 1.787553e-04 [5,] -9.656287e-06 1.787553e-04 [6,] -9.656287e-06 1.787553e-04 [7,] -8.382297e-06 1.786396e-04 [8,] -8.382297e-06 1.786396e-04 [9,] -1.002028e-05 1.784906e-04 [10,] -9.693737e-06 1.784611e-04 [11,] -8.271365e-06 1.783334e-04 [12,] -1.021490e-05 1.781581e-04 [13,] -1.021490e-05 1.781581e-04 [14,] -1.021490e-05 1.781581e-04 [15,] -1.021490e-05 1.781581e-04 [16,] -1.021490e-05 1.781581e-04 [17,] -1.021490e-05 1.781581e-04 [18,] -1.021490e-05 1.780548e-04 [19,] -1.021490e-05 1.780548e-04 [20,] -1.021490e-05 1.780548e-04 [21,] -1.021490e-05 1.780548e-04 [22,] -1.021490e-05 1.780548e-04 [23,] -1.021490e-05 1.780548e-04 [24,] -1.021490e-05 1.780548e-04 [25,] -1.021490e-05 1.780548e-04 [26,] -1.021490e-05 1.780548e-04 [27,] -9.354297e-06 1.779779e-04 [28,] -1.024678e-05 1.778979e-04 [29,] -1.024678e-05 1.778979e-04 [30,] -1.024678e-05 1.778979e-04 [31,] -1.024678e-05 1.778979e-04 [32,] -1.024678e-05 1.778979e-04 [33,] -1.024678e-05 1.778979e-04 [34,] -5.954929e-06 1.775154e-04 [35,] 3.241582e-05 1.734109e-04 [36,] -1.159554e-05 1.695199e-04 [37,] -1.159554e-05 1.695199e-04 [38,] -1.159554e-05 1.695199e-04 [39,] -1.159554e-05 1.695199e-04 [40,] -5.613405e-07 1.685790e-04 [41,] -1.187140e-05 1.676077e-04 [42,] -1.187140e-05 1.676077e-04 [43,] -1.187140e-05 1.676077e-04 [44,] -1.187140e-05 1.584585e-04 [45,] -1.247190e-05 1.584090e-04 [46,] -1.247064e-05 1.583081e-04 [47,] -1.184474e-05 1.582567e-04 [48,] -1.184474e-05 1.582567e-04 [49,] -1.184474e-05 1.582567e-04 [50,] -1.184474e-05 1.582567e-04 [51,] -1.589078e-05 1.579240e-04 [52,] -1.176540e-05 1.575852e-04 [53,] -1.176540e-05 1.575852e-04 [54,] -1.176540e-05 1.575852e-04 [55,] -1.176540e-05 1.575852e-04 [56,] -1.615421e-05 1.572249e-04 [57,] -1.615421e-05 1.572249e-04 [58,] -1.160866e-05 1.568521e-04 [59,] -2.514980e-05 1.557454e-04 [60,] -2.514980e-05 1.557454e-04 [61,] -1.796369e-05 1.540418e-04 [62,] -9.242162e-05 1.481256e-04 [63,] -9.725293e-06 1.270181e-04 [64,] -9.725293e-06 1.270181e-04 [65,] -9.725293e-06 1.270181e-04 [66,] -9.725293e-06 1.270181e-04 [67,] -9.725293e-06 1.270181e-04 [68,] -9.725293e-06 1.270181e-04 [69,] -9.725293e-06 1.270181e-04 [70,] -9.725293e-06 1.270181e-04 [71,] -9.725293e-06 1.261247e-04 [72,] -9.725293e-06 1.257890e-04 [73,] -9.725293e-06 1.257890e-04 [74,] -1.846011e-05 1.249973e-04 [75,] -9.607254e-06 1.243104e-04 [76,] -1.401164e-05 1.239689e-04 [77,] -1.401164e-05 1.239689e-04 [78,] -2.489232e-05 1.224288e-04 [79,] -2.489232e-05 1.224288e-04 [80,] -1.087309e-04 1.136961e-04 [81,] -1.167066e-04 1.131151e-04 [82,] -2.548110e-04 1.035736e-04 [83,] -2.797257e-04 8.661489e-05 [84,] 1.255363e-04 5.387612e-05 [85,] 1.253905e-04 5.352359e-05 [86,] -2.857606e-06 4.303893e-05 [87,] -1.345171e-04 3.289716e-05 [88,] -1.780308e-06 2.180362e-05 [89,] -4.795056e-06 2.155802e-05 [90,] -1.453895e-04 1.410852e-05 [91,] -9.671276e-15 1.261676e-15 [92,] -9.606425e-15 1.258595e-15 [93,] -9.606425e-15 1.258595e-15 [94,] -9.606425e-15 1.258595e-15 [95,] -9.606425e-15 1.258595e-15 [96,] -1.102751e-14 1.178112e-15 [97,] -1.082238e-14 1.168249e-15 [98,] -1.082238e-14 1.168249e-15 [99,] -1.082238e-14 1.168249e-15 [100,] -1.082238e-14 1.168249e-15 [101,] -1.082238e-14 1.168249e-15 [102,] -1.082238e-14 1.168249e-15 [103,] -1.082238e-14 1.168249e-15 [104,] -1.082238e-14 1.168249e-15 [105,] -1.082238e-14 1.168249e-15 [106,] -1.082238e-14 1.168249e-15 [107,] -1.082238e-14 1.168249e-15 [108,] -1.082238e-14 1.168249e-15 [109,] -1.082238e-14 1.168249e-15 [110,] -1.082238e-14 1.168249e-15 [111,] -1.129185e-14 1.143619e-15 [112,] -1.129185e-14 1.143619e-15 [113,] -1.129185e-14 1.143619e-15 [114,] -1.137221e-14 1.139493e-15 [115,] -1.137221e-14 1.139493e-15 [116,] -1.137221e-14 1.139493e-15 [117,] -1.137221e-14 1.139493e-15 [118,] -1.137221e-14 1.139493e-15 [119,] -1.137221e-14 1.139493e-15 [120,] -1.137221e-14 1.139493e-15 [121,] -1.137221e-14 1.139493e-15 [122,] -1.137221e-14 1.139493e-15 [123,] -1.137221e-14 1.139493e-15 [124,] -1.137221e-14 1.139493e-15 [125,] -1.234145e-14 1.091521e-15 [126,] -1.234145e-14 1.091521e-15 [127,] -1.207288e-14 1.078587e-15 [128,] -1.207288e-14 1.078587e-15 [129,] -1.234568e-14 1.065654e-15 [130,] -1.234568e-14 1.065654e-15 [131,] -1.234568e-14 1.065654e-15 [132,] -1.234568e-14 1.065654e-15 [133,] -1.234568e-14 1.065654e-15 [134,] -1.234568e-14 1.065654e-15 [135,] -1.234568e-14 1.065654e-15 [136,] -1.234568e-14 1.065654e-15 [137,] -1.234568e-14 1.065654e-15 [138,] -1.234568e-14 1.065654e-15 [139,] -1.234568e-14 1.065654e-15 [140,] -1.234568e-14 1.065654e-15 [141,] -1.234568e-14 1.065654e-15 [142,] -1.244578e-14 1.060975e-15 [143,] -1.244578e-14 1.060975e-15 [144,] -1.244578e-14 1.060975e-15 [145,] -1.244578e-14 1.060975e-15 [146,] -1.244578e-14 1.060975e-15 [147,] -1.244578e-14 1.060975e-15 [148,] -1.244578e-14 1.060975e-15 [149,] -1.244578e-14 1.060975e-15 [150,] -1.244578e-14 1.060975e-15 [151,] -1.244578e-14 1.060975e-15 [152,] -1.244578e-14 1.060975e-15 [153,] -1.244578e-14 1.060975e-15 [154,] -1.244578e-14 1.060975e-15 [155,] -1.244578e-14 1.060975e-15 [156,] -1.244578e-14 1.060975e-15 [157,] -1.244578e-14 1.060975e-15 [158,] -1.233440e-14 1.055655e-15 [159,] -1.233440e-14 1.055655e-15 [160,] -1.233440e-14 1.055655e-15 [161,] -1.233440e-14 1.055655e-15 [162,] -1.233440e-14 1.055655e-15 [163,] -1.233440e-14 1.055655e-15 [164,] -1.233440e-14 1.055655e-15 [165,] -1.233440e-14 1.055655e-15 [166,] -1.233440e-14 1.055655e-15 [167,] -1.233440e-14 1.055655e-15 [168,] -1.233440e-14 1.055655e-15 [169,] -1.233440e-14 1.055655e-15 [170,] -1.233440e-14 1.055655e-15 [171,] -1.233440e-14 1.055655e-15 [172,] -1.233440e-14 1.055655e-15 [173,] -1.233440e-14 1.055655e-15 [174,] -1.233440e-14 1.055655e-15 [175,] -1.233440e-14 1.055655e-15 [176,] -1.233440e-14 1.055655e-15 [177,] -1.233440e-14 1.055655e-15 [178,] -1.271082e-14 1.038216e-15 [179,] -1.258464e-14 1.032190e-15 [180,] -1.068140e-14 9.447750e-16 [181,] -1.004347e-14 9.171176e-16 [182,] -1.004347e-14 9.171176e-16 [183,] -1.004347e-14 9.171176e-16 [184,] -1.004347e-14 9.171176e-16 [185,] -1.004347e-14 9.171176e-16 [186,] -1.004347e-14 9.171176e-16 [187,] -1.004347e-14 9.171176e-16 [188,] -1.004347e-14 9.171176e-16 [189,] -1.004347e-14 9.171176e-16 [190,] -1.004347e-14 9.171176e-16 [191,] -1.004347e-14 9.171176e-16 [192,] -1.004347e-14 9.171176e-16 [193,] -1.004347e-14 9.171176e-16 [194,] -1.004347e-14 9.171176e-16 [195,] -1.004347e-14 9.171176e-16 [196,] -1.018163e-14 9.105775e-16 [197,] -1.018163e-14 9.105775e-16 [198,] -1.018163e-14 9.105775e-16 [199,] -1.018163e-14 9.105775e-16 [200,] -1.018163e-14 9.105775e-16 [201,] -1.018163e-14 9.105775e-16 [202,] -1.018163e-14 9.105775e-16 [203,] -1.018163e-14 9.105775e-16 [204,] -1.018163e-14 9.105775e-16 [205,] -1.018163e-14 9.105775e-16 [206,] -1.018163e-14 9.105775e-16 [207,] -1.018163e-14 9.105775e-16 [208,] -1.018163e-14 9.105775e-16 [209,] -1.018163e-14 9.105775e-16 [210,] -1.018163e-14 9.105775e-16 > (tri <- trimean(x)) [,1] [,2] [1,] 7.965976e-06 1.820289e-04 [2,] -2.072721e-06 1.777704e-04 [3,] -9.647150e-06 1.759200e-04 [4,] -9.780484e-06 1.749232e-04 [5,] -9.812034e-06 1.739157e-04 [6,] -9.843788e-06 1.728873e-04 [7,] -9.875748e-06 1.718374e-04 [8,] -1.009466e-05 1.707835e-04 [9,] -1.031500e-05 1.697073e-04 [10,] -1.034882e-05 1.686268e-04 [11,] -1.041670e-05 1.675264e-04 [12,] -1.061945e-05 1.664156e-04 [13,] -1.065462e-05 1.652977e-04 [14,] -1.069001e-05 1.641553e-04 [15,] -1.072565e-05 1.629876e-04 [16,] -1.076152e-05 1.617938e-04 [17,] -1.079763e-05 1.605730e-04 [18,] -1.083399e-05 1.593243e-04 [19,] -1.087059e-05 1.580541e-04 [20,] -1.090744e-05 1.567542e-04 [21,] -1.094454e-05 1.554237e-04 [22,] -1.098189e-05 1.540613e-04 [23,] -1.101950e-05 1.526660e-04 [24,] -1.105737e-05 1.512364e-04 [25,] -1.109550e-05 1.497713e-04 [26,] -1.113389e-05 1.482693e-04 [27,] -1.117255e-05 1.467287e-04 [28,] -1.124646e-05 1.451524e-04 [29,] -1.128578e-05 1.435387e-04 [30,] -1.132538e-05 1.418817e-04 [31,] -1.136526e-05 1.401795e-04 [32,] -1.140542e-05 1.384299e-04 [33,] -1.144587e-05 1.366308e-04 [34,] -1.148660e-05 1.347797e-04 [35,] -1.166963e-05 1.328930e-04 [36,] -1.309174e-05 1.311494e-04 [37,] -1.313884e-05 1.295391e-04 [38,] -1.318627e-05 1.278830e-04 [39,] -1.323404e-05 1.261789e-04 [40,] -1.328217e-05 1.244243e-04 [41,] -1.364778e-05 1.226595e-04 [42,] -1.369777e-05 1.208846e-04 [43,] -1.369777e-05 1.190548e-04 [44,] -1.379886e-05 1.171670e-04 [45,] -1.384997e-05 1.156177e-04 [46,] -1.388583e-05 1.140242e-04 [47,] -1.392199e-05 1.123864e-04 [48,] -1.397413e-05 1.106998e-04 [49,] -1.402666e-05 1.089593e-04 [50,] -1.407959e-05 1.071620e-04 [51,] -1.413292e-05 1.053044e-04 [52,] -1.409164e-05 1.033970e-04 [53,] -1.414543e-05 1.014368e-04 [54,] -1.419962e-05 9.940528e-05 [55,] -1.425424e-05 9.729750e-05 [56,] -1.430927e-05 9.510772e-05 [57,] -1.426905e-05 9.284577e-05 [58,] -1.422851e-05 9.048915e-05 [59,] -1.428409e-05 8.804691e-05 [60,] -1.405660e-05 8.554568e-05 [61,] -1.382731e-05 8.292630e-05 [62,] -1.374288e-05 8.026003e-05 [63,] -1.215662e-05 7.774475e-05 [64,] -1.220505e-05 7.606420e-05 [65,] -1.225387e-05 7.431646e-05 [66,] -1.230308e-05 7.249612e-05 [67,] -1.235269e-05 7.059699e-05 [68,] -1.240270e-05 6.861193e-05 [69,] -1.245312e-05 6.653259e-05 [70,] -1.250395e-05 6.434917e-05 [71,] -1.255519e-05 6.204996e-05 [72,] -1.260686e-05 5.966570e-05 [73,] -1.265896e-05 5.715628e-05 [74,] -1.271148e-05 5.448561e-05 [75,] -1.260952e-05 5.167251e-05 [76,] -1.266228e-05 4.868538e-05 [77,] -1.263878e-05 4.546746e-05 [78,] -1.261509e-05 4.194163e-05 [79,] -1.240500e-05 3.813479e-05 [80,] -1.219312e-05 3.382535e-05 [81,] -1.219312e-05 2.958165e-05 [82,] -8.797188e-06 2.458090e-05 [83,] -4.723683e-06 1.932273e-05 [84,] -2.055770e-07 1.449953e-05 [85,] -2.255716e-06 1.221355e-05 [86,] -2.255716e-06 9.372997e-06 [87,] -4.344916e-06 6.965564e-06 [88,] -2.268646e-06 4.993626e-06 [89,] -2.276381e-06 3.865501e-06 [90,] -2.236761e-06 2.236761e-06 [91,] -2.236761e-06 1.211946e-15 [92,] -1.164589e-14 1.209651e-15 [93,] -1.167735e-14 1.207351e-15 [94,] -1.170909e-14 1.204977e-15 [95,] -1.174111e-14 1.202529e-15 [96,] -1.177343e-14 1.200003e-15 [97,] -1.178466e-14 1.199239e-15 [98,] -1.179906e-14 1.198640e-15 [99,] -1.181360e-14 1.198002e-15 [100,] -1.182826e-14 1.197327e-15 [101,] -1.184307e-14 1.196611e-15 [102,] -1.185802e-14 1.195855e-15 [103,] -1.187310e-14 1.195056e-15 [104,] -1.188833e-14 1.194215e-15 [105,] -1.190371e-14 1.193329e-15 [106,] -1.191923e-14 1.192397e-15 [107,] -1.193490e-14 1.191419e-15 [108,] -1.195072e-14 1.190392e-15 [109,] -1.196670e-14 1.189316e-15 [110,] -1.198283e-14 1.188188e-15 [111,] -1.199912e-14 1.187008e-15 [112,] -1.200900e-14 1.186308e-15 [113,] -1.201899e-14 1.185565e-15 [114,] -1.202907e-14 1.184776e-15 [115,] -1.203815e-14 1.184030e-15 [116,] -1.204731e-14 1.183237e-15 [117,] -1.205657e-14 1.182398e-15 [118,] -1.206593e-14 1.181510e-15 [119,] -1.207537e-14 1.180572e-15 [120,] -1.208492e-14 1.179583e-15 [121,] -1.209456e-14 1.178540e-15 [122,] -1.210431e-14 1.177443e-15 [123,] -1.211415e-14 1.176289e-15 [124,] -1.212410e-14 1.175078e-15 [125,] -1.213415e-14 1.173806e-15 [126,] -1.213139e-14 1.173496e-15 [127,] -1.212860e-14 1.173146e-15 [128,] -1.212933e-14 1.173023e-15 [129,] -1.213008e-14 1.172863e-15 [130,] -1.212724e-14 1.172931e-15 [131,] -1.212436e-14 1.172965e-15 [132,] -1.212145e-14 1.172965e-15 [133,] -1.211851e-14 1.172930e-15 [134,] -1.211554e-14 1.172858e-15 [135,] -1.211253e-14 1.172747e-15 [136,] -1.210949e-14 1.172598e-15 [137,] -1.210642e-14 1.172408e-15 [138,] -1.210331e-14 1.172176e-15 [139,] -1.210017e-14 1.171901e-15 [140,] -1.209699e-14 1.171581e-15 [141,] -1.209377e-14 1.171214e-15 [142,] -1.209052e-14 1.170800e-15 [143,] -1.208594e-14 1.170436e-15 [144,] -1.208130e-14 1.170023e-15 [145,] -1.207661e-14 1.169559e-15 [146,] -1.207187e-14 1.169042e-15 [147,] -1.206707e-14 1.168472e-15 [148,] -1.206221e-14 1.167845e-15 [149,] -1.205729e-14 1.167160e-15 [150,] -1.205231e-14 1.166415e-15 [151,] -1.204727e-14 1.165608e-15 [152,] -1.204217e-14 1.164736e-15 [153,] -1.203701e-14 1.163797e-15 [154,] -1.203178e-14 1.162788e-15 [155,] -1.202649e-14 1.161707e-15 [156,] -1.202113e-14 1.160552e-15 [157,] -1.201570e-14 1.159319e-15 [158,] -1.201021e-14 1.158006e-15 [159,] -1.200607e-14 1.156737e-15 [160,] -1.200187e-14 1.155385e-15 [161,] -1.199762e-14 1.153947e-15 [162,] -1.199762e-14 1.152420e-15 [163,] -1.198895e-14 1.150800e-15 [164,] -1.198453e-14 1.149084e-15 [165,] -1.198005e-14 1.147267e-15 [166,] -1.197551e-14 1.145346e-15 [167,] -1.197090e-14 1.143317e-15 [168,] -1.196624e-14 1.141175e-15 [169,] -1.196151e-14 1.138915e-15 [170,] -1.195672e-14 1.136534e-15 [171,] -1.195186e-14 1.134026e-15 [172,] -1.195186e-14 1.131385e-15 [173,] -1.194193e-14 1.128606e-15 [174,] -1.193687e-14 1.125683e-15 [175,] -1.193173e-14 1.122610e-15 [176,] -1.192651e-14 1.119381e-15 [177,] -1.192122e-14 1.115988e-15 [178,] -1.191585e-14 1.112424e-15 [179,] -1.190551e-14 1.109169e-15 [180,] -1.189666e-14 1.105922e-15 [181,] -1.191253e-14 1.104890e-15 [182,] -1.191253e-14 1.104442e-15 [183,] -1.196181e-14 1.103914e-15 [184,] -1.198702e-14 1.103303e-15 [185,] -1.201261e-14 1.102604e-15 [186,] -1.203860e-14 1.101814e-15 [187,] -1.206500e-14 1.100929e-15 [188,] -1.209181e-14 1.099944e-15 [189,] -1.211905e-14 1.098854e-15 [190,] -1.214673e-14 1.097654e-15 [191,] -1.217485e-14 1.096339e-15 [192,] -1.220343e-14 1.094903e-15 [193,] -1.223247e-14 1.093341e-15 [194,] -1.226200e-14 1.091645e-15 [195,] -1.229202e-14 1.089810e-15 [196,] -1.232254e-14 1.087829e-15 [197,] -1.235170e-14 1.085910e-15 [198,] -1.238136e-14 1.083839e-15 [199,] -1.241153e-14 1.081606e-15 [200,] -1.244222e-14 1.079204e-15 [201,] -1.244222e-14 1.076624e-15 [202,] -1.250524e-14 1.073855e-15 [203,] -1.253759e-14 1.070889e-15 [204,] -1.257053e-14 1.067713e-15 [205,] -1.260406e-14 1.064317e-15 [206,] -1.263821e-14 1.060689e-15 [207,] -1.267299e-14 1.056814e-15 [208,] -1.270842e-14 1.052679e-15 [209,] -1.274452e-14 1.048268e-15 [210,] -1.278131e-14 1.043565e-15 > (midr <- midrange(x)) [1] 0.003462969 > midm <- array(NA,dim=8) > for (j in 1:8) midm[j] <- midmean(x,j) > midm [1] -1.343339e-14 -1.343339e-14 -1.343339e-14 -1.343339e-14 -1.052357e-14 [6] -1.343339e-14 -1.343339e-14 -1.343339e-14 > postscript(file="/var/www/html/freestat/rcomp/tmp/1u4p41287598601.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/2u4p41287598601.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/3pgx61287598601.tab") > > try(system("convert tmp/1u4p41287598601.ps tmp/1u4p41287598601.png",intern=TRUE)) character(0) > try(system("convert tmp/2u4p41287598601.ps tmp/2u4p41287598601.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 5.111 0.589 5.206