R version 2.9.0 (2009-04-17) Copyright (C) 2009 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(15086 + ,20559 + ,39383 + ,68032 + ,42346 + ,40184 + ,29183 + ,18049 + ,24377 + ,48033 + ,61410 + ,40637 + ,37914 + ,22481 + ,17061 + ,21921 + ,34043 + ,36152 + ,38929 + ,36713 + ,27928 + ,14178 + ,19331 + ,35030 + ,37834 + ,45257 + ,39650 + ,33001 + ,22268 + ,34283 + ,63733 + ,54548 + ,31666 + ,32654 + ,28409 + ,16234 + ,19544 + ,33375 + ,30839 + ,27581 + ,28596 + ,28542 + ,14285 + ,20799 + ,28729 + ,30278 + ,31266 + ,38341 + ,24938 + ,15700 + ,17382 + ,34229 + ,29584 + ,33909 + ,29397 + ,24644 + ,14124 + ,18770 + ,32681 + ,30518 + ,34310 + ,24217 + ,15806 + ,12362 + ,13751 + ,25312 + ,23549 + ,30305 + ,30411 + ,15593 + ,17382 + ,15806 + ,29183 + ,33348 + ,33589 + ,30198 + ,23469 + ,15593 + ,16474 + ,27341 + ,26673 + ,33268 + ,26406 + ,23656 + ,15086 + ,18316 + ,27528 + ,30305 + ,34550 + ,34710 + ,23843 + ,12362 + ,21200 + ,27261 + ,31586 + ,33322 + ,32814 + ,21974 + ,17115 + ,14445 + ,26673 + ,31746 + ,28249 + ,29103 + ,21974 + ,17622 + ,21066 + ,28382 + ,32147 + ,35885 + ,40824 + ,28622 + ,24297 + ,30465 + ,47927 + ,51665 + ,57058 + ,48407 + ,30011 + ,21413 + ,30518 + ,41171 + ,37487 + ,44456 + ,41465 + ,28409 + ,34870 + ,29557 + ,44936 + ,44028 + ,45310 + ,39302 + ,28516 + ,19865 + ,30091 + ,38101 + ,44402 + ,52359 + ,71182 + ,60796 + ,43735 + ,63947 + ,89071 + ,86321 + ,90967 + ,54441 + ,36339 + ,29290 + ,65842 + ,99404 + ,83865 + ,84452 + ,76309 + ,69554 + ,44616 + ,67337 + ,94144 + ,84425 + ,76763 + ,73559 + ,60956 + ,46271 + ,74066 + ,91634 + ,63947 + ,56017 + ,53293 + ,54895 + ,48781 + ,60449 + ,99965 + ,84746 + ,81088 + ,77777 + ,60315 + ,46218 + ,62985 + ,92596 + ,78525 + ,72811 + ,62318 + ,48140 + ,45443 + ,66910 + ,90032 + ,91661 + ,103516 + ,100285 + ,72384 + ,51985 + ,72224 + ,99351 + ,108616 + ,102128 + ,97669 + ,71716 + ,51584 + ,63413 + ,88137 + ,86535 + ,91955 + ,84853 + ,54388 + ,40931 + ,58900 + ,76576 + ,80447 + ,60315 + ,52012 + ,37487 + ,19972 + ,25338 + ,23549 + ,19785 + ,29637 + ,34416 + ,26914 + ,25071 + ,24003 + ,19358 + ,28783 + ,46992 + ,49796 + ,40904 + ,33722 + ,42640 + ,56097 + ,62184 + ,63039 + ,61730 + ,44349 + ,32307 + ,45123 + ,63466 + ,71903 + ,73211 + ,67845 + ,51237 + ,31506 + ,45043 + ,74386 + ,76602 + ,72704 + ,71770 + ,52225 + ,40958 + ,55856 + ,79700 + ,81889 + ,82049 + ,66643 + ,54655 + ,36739 + ,53026 + ,78098 + ,81248 + ,84372 + ,75614 + ,63012 + ,38795 + ,55670 + ,81168 + ,85787 + ,88591 + ,75401 + ,63413 + ,42907 + ,62291 + ,86241 + ,85200 + ,84559 + ,75801 + ,72304 + ,50944 + ,70034 + ,104691 + ,88404 + ,93397 + ,88030 + ,60582 + ,50890 + ,69847 + ,90166 + ,84158 + ,82503 + ,83811 + ,66296 + ,42640 + ,55563 + ,82450 + ,84799 + ,82370 + ,83384 + ,60075 + ,41091 + ,55589 + ,79513 + ,87656 + ,74760 + ,64214 + ,47019 + ,34229 + ,45150 + ,64187 + ,73879 + ,76175 + ,83785 + ,52786 + ,43681 + ,61864 + ,74573 + ,79913 + ,78765 + ,79112 + ,62718 + ,40557 + ,57058 + ,88804 + ,89205 + ,82904 + ,78445 + ,62825 + ,57832 + ,70675 + ,102795 + ,107094 + ,105625 + ,98390 + ,74360 + ,64641 + ,82423 + ,116759 + ,101941 + ,101140 + ,98790 + ,65949 + ,49635 + ,77323 + ,109203 + ,102902 + ,93637 + ,89365 + ,72250 + ,53267 + ,73238 + ,99725 + ,69100 + ,66323 + ,80447 + ,80954 + ,48194 + ,52546 + ,91821 + ,84372 + ,90994 + ,72410 + ,62745 + ,47366 + ,63947 + ,82690 + ,89712 + ,79085 + ,80527 + ,61303 + ,40691 + ,58340 + ,79166 + ,82850 + ,82076 + ,73292 + ,60823 + ,51157 + ,61410 + ,64294 + ,67658 + ,66323 + ,66563 + ,52145 + ,36980 + ,50703 + ,68112 + ,68379 + ,65735 + ,60716 + ,50463 + ,33749 + ,39436 + ,66990 + ,70541 + ,75027 + ,69233 + ,53213 + ,31186 + ,41145 + ,67604 + ,68699 + ,71690 + ,60022 + ,51718 + ,32814 + ,42400 + ,71316 + ,70141 + ,64694 + ,60636 + ,37620 + ,35431 + ,41145 + ,69393 + ,72891 + ,62638 + ,63439 + ,54014 + ,37567 + ,38875 + ,64854 + ,75214 + ,72944 + ,71796 + ,60689 + ,53000 + ,50757 + ,71289 + ,75080 + ,77377 + ,68779 + ,49155 + ,30678 + ,49342 + ,67604 + ,63733 + ,71209 + ,59461 + ,49609 + ,30839 + ,39302 + ,67685 + ,67925 + ,66056 + ,60743 + ,49422 + ,30972 + ,34924 + ,56818 + ,62024 + ,59274 + ,59007 + ,43494 + ,26486 + ,41252 + ,63146 + ,70114 + ,57779 + ,61944 + ,48728 + ,32708 + ,40237 + ,67845 + ,68913 + ,66697 + ,68272 + ,51317 + ,30812 + ,37300 + ,63599 + ,72918 + ,74386 + ,73292 + ,56577 + ,41652 + ,58900 + ,86188 + ,88350 + ,85787 + ,86535 + ,67978 + ,42720 + ,64187 + ,98897 + ,108429 + ,92249 + ,111820 + ,77163 + ,59648 + ,83651 + ,103569 + ,108749 + ,96948 + ,90113 + ,75614 + ,52733 + ,78178 + ,103302 + ,96040 + ,98843 + ,93183 + ,75267 + ,54815 + ,93904 + ,110725 + ,107735 + ,100793 + ,102795 + ,73772 + ,63092 + ,82476 + ,162042 + ,108215 + ,100846 + ,104183 + ,82156 + ,59728 + ,97028 + ,148746 + ,101380 + ,100953 + ,100205 + ,73292 + ,57138 + ,63039 + ,135236 + ,104370 + ,109844 + ,104290 + ,78231 + ,58153 + ,88938 + ,132459 + ,107654 + ,113021 + ,117587 + ,85387 + ,64347 + ,88778 + ,130777 + ,159185 + ,139988 + ,136731 + ,90273 + ,61277 + ,81275 + ,82503 + ,74333 + ,77350 + ,102662 + ,76923 + ,62291 + ,97455 + ,127733 + ,101861 + ,87763 + ,80954 + ,72918 + ,67952 + ,96574 + ,148105 + ,151576 + ,132592 + ,117213 + ,102982 + ,92622 + ,114196 + ,145515 + ,157664 + ,129522 + ,121886 + ,99271 + ,76923 + ,103649 + ,111232 + ,123621 + ,104077 + ,112354 + ,69260 + ,47286 + ,66910 + ,78738 + ,62905 + ,44082 + ,25392 + ,31666 + ,30732 + ,35751 + ,47579 + ,44002 + ,30064 + ,28809 + ,44375 + ,39249 + ,48007 + ,73452 + ,91154 + ,71342 + ,63226 + ,58820 + ,50623 + ,63573 + ,75347 + ,80714 + ,88644 + ,82370 + ,63439 + ,44776 + ,62024 + ,86081 + ,95212 + ,88244 + ,80287 + ,61330 + ,47206 + ,64400 + ,89579 + ,97722 + ,109203 + ,120551 + ,80234 + ,52679 + ,76923 + ,99164 + ,100499 + ,101620 + ,94385 + ,75241 + ,58420 + ,72677 + ,76148 + ,70568 + ,76469 + ,71449 + ,49368 + ,33135 + ,62131 + ,77003 + ,81836 + ,75107 + ,112968 + ,84399 + ,65015 + ,83758 + ,116225 + ,104984 + ,115451 + ,115718 + ,82930 + ,61170 + ,88377 + ,118254 + ,124208 + ,124502 + ,109871 + ,81782 + ,59034 + ,79593 + ,117667 + ,107788 + ,108135 + ,104958 + ,78685 + ,57378 + ,76496 + ,104370 + ,102875 + ,99591 + ,97802 + ,79299 + ,57138 + ,78738 + ,123007 + ,117293 + ,147731 + ,156409 + ,85173 + ,66430 + ,85387 + ,123728 + ,114703 + ,113769 + ,96307 + ,79433 + ,58607 + ,75828 + ,108429 + ,104157 + ,99538 + ,90353 + ,69874 + ,53694 + ,66910 + ,102261 + ,104584 + ,109550 + ,99858 + ,83598 + ,64748 + ,86642 + ,84719 + ,79566 + ,93236 + ,81649 + ,61891 + ,44723 + ,62291 + ,78365 + ,127493 + ,123568 + ,114356 + ,88911 + ,63092 + ,85574 + ,115104 + ,115638 + ,107975 + ,117533 + ,80474 + ,58259 + ,78231 + ,70408 + ,73986 + ,78471 + ,83037 + ,52813 + ,34123 + ,45791 + ,64454 + ,61196 + ,63039 + ,59781 + ,50997 + ,31346 + ,37861 + ,53614 + ,66857 + ,73478 + ,82183 + ,51531 + ,35564 + ,42346 + ,75054 + ,81088 + ,76255 + ,71503 + ,57992 + ,46298 + ,58313 + ,92943 + ,97055 + ,93664 + ,79993 + ,67364 + ,53106 + ,63813 + ,87843 + ,95239 + ,91661 + ,85947 + ,70408 + ,47152 + ,56043 + ,79753 + ,82530 + ,76148 + ,71583 + ,59060 + ,38582 + ,54521 + ,81702 + ,81755 + ,77136 + ,64748 + ,51932 + ,34683 + ,39730 + ,68325 + ,71049 + ,68779 + ,72758 + ,51878 + ,35484 + ,47499 + ,68405 + ,77350 + ,59461 + ,67578 + ,56524 + ,35030 + ,41839 + ,68112 + ,61570 + ,76335 + ,68432 + ,46832 + ,36766 + ,42320 + ,66563 + ,67498 + ,71903 + ,63813 + ,53000 + ,32467 + ,41065 + ,64187 + ,69740 + ,68539 + ,59327 + ,49876 + ,35111 + ,48487 + ,71503 + ,68726 + ,66803 + ,65735 + ,54068 + ,30732 + ,47606 + ,72117 + ,80207 + ,77136 + ,67231 + ,54468 + ,39169 + ,53000 + ,78551 + ,76148 + ,78738 + ,74119 + ,54602 + ,41679 + ,52279 + ,79379 + ,90219 + ,83571 + ,86081 + ,58019 + ,41465 + ,56871 + ,83090 + ,96280 + ,94091 + ,97161 + ,70141 + ,50303 + ,66002 + ,84746 + ,88404 + ,119936 + ,97669 + ,78712 + ,59114 + ,85200 + ,101380 + ,103382 + ,111713 + ,106373 + ,75134 + ,63012 + ,77056 + ,100712 + ,117800 + ,118254 + ,102448 + ,79272 + ,68592 + ,93877 + ,133206 + ,147571 + ,144687 + ,111846 + ,87069 + ,65976 + ,95880 + ,150828 + ,128614 + ,122526 + ,105946 + ,82984 + ,66216 + ,98897 + ,128080 + ,140122 + ,125223 + ,122473 + ,85467 + ,66350 + ,86722 + ,106266 + ,123381 + ,128294 + ,117080 + ,84292 + ,70007 + ,101647 + ,151336 + ,149200 + ,143646 + ,106319 + ,78284 + ,58553 + ,88564 + ,114303 + ,120337 + ,123354 + ,117133 + ,101674 + ,60876 + ,87816 + ,118495 + ,114169 + ,117720 + ,120871 + ,84826 + ,64774 + ,91955 + ,115184 + ,103783 + ,96948 + ,89979 + ,84933 + ,72143 + ,87896 + ,122286 + ,131631 + ,132806 + ,139508 + ,102074 + ,74466 + ,104717 + ,132272 + ,130323 + ,138066 + ,119990 + ,102261 + ,78605 + ,104584 + ,127599 + ,111579 + ,123381 + ,103890 + ,61811 + ,47366 + ,53507 + ,75614 + ,60262 + ,47980 + ,28703 + ,20986 + ,29290 + ,37594 + ,48888 + ,50757 + ,56417 + ,30465 + ,30411 + ,38822 + ,50383 + ,73051 + ,68325 + ,84078 + ,93477 + ,72090 + ,56177 + ,77243 + ,102181 + ,93877 + ,90433 + ,88617 + ,65655 + ,67177 + ,62665 + ,95212 + ,103543 + ,109871 + ,105305 + ,78738 + ,56898 + ,109497 + ,168344 + ,123915 + ,120711 + ,113128 + ,83838 + ,67151 + ,85574 + ,116999 + ,136117 + ,127065 + ,111473 + ,80100 + ,74093 + ,102021 + ,146690 + ,150508 + ,135289 + ,116012 + ,78498 + ,59701 + ,86294 + ,107468 + ,115371 + ,118388 + ,110485 + ,100659 + ,79139 + ,109310 + ,139935 + ,126051 + ,141270 + ,118949 + ,97455 + ,77777 + ,107040 + ,149760 + ,147731 + ,159132 + ,124369 + ,98016 + ,73185 + ,102528 + ,242516 + ,293460 + ,320987 + ,268736 + ,239873 + ,189143 + ,248550 + ,269670 + ,278801 + ,272634 + ,181240 + ,90566 + ,67017 + ,92809 + ,134755 + ,137291 + ,138199 + ,125704 + ,107414 + ,75588 + ,100606 + ,120791 + ,124075 + ,101460 + ,105198 + ,70007 + ,58446 + ,70808 + ,101220 + ,127599 + ,120497 + ,110752 + ,90113 + ,70541 + ,84559 + ,120524 + ,118922 + ,123648 + ,113448 + ,84132 + ,65228 + ,101460 + ,125303 + ,141029 + ,130136 + ,124502 + ,97241 + ,71556 + ,97188 + ,141430 + ,134621 + ,134354 + ,113261 + ,74520 + ,70835 + ,87736 + ,128908 + ,118575 + ,118014 + ,111232 + ,79593 + ,56043 + ,81221 + ,115024 + ,91661 + ,110298 + ,99591 + ,70488 + ,50757 + ,69473 + ,117186 + ,82637 + ,91207 + ,76682 + ,66964 + ,44536 + ,64267 + ,99458 + ,89338 + ,104397 + ,90486 + ,70702 + ,58633 + ,68672 + ,93477 + ,110485 + ,155875 + ,86828 + ,74653 + ,57512 + ,88484 + ,78979 + ,84158 + ,76389 + ,77136 + ,55189 + ,36686 + ,50703 + ,75428 + ,78979 + ,78151 + ,69260 + ,56844 + ,39436 + ,51985 + ,73585 + ,67177 + ,65442 + ,61757 + ,50890 + ,29050 + ,48647 + ,70488 + ,63199 + ,62825 + ,66803 + ,49555 + ,30892 + ,40103 + ,59434 + ,67604 + ,70194 + ,63760 + ,53106 + ,34310 + ,44883 + ,68058 + ,74333 + ,71156 + ,67337 + ,52279 + ,37727 + ,43254 + ,74493 + ,77991 + ,68699 + ,66483 + ,62318 + ,35885 + ,49902 + ,78231 + ,67044 + ,71476 + ,68112 + ,52492 + ,41492 + ,46805 + ,69847 + ,73478 + ,78632 + ,74813 + ,60182 + ,41305 + ,39116 + ,72010 + ,79246 + ,71690 + ,69420 + ,57939 + ,46298 + ,60422 + ,81035 + ,82637 + ,79726 + ,86348 + ,66296 + ,33695 + ,52199 + ,72517 + ,72277 + ,73078 + ,76763 + ,65655 + ,34657 + ,56417 + ,68699 + ,77163 + ,79700 + ,78151 + ,59781 + ,44829 + ,46992 + ,77030 + ,92596 + ,84933 + ,83918 + ,63306 + ,47286 + ,77403 + ,102768 + ,118361 + ,107708 + ,94705 + ,76122 + ,45417 + ,79806 + ,106346 + ,95826 + ,93797 + ,99351 + ,83838 + ,54174 + ,84559 + ,116492 + ,114543 + ,123407 + ,114677 + ,110378 + ,96040 + ,141857 + ,233465 + ,146022 + ,142124 + ,118067 + ,83010 + ,75695 + ,104157 + ,150161 + ,174832 + ,139935 + ,108749 + ,109977 + ,78872 + ,126451 + ,161295) > par1 = '4' > par1 <- 7 > (n <- length(x)) [1] 1242 > (np <- floor(n / par1)) [1] 177 > arr <- array(NA,dim=c(par1,np)) > j <- 0 > k <- 1 > for (i in 1:(np*par1)) + { + j = j + 1 + arr[j,k] <- x[i] + if (j == par1) { + j = 0 + k=k+1 + } + } > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 15086 18049 17061 14178 22268 16234 14285 15700 14124 12362 17382 15593 [2,] 20559 24377 21921 19331 34283 19544 20799 17382 18770 13751 15806 16474 [3,] 39383 48033 34043 35030 63733 33375 28729 34229 32681 25312 29183 27341 [4,] 68032 61410 36152 37834 54548 30839 30278 29584 30518 23549 33348 26673 [5,] 42346 40637 38929 45257 31666 27581 31266 33909 34310 30305 33589 33268 [6,] 40184 37914 36713 39650 32654 28596 38341 29397 24217 30411 30198 26406 [7,] 29183 22481 27928 33001 28409 28542 24938 24644 15806 15593 23469 23656 [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 15086 12362 17115 17622 24297 21413 34870 19865 43735 29290 44616 46271 [2,] 18316 21200 14445 21066 30465 30518 29557 30091 63947 65842 67337 74066 [3,] 27528 27261 26673 28382 47927 41171 44936 38101 89071 99404 94144 91634 [4,] 30305 31586 31746 32147 51665 37487 44028 44402 86321 83865 84425 63947 [5,] 34550 33322 28249 35885 57058 44456 45310 52359 90967 84452 76763 56017 [6,] 34710 32814 29103 40824 48407 41465 39302 71182 54441 76309 73559 53293 [7,] 23843 21974 21974 28622 30011 28409 28516 60796 36339 69554 60956 54895 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [1,] 48781 46218 45443 51985 51584 40931 19972 25071 33722 32307 31506 40958 [2,] 60449 62985 66910 72224 63413 58900 25338 24003 42640 45123 45043 55856 [3,] 99965 92596 90032 99351 88137 76576 23549 19358 56097 63466 74386 79700 [4,] 84746 78525 91661 108616 86535 80447 19785 28783 62184 71903 76602 81889 [5,] 81088 72811 103516 102128 91955 60315 29637 46992 63039 73211 72704 82049 [6,] 77777 62318 100285 97669 84853 52012 34416 49796 61730 67845 71770 66643 [7,] 60315 48140 72384 71716 54388 37487 26914 40904 44349 51237 52225 54655 [,37] [,38] [,39] [,40] [,41] [,42] [,43] [,44] [,45] [,46] [,47] [,48] [1,] 36739 38795 42907 50944 50890 42640 41091 34229 43681 40557 57832 64641 [2,] 53026 55670 62291 70034 69847 55563 55589 45150 61864 57058 70675 82423 [3,] 78098 81168 86241 104691 90166 82450 79513 64187 74573 88804 102795 116759 [4,] 81248 85787 85200 88404 84158 84799 87656 73879 79913 89205 107094 101941 [5,] 84372 88591 84559 93397 82503 82370 74760 76175 78765 82904 105625 101140 [6,] 75614 75401 75801 88030 83811 83384 64214 83785 79112 78445 98390 98790 [7,] 63012 63413 72304 60582 66296 60075 47019 52786 62718 62825 74360 65949 [,49] [,50] [,51] [,52] [,53] [,54] [,55] [,56] [,57] [,58] [,59] [,60] [1,] 49635 53267 48194 47366 40691 51157 36980 33749 31186 32814 35431 37567 [2,] 77323 73238 52546 63947 58340 61410 50703 39436 41145 42400 41145 38875 [3,] 109203 99725 91821 82690 79166 64294 68112 66990 67604 71316 69393 64854 [4,] 102902 69100 84372 89712 82850 67658 68379 70541 68699 70141 72891 75214 [5,] 93637 66323 90994 79085 82076 66323 65735 75027 71690 64694 62638 72944 [6,] 89365 80447 72410 80527 73292 66563 60716 69233 60022 60636 63439 71796 [7,] 72250 80954 62745 61303 60823 52145 50463 53213 51718 37620 54014 60689 [,61] [,62] [,63] [,64] [,65] [,66] [,67] [,68] [,69] [,70] [,71] [1,] 53000 30678 30839 30972 26486 32708 30812 41652 42720 59648 52733 [2,] 50757 49342 39302 34924 41252 40237 37300 58900 64187 83651 78178 [3,] 71289 67604 67685 56818 63146 67845 63599 86188 98897 103569 103302 [4,] 75080 63733 67925 62024 70114 68913 72918 88350 108429 108749 96040 [5,] 77377 71209 66056 59274 57779 66697 74386 85787 92249 96948 98843 [6,] 68779 59461 60743 59007 61944 68272 73292 86535 111820 90113 93183 [7,] 49155 49609 49422 43494 48728 51317 56577 67978 77163 75614 75267 [,72] [,73] [,74] [,75] [,76] [,77] [,78] [,79] [,80] [,81] [1,] 54815 63092 59728 57138 58153 64347 61277 62291 67952 92622 [2,] 93904 82476 97028 63039 88938 88778 81275 97455 96574 114196 [3,] 110725 162042 148746 135236 132459 130777 82503 127733 148105 145515 [4,] 107735 108215 101380 104370 107654 159185 74333 101861 151576 157664 [5,] 100793 100846 100953 109844 113021 139988 77350 87763 132592 129522 [6,] 102795 104183 100205 104290 117587 136731 102662 80954 117213 121886 [7,] 73772 82156 73292 78231 85387 90273 76923 72918 102982 99271 [,82] [,83] [,84] [,85] [,86] [,87] [,88] [,89] [,90] [,91] [,92] [1,] 76923 47286 30732 39249 50623 44776 47206 52679 58420 33135 65015 [2,] 103649 66910 35751 48007 63573 62024 64400 76923 72677 62131 83758 [3,] 111232 78738 47579 73452 75347 86081 89579 99164 76148 77003 116225 [4,] 123621 62905 44002 91154 80714 95212 97722 100499 70568 81836 104984 [5,] 104077 44082 30064 71342 88644 88244 109203 101620 76469 75107 115451 [6,] 112354 25392 28809 63226 82370 80287 120551 94385 71449 112968 115718 [7,] 69260 31666 44375 58820 63439 61330 80234 75241 49368 84399 82930 [,93] [,94] [,95] [,96] [,97] [,98] [,99] [,100] [,101] [,102] [1,] 61170 59034 57378 57138 66430 58607 53694 64748 44723 63092 [2,] 88377 79593 76496 78738 85387 75828 66910 86642 62291 85574 [3,] 118254 117667 104370 123007 123728 108429 102261 84719 78365 115104 [4,] 124208 107788 102875 117293 114703 104157 104584 79566 127493 115638 [5,] 124502 108135 99591 147731 113769 99538 109550 93236 123568 107975 [6,] 109871 104958 97802 156409 96307 90353 99858 81649 114356 117533 [7,] 81782 78685 79299 85173 79433 69874 83598 61891 88911 80474 [,103] [,104] [,105] [,106] [,107] [,108] [,109] [,110] [,111] [,112] [1,] 58259 34123 31346 35564 46298 53106 47152 38582 34683 35484 [2,] 78231 45791 37861 42346 58313 63813 56043 54521 39730 47499 [3,] 70408 64454 53614 75054 92943 87843 79753 81702 68325 68405 [4,] 73986 61196 66857 81088 97055 95239 82530 81755 71049 77350 [5,] 78471 63039 73478 76255 93664 91661 76148 77136 68779 59461 [6,] 83037 59781 82183 71503 79993 85947 71583 64748 72758 67578 [7,] 52813 50997 51531 57992 67364 70408 59060 51932 51878 56524 [,113] [,114] [,115] [,116] [,117] [,118] [,119] [,120] [,121] [,122] [1,] 35030 36766 32467 35111 30732 39169 41679 41465 50303 59114 [2,] 41839 42320 41065 48487 47606 53000 52279 56871 66002 85200 [3,] 68112 66563 64187 71503 72117 78551 79379 83090 84746 101380 [4,] 61570 67498 69740 68726 80207 76148 90219 96280 88404 103382 [5,] 76335 71903 68539 66803 77136 78738 83571 94091 119936 111713 [6,] 68432 63813 59327 65735 67231 74119 86081 97161 97669 106373 [7,] 46832 53000 49876 54068 54468 54602 58019 70141 78712 75134 [,123] [,124] [,125] [,126] [,127] [,128] [,129] [,130] [,131] [,132] [1,] 63012 68592 65976 66216 66350 70007 58553 60876 64774 72143 [2,] 77056 93877 95880 98897 86722 101647 88564 87816 91955 87896 [3,] 100712 133206 150828 128080 106266 151336 114303 118495 115184 122286 [4,] 117800 147571 128614 140122 123381 149200 120337 114169 103783 131631 [5,] 118254 144687 122526 125223 128294 143646 123354 117720 96948 132806 [6,] 102448 111846 105946 122473 117080 106319 117133 120871 89979 139508 [7,] 79272 87069 82984 85467 84292 78284 101674 84826 84933 102074 [,133] [,134] [,135] [,136] [,137] [,138] [,139] [,140] [,141] [,142] [1,] 74466 78605 47366 29290 38822 56177 67177 56898 67151 74093 [2,] 104717 104584 53507 37594 50383 77243 62665 109497 85574 102021 [3,] 132272 127599 75614 48888 73051 102181 95212 168344 116999 146690 [4,] 130323 111579 60262 50757 68325 93877 103543 123915 136117 150508 [5,] 138066 123381 47980 56417 84078 90433 109871 120711 127065 135289 [6,] 119990 103890 28703 30465 93477 88617 105305 113128 111473 116012 [7,] 102261 61811 20986 30411 72090 65655 78738 83838 80100 78498 [,143] [,144] [,145] [,146] [,147] [,148] [,149] [,150] [,151] [,152] [1,] 59701 79139 77777 73185 189143 67017 75588 58446 70541 65228 [2,] 86294 109310 107040 102528 248550 92809 100606 70808 84559 101460 [3,] 107468 139935 149760 242516 269670 134755 120791 101220 120524 125303 [4,] 115371 126051 147731 293460 278801 137291 124075 127599 118922 141029 [5,] 118388 141270 159132 320987 272634 138199 101460 120497 123648 130136 [6,] 110485 118949 124369 268736 181240 125704 105198 110752 113448 124502 [7,] 100659 97455 98016 239873 90566 107414 70007 90113 84132 97241 [,153] [,154] [,155] [,156] [,157] [,158] [,159] [,160] [,161] [,162] [1,] 71556 70835 56043 50757 44536 58633 57512 36686 39436 29050 [2,] 97188 87736 81221 69473 64267 68672 88484 50703 51985 48647 [3,] 141430 128908 115024 117186 99458 93477 78979 75428 73585 70488 [4,] 134621 118575 91661 82637 89338 110485 84158 78979 67177 63199 [5,] 134354 118014 110298 91207 104397 155875 76389 78151 65442 62825 [6,] 113261 111232 99591 76682 90486 86828 77136 69260 61757 66803 [7,] 74520 79593 70488 66964 70702 74653 55189 56844 50890 49555 [,163] [,164] [,165] [,166] [,167] [,168] [,169] [,170] [,171] [,172] [1,] 30892 34310 37727 35885 41492 41305 46298 33695 34657 44829 [2,] 40103 44883 43254 49902 46805 39116 60422 52199 56417 46992 [3,] 59434 68058 74493 78231 69847 72010 81035 72517 68699 77030 [4,] 67604 74333 77991 67044 73478 79246 82637 72277 77163 92596 [5,] 70194 71156 68699 71476 78632 71690 79726 73078 79700 84933 [6,] 63760 67337 66483 68112 74813 69420 86348 76763 78151 83918 [7,] 53106 52279 62318 52492 60182 57939 66296 65655 59781 63306 [,173] [,174] [,175] [,176] [,177] [1,] 47286 45417 54174 96040 75695 [2,] 77403 79806 84559 141857 104157 [3,] 102768 106346 116492 233465 150161 [4,] 118361 95826 114543 146022 174832 [5,] 107708 93797 123407 142124 139935 [6,] 94705 99351 114677 118067 108749 [7,] 76122 83838 110378 83010 109977 > arr.mean <- array(NA,dim=np) > arr.sd <- array(NA,dim=np) > arr.range <- array(NA,dim=np) > for (j in 1:np) + { + arr.mean[j] <- mean(arr[,j],na.rm=TRUE) + arr.sd[j] <- sd(arr[,j],na.rm=TRUE) + arr.range[j] <- max(arr[,j],na.rm=TRUE) - min(arr[,j],na.rm=TRUE) + } > arr.mean [1] 36396.14 36128.71 30392.43 32040.14 38223.00 26387.29 26948.00 [8] 26406.43 24346.57 21611.86 26139.29 24201.57 26334.00 25788.43 [15] 24186.43 29221.14 41404.29 34988.43 38074.14 45256.57 66403.00 [22] 72673.71 71685.71 62874.71 73303.00 66227.57 81461.57 86241.29 [29] 74409.29 58095.43 25658.71 33558.14 51965.86 57870.29 60605.14 [36] 65964.29 67444.14 69832.14 72757.57 79440.29 75381.57 70183.00 [43] 64263.14 61455.86 68660.86 71399.71 88110.14 90234.71 84902.14 [50] 74722.00 71868.86 72090.00 68176.86 61364.29 57298.29 58312.71 [57] 56009.14 54231.57 56993.00 60277.00 63633.86 55948.00 54567.43 [64] 49501.86 52778.43 56569.86 58412.00 73627.14 85066.43 88327.43 [71] 85363.71 92077.00 100430.00 97333.14 93164.00 100457.00 115725.57 [78] 79474.71 90139.29 116713.43 122953.71 100159.43 50997.00 37330.29 [85] 63607.14 72101.43 73993.43 86985.00 85787.29 67871.29 75225.57 [92] 97725.86 101166.29 93694.29 88258.71 109355.57 97108.14 86683.71 [99] 88636.43 78921.57 91386.71 97912.86 70743.57 54197.29 56695.71 [106] 62828.86 76518.57 78288.14 67467.00 64339.43 58171.71 58900.14 [113] 56878.57 57409.00 55028.71 58633.29 61356.71 64903.86 70175.29 [120] 77014.14 83681.71 91756.57 94079.14 112406.86 107536.29 109496.86 [127] 101769.29 114348.43 103416.86 100681.86 92508.00 112620.57 114585.00 [134] 101635.57 47774.00 40546.00 68603.71 82026.14 88930.14 110904.43 [141] 103497.00 114730.14 99766.57 116015.57 123403.57 220183.57 218657.71 [148] 114741.29 99675.00 97062.14 102253.43 112128.43 109561.43 102127.57 [155] 89189.43 79272.29 80454.86 92660.43 73978.14 63721.57 58610.29 [162] 55795.29 55013.29 58908.00 61566.43 60448.86 63607.00 61532.29 [169] 71823.14 63740.57 64938.29 70514.86 89193.29 86340.14 102604.29 [176] 137226.43 123358.00 > arr.sd [1] 17399.541 15575.656 8315.523 11229.463 15042.573 6184.818 7798.593 [8] 7479.926 8320.421 7681.350 7346.183 6286.813 7657.952 7721.897 [15] 6503.880 8059.193 12806.502 8402.518 7194.501 17751.550 22659.651 [22] 22099.200 16138.771 15447.156 17595.782 16548.940 20855.407 20962.006 [29] 17297.255 16345.393 5251.754 12150.971 11670.877 15400.416 17682.454 [36] 16108.159 17453.188 18124.762 15737.241 19333.353 13725.835 17133.928 [43] 17318.197 18071.094 13406.418 18404.735 19992.891 19742.298 20297.918 [50] 14499.614 17970.332 14927.334 15608.377 6940.167 11730.510 16376.608 [57] 15346.067 16176.678 14170.843 15870.111 12203.084 13929.004 14961.993 [64] 12833.137 15042.043 15180.712 17910.685 18074.954 25153.105 16994.403 [71] 17780.791 20455.776 31414.737 27908.261 28065.461 24754.247 34542.060 [78] 12373.907 21472.158 30069.458 23507.727 19765.228 19388.482 7859.375 [85] 17201.377 13372.266 18240.581 25430.805 18284.442 10143.062 24147.456 [92] 20480.948 24434.461 21347.090 17630.281 36934.560 21163.515 18910.246 [99] 21298.090 11524.631 31797.914 21466.121 11216.795 11161.096 18565.309 [106] 17910.274 19689.575 15926.315 13441.389 16745.503 15945.821 14084.904 [113] 15624.319 13600.990 14328.677 13302.618 17930.513 15806.594 19144.130 [120] 21688.446 22304.181 19254.852 21351.587 30595.609 28900.992 26651.917 [127] 23136.073 33994.273 23190.631 23088.798 15781.303 25645.567 22360.227 [134] 23720.647 18463.441 11295.270 18748.729 16430.899 19263.911 34691.399 [141] 25997.182 31228.900 20623.648 22614.822 30435.771 95066.581 69167.848 [148] 27054.900 20559.571 25544.632 21763.342 25888.515 29147.220 22424.420 [155] 21366.731 21049.780 21452.307 32675.321 12773.723 16105.995 11728.564 [162] 14414.269 14678.889 15211.950 15363.036 14820.892 14553.409 15882.370 [169] 14647.484 15525.660 16187.633 19072.697 24063.430 20160.863 24649.072 [176] 48968.493 33355.870 > arr.range [1] 52946 43361 21868 31079 41465 17141 24056 18529 20186 18049 [11] 17783 17675 19624 20960 17301 23202 32761 23043 16794 51317 [21] 54628 70114 49528 45363 51184 46378 58073 56631 40371 42960 [31] 14631 30438 29317 40904 45096 41091 47633 49796 43334 53747 [41] 39276 42159 46565 49556 36232 48648 49262 52118 59568 46458 [51] 43627 42346 42159 16501 31399 41278 40504 38502 37460 37647 [61] 28222 40531 37086 31052 43628 36205 43574 46698 69100 49101 [71] 50569 55910 98950 89018 78098 74306 94838 41385 65442 83624 [81] 65042 54361 53346 18770 51905 38021 50436 73345 48941 27101 [91] 79833 51210 63332 58633 46992 99271 57298 49822 55856 31345 [101] 82770 54441 30224 30331 50837 45524 50757 42133 35378 43173 [111] 38075 41866 41305 35137 37273 36392 49475 39569 48540 55696 [121] 69633 52599 55242 78979 84852 73906 61944 81329 64801 59995 [131] 50410 67365 63600 65788 54628 27127 54655 46004 47206 111446 [141] 68966 76415 58687 62131 81355 247802 188235 71182 54068 69153 [151] 53107 75801 69874 58073 58981 66429 59861 97242 33295 42293 [161] 34149 41438 39302 40023 40264 42346 37140 40130 40050 43068 [171] 45043 47767 71075 60929 69233 150455 99137 > (lm1 <- lm(arr.sd~arr.mean)) Call: lm(formula = arr.sd ~ arr.mean) Coefficients: (Intercept) arr.mean -2567.5247 0.2882 > (lnlm1 <- lm(log(arr.sd)~log(arr.mean))) Call: lm(formula = log(arr.sd) ~ log(arr.mean)) Coefficients: (Intercept) log(arr.mean) -0.1798 0.8920 > (lm2 <- lm(arr.range~arr.mean)) Call: lm(formula = arr.range ~ arr.mean) Coefficients: (Intercept) arr.mean -7738.8970 0.7964 > postscript(file="/var/www/html/rcomp/tmp/1mshz1291029052.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.mean,arr.sd,main='Standard Deviation-Mean Plot',xlab='mean',ylab='standard deviation') > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/2mshz1291029052.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.mean,arr.range,main='Range-Mean Plot',xlab='mean',ylab='range') > dev.off() null device 1 > > #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'Standard Deviation-Mean Plot',4,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Section',header=TRUE) > a<-table.element(a,'Mean',header=TRUE) > a<-table.element(a,'Standard Deviation',header=TRUE) > a<-table.element(a,'Range',header=TRUE) > a<-table.row.end(a) > for (j in 1:np) { + a<-table.row.start(a) + a<-table.element(a,j,header=TRUE) + a<-table.element(a,arr.mean[j]) + a<-table.element(a,arr.sd[j] ) + a<-table.element(a,arr.range[j] ) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/3qbg51291029052.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'Regression: S.E.(k) = alpha + beta * Mean(k)',2,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'alpha',header=TRUE) > a<-table.element(a,lm1$coefficients[[1]]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'beta',header=TRUE) > a<-table.element(a,lm1$coefficients[[2]]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'S.D.',header=TRUE) > a<-table.element(a,summary(lm1)$coefficients[2,2]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'T-STAT',header=TRUE) > a<-table.element(a,summary(lm1)$coefficients[2,3]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'p-value',header=TRUE) > a<-table.element(a,summary(lm1)$coefficients[2,4]) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/4bbea1291029052.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,'Regression: ln S.E.(k) = alpha + beta * ln Mean(k)',2,TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'alpha',header=TRUE) > a<-table.element(a,lnlm1$coefficients[[1]]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'beta',header=TRUE) > a<-table.element(a,lnlm1$coefficients[[2]]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'S.D.',header=TRUE) > a<-table.element(a,summary(lnlm1)$coefficients[2,2]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'T-STAT',header=TRUE) > a<-table.element(a,summary(lnlm1)$coefficients[2,3]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'p-value',header=TRUE) > a<-table.element(a,summary(lnlm1)$coefficients[2,4]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Lambda',header=TRUE) > a<-table.element(a,1-lnlm1$coefficients[[2]]) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/573uj1291029052.tab") > > try(system("convert tmp/1mshz1291029052.ps tmp/1mshz1291029052.png",intern=TRUE)) character(0) > try(system("convert tmp/2mshz1291029052.ps tmp/2mshz1291029052.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.240 0.342 3.677