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 <- array(list(1 + ,162556 + ,162556 + ,807 + ,807 + ,213118 + ,213118 + ,6282154 + ,1 + ,29790 + ,29790 + ,444 + ,444 + ,81767 + ,81767 + ,4321023 + ,1 + ,87550 + ,87550 + ,412 + ,412 + ,153198 + ,153198 + ,4111912 + ,0 + ,84738 + ,0 + ,428 + ,0 + ,-26007 + ,0 + ,223193 + ,1 + ,54660 + ,54660 + ,315 + ,315 + ,126942 + ,126942 + ,1491348 + ,1 + ,42634 + ,42634 + ,168 + ,168 + ,157214 + ,157214 + ,1629616 + ,0 + ,40949 + ,0 + ,263 + ,0 + ,129352 + ,0 + ,1398893 + ,1 + ,45187 + ,45187 + ,267 + ,267 + ,234817 + ,234817 + ,1926517 + ,1 + ,37704 + ,37704 + ,228 + ,228 + ,60448 + ,60448 + ,983660 + ,1 + ,16275 + ,16275 + ,129 + ,129 + ,47818 + ,47818 + ,1443586 + ,0 + ,25830 + ,0 + ,104 + ,0 + ,245546 + ,0 + ,1073089 + ,0 + ,12679 + ,0 + ,122 + ,0 + ,48020 + ,0 + ,984885 + ,1 + ,18014 + ,18014 + ,393 + ,393 + ,-1710 + ,-1710 + ,1405225 + ,0 + ,43556 + ,0 + ,190 + ,0 + ,32648 + ,0 + ,227132 + ,1 + ,24811 + ,24811 + ,280 + ,280 + ,95350 + ,95350 + ,929118 + ,0 + ,6575 + ,0 + ,63 + ,0 + ,151352 + ,0 + ,1071292 + ,0 + ,7123 + ,0 + ,102 + ,0 + ,288170 + ,0 + ,638830 + ,1 + ,21950 + ,21950 + ,265 + ,265 + ,114337 + ,114337 + ,856956 + ,1 + ,37597 + ,37597 + ,234 + ,234 + ,37884 + ,37884 + ,992426 + ,0 + ,17821 + ,0 + ,277 + ,0 + ,122844 + ,0 + ,444477 + ,1 + ,12988 + ,12988 + ,73 + ,73 + ,82340 + ,82340 + ,857217 + ,1 + ,22330 + ,22330 + ,67 + ,67 + ,79801 + ,79801 + ,711969 + ,0 + ,13326 + ,0 + ,103 + ,0 + ,165548 + ,0 + ,702380 + ,0 + ,16189 + ,0 + ,290 + ,0 + ,116384 + ,0 + ,358589 + ,0 + ,7146 + ,0 + ,83 + ,0 + ,134028 + ,0 + ,297978 + ,0 + ,15824 + ,0 + ,56 + ,0 + ,63838 + ,0 + ,585715 + ,1 + ,27664 + ,27664 + ,236 + ,236 + ,74996 + ,74996 + ,657954 + ,0 + ,11920 + ,0 + ,73 + ,0 + ,31080 + ,0 + ,209458 + ,0 + ,8568 + ,0 + ,34 + ,0 + ,32168 + ,0 + ,786690 + ,0 + ,14416 + ,0 + ,139 + ,0 + ,49857 + ,0 + ,439798 + ,1 + ,3369 + ,3369 + ,26 + ,26 + ,87161 + ,87161 + ,688779 + ,1 + ,11819 + ,11819 + ,70 + ,70 + ,106113 + ,106113 + ,574339 + ,1 + ,6984 + ,6984 + ,40 + ,40 + ,80570 + ,80570 + ,741409 + ,1 + ,4519 + ,4519 + ,42 + ,42 + ,102129 + ,102129 + ,597793 + ,0 + ,2220 + ,0 + ,12 + ,0 + ,301670 + ,0 + ,644190 + ,0 + ,18562 + ,0 + ,211 + ,0 + ,102313 + ,0 + ,377934 + ,0 + ,10327 + ,0 + ,74 + ,0 + ,88577 + ,0 + ,640273 + ,1 + ,5336 + ,5336 + ,80 + ,80 + ,112477 + ,112477 + ,697458 + ,1 + ,2365 + ,2365 + ,83 + ,83 + ,191778 + ,191778 + ,550608 + ,0 + ,4069 + ,0 + ,131 + ,0 + ,79804 + ,0 + ,207393 + ,0 + ,8636 + ,0 + ,203 + ,0 + ,128294 + ,0 + ,301607 + ,0 + ,13718 + ,0 + ,56 + ,0 + ,96448 + ,0 + ,345783 + ,0 + ,4525 + ,0 + ,89 + ,0 + ,93811 + ,0 + ,501749 + ,0 + ,6869 + ,0 + ,88 + ,0 + ,117520 + ,0 + ,379983 + ,0 + ,4628 + ,0 + ,39 + ,0 + ,69159 + ,0 + ,387475 + ,1 + ,3689 + ,3689 + ,25 + ,25 + ,101792 + ,101792 + ,377305 + ,1 + ,4891 + ,4891 + ,49 + ,49 + ,210568 + ,210568 + ,370837 + ,1 + ,7489 + ,7489 + ,149 + ,149 + ,136996 + ,136996 + ,430866 + ,0 + ,4901 + ,0 + ,58 + ,0 + ,121920 + ,0 + ,469107 + ,0 + ,2284 + ,0 + ,41 + ,0 + ,76403 + ,0 + ,194493 + ,1 + ,3160 + ,3160 + ,90 + ,90 + ,108094 + ,108094 + ,530670 + ,1 + ,4150 + ,4150 + ,136 + ,136 + ,134759 + ,134759 + ,518365 + ,1 + ,7285 + ,7285 + ,97 + ,97 + ,188873 + ,188873 + ,491303 + ,1 + ,1134 + ,1134 + ,63 + ,63 + ,146216 + ,146216 + ,527021 + ,1 + ,4658 + ,4658 + ,114 + ,114 + ,156608 + ,156608 + ,233773 + ,0 + ,2384 + ,0 + ,77 + ,0 + ,61348 + ,0 + ,405972 + ,0 + ,3748 + ,0 + ,6 + ,0 + ,50350 + ,0 + ,652925 + ,0 + ,5371 + ,0 + ,47 + ,0 + ,87720 + ,0 + ,446211 + ,0 + ,1285 + ,0 + ,51 + ,0 + ,99489 + ,0 + ,341340 + ,1 + ,9327 + ,9327 + ,85 + ,85 + ,87419 + ,87419 + ,387699 + ,1 + ,5565 + ,5565 + ,43 + ,43 + ,94355 + ,94355 + ,493408 + ,0 + ,1528 + ,0 + ,32 + ,0 + ,60326 + ,0 + ,146494 + ,1 + ,3122 + ,3122 + ,25 + ,25 + ,94670 + ,94670 + ,414462 + ,1 + ,7561 + ,7561 + ,77 + ,77 + ,82425 + ,82425 + ,364304 + ,0 + ,2675 + ,0 + ,54 + ,0 + ,59017 + ,0 + ,355178 + ,0 + ,13253 + ,0 + ,251 + ,0 + ,90829 + ,0 + ,357760 + ,0 + ,880 + ,0 + ,15 + ,0 + ,80791 + ,0 + ,261216 + ,1 + ,2053 + ,2053 + ,44 + ,44 + ,100423 + ,100423 + ,397144 + ,0 + ,1424 + ,0 + ,73 + ,0 + ,131116 + ,0 + ,374943 + ,1 + ,4036 + ,4036 + ,85 + ,85 + ,100269 + ,100269 + ,424898 + ,1 + ,3045 + ,3045 + ,49 + ,49 + ,27330 + ,27330 + ,202055 + ,0 + ,5119 + ,0 + ,38 + ,0 + ,39039 + ,0 + ,378525 + ,0 + ,1431 + ,0 + ,35 + ,0 + ,106885 + ,0 + ,310768 + ,0 + ,554 + ,0 + ,9 + ,0 + ,79285 + ,0 + ,325738 + ,0 + ,1975 + ,0 + ,34 + ,0 + ,118881 + ,0 + ,394510 + ,1 + ,1765 + ,1765 + ,20 + ,20 + ,77623 + ,77623 + ,247060 + ,0 + ,1012 + ,0 + ,29 + ,0 + ,114768 + ,0 + ,368078 + ,0 + ,810 + ,0 + ,11 + ,0 + ,74015 + ,0 + ,236761 + ,0 + ,1280 + ,0 + ,52 + ,0 + ,69465 + ,0 + ,312378 + ,1 + ,666 + ,666 + ,13 + ,13 + ,117869 + ,117869 + ,339836 + ,0 + ,1380 + ,0 + ,29 + ,0 + ,60982 + ,0 + ,347385 + ,1 + ,4677 + ,4677 + ,66 + ,66 + ,90131 + ,90131 + ,426280 + ,0 + ,876 + ,0 + ,33 + ,0 + ,138971 + ,0 + ,352850 + ,0 + ,814 + ,0 + ,15 + ,0 + ,39625 + ,0 + ,301881 + ,0 + ,514 + ,0 + ,15 + ,0 + ,102725 + ,0 + ,377516 + ,1 + ,5692 + ,5692 + ,68 + ,68 + ,64239 + ,64239 + ,357312 + ,0 + ,3642 + ,0 + ,100 + ,0 + ,90262 + ,0 + ,458343 + ,0 + ,540 + ,0 + ,13 + ,0 + ,103960 + ,0 + ,354228 + ,0 + ,2099 + ,0 + ,45 + ,0 + ,106611 + ,0 + ,308636 + ,0 + ,567 + ,0 + ,14 + ,0 + ,103345 + ,0 + ,386212 + ,0 + ,2001 + ,0 + ,36 + ,0 + ,95551 + ,0 + ,393343 + ,1 + ,2949 + ,2949 + ,40 + ,40 + ,82903 + ,82903 + ,378509 + ,0 + ,2253 + ,0 + ,68 + ,0 + ,63593 + ,0 + ,452469 + ,1 + ,6533 + ,6533 + ,29 + ,29 + ,126910 + ,126910 + ,364839 + ,0 + ,1889 + ,0 + ,43 + ,0 + ,37527 + ,0 + ,358649 + ,1 + ,3055 + ,3055 + ,30 + ,30 + ,60247 + ,60247 + ,376641 + ,0 + ,272 + ,0 + ,9 + ,0 + ,112995 + ,0 + ,429112 + ,1 + ,1414 + ,1414 + ,22 + ,22 + ,70184 + ,70184 + ,330546 + ,0 + ,2564 + ,0 + ,19 + ,0 + ,130140 + ,0 + ,403560 + ,1 + ,1383 + ,1383 + ,9 + ,9 + ,73221 + ,73221 + ,317892 + ,1 + ,1261 + ,1261 + ,31 + ,31 + ,76114 + ,76114 + ,307528 + ,0 + ,975 + ,0 + ,19 + ,0 + ,90534 + ,0 + ,235133 + ,0 + ,3366 + ,0 + ,55 + ,0 + ,108479 + ,0 + ,299243 + ,0 + ,576 + ,0 + ,8 + ,0 + ,113761 + ,0 + ,314073 + ,0 + ,1686 + ,0 + ,28 + ,0 + ,68696 + ,0 + ,368186 + ,0 + ,746 + ,0 + ,29 + ,0 + ,71561 + ,0 + ,269661 + ,1 + ,3192 + ,3192 + ,48 + ,48 + ,59831 + ,59831 + ,125390 + ,1 + ,2045 + ,2045 + ,16 + ,16 + ,97890 + ,97890 + ,510834 + ,0 + ,5702 + ,0 + ,47 + ,0 + ,101481 + ,0 + ,321896 + ,1 + ,1932 + ,1932 + ,20 + ,20 + ,72954 + ,72954 + ,249898 + ,0 + ,936 + ,0 + ,22 + ,0 + ,67939 + ,0 + ,408881 + ,1 + ,3437 + ,3437 + ,33 + ,33 + ,48022 + ,48022 + ,158492 + ,0 + ,5131 + ,0 + ,44 + ,0 + ,86111 + ,0 + ,292154 + ,1 + ,2397 + ,2397 + ,13 + ,13 + ,74020 + ,74020 + ,289513 + ,1 + ,1389 + ,1389 + ,6 + ,6 + ,57530 + ,57530 + ,378049 + ,0 + ,1503 + ,0 + ,35 + ,0 + ,56364 + ,0 + ,343466 + ,0 + ,402 + ,0 + ,8 + ,0 + ,84990 + ,0 + ,332743 + ,0 + ,2239 + ,0 + ,17 + ,0 + ,88590 + ,0 + ,442882 + ,1 + ,2234 + ,2234 + ,11 + ,11 + ,77200 + ,77200 + ,214215 + ,0 + ,837 + ,0 + ,21 + ,0 + ,61262 + ,0 + ,315688 + ,0 + ,10579 + ,0 + ,92 + ,0 + ,110309 + ,0 + ,375195 + ,0 + ,875 + ,0 + ,12 + ,0 + ,67000 + ,0 + ,334280 + ,0 + ,1585 + ,0 + ,112 + ,0 + ,93099 + ,0 + ,355864 + ,1 + ,1659 + ,1659 + ,25 + ,25 + ,107577 + ,107577 + ,480382 + ,1 + ,2647 + ,2647 + ,17 + ,17 + ,62920 + ,62920 + ,353058 + ,1 + ,3294 + ,3294 + ,23 + ,23 + ,75832 + ,75832 + ,217193 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,94 + ,0 + ,10 + ,0 + ,60793 + ,0 + ,314533 + ,0 + ,422 + ,0 + ,23 + ,0 + ,57935 + ,0 + ,318056 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,34 + ,0 + ,7 + ,0 + ,60630 + ,0 + ,314353 + ,0 + ,1558 + ,0 + ,25 + ,0 + ,55637 + ,0 + ,369448 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,43 + ,0 + ,20 + ,0 + ,60887 + ,0 + ,312846 + ,0 + ,645 + ,0 + ,4 + ,0 + ,60720 + ,0 + ,312075 + ,0 + ,316 + ,0 + ,4 + ,0 + ,60505 + ,0 + ,315009 + ,0 + ,115 + ,0 + ,10 + ,0 + ,60945 + ,0 + ,318903 + ,0 + ,5 + ,0 + ,1 + ,0 + ,60720 + ,0 + ,314887 + ,0 + ,897 + ,0 + ,4 + ,0 + ,60720 + ,0 + ,314913 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,389 + ,0 + ,8 + ,0 + ,58990 + ,0 + ,325506 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,1002 + ,0 + ,11 + ,0 + ,56750 + ,0 + ,298568 + ,0 + ,36 + ,0 + ,4 + ,0 + ,60894 + ,0 + ,315834 + ,0 + ,460 + ,0 + ,15 + ,0 + ,63346 + ,0 + ,329784 + ,0 + ,309 + ,0 + ,9 + ,0 + ,56535 + ,0 + ,312878 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,9 + ,0 + ,7 + ,0 + ,60835 + ,0 + ,314987 + ,0 + ,271 + ,0 + ,2 + ,0 + ,60720 + ,0 + ,325249 + ,0 + ,14 + ,0 + ,0 + ,0 + ,61016 + ,0 + ,315877 + ,0 + ,520 + ,0 + ,7 + ,0 + ,58650 + ,0 + ,291650 + ,0 + ,1766 + ,0 + ,46 + ,0 + ,60438 + ,0 + ,305959 + ,0 + ,0 + ,0 + ,5 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,458 + ,0 + ,7 + ,0 + ,58625 + ,0 + ,297765 + ,0 + ,20 + ,0 + ,2 + ,0 + ,60938 + ,0 + ,315245 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,98 + ,0 + ,2 + ,0 + ,61490 + ,0 + ,315236 + ,0 + ,405 + ,0 + ,5 + ,0 + ,60845 + ,0 + ,336425 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,483 + ,0 + ,7 + ,0 + ,60830 + ,0 + ,306268 + ,0 + ,454 + ,0 + ,24 + ,0 + ,63261 + ,0 + ,302187 + ,0 + ,47 + ,0 + ,1 + ,0 + ,60720 + ,0 + ,314882 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,757 + ,0 + ,18 + ,0 + ,45689 + ,0 + ,382712 + ,0 + ,4655 + ,0 + ,55 + ,0 + ,60720 + ,0 + ,341570 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,36 + ,0 + ,3 + ,0 + ,61564 + ,0 + ,312412 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,203 + ,0 + ,9 + ,0 + ,61938 + ,0 + ,309596 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,126 + ,0 + ,8 + ,0 + ,60951 + ,0 + ,315547 + ,0 + ,400 + ,0 + ,113 + ,0 + ,60720 + ,0 + ,313267 + ,1 + ,71 + ,71 + ,0 + ,0 + ,60745 + ,60745 + ,316176 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,972 + ,0 + ,19 + ,0 + ,71642 + ,0 + ,359335 + ,1 + ,531 + ,531 + ,11 + ,11 + ,71641 + ,71641 + ,330068 + ,0 + ,2461 + ,0 + ,25 + ,0 + ,55792 + ,0 + ,314289 + ,1 + ,378 + ,378 + ,16 + ,16 + ,71873 + ,71873 + ,297413 + ,1 + ,23 + ,23 + ,5 + ,5 + ,62555 + ,62555 + ,314806 + ,1 + ,638 + ,638 + ,11 + ,11 + ,60370 + ,60370 + ,333210 + ,1 + ,2300 + ,2300 + ,23 + ,23 + ,64873 + ,64873 + ,352108 + ,0 + ,149 + ,0 + ,6 + ,0 + ,62041 + ,0 + ,313332 + ,0 + ,226 + ,0 + ,5 + ,0 + ,65745 + ,0 + ,291787 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,275 + ,0 + ,7 + ,0 + ,59500 + ,0 + ,318745 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,141 + ,0 + ,7 + ,0 + ,61630 + ,0 + ,315366 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,28 + ,0 + ,3 + ,0 + ,60890 + ,0 + ,315688 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,1 + ,4980 + ,4980 + ,89 + ,89 + ,113521 + ,113521 + ,409642 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,1 + ,472 + ,472 + ,19 + ,19 + ,80045 + ,80045 + ,269587 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,1 + ,203 + ,203 + ,12 + ,12 + ,50804 + ,50804 + ,300962 + ,1 + ,496 + ,496 + ,12 + ,12 + ,87390 + ,87390 + ,325479 + ,1 + ,10 + ,10 + ,5 + ,5 + ,61656 + ,61656 + ,316155 + ,1 + ,63 + ,63 + ,2 + ,2 + ,65688 + ,65688 + ,318574 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,1 + ,1136 + ,1136 + ,26 + ,26 + ,48522 + ,48522 + ,343613 + ,1 + ,265 + ,265 + ,3 + ,3 + ,60720 + ,60720 + ,306948 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,267 + ,0 + ,11 + ,0 + ,57640 + ,0 + ,330059 + ,0 + ,474 + ,0 + ,10 + ,0 + ,61977 + ,0 + ,288985 + ,0 + ,534 + ,0 + ,5 + ,0 + ,62620 + ,0 + ,304485 + ,0 + ,0 + ,0 + ,2 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,15 + ,0 + ,6 + ,0 + ,60831 + ,0 + ,315688 + ,0 + ,397 + ,0 + ,7 + ,0 + ,60646 + ,0 + ,317736 + ,0 + ,0 + ,0 + ,2 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,1866 + ,0 + ,28 + ,0 + ,56225 + ,0 + ,322331 + ,0 + ,288 + ,0 + ,3 + ,0 + ,60510 + ,0 + ,296656 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,3 + ,0 + ,1 + ,0 + ,60698 + ,0 + ,315354 + ,0 + ,468 + ,0 + ,20 + ,0 + ,60720 + ,0 + ,312161 + ,0 + ,20 + ,0 + ,1 + ,0 + ,60805 + ,0 + ,315576 + ,0 + ,278 + ,0 + ,22 + ,0 + ,61404 + ,0 + ,314922 + ,0 + ,61 + ,0 + ,9 + ,0 + ,60720 + ,0 + ,314551 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,192 + ,0 + ,2 + ,0 + ,65276 + ,0 + ,312339 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,317 + ,0 + ,7 + ,0 + ,63915 + ,0 + ,298700 + ,0 + ,738 + ,0 + ,9 + ,0 + ,60720 + ,0 + ,321376 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,368 + ,0 + ,13 + ,0 + ,61686 + ,0 + ,303230 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,2 + ,0 + ,0 + ,0 + ,60743 + ,0 + ,315487 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,53 + ,0 + ,6 + ,0 + ,60349 + ,0 + ,315793 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,94 + ,0 + ,3 + ,0 + ,61360 + ,0 + ,312887 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,24 + ,0 + ,7 + ,0 + ,59818 + ,0 + ,315637 + ,0 + ,2332 + ,0 + ,2 + ,0 + ,72680 + ,0 + ,324385 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,0 + ,131 + ,0 + ,15 + ,0 + ,61808 + ,0 + ,308989 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,0 + ,206 + ,0 + ,9 + ,0 + ,53110 + ,0 + ,296702 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,0 + ,167 + ,0 + ,1 + ,0 + ,64245 + ,0 + ,307322 + ,0 + ,622 + ,0 + ,38 + ,0 + ,73007 + ,0 + ,304376 + ,0 + ,2328 + ,0 + ,57 + ,0 + ,82732 + ,0 + ,253588 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,365 + ,0 + ,7 + ,0 + ,54820 + ,0 + ,309560 + ,0 + ,364 + ,0 + ,26 + ,0 + ,47705 + ,0 + ,298466 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,0 + ,226 + ,0 + ,13 + ,0 + ,72835 + ,0 + ,343929 + ,0 + ,307 + ,0 + ,10 + ,0 + ,58856 + ,0 + ,331955 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,0 + ,188 + ,0 + ,9 + ,0 + ,77655 + ,0 + ,381180 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,0 + ,138 + ,0 + ,26 + ,0 + ,69817 + ,0 + ,331420 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,0 + ,125 + ,0 + ,19 + ,0 + ,60798 + ,0 + ,310201 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,282 + ,0 + ,12 + ,0 + ,62452 + ,0 + ,320016 + ,0 + ,335 + ,0 + ,23 + ,0 + ,64175 + ,0 + ,320398 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,1324 + ,1324 + ,29 + ,29 + ,67440 + ,67440 + ,291841 + ,0 + ,176 + ,0 + ,8 + ,0 + ,68136 + ,0 + ,310670 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,249 + ,0 + ,26 + ,0 + ,56726 + ,0 + ,313491 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,333 + ,0 + ,9 + ,0 + ,70811 + ,0 + ,331323 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,601 + ,601 + ,5 + ,5 + ,60720 + ,60720 + ,319210 + ,0 + ,30 + ,0 + ,3 + ,0 + ,62045 + ,0 + ,318098 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,0 + ,249 + ,0 + ,13 + ,0 + ,54323 + ,0 + ,292754 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,0 + ,165 + ,0 + ,12 + ,0 + ,62841 + ,0 + ,325176 + ,0 + ,453 + ,0 + ,19 + ,0 + ,81125 + ,0 + ,365959 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,53 + ,0 + ,10 + ,0 + ,59506 + ,0 + ,302409 + ,1 + ,382 + ,382 + ,9 + ,9 + ,59365 + ,59365 + ,340968 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,9 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,30 + ,30 + ,4 + ,4 + ,60798 + ,60798 + ,313164 + ,0 + ,290 + ,0 + ,1 + ,0 + ,58790 + ,0 + ,301164 + ,0 + ,0 + ,0 + ,1 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,0 + ,366 + ,0 + ,14 + ,0 + ,61808 + ,0 + ,344425 + ,0 + ,2 + ,0 + ,12 + ,0 + ,60735 + ,0 + ,315394 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,1 + ,209 + ,209 + ,19 + ,19 + ,64016 + ,64016 + ,316647 + ,0 + ,384 + ,0 + ,17 + ,0 + ,54683 + ,0 + ,309836 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,0 + ,365 + ,0 + ,32 + ,0 + ,87192 + ,0 + ,346611 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,1 + ,49 + ,49 + ,14 + ,14 + ,64107 + ,64107 + ,322031 + ,0 + ,3 + ,0 + ,8 + ,0 + ,60761 + ,0 + ,315656 + ,0 + ,133 + ,0 + ,4 + ,0 + ,65990 + ,0 + ,339445 + ,0 + ,32 + ,0 + ,0 + ,0 + ,59988 + ,0 + ,314964 + ,0 + ,368 + ,0 + ,20 + ,0 + ,61167 + ,0 + ,297141 + ,0 + ,1 + ,0 + ,5 + ,0 + ,60719 + ,0 + ,315372 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,22 + ,0 + ,1 + ,0 + ,60722 + ,0 + ,312502 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,96 + ,0 + ,4 + ,0 + ,60379 + ,0 + ,313729 + ,0 + ,1 + ,0 + ,1 + ,0 + ,60727 + ,0 + ,315388 + ,0 + ,314 + ,0 + ,4 + ,0 + ,60720 + ,0 + ,315371 + ,0 + ,844 + ,0 + ,20 + ,0 + ,60925 + ,0 + ,296139 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,26 + ,0 + ,1 + ,0 + ,60896 + ,0 + ,313880 + ,0 + ,125 + ,0 + ,10 + ,0 + ,59734 + ,0 + ,317698 + ,0 + ,304 + ,0 + ,12 + ,0 + ,62969 + ,0 + ,295580 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,621 + ,0 + ,13 + ,0 + ,60720 + ,0 + ,308256 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,119 + ,0 + ,3 + ,0 + ,59118 + ,0 + ,303677 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,1595 + ,0 + ,10 + ,0 + ,60720 + ,0 + ,319369 + ,0 + ,312 + ,0 + ,3 + ,0 + ,58598 + ,0 + ,318690 + ,0 + ,60 + ,0 + ,7 + ,0 + ,61124 + ,0 + ,314049 + ,0 + ,587 + ,0 + ,10 + ,0 + ,59595 + ,0 + ,325699 + ,0 + ,135 + ,0 + ,1 + ,0 + ,62065 + ,0 + ,314210 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,514 + ,0 + ,15 + ,0 + ,78780 + ,0 + ,322378 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,1 + ,0 + ,4 + ,0 + ,60722 + ,0 + ,315398 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,1763 + ,1763 + ,28 + ,28 + ,61600 + ,61600 + ,308336 + ,0 + ,180 + ,0 + ,9 + ,0 + ,59635 + ,0 + ,316386 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,218 + ,0 + ,7 + ,0 + ,60720 + ,0 + ,315553 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,448 + ,0 + ,7 + ,0 + ,59781 + ,0 + ,323361 + ,0 + ,227 + ,0 + ,7 + ,0 + ,76644 + ,0 + ,336639 + ,0 + ,174 + ,0 + ,3 + ,0 + ,64820 + ,0 + ,307424 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,121 + ,0 + ,11 + ,0 + ,56178 + ,0 + ,295370 + ,0 + ,607 + ,0 + ,7 + ,0 + ,60436 + ,0 + ,322340 + ,0 + ,2212 + ,0 + ,10 + ,0 + ,60720 + ,0 + ,319864 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,530 + ,0 + ,18 + ,0 + ,73433 + ,0 + ,317291 + ,0 + ,571 + ,0 + ,14 + ,0 + ,41477 + ,0 + ,280398 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,78 + ,0 + ,12 + ,0 + ,62700 + ,0 + ,317330 + ,0 + ,2489 + ,0 + ,29 + ,0 + ,67804 + ,0 + ,238125 + ,0 + ,131 + ,0 + ,3 + ,0 + ,59661 + ,0 + ,327071 + ,0 + ,923 + ,0 + ,6 + ,0 + ,58620 + ,0 + ,309038 + ,0 + ,72 + ,0 + ,3 + ,0 + ,60398 + ,0 + ,314210 + ,0 + ,572 + ,0 + ,8 + ,0 + ,58580 + ,0 + ,307930 + ,0 + ,397 + ,0 + ,10 + ,0 + ,62710 + ,0 + ,322327 + ,0 + ,450 + ,0 + ,6 + ,0 + ,59325 + ,0 + ,292136 + ,0 + ,622 + ,0 + ,8 + ,0 + ,60950 + ,0 + ,263276 + ,0 + ,694 + ,0 + ,6 + ,0 + ,68060 + ,0 + ,367655 + ,1 + ,3425 + ,3425 + ,9 + ,9 + ,83620 + ,83620 + ,283910 + ,0 + ,562 + ,0 + ,8 + ,0 + ,58456 + ,0 + ,283587 + ,0 + ,4917 + ,0 + ,26 + ,0 + ,52811 + ,0 + ,243650 + ,1 + ,1442 + ,1442 + ,239 + ,239 + ,121173 + ,121173 + ,438493 + ,0 + ,529 + ,0 + ,7 + ,0 + ,63870 + ,0 + ,296261 + ,1 + ,2126 + ,2126 + ,41 + ,41 + ,21001 + ,21001 + ,230621 + ,0 + ,1061 + ,0 + ,3 + ,0 + ,70415 + ,0 + ,304252 + ,0 + ,776 + ,0 + ,8 + ,0 + ,64230 + ,0 + ,333505 + ,0 + ,611 + ,0 + ,6 + ,0 + ,59190 + ,0 + ,296919 + ,1 + ,1526 + ,1526 + ,21 + ,21 + ,69351 + ,69351 + ,278990 + ,0 + ,592 + ,0 + ,7 + ,0 + ,64270 + ,0 + ,276898 + ,0 + ,1182 + ,0 + ,11 + ,0 + ,70694 + ,0 + ,327007 + ,0 + ,621 + ,0 + ,11 + ,0 + ,68005 + ,0 + ,317046 + ,0 + ,989 + ,0 + ,12 + ,0 + ,58930 + ,0 + ,304555 + ,0 + ,438 + ,0 + ,9 + ,0 + ,58320 + ,0 + ,298096 + ,0 + ,726 + ,0 + ,3 + ,0 + ,69980 + ,0 + ,231861 + ,0 + ,1303 + ,0 + ,57 + ,0 + ,69863 + ,0 + ,309422 + ,1 + ,7419 + ,7419 + ,21 + ,21 + ,63255 + ,63255 + ,286963 + ,1 + ,1164 + ,1164 + ,15 + ,15 + ,57320 + ,57320 + ,269753 + ,1 + ,3310 + ,3310 + ,32 + ,32 + ,75230 + ,75230 + ,448243 + ,0 + ,1920 + ,0 + ,11 + ,0 + ,79420 + ,0 + ,165404 + ,0 + ,965 + ,0 + ,2 + ,0 + ,73490 + ,0 + ,204325 + ,0 + ,3256 + ,0 + ,23 + ,0 + ,35250 + ,0 + ,407159 + ,1 + ,1135 + ,1135 + ,20 + ,20 + ,62285 + ,62285 + ,290476 + ,0 + ,1270 + ,0 + ,24 + ,0 + ,69206 + ,0 + ,275311 + ,0 + ,661 + ,0 + ,1 + ,0 + ,65920 + ,0 + ,246541 + ,0 + ,1013 + ,0 + ,1 + ,0 + ,69770 + ,0 + ,253468 + ,0 + ,2844 + ,0 + ,74 + ,0 + ,72683 + ,0 + ,240897 + ,1 + ,11528 + ,11528 + ,68 + ,68 + ,-14545 + ,-14545 + ,-83265 + ,0 + ,6526 + ,0 + ,20 + ,0 + ,55830 + ,0 + ,-42143 + ,0 + ,2264 + ,0 + ,20 + ,0 + ,55174 + ,0 + ,272713 + ,1 + ,5109 + ,5109 + ,82 + ,82 + ,67038 + ,67038 + ,215362 + ,0 + ,3999 + ,0 + ,21 + ,0 + ,51252 + ,0 + ,42754 + ,0 + ,35624 + ,0 + ,244 + ,0 + ,157278 + ,0 + ,306275 + ,0 + ,9252 + ,0 + ,32 + ,0 + ,79510 + ,0 + ,253537 + ,0 + ,15236 + ,0 + ,86 + ,0 + ,77440 + ,0 + ,372631 + ,0 + ,18073 + ,0 + ,69 + ,0 + ,27284 + ,0 + ,-7170) + ,dim=c(8 + ,431) + ,dimnames=list(c('Group' + ,'Costs' + ,'Costs_g' + ,'Orders' + ,'Orders_g' + ,'Dividends' + ,'Dividends_g' + ,'Wealth') + ,1:431)) > y <- array(NA,dim=c(8,431),dimnames=list(c('Group','Costs','Costs_g','Orders','Orders_g','Dividends','Dividends_g','Wealth'),1:431)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par3 = 'No Linear Trend' > par2 = 'Do not include Seasonal Dummies' > par1 = '8' > #'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!) > library(lattice) > library(lmtest) Loading required package: zoo Attaching package: 'zoo' The following object(s) are masked from package:base : as.Date.numeric > n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test > par1 <- as.numeric(par1) > x <- t(y) > k <- length(x[1,]) > n <- length(x[,1]) > x1 <- cbind(x[,par1], x[,1:k!=par1]) > mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1]) > colnames(x1) <- mycolnames #colnames(x)[par1] > x <- x1 > if (par3 == 'First Differences'){ + x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep=''))) + for (i in 1:n-1) { + for (j in 1:k) { + x2[i,j] <- x[i+1,j] - x[i,j] + } + } + x <- x2 + } > if (par2 == 'Include Monthly Dummies'){ + x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep =''))) + for (i in 1:11){ + x2[seq(i,n,12),i] <- 1 + } + x <- cbind(x, x2) + } > if (par2 == 'Include Quarterly Dummies'){ + x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep =''))) + for (i in 1:3){ + x2[seq(i,n,4),i] <- 1 + } + x <- cbind(x, x2) + } > k <- length(x[1,]) > if (par3 == 'Linear Trend'){ + x <- cbind(x, c(1:n)) + colnames(x)[k+1] <- 't' + } > x Wealth Group Costs Costs_g Orders Orders_g Dividends Dividends_g 1 6282154 1 162556 162556 807 807 213118 213118 2 4321023 1 29790 29790 444 444 81767 81767 3 4111912 1 87550 87550 412 412 153198 153198 4 223193 0 84738 0 428 0 -26007 0 5 1491348 1 54660 54660 315 315 126942 126942 6 1629616 1 42634 42634 168 168 157214 157214 7 1398893 0 40949 0 263 0 129352 0 8 1926517 1 45187 45187 267 267 234817 234817 9 983660 1 37704 37704 228 228 60448 60448 10 1443586 1 16275 16275 129 129 47818 47818 11 1073089 0 25830 0 104 0 245546 0 12 984885 0 12679 0 122 0 48020 0 13 1405225 1 18014 18014 393 393 -1710 -1710 14 227132 0 43556 0 190 0 32648 0 15 929118 1 24811 24811 280 280 95350 95350 16 1071292 0 6575 0 63 0 151352 0 17 638830 0 7123 0 102 0 288170 0 18 856956 1 21950 21950 265 265 114337 114337 19 992426 1 37597 37597 234 234 37884 37884 20 444477 0 17821 0 277 0 122844 0 21 857217 1 12988 12988 73 73 82340 82340 22 711969 1 22330 22330 67 67 79801 79801 23 702380 0 13326 0 103 0 165548 0 24 358589 0 16189 0 290 0 116384 0 25 297978 0 7146 0 83 0 134028 0 26 585715 0 15824 0 56 0 63838 0 27 657954 1 27664 27664 236 236 74996 74996 28 209458 0 11920 0 73 0 31080 0 29 786690 0 8568 0 34 0 32168 0 30 439798 0 14416 0 139 0 49857 0 31 688779 1 3369 3369 26 26 87161 87161 32 574339 1 11819 11819 70 70 106113 106113 33 741409 1 6984 6984 40 40 80570 80570 34 597793 1 4519 4519 42 42 102129 102129 35 644190 0 2220 0 12 0 301670 0 36 377934 0 18562 0 211 0 102313 0 37 640273 0 10327 0 74 0 88577 0 38 697458 1 5336 5336 80 80 112477 112477 39 550608 1 2365 2365 83 83 191778 191778 40 207393 0 4069 0 131 0 79804 0 41 301607 0 8636 0 203 0 128294 0 42 345783 0 13718 0 56 0 96448 0 43 501749 0 4525 0 89 0 93811 0 44 379983 0 6869 0 88 0 117520 0 45 387475 0 4628 0 39 0 69159 0 46 377305 1 3689 3689 25 25 101792 101792 47 370837 1 4891 4891 49 49 210568 210568 48 430866 1 7489 7489 149 149 136996 136996 49 469107 0 4901 0 58 0 121920 0 50 194493 0 2284 0 41 0 76403 0 51 530670 1 3160 3160 90 90 108094 108094 52 518365 1 4150 4150 136 136 134759 134759 53 491303 1 7285 7285 97 97 188873 188873 54 527021 1 1134 1134 63 63 146216 146216 55 233773 1 4658 4658 114 114 156608 156608 56 405972 0 2384 0 77 0 61348 0 57 652925 0 3748 0 6 0 50350 0 58 446211 0 5371 0 47 0 87720 0 59 341340 0 1285 0 51 0 99489 0 60 387699 1 9327 9327 85 85 87419 87419 61 493408 1 5565 5565 43 43 94355 94355 62 146494 0 1528 0 32 0 60326 0 63 414462 1 3122 3122 25 25 94670 94670 64 364304 1 7561 7561 77 77 82425 82425 65 355178 0 2675 0 54 0 59017 0 66 357760 0 13253 0 251 0 90829 0 67 261216 0 880 0 15 0 80791 0 68 397144 1 2053 2053 44 44 100423 100423 69 374943 0 1424 0 73 0 131116 0 70 424898 1 4036 4036 85 85 100269 100269 71 202055 1 3045 3045 49 49 27330 27330 72 378525 0 5119 0 38 0 39039 0 73 310768 0 1431 0 35 0 106885 0 74 325738 0 554 0 9 0 79285 0 75 394510 0 1975 0 34 0 118881 0 76 247060 1 1765 1765 20 20 77623 77623 77 368078 0 1012 0 29 0 114768 0 78 236761 0 810 0 11 0 74015 0 79 312378 0 1280 0 52 0 69465 0 80 339836 1 666 666 13 13 117869 117869 81 347385 0 1380 0 29 0 60982 0 82 426280 1 4677 4677 66 66 90131 90131 83 352850 0 876 0 33 0 138971 0 84 301881 0 814 0 15 0 39625 0 85 377516 0 514 0 15 0 102725 0 86 357312 1 5692 5692 68 68 64239 64239 87 458343 0 3642 0 100 0 90262 0 88 354228 0 540 0 13 0 103960 0 89 308636 0 2099 0 45 0 106611 0 90 386212 0 567 0 14 0 103345 0 91 393343 0 2001 0 36 0 95551 0 92 378509 1 2949 2949 40 40 82903 82903 93 452469 0 2253 0 68 0 63593 0 94 364839 1 6533 6533 29 29 126910 126910 95 358649 0 1889 0 43 0 37527 0 96 376641 1 3055 3055 30 30 60247 60247 97 429112 0 272 0 9 0 112995 0 98 330546 1 1414 1414 22 22 70184 70184 99 403560 0 2564 0 19 0 130140 0 100 317892 1 1383 1383 9 9 73221 73221 101 307528 1 1261 1261 31 31 76114 76114 102 235133 0 975 0 19 0 90534 0 103 299243 0 3366 0 55 0 108479 0 104 314073 0 576 0 8 0 113761 0 105 368186 0 1686 0 28 0 68696 0 106 269661 0 746 0 29 0 71561 0 107 125390 1 3192 3192 48 48 59831 59831 108 510834 1 2045 2045 16 16 97890 97890 109 321896 0 5702 0 47 0 101481 0 110 249898 1 1932 1932 20 20 72954 72954 111 408881 0 936 0 22 0 67939 0 112 158492 1 3437 3437 33 33 48022 48022 113 292154 0 5131 0 44 0 86111 0 114 289513 1 2397 2397 13 13 74020 74020 115 378049 1 1389 1389 6 6 57530 57530 116 343466 0 1503 0 35 0 56364 0 117 332743 0 402 0 8 0 84990 0 118 442882 0 2239 0 17 0 88590 0 119 214215 1 2234 2234 11 11 77200 77200 120 315688 0 837 0 21 0 61262 0 121 375195 0 10579 0 92 0 110309 0 122 334280 0 875 0 12 0 67000 0 123 355864 0 1585 0 112 0 93099 0 124 480382 1 1659 1659 25 25 107577 107577 125 353058 1 2647 2647 17 17 62920 62920 126 217193 1 3294 3294 23 23 75832 75832 127 315380 0 0 0 0 0 60720 0 128 314533 0 94 0 10 0 60793 0 129 318056 0 422 0 23 0 57935 0 130 315380 0 0 0 0 0 60720 0 131 314353 0 34 0 7 0 60630 0 132 369448 0 1558 0 25 0 55637 0 133 315380 0 0 0 0 0 60720 0 134 312846 0 43 0 20 0 60887 0 135 312075 0 645 0 4 0 60720 0 136 315009 0 316 0 4 0 60505 0 137 318903 0 115 0 10 0 60945 0 138 314887 0 5 0 1 0 60720 0 139 314913 0 897 0 4 0 60720 0 140 315380 0 0 0 0 0 60720 0 141 325506 0 389 0 8 0 58990 0 142 315380 0 0 0 0 0 60720 0 143 298568 0 1002 0 11 0 56750 0 144 315834 0 36 0 4 0 60894 0 145 329784 0 460 0 15 0 63346 0 146 312878 0 309 0 9 0 56535 0 147 315380 0 0 0 0 0 60720 0 148 314987 0 9 0 7 0 60835 0 149 325249 0 271 0 2 0 60720 0 150 315877 0 14 0 0 0 61016 0 151 291650 0 520 0 7 0 58650 0 152 305959 0 1766 0 46 0 60438 0 153 315380 0 0 0 5 0 60720 0 154 297765 0 458 0 7 0 58625 0 155 315245 0 20 0 2 0 60938 0 156 315380 0 0 0 0 0 60720 0 157 315380 0 0 0 0 0 60720 0 158 315236 0 98 0 2 0 61490 0 159 336425 0 405 0 5 0 60845 0 160 315380 0 0 0 0 0 60720 0 161 315380 0 0 0 0 0 60720 0 162 315380 0 0 0 0 0 60720 0 163 315380 0 0 0 0 0 60720 0 164 306268 0 483 0 7 0 60830 0 165 302187 0 454 0 24 0 63261 0 166 314882 0 47 0 1 0 60720 0 167 315380 0 0 0 0 0 60720 0 168 382712 0 757 0 18 0 45689 0 169 341570 0 4655 0 55 0 60720 0 170 315380 0 0 0 0 0 60720 0 171 315380 0 0 0 0 0 60720 0 172 312412 0 36 0 3 0 61564 0 173 315380 0 0 0 0 0 60720 0 174 309596 0 203 0 9 0 61938 0 175 315380 0 0 0 0 0 60720 0 176 315547 0 126 0 8 0 60951 0 177 313267 0 400 0 113 0 60720 0 178 316176 1 71 71 0 0 60745 60745 179 315380 0 0 0 0 0 60720 0 180 315380 0 0 0 0 0 60720 0 181 359335 0 972 0 19 0 71642 0 182 330068 1 531 531 11 11 71641 71641 183 314289 0 2461 0 25 0 55792 0 184 297413 1 378 378 16 16 71873 71873 185 314806 1 23 23 5 5 62555 62555 186 333210 1 638 638 11 11 60370 60370 187 352108 1 2300 2300 23 23 64873 64873 188 313332 0 149 0 6 0 62041 0 189 291787 0 226 0 5 0 65745 0 190 315380 0 0 0 0 0 60720 0 191 318745 0 275 0 7 0 59500 0 192 315380 0 0 0 0 0 60720 0 193 315366 0 141 0 7 0 61630 0 194 315380 0 0 0 0 0 60720 0 195 315688 0 28 0 3 0 60890 0 196 315380 1 0 0 0 0 60720 60720 197 409642 1 4980 4980 89 89 113521 113521 198 315380 1 0 0 0 0 60720 60720 199 315380 1 0 0 0 0 60720 60720 200 269587 1 472 472 19 19 80045 80045 201 315380 1 0 0 0 0 60720 60720 202 315380 1 0 0 0 0 60720 60720 203 315380 1 0 0 0 0 60720 60720 204 300962 1 203 203 12 12 50804 50804 205 325479 1 496 496 12 12 87390 87390 206 316155 1 10 10 5 5 61656 61656 207 318574 1 63 63 2 2 65688 65688 208 315380 1 0 0 0 0 60720 60720 209 343613 1 1136 1136 26 26 48522 48522 210 306948 1 265 265 3 3 60720 60720 211 315380 0 0 0 0 0 60720 0 212 315380 0 0 0 0 0 60720 0 213 330059 0 267 0 11 0 57640 0 214 288985 0 474 0 10 0 61977 0 215 304485 0 534 0 5 0 62620 0 216 315380 0 0 0 2 0 60720 0 217 315688 0 15 0 6 0 60831 0 218 317736 0 397 0 7 0 60646 0 219 315380 0 0 0 2 0 60720 0 220 322331 0 1866 0 28 0 56225 0 221 296656 0 288 0 3 0 60510 0 222 315380 0 0 0 0 0 60720 0 223 315354 0 3 0 1 0 60698 0 224 312161 0 468 0 20 0 60720 0 225 315576 0 20 0 1 0 60805 0 226 314922 0 278 0 22 0 61404 0 227 314551 0 61 0 9 0 60720 0 228 315380 0 0 0 0 0 60720 0 229 312339 0 192 0 2 0 65276 0 230 315380 0 0 0 0 0 60720 0 231 298700 0 317 0 7 0 63915 0 232 321376 0 738 0 9 0 60720 0 233 315380 0 0 0 0 0 60720 0 234 303230 0 368 0 13 0 61686 0 235 315380 0 0 0 0 0 60720 0 236 315487 0 2 0 0 0 60743 0 237 315380 0 0 0 0 0 60720 0 238 315793 0 53 0 6 0 60349 0 239 315380 0 0 0 0 0 60720 0 240 315380 0 0 0 0 0 60720 0 241 315380 0 0 0 0 0 60720 0 242 312887 0 94 0 3 0 61360 0 243 315380 0 0 0 0 0 60720 0 244 315637 0 24 0 7 0 59818 0 245 324385 0 2332 0 2 0 72680 0 246 315380 1 0 0 0 0 60720 60720 247 315380 1 0 0 0 0 60720 60720 248 308989 0 131 0 15 0 61808 0 249 315380 1 0 0 0 0 60720 60720 250 315380 1 0 0 0 0 60720 60720 251 296702 0 206 0 9 0 53110 0 252 315380 1 0 0 0 0 60720 60720 253 307322 0 167 0 1 0 64245 0 254 304376 0 622 0 38 0 73007 0 255 253588 0 2328 0 57 0 82732 0 256 315380 0 0 0 0 0 60720 0 257 309560 0 365 0 7 0 54820 0 258 298466 0 364 0 26 0 47705 0 259 315380 1 0 0 0 0 60720 60720 260 315380 0 0 0 0 0 60720 0 261 315380 1 0 0 0 0 60720 60720 262 315380 1 0 0 0 0 60720 60720 263 343929 0 226 0 13 0 72835 0 264 331955 0 307 0 10 0 58856 0 265 315380 1 0 0 0 0 60720 60720 266 315380 0 0 0 0 0 60720 0 267 315380 1 0 0 0 0 60720 60720 268 381180 0 188 0 9 0 77655 0 269 315380 1 0 0 0 0 60720 60720 270 331420 0 138 0 26 0 69817 0 271 315380 1 0 0 0 0 60720 60720 272 315380 1 0 0 0 0 60720 60720 273 315380 1 0 0 0 0 60720 60720 274 310201 0 125 0 19 0 60798 0 275 315380 0 0 0 0 0 60720 0 276 320016 0 282 0 12 0 62452 0 277 320398 0 335 0 23 0 64175 0 278 315380 0 0 0 0 0 60720 0 279 291841 1 1324 1324 29 29 67440 67440 280 310670 0 176 0 8 0 68136 0 281 315380 1 0 0 0 0 60720 60720 282 315380 0 0 0 0 0 60720 0 283 313491 0 249 0 26 0 56726 0 284 315380 0 0 0 0 0 60720 0 285 331323 0 333 0 9 0 70811 0 286 315380 0 0 0 0 0 60720 0 287 319210 1 601 601 5 5 60720 60720 288 318098 0 30 0 3 0 62045 0 289 315380 1 0 0 0 0 60720 60720 290 292754 0 249 0 13 0 54323 0 291 315380 1 0 0 0 0 60720 60720 292 325176 0 165 0 12 0 62841 0 293 365959 0 453 0 19 0 81125 0 294 315380 0 0 0 0 0 60720 0 295 302409 0 53 0 10 0 59506 0 296 340968 1 382 382 9 9 59365 59365 297 315380 0 0 0 0 0 60720 0 298 315380 1 0 0 0 0 60720 60720 299 315380 0 0 0 0 0 60720 0 300 315380 0 0 0 9 0 60720 0 301 313164 1 30 30 4 4 60798 60798 302 301164 0 290 0 1 0 58790 0 303 315380 0 0 0 1 0 60720 0 304 315380 1 0 0 0 0 60720 60720 305 344425 0 366 0 14 0 61808 0 306 315394 0 2 0 12 0 60735 0 307 315380 1 0 0 0 0 60720 60720 308 316647 1 209 209 19 19 64016 64016 309 309836 0 384 0 17 0 54683 0 310 315380 1 0 0 0 0 60720 60720 311 315380 1 0 0 0 0 60720 60720 312 346611 0 365 0 32 0 87192 0 313 315380 1 0 0 0 0 60720 60720 314 322031 1 49 49 14 14 64107 64107 315 315656 0 3 0 8 0 60761 0 316 339445 0 133 0 4 0 65990 0 317 314964 0 32 0 0 0 59988 0 318 297141 0 368 0 20 0 61167 0 319 315372 0 1 0 5 0 60719 0 320 315380 0 0 0 0 0 60720 0 321 315380 0 0 0 0 0 60720 0 322 315380 0 0 0 0 0 60720 0 323 315380 0 0 0 0 0 60720 0 324 315380 1 0 0 0 0 60720 60720 325 315380 0 0 0 0 0 60720 0 326 312502 0 22 0 1 0 60722 0 327 315380 1 0 0 0 0 60720 60720 328 315380 1 0 0 0 0 60720 60720 329 315380 0 0 0 0 0 60720 0 330 315380 0 0 0 0 0 60720 0 331 315380 0 0 0 0 0 60720 0 332 315380 0 0 0 0 0 60720 0 333 315380 0 0 0 0 0 60720 0 334 313729 0 96 0 4 0 60379 0 335 315388 0 1 0 1 0 60727 0 336 315371 0 314 0 4 0 60720 0 337 296139 0 844 0 20 0 60925 0 338 315380 0 0 0 0 0 60720 0 339 313880 0 26 0 1 0 60896 0 340 317698 0 125 0 10 0 59734 0 341 295580 0 304 0 12 0 62969 0 342 315380 0 0 0 0 0 60720 0 343 315380 0 0 0 0 0 60720 0 344 315380 0 0 0 0 0 60720 0 345 308256 0 621 0 13 0 60720 0 346 315380 0 0 0 0 0 60720 0 347 303677 0 119 0 3 0 59118 0 348 315380 0 0 0 0 0 60720 0 349 315380 0 0 0 0 0 60720 0 350 319369 0 1595 0 10 0 60720 0 351 318690 0 312 0 3 0 58598 0 352 314049 0 60 0 7 0 61124 0 353 325699 0 587 0 10 0 59595 0 354 314210 0 135 0 1 0 62065 0 355 315380 0 0 0 0 0 60720 0 356 315380 0 0 0 0 0 60720 0 357 322378 0 514 0 15 0 78780 0 358 315380 0 0 0 0 0 60720 0 359 315380 0 0 0 0 0 60720 0 360 315380 0 0 0 0 0 60720 0 361 315398 0 1 0 4 0 60722 0 362 315380 0 0 0 0 0 60720 0 363 315380 0 0 0 0 0 60720 0 364 308336 1 1763 1763 28 28 61600 61600 365 316386 0 180 0 9 0 59635 0 366 315380 0 0 0 0 0 60720 0 367 315380 0 0 0 0 0 60720 0 368 315380 0 0 0 0 0 60720 0 369 315380 0 0 0 0 0 60720 0 370 315553 0 218 0 7 0 60720 0 371 315380 0 0 0 0 0 60720 0 372 323361 0 448 0 7 0 59781 0 373 336639 0 227 0 7 0 76644 0 374 307424 0 174 0 3 0 64820 0 375 315380 0 0 0 0 0 60720 0 376 315380 0 0 0 0 0 60720 0 377 295370 0 121 0 11 0 56178 0 378 322340 0 607 0 7 0 60436 0 379 319864 0 2212 0 10 0 60720 0 380 315380 0 0 0 0 0 60720 0 381 315380 0 0 0 0 0 60720 0 382 317291 0 530 0 18 0 73433 0 383 280398 0 571 0 14 0 41477 0 384 315380 0 0 0 0 0 60720 0 385 317330 0 78 0 12 0 62700 0 386 238125 0 2489 0 29 0 67804 0 387 327071 0 131 0 3 0 59661 0 388 309038 0 923 0 6 0 58620 0 389 314210 0 72 0 3 0 60398 0 390 307930 0 572 0 8 0 58580 0 391 322327 0 397 0 10 0 62710 0 392 292136 0 450 0 6 0 59325 0 393 263276 0 622 0 8 0 60950 0 394 367655 0 694 0 6 0 68060 0 395 283910 1 3425 3425 9 9 83620 83620 396 283587 0 562 0 8 0 58456 0 397 243650 0 4917 0 26 0 52811 0 398 438493 1 1442 1442 239 239 121173 121173 399 296261 0 529 0 7 0 63870 0 400 230621 1 2126 2126 41 41 21001 21001 401 304252 0 1061 0 3 0 70415 0 402 333505 0 776 0 8 0 64230 0 403 296919 0 611 0 6 0 59190 0 404 278990 1 1526 1526 21 21 69351 69351 405 276898 0 592 0 7 0 64270 0 406 327007 0 1182 0 11 0 70694 0 407 317046 0 621 0 11 0 68005 0 408 304555 0 989 0 12 0 58930 0 409 298096 0 438 0 9 0 58320 0 410 231861 0 726 0 3 0 69980 0 411 309422 0 1303 0 57 0 69863 0 412 286963 1 7419 7419 21 21 63255 63255 413 269753 1 1164 1164 15 15 57320 57320 414 448243 1 3310 3310 32 32 75230 75230 415 165404 0 1920 0 11 0 79420 0 416 204325 0 965 0 2 0 73490 0 417 407159 0 3256 0 23 0 35250 0 418 290476 1 1135 1135 20 20 62285 62285 419 275311 0 1270 0 24 0 69206 0 420 246541 0 661 0 1 0 65920 0 421 253468 0 1013 0 1 0 69770 0 422 240897 0 2844 0 74 0 72683 0 423 -83265 1 11528 11528 68 68 -14545 -14545 424 -42143 0 6526 0 20 0 55830 0 425 272713 0 2264 0 20 0 55174 0 426 215362 1 5109 5109 82 82 67038 67038 427 42754 0 3999 0 21 0 51252 0 428 306275 0 35624 0 244 0 157278 0 429 253537 0 9252 0 32 0 79510 0 430 372631 0 15236 0 86 0 77440 0 431 -7170 0 18073 0 69 0 27284 0 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Group Costs Costs_g Orders Orders_g 195390.250 13727.221 5.051 18.168 -186.494 2735.481 Dividends Dividends_g 1.899 -1.742 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -809715 -16639 4696 20133 2275663 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 1.954e+05 2.778e+04 7.033 8.19e-12 *** Group 1.373e+04 4.962e+04 0.277 0.78218 Costs 5.051e+00 2.855e+00 1.769 0.07756 . Costs_g 1.817e+01 3.481e+00 5.219 2.83e-07 *** Orders -1.865e+02 4.150e+02 -0.449 0.65339 Orders_g 2.735e+03 5.325e+02 5.137 4.26e-07 *** Dividends 1.899e+00 3.894e-01 4.876 1.54e-06 *** Dividends_g -1.742e+00 6.362e-01 -2.738 0.00644 ** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 186000 on 423 degrees of freedom Multiple R-squared: 0.8165, Adjusted R-squared: 0.8134 F-statistic: 268.8 on 7 and 423 DF, p-value: < 2.2e-16 > if (n > n25) { + kp3 <- k + 3 + nmkm3 <- n - k - 3 + gqarr <- array(NA, dim=c(nmkm3-kp3+1,3)) + numgqtests <- 0 + numsignificant1 <- 0 + numsignificant5 <- 0 + numsignificant10 <- 0 + for (mypoint in kp3:nmkm3) { + j <- 0 + numgqtests <- numgqtests + 1 + for (myalt in c('greater', 'two.sided', 'less')) { + j <- j + 1 + gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value + } + if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1 + } + gqarr + } [,1] [,2] [,3] [1,] 1.0000000 1.859270e-53 9.296352e-54 [2,] 1.0000000 1.002634e-57 5.013170e-58 [3,] 1.0000000 1.106020e-88 5.530099e-89 [4,] 1.0000000 8.162639e-88 4.081319e-88 [5,] 1.0000000 1.057709e-96 5.288546e-97 [6,] 1.0000000 5.741315e-105 2.870658e-105 [7,] 1.0000000 1.546406e-111 7.732028e-112 [8,] 1.0000000 5.957946e-117 2.978973e-117 [9,] 1.0000000 3.571888e-118 1.785944e-118 [10,] 1.0000000 4.811510e-119 2.405755e-119 [11,] 1.0000000 4.323433e-125 2.161716e-125 [12,] 1.0000000 1.338208e-125 6.691038e-126 [13,] 1.0000000 3.738779e-127 1.869390e-127 [14,] 1.0000000 1.765264e-126 8.826319e-127 [15,] 1.0000000 1.119176e-127 5.595880e-128 [16,] 1.0000000 3.845026e-129 1.922513e-129 [17,] 1.0000000 1.192014e-133 5.960070e-134 [18,] 1.0000000 1.646223e-133 8.231116e-134 [19,] 1.0000000 4.150782e-143 2.075391e-143 [20,] 1.0000000 3.023552e-143 1.511776e-143 [21,] 1.0000000 7.806372e-150 3.903186e-150 [22,] 1.0000000 2.756734e-150 1.378367e-150 [23,] 1.0000000 1.602737e-159 8.013686e-160 [24,] 1.0000000 7.728812e-162 3.864406e-162 [25,] 1.0000000 6.560252e-162 3.280126e-162 [26,] 1.0000000 3.388429e-161 1.694214e-161 [27,] 1.0000000 1.880456e-166 9.402279e-167 [28,] 1.0000000 8.574646e-172 4.287323e-172 [29,] 1.0000000 3.998124e-171 1.999062e-171 [30,] 1.0000000 4.767691e-172 2.383845e-172 [31,] 1.0000000 4.083583e-172 2.041791e-172 [32,] 1.0000000 5.363968e-172 2.681984e-172 [33,] 1.0000000 2.233137e-172 1.116568e-172 [34,] 1.0000000 1.304259e-171 6.521295e-172 [35,] 1.0000000 3.583430e-171 1.791715e-171 [36,] 1.0000000 2.672162e-170 1.336081e-170 [37,] 1.0000000 1.521574e-172 7.607869e-173 [38,] 1.0000000 4.469866e-173 2.234933e-173 [39,] 1.0000000 5.129913e-173 2.564956e-173 [40,] 1.0000000 7.752304e-174 3.876152e-174 [41,] 1.0000000 3.050100e-174 1.525050e-174 [42,] 1.0000000 3.621719e-174 1.810860e-174 [43,] 1.0000000 1.515167e-173 7.575833e-174 [44,] 1.0000000 8.020644e-173 4.010322e-173 [45,] 1.0000000 1.118021e-178 5.590103e-179 [46,] 1.0000000 4.265351e-178 2.132676e-178 [47,] 1.0000000 3.039374e-190 1.519687e-190 [48,] 1.0000000 6.006216e-191 3.003108e-191 [49,] 1.0000000 4.021584e-190 2.010792e-190 [50,] 1.0000000 7.905830e-190 3.952915e-190 [51,] 1.0000000 3.546664e-191 1.773332e-191 [52,] 1.0000000 1.171620e-193 5.858101e-194 [53,] 1.0000000 3.564972e-193 1.782486e-193 [54,] 1.0000000 6.993212e-193 3.496606e-193 [55,] 1.0000000 5.066030e-192 2.533015e-192 [56,] 1.0000000 4.541343e-191 2.270672e-191 [57,] 1.0000000 8.524810e-191 4.262405e-191 [58,] 1.0000000 6.712301e-190 3.356151e-190 [59,] 1.0000000 3.767170e-189 1.883585e-189 [60,] 1.0000000 1.544637e-188 7.723184e-189 [61,] 1.0000000 1.513599e-187 7.567996e-188 [62,] 1.0000000 5.227659e-188 2.613829e-188 [63,] 1.0000000 1.500655e-187 7.503275e-188 [64,] 1.0000000 1.260969e-186 6.304844e-187 [65,] 1.0000000 9.862948e-186 4.931474e-186 [66,] 1.0000000 1.174082e-185 5.870411e-186 [67,] 1.0000000 9.393846e-185 4.696923e-185 [68,] 1.0000000 9.110194e-185 4.555097e-185 [69,] 1.0000000 6.653659e-184 3.326830e-184 [70,] 1.0000000 1.162629e-184 5.813146e-185 [71,] 1.0000000 8.523358e-184 4.261679e-184 [72,] 1.0000000 8.743694e-184 4.371847e-184 [73,] 1.0000000 2.143346e-183 1.071673e-183 [74,] 1.0000000 1.858758e-182 9.293790e-183 [75,] 1.0000000 1.588619e-181 7.943096e-182 [76,] 1.0000000 4.934589e-182 2.467294e-182 [77,] 1.0000000 6.133660e-182 3.066830e-182 [78,] 1.0000000 4.970145e-181 2.485072e-181 [79,] 1.0000000 1.150320e-180 5.751601e-181 [80,] 1.0000000 9.290865e-180 4.645433e-180 [81,] 1.0000000 5.684153e-179 2.842077e-179 [82,] 1.0000000 2.290821e-178 1.145411e-178 [83,] 1.0000000 2.400047e-179 1.200023e-179 [84,] 1.0000000 1.377684e-178 6.888419e-179 [85,] 1.0000000 2.030862e-178 1.015431e-178 [86,] 1.0000000 4.713138e-179 2.356569e-179 [87,] 1.0000000 2.226320e-178 1.113160e-178 [88,] 1.0000000 1.384287e-177 6.921436e-178 [89,] 1.0000000 9.860386e-177 4.930193e-177 [90,] 1.0000000 5.385393e-176 2.692696e-176 [91,] 1.0000000 3.414243e-175 1.707121e-175 [92,] 1.0000000 5.073707e-176 2.536854e-176 [93,] 1.0000000 7.426381e-176 3.713191e-176 [94,] 1.0000000 7.198109e-176 3.599054e-176 [95,] 1.0000000 3.065713e-175 1.532857e-175 [96,] 1.0000000 7.386030e-175 3.693015e-175 [97,] 1.0000000 1.174400e-176 5.872002e-177 [98,] 1.0000000 1.926782e-178 9.633909e-179 [99,] 1.0000000 1.175430e-177 5.877149e-178 [100,] 1.0000000 2.289161e-177 1.144580e-177 [101,] 1.0000000 1.741724e-177 8.708619e-178 [102,] 1.0000000 1.886852e-177 9.434261e-178 [103,] 1.0000000 1.023256e-176 5.116282e-177 [104,] 1.0000000 6.174353e-176 3.087177e-176 [105,] 1.0000000 1.907022e-176 9.535111e-177 [106,] 1.0000000 8.889662e-176 4.444831e-176 [107,] 1.0000000 6.432145e-175 3.216072e-175 [108,] 1.0000000 2.463997e-175 1.231998e-175 [109,] 1.0000000 2.196189e-176 1.098094e-176 [110,] 1.0000000 1.951773e-175 9.758867e-176 [111,] 1.0000000 7.608143e-175 3.804071e-175 [112,] 1.0000000 6.038643e-174 3.019321e-174 [113,] 1.0000000 3.227808e-173 1.613904e-173 [114,] 1.0000000 3.111016e-173 1.555508e-173 [115,] 1.0000000 4.216492e-173 2.108246e-173 [116,] 1.0000000 1.741588e-173 8.707942e-174 [117,] 1.0000000 1.522130e-172 7.610649e-173 [118,] 1.0000000 1.342996e-171 6.714982e-172 [119,] 1.0000000 1.158537e-170 5.792685e-171 [120,] 1.0000000 9.955929e-170 4.977965e-170 [121,] 1.0000000 8.620859e-169 4.310429e-169 [122,] 1.0000000 1.013482e-168 5.067411e-169 [123,] 1.0000000 8.614615e-168 4.307308e-168 [124,] 1.0000000 7.369994e-167 3.684997e-167 [125,] 1.0000000 6.130761e-166 3.065381e-166 [126,] 1.0000000 5.085850e-165 2.542925e-165 [127,] 1.0000000 4.237830e-164 2.118915e-164 [128,] 1.0000000 3.518238e-163 1.759119e-163 [129,] 1.0000000 2.750148e-162 1.375074e-162 [130,] 1.0000000 2.249791e-161 1.124896e-161 [131,] 1.0000000 1.608721e-160 8.043607e-161 [132,] 1.0000000 1.301862e-159 6.509309e-160 [133,] 1.0000000 1.057447e-158 5.287235e-159 [134,] 1.0000000 8.516373e-158 4.258186e-158 [135,] 1.0000000 6.398076e-157 3.199038e-157 [136,] 1.0000000 4.941617e-156 2.470809e-156 [137,] 1.0000000 3.893192e-155 1.946596e-155 [138,] 1.0000000 3.092629e-154 1.546315e-154 [139,] 1.0000000 2.174795e-153 1.087398e-153 [140,] 1.0000000 1.685108e-152 8.425541e-153 [141,] 1.0000000 1.260904e-151 6.304518e-152 [142,] 1.0000000 9.802149e-151 4.901075e-151 [143,] 1.0000000 7.557074e-150 3.778537e-150 [144,] 1.0000000 5.760767e-149 2.880383e-149 [145,] 1.0000000 4.377791e-148 2.188896e-148 [146,] 1.0000000 3.291658e-147 1.645829e-147 [147,] 1.0000000 2.461976e-146 1.230988e-146 [148,] 1.0000000 1.842507e-145 9.212537e-146 [149,] 1.0000000 1.029251e-144 5.146256e-145 [150,] 1.0000000 7.581243e-144 3.790622e-144 [151,] 1.0000000 5.555053e-143 2.777526e-143 [152,] 1.0000000 4.049198e-142 2.024599e-142 [153,] 1.0000000 2.936212e-141 1.468106e-141 [154,] 1.0000000 2.153285e-140 1.076642e-140 [155,] 1.0000000 1.444617e-139 7.223084e-140 [156,] 1.0000000 1.034609e-138 5.173043e-139 [157,] 1.0000000 7.355165e-138 3.677582e-138 [158,] 1.0000000 1.576302e-138 7.881511e-139 [159,] 1.0000000 3.480937e-138 1.740469e-138 [160,] 1.0000000 2.480282e-137 1.240141e-137 [161,] 1.0000000 1.758249e-136 8.791243e-137 [162,] 1.0000000 1.262759e-135 6.313797e-136 [163,] 1.0000000 8.861637e-135 4.430818e-135 [164,] 1.0000000 6.311552e-134 3.155776e-134 [165,] 1.0000000 4.385392e-133 2.192696e-133 [166,] 1.0000000 3.058349e-132 1.529174e-132 [167,] 1.0000000 1.363451e-131 6.817257e-132 [168,] 1.0000000 6.326054e-131 3.163027e-131 [169,] 1.0000000 4.306306e-130 2.153153e-130 [170,] 1.0000000 2.916511e-129 1.458256e-129 [171,] 1.0000000 1.434010e-128 7.170048e-129 [172,] 1.0000000 7.659527e-128 3.829763e-128 [173,] 1.0000000 3.898096e-127 1.949048e-127 [174,] 1.0000000 1.373586e-126 6.867931e-127 [175,] 1.0000000 7.052713e-126 3.526356e-126 [176,] 1.0000000 3.125534e-125 1.562767e-125 [177,] 1.0000000 7.533314e-125 3.766657e-125 [178,] 1.0000000 5.089730e-124 2.544865e-124 [179,] 1.0000000 2.823128e-123 1.411564e-123 [180,] 1.0000000 1.848844e-122 9.244221e-123 [181,] 1.0000000 1.161433e-121 5.807166e-122 [182,] 1.0000000 7.530845e-121 3.765423e-121 [183,] 1.0000000 4.945832e-120 2.472916e-120 [184,] 1.0000000 3.175377e-119 1.587689e-119 [185,] 1.0000000 2.042847e-118 1.021424e-118 [186,] 1.0000000 1.013895e-117 5.069476e-118 [187,] 1.0000000 3.938831e-117 1.969415e-117 [188,] 1.0000000 1.977995e-116 9.889977e-117 [189,] 1.0000000 1.007420e-115 5.037100e-116 [190,] 1.0000000 6.375272e-116 3.187636e-116 [191,] 1.0000000 3.314475e-115 1.657238e-115 [192,] 1.0000000 1.741010e-114 8.705049e-115 [193,] 1.0000000 9.227527e-114 4.613763e-114 [194,] 1.0000000 5.265745e-113 2.632873e-113 [195,] 1.0000000 1.187410e-112 5.937050e-113 [196,] 1.0000000 6.599407e-112 3.299704e-112 [197,] 1.0000000 3.442296e-111 1.721148e-111 [198,] 1.0000000 1.867705e-110 9.338525e-111 [199,] 1.0000000 1.851183e-110 9.255914e-111 [200,] 1.0000000 1.052821e-109 5.264107e-110 [201,] 1.0000000 6.465497e-109 3.232749e-109 [202,] 1.0000000 3.950933e-108 1.975467e-108 [203,] 1.0000000 1.971419e-107 9.857097e-108 [204,] 1.0000000 1.091789e-106 5.458946e-107 [205,] 1.0000000 6.739202e-106 3.369601e-106 [206,] 1.0000000 4.070813e-105 2.035406e-105 [207,] 1.0000000 2.467736e-104 1.233868e-104 [208,] 1.0000000 1.437769e-103 7.188847e-104 [209,] 1.0000000 8.559181e-103 4.279590e-103 [210,] 1.0000000 3.890513e-102 1.945257e-102 [211,] 1.0000000 2.307902e-101 1.153951e-101 [212,] 1.0000000 1.349082e-100 6.745409e-101 [213,] 1.0000000 7.866339e-100 3.933170e-100 [214,] 1.0000000 4.680174e-99 2.340087e-99 [215,] 1.0000000 2.699914e-98 1.349957e-98 [216,] 1.0000000 1.593521e-97 7.967605e-98 [217,] 1.0000000 9.277306e-97 4.638653e-97 [218,] 1.0000000 5.267942e-96 2.633971e-96 [219,] 1.0000000 3.056520e-95 1.528260e-95 [220,] 1.0000000 1.719012e-94 8.595062e-95 [221,] 1.0000000 9.423447e-94 4.711724e-94 [222,] 1.0000000 4.864129e-93 2.432065e-93 [223,] 1.0000000 2.698320e-92 1.349160e-92 [224,] 1.0000000 1.515895e-91 7.579474e-92 [225,] 1.0000000 8.330060e-91 4.165030e-91 [226,] 1.0000000 4.552405e-90 2.276202e-90 [227,] 1.0000000 2.476545e-89 1.238272e-89 [228,] 1.0000000 1.350140e-88 6.750698e-89 [229,] 1.0000000 7.271984e-88 3.635992e-88 [230,] 1.0000000 3.896967e-87 1.948484e-87 [231,] 1.0000000 2.077760e-86 1.038880e-86 [232,] 1.0000000 1.123832e-85 5.619158e-86 [233,] 1.0000000 5.932133e-85 2.966066e-85 [234,] 1.0000000 3.138790e-84 1.569395e-84 [235,] 1.0000000 1.498153e-83 7.490766e-84 [236,] 1.0000000 7.200266e-83 3.600133e-83 [237,] 1.0000000 3.465835e-82 1.732918e-82 [238,] 1.0000000 1.840094e-81 9.200468e-82 [239,] 1.0000000 8.828882e-81 4.414441e-81 [240,] 1.0000000 4.238824e-80 2.119412e-80 [241,] 1.0000000 2.215489e-79 1.107744e-79 [242,] 1.0000000 1.059322e-78 5.296609e-79 [243,] 1.0000000 5.498331e-78 2.749165e-78 [244,] 1.0000000 1.987154e-77 9.935772e-78 [245,] 1.0000000 5.802622e-78 2.901311e-78 [246,] 1.0000000 2.928944e-77 1.464472e-77 [247,] 1.0000000 1.427303e-76 7.136516e-77 [248,] 1.0000000 7.388195e-76 3.694097e-76 [249,] 1.0000000 3.511975e-75 1.755987e-75 [250,] 1.0000000 1.738928e-74 8.694640e-75 [251,] 1.0000000 8.225467e-74 4.112734e-74 [252,] 1.0000000 3.887067e-73 1.943533e-73 [253,] 1.0000000 1.936282e-72 9.681411e-73 [254,] 1.0000000 8.101872e-72 4.050936e-72 [255,] 1.0000000 3.794102e-71 1.897051e-71 [256,] 1.0000000 1.826726e-70 9.133631e-71 [257,] 1.0000000 8.503974e-70 4.251987e-70 [258,] 1.0000000 2.881226e-69 1.440613e-69 [259,] 1.0000000 1.336075e-68 6.680375e-69 [260,] 1.0000000 6.434815e-68 3.217408e-68 [261,] 1.0000000 2.964194e-67 1.482097e-67 [262,] 1.0000000 1.362323e-66 6.811614e-67 [263,] 1.0000000 6.245535e-66 3.122767e-66 [264,] 1.0000000 2.973588e-65 1.486794e-65 [265,] 1.0000000 1.372955e-64 6.864776e-65 [266,] 1.0000000 6.483862e-64 3.241931e-64 [267,] 1.0000000 3.065062e-63 1.532531e-63 [268,] 1.0000000 1.393436e-62 6.967182e-63 [269,] 1.0000000 5.613624e-62 2.806812e-62 [270,] 1.0000000 2.563572e-61 1.281786e-61 [271,] 1.0000000 1.133035e-60 5.665176e-61 [272,] 1.0000000 5.057098e-60 2.528549e-60 [273,] 1.0000000 2.331364e-59 1.165682e-59 [274,] 1.0000000 1.029533e-58 5.147665e-59 [275,] 1.0000000 4.661273e-58 2.330636e-58 [276,] 1.0000000 2.035881e-57 1.017940e-57 [277,] 1.0000000 8.859554e-57 4.429777e-57 [278,] 1.0000000 3.874733e-56 1.937366e-56 [279,] 1.0000000 1.653626e-55 8.268128e-56 [280,] 1.0000000 7.306698e-55 3.653349e-55 [281,] 1.0000000 3.091235e-54 1.545617e-54 [282,] 1.0000000 1.332340e-53 6.661700e-54 [283,] 1.0000000 5.603748e-53 2.801874e-53 [284,] 1.0000000 2.347421e-52 1.173710e-52 [285,] 1.0000000 1.011522e-51 5.057611e-52 [286,] 1.0000000 3.470761e-51 1.735380e-51 [287,] 1.0000000 1.433157e-50 7.165787e-51 [288,] 1.0000000 5.892636e-50 2.946318e-50 [289,] 1.0000000 2.406769e-49 1.203385e-49 [290,] 1.0000000 1.010638e-48 5.053188e-49 [291,] 1.0000000 4.131765e-48 2.065883e-48 [292,] 1.0000000 1.709760e-47 8.548800e-48 [293,] 1.0000000 6.862904e-47 3.431452e-47 [294,] 1.0000000 2.739420e-46 1.369710e-46 [295,] 1.0000000 9.243566e-46 4.621783e-46 [296,] 1.0000000 3.773930e-45 1.886965e-45 [297,] 1.0000000 1.485752e-44 7.428762e-45 [298,] 1.0000000 5.953781e-44 2.976891e-44 [299,] 1.0000000 2.355398e-43 1.177699e-43 [300,] 1.0000000 9.123191e-43 4.561596e-43 [301,] 1.0000000 3.515290e-42 1.757645e-42 [302,] 1.0000000 1.202977e-41 6.014883e-42 [303,] 1.0000000 4.589251e-41 2.294626e-41 [304,] 1.0000000 1.769288e-40 8.846439e-41 [305,] 1.0000000 6.816804e-40 3.408402e-40 [306,] 1.0000000 2.343019e-39 1.171510e-39 [307,] 1.0000000 8.608296e-39 4.304148e-39 [308,] 1.0000000 3.114822e-38 1.557411e-38 [309,] 1.0000000 1.163686e-37 5.818430e-38 [310,] 1.0000000 4.230651e-37 2.115326e-37 [311,] 1.0000000 1.527991e-36 7.639955e-37 [312,] 1.0000000 5.482204e-36 2.741102e-36 [313,] 1.0000000 1.953839e-35 9.769197e-36 [314,] 1.0000000 6.994808e-35 3.497404e-35 [315,] 1.0000000 2.461911e-34 1.230956e-34 [316,] 1.0000000 8.743055e-34 4.371528e-34 [317,] 1.0000000 3.067151e-33 1.533576e-33 [318,] 1.0000000 1.063840e-32 5.319201e-33 [319,] 1.0000000 3.650664e-32 1.825332e-32 [320,] 1.0000000 1.244026e-31 6.220131e-32 [321,] 1.0000000 4.209433e-31 2.104717e-31 [322,] 1.0000000 1.414255e-30 7.071274e-31 [323,] 1.0000000 4.717527e-30 2.358763e-30 [324,] 1.0000000 1.590148e-29 7.950739e-30 [325,] 1.0000000 5.255731e-29 2.627866e-29 [326,] 1.0000000 1.716936e-28 8.584679e-29 [327,] 1.0000000 5.720107e-28 2.860054e-28 [328,] 1.0000000 1.842462e-27 9.212309e-28 [329,] 1.0000000 5.958614e-27 2.979307e-27 [330,] 1.0000000 1.933843e-26 9.669214e-27 [331,] 1.0000000 6.223739e-26 3.111869e-26 [332,] 1.0000000 1.949054e-25 9.745268e-26 [333,] 1.0000000 6.056215e-25 3.028108e-25 [334,] 1.0000000 1.867015e-24 9.335076e-25 [335,] 1.0000000 5.958287e-24 2.979144e-24 [336,] 1.0000000 1.808832e-23 9.044158e-24 [337,] 1.0000000 5.651498e-23 2.825749e-23 [338,] 1.0000000 1.688598e-22 8.442991e-23 [339,] 1.0000000 5.003528e-22 2.501764e-22 [340,] 1.0000000 1.345220e-21 6.726102e-22 [341,] 1.0000000 3.750855e-21 1.875428e-21 [342,] 1.0000000 1.124412e-20 5.622058e-21 [343,] 1.0000000 3.069099e-20 1.534550e-20 [344,] 1.0000000 8.841589e-20 4.420794e-20 [345,] 1.0000000 2.491506e-19 1.245753e-19 [346,] 1.0000000 6.957594e-19 3.478797e-19 [347,] 1.0000000 2.010850e-18 1.005425e-18 [348,] 1.0000000 5.528211e-18 2.764105e-18 [349,] 1.0000000 1.505568e-17 7.527842e-18 [350,] 1.0000000 4.061301e-17 2.030651e-17 [351,] 1.0000000 1.108940e-16 5.544701e-17 [352,] 1.0000000 2.934615e-16 1.467308e-16 [353,] 1.0000000 7.688602e-16 3.844301e-16 [354,] 1.0000000 2.114519e-15 1.057260e-15 [355,] 1.0000000 5.581670e-15 2.790835e-15 [356,] 1.0000000 1.422068e-14 7.110342e-15 [357,] 1.0000000 3.584977e-14 1.792488e-14 [358,] 1.0000000 8.940791e-14 4.470396e-14 [359,] 1.0000000 2.205456e-13 1.102728e-13 [360,] 1.0000000 5.527686e-13 2.763843e-13 [361,] 1.0000000 1.334095e-12 6.670476e-13 [362,] 1.0000000 3.072227e-12 1.536114e-12 [363,] 1.0000000 7.622041e-12 3.811021e-12 [364,] 1.0000000 1.890643e-11 9.453213e-12 [365,] 1.0000000 4.358722e-11 2.179361e-11 [366,] 1.0000000 9.921869e-11 4.960935e-11 [367,] 1.0000000 2.427350e-10 1.213675e-10 [368,] 1.0000000 5.225316e-10 2.612658e-10 [369,] 1.0000000 9.975512e-10 4.987756e-10 [370,] 1.0000000 2.160301e-09 1.080151e-09 [371,] 1.0000000 4.611232e-09 2.305616e-09 [372,] 1.0000000 1.084087e-08 5.420436e-09 [373,] 1.0000000 2.416476e-08 1.208238e-08 [374,] 1.0000000 4.973695e-08 2.486848e-08 [375,] 0.9999999 1.086228e-07 5.431139e-08 [376,] 0.9999999 2.072455e-07 1.036227e-07 [377,] 0.9999998 3.905439e-07 1.952720e-07 [378,] 0.9999996 7.636398e-07 3.818199e-07 [379,] 0.9999992 1.508221e-06 7.541106e-07 [380,] 0.9999985 2.981606e-06 1.490803e-06 [381,] 0.9999971 5.748432e-06 2.874216e-06 [382,] 0.9999942 1.162291e-05 5.811454e-06 [383,] 0.9999880 2.408334e-05 1.204167e-05 [384,] 0.9999838 3.237006e-05 1.618503e-05 [385,] 0.9999746 5.085836e-05 2.542918e-05 [386,] 0.9999498 1.004109e-04 5.020546e-05 [387,] 0.9999033 1.934453e-04 9.672267e-05 [388,] 0.9998421 3.158941e-04 1.579471e-04 [389,] 0.9997056 5.888583e-04 2.944291e-04 [390,] 0.9996335 7.330176e-04 3.665088e-04 [391,] 0.9993539 1.292101e-03 6.460506e-04 [392,] 0.9990169 1.966253e-03 9.831267e-04 [393,] 0.9983338 3.332432e-03 1.666216e-03 [394,] 0.9971950 5.610010e-03 2.805005e-03 [395,] 0.9950441 9.911888e-03 4.955944e-03 [396,] 0.9925095 1.498099e-02 7.490493e-03 [397,] 0.9883788 2.324236e-02 1.162118e-02 [398,] 0.9825882 3.482356e-02 1.741178e-02 [399,] 0.9739988 5.200233e-02 2.600117e-02 [400,] 0.9571630 8.567394e-02 4.283697e-02 [401,] 0.9313077 1.373846e-01 6.869229e-02 [402,] 0.9152946 1.694107e-01 8.470536e-02 [403,] 0.8692038 2.615924e-01 1.307962e-01 [404,] 0.8056023 3.887953e-01 1.943977e-01 [405,] 0.7418097 5.163806e-01 2.581903e-01 [406,] 0.6467826 7.064347e-01 3.532174e-01 [407,] 0.7987507 4.024987e-01 2.012493e-01 [408,] 0.6871032 6.257936e-01 3.128968e-01 [409,] 0.5494810 9.010379e-01 4.505190e-01 [410,] 0.3868022 7.736045e-01 6.131978e-01 > postscript(file="/var/www/html/rcomp/tmp/1ph0j1291225443.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index') > points(x[,1]-mysum$resid) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/2ph0j1291225443.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/3ph0j1291225443.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/40qhm1291225443.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals') > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/50qhm1291225443.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > qqnorm(mysum$resid, main='Residual Normal Q-Q Plot') > qqline(mysum$resid) > grid() > dev.off() null device 1 > (myerror <- as.ts(mysum$resid)) Time Series: Start = 1 End = 431 Frequency = 1 1 2 3 4 5 208308.99237 2275662.88080 795828.19898 -270987.19826 -809715.40781 6 7 8 9 10 -22266.82622 800116.23761 -49155.83317 -691527.85593 520274.61695 11 12 13 14 15 300395.15083 657029.06624 -223632.43021 -214806.51231 -584732.95539 16 17 18 19 20 567057.94638 -120685.92990 -555208.39746 -692034.64645 -22517.65388 21 22 23 24 25 147556.65102 -198905.77113 144552.96763 -85472.26929 -172515.67738 26 27 28 29 30 199630.89796 -806793.29430 -91537.29089 493285.69665 102851.24420 31 32 33 34 35 321503.16649 -104258.33967 255545.92790 160686.08391 -132980.09112 36 37 38 39 40 -66128.03087 238335.88236 142898.18103 44953.18110 -135648.29772 41 42 43 44 45 -143145.37863 -91583.43595 121975.57705 -56834.07852 44665.36706 46 47 48 49 50 2854.95073 -109747.17134 -353406.30019 28280.34861 -149859.43355 51 52 53 54 55 1830.56568 -154893.59582 -163816.82868 108069.30970 -398627.52899 56 57 58 59 60 96414.42623 344119.88167 65897.34530 -39937.02809 -268342.71677 61 62 63 64 65 30684.17243 -165191.81686 54293.11167 -229559.34975 44287.34309 66 67 68 69 70 -30222.19008 -89225.87105 12463.09493 -62985.93646 -110309.30452 71 72 73 74 75 -206946.64271 90240.21711 -88273.72819 -21316.72675 -30243.60472 76 77 78 79 80 -66184.57835 -44934.51284 -101207.15920 -11678.29645 63643.05378 81 82 83 84 85 34641.61231 -73790.49159 -104685.81508 29937.67397 -12725.14358 86 87 88 89 90 -167364.74107 91819.60560 -38862.44579 -91394.42511 -5660.57195 91 92 93 94 95 13129.61159 -14033.50196 137631.86438 -89777.21321 90481.59227 96 97 98 99 100 10678.39226 19473.69877 21519.12241 -48344.28494 42245.68748 101 102 103 104 105 -21816.82929 -133542.53336 -108868.71359 -98741.69172 39063.39389 106 107 108 109 110 -59967.40542 -289570.08055 198107.49533 -86218.57294 -66492.21629 111 112 113 114 115 83864.87293 -222070.80804 -84451.81816 -19997.95631 112369.79821 116 117 118 119 120 39988.90315 -24563.04824 71140.54934 -86911.82797 3663.56051 121 122 123 124 125 -65922.30668 9489.98057 -3418.89374 152158.54615 29286.34957 126 127 128 129 130 -138918.74718 4695.80258 5100.36378 14817.85613 4695.80258 131 132 133 134 135 4973.42744 65208.58145 4695.80258 5357.41047 -1120.95879 136 137 138 139 140 3882.97829 9075.68325 4364.04313 444.25026 4695.80258 141 142 143 144 145 17633.90268 4695.80258 -7587.46063 5383.56487 14587.67173 146 147 148 149 150 10257.96358 4695.80258 5344.44639 13569.03592 4560.05297 151 152 153 154 155 -16424.65558 -4530.63985 5628.27428 -9949.03909 4418.84197 156 157 158 159 160 4695.80258 4695.80258 2967.75607 24390.36985 4695.80258 161 162 163 164 165 4695.80258 4695.80258 4695.80258 -5759.11887 -11139.17833 166 167 168 169 170 4146.91131 4695.80258 100101.84568 17631.71406 4695.80258 171 172 173 174 175 4695.80258 502.88764 4695.80258 -2747.76674 4695.80258 176 177 178 179 180 5279.74355 21636.35992 95888.64101 4695.80258 4695.80258 181 182 183 184 185 26546.37872 69353.44777 5194.43625 27469.56632 82604.48615 186 187 188 189 190 71777.73429 20793.09613 505.91678 -28647.56861 4695.80258 191 192 193 194 195 10293.81670 4695.80258 3547.21593 4695.80258 5099.07171 196 197 198 199 200 96745.06614 -159756.58648 96745.06614 96745.06614 -11466.83054 201 202 203 204 205 96745.06614 96745.06614 96745.06614 48580.16067 60559.55258 206 207 208 209 210 84396.23774 92599.63230 96745.06614 34240.26905 74513.23008 211 212 213 214 215 4695.80258 4695.80258 25925.93172 -24615.08018 -11571.51316 216 217 218 219 220 5068.79126 5836.24263 6492.62189 5068.79126 15978.94819 221 222 223 224 225 -14524.58874 4695.80258 4882.91782 2842.93476 4815.88543 226 227 228 229 230 5637.80160 5237.15542 4695.80258 -7592.79784 4695.80258 231 232 233 234 235 -18346.43070 8642.79244 4695.80258 -8722.67173 4695.80258 236 237 238 239 240 4749.02911 4695.80258 6664.52600 4695.80258 4695.80258 241 242 243 244 245 4695.80258 1072.29490 4695.80258 7849.74472 -20413.98966 246 247 248 249 250 96745.06614 96745.06614 -1625.30469 96745.06614 96745.06614 251 252 253 254 255 1105.51459 96745.06614 -10712.38062 -25693.29890 -100020.13876 256 257 258 259 260 4695.80258 9540.54090 15504.80685 96745.06614 4695.80258 261 262 263 264 265 96745.06614 96745.06614 11524.03289 25124.49007 96745.06614 266 267 268 269 270 4695.80258 96745.06614 39068.86267 96745.06614 7614.44532 271 272 273 274 275 96745.06614 96745.06614 96745.06614 2280.74545 4695.80258 276 277 278 279 280 6856.73328 5750.88219 4695.80258 -32509.02284 -13492.53170 281 282 283 284 285 96745.06614 4695.80258 13981.74577 4695.80258 1474.75604 286 287 288 289 290 4695.80258 73875.87396 5305.87880 96745.06614 -4616.91129 291 292 293 294 295 96745.06614 11869.04635 17785.58720 4695.80258 -4372.82474 296 297 298 299 300 90735.14542 4695.80258 96745.06614 4695.80258 6374.25164 301 302 303 304 305 83624.34209 -7133.77660 4882.29692 96745.06614 32437.27291 306 307 308 309 310 6909.15144 96745.06614 44212.05015 11845.65282 96745.06614 311 312 313 314 315 96745.06614 -10213.42054 96745.06614 66041.66334 6370.75503 316 317 318 319 320 18828.45632 5508.08560 -12520.74431 5617.12230 4695.80258 321 322 323 324 325 4695.80258 4695.80258 4695.80258 96745.06614 4695.80258 326 327 328 329 330 1889.38269 96745.06614 96745.06614 4695.80258 4695.80258 331 332 333 334 335 4695.80258 4695.80258 4695.80258 3953.39132 4871.95470 336 337 338 339 340 3846.84201 -15467.40013 4695.80258 2916.79186 10119.59876 341 342 343 344 345 -18672.05287 4695.80258 4695.80258 4695.80258 -3140.29155 346 347 348 349 350 4695.80258 -4006.90833 4695.80258 4695.80258 2493.78740 351 352 353 354 355 11018.66125 3600.11020 16051.07917 476.58495 4695.80258 356 357 358 359 360 4695.80258 -22396.84628 4695.80258 4695.80258 4695.80258 361 362 363 364 365 5450.93162 4695.80258 4695.80258 -22742.53307 8531.29201 366 367 368 369 370 4695.80258 4695.80258 4695.80258 4695.80258 5073.19778 371 372 373 374 375 4695.80258 13502.47832 -4122.43855 -11364.54599 4695.80258 376 377 378 379 380 4695.80258 -5249.64067 10434.70666 -127.53011 4695.80258 381 382 383 384 385 4695.80258 -16852.39645 5978.97241 4695.80258 4730.19031 386 387 388 389 390 -93173.15808 18295.44481 -1201.64187 4333.03834 -87.88601 391 392 393 394 395 7724.02211 -17053.27366 -49494.53350 40647.49443 -40778.26494 396 397 398 399 400 -24144.92966 -72002.46594 -432306.25565 -21770.74622 -135659.00978 401 402 403 404 405 -29640.24458 13728.65001 -12827.11030 -29957.78483 -42211.45613 406 407 408 409 410 -6534.19128 -8555.89563 -5487.64777 -8564.90722 -99513.27125 411 412 413 414 415 -14573.70679 -157855.32051 -13610.02368 68913.23945 -188433.40848 416 417 418 419 420 -135107.61599 132680.84269 -5736.85380 -53424.84627 -77168.91215 421 422 423 424 425 -79330.08351 -93066.08147 -731095.45915 -372773.51927 -35145.58987 426 427 428 429 430 -331902.94857 -266234.14336 -322175.21212 -133587.07324 -30715.63764 431 -332780.81006 > postscript(file="/var/www/html/rcomp/tmp/6shy71291225443.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > dum <- cbind(lag(myerror,k=1),myerror) > dum Time Series: Start = 0 End = 431 Frequency = 1 lag(myerror, k = 1) myerror 0 208308.99237 NA 1 2275662.88080 208308.99237 2 795828.19898 2275662.88080 3 -270987.19826 795828.19898 4 -809715.40781 -270987.19826 5 -22266.82622 -809715.40781 6 800116.23761 -22266.82622 7 -49155.83317 800116.23761 8 -691527.85593 -49155.83317 9 520274.61695 -691527.85593 10 300395.15083 520274.61695 11 657029.06624 300395.15083 12 -223632.43021 657029.06624 13 -214806.51231 -223632.43021 14 -584732.95539 -214806.51231 15 567057.94638 -584732.95539 16 -120685.92990 567057.94638 17 -555208.39746 -120685.92990 18 -692034.64645 -555208.39746 19 -22517.65388 -692034.64645 20 147556.65102 -22517.65388 21 -198905.77113 147556.65102 22 144552.96763 -198905.77113 23 -85472.26929 144552.96763 24 -172515.67738 -85472.26929 25 199630.89796 -172515.67738 26 -806793.29430 199630.89796 27 -91537.29089 -806793.29430 28 493285.69665 -91537.29089 29 102851.24420 493285.69665 30 321503.16649 102851.24420 31 -104258.33967 321503.16649 32 255545.92790 -104258.33967 33 160686.08391 255545.92790 34 -132980.09112 160686.08391 35 -66128.03087 -132980.09112 36 238335.88236 -66128.03087 37 142898.18103 238335.88236 38 44953.18110 142898.18103 39 -135648.29772 44953.18110 40 -143145.37863 -135648.29772 41 -91583.43595 -143145.37863 42 121975.57705 -91583.43595 43 -56834.07852 121975.57705 44 44665.36706 -56834.07852 45 2854.95073 44665.36706 46 -109747.17134 2854.95073 47 -353406.30019 -109747.17134 48 28280.34861 -353406.30019 49 -149859.43355 28280.34861 50 1830.56568 -149859.43355 51 -154893.59582 1830.56568 52 -163816.82868 -154893.59582 53 108069.30970 -163816.82868 54 -398627.52899 108069.30970 55 96414.42623 -398627.52899 56 344119.88167 96414.42623 57 65897.34530 344119.88167 58 -39937.02809 65897.34530 59 -268342.71677 -39937.02809 60 30684.17243 -268342.71677 61 -165191.81686 30684.17243 62 54293.11167 -165191.81686 63 -229559.34975 54293.11167 64 44287.34309 -229559.34975 65 -30222.19008 44287.34309 66 -89225.87105 -30222.19008 67 12463.09493 -89225.87105 68 -62985.93646 12463.09493 69 -110309.30452 -62985.93646 70 -206946.64271 -110309.30452 71 90240.21711 -206946.64271 72 -88273.72819 90240.21711 73 -21316.72675 -88273.72819 74 -30243.60472 -21316.72675 75 -66184.57835 -30243.60472 76 -44934.51284 -66184.57835 77 -101207.15920 -44934.51284 78 -11678.29645 -101207.15920 79 63643.05378 -11678.29645 80 34641.61231 63643.05378 81 -73790.49159 34641.61231 82 -104685.81508 -73790.49159 83 29937.67397 -104685.81508 84 -12725.14358 29937.67397 85 -167364.74107 -12725.14358 86 91819.60560 -167364.74107 87 -38862.44579 91819.60560 88 -91394.42511 -38862.44579 89 -5660.57195 -91394.42511 90 13129.61159 -5660.57195 91 -14033.50196 13129.61159 92 137631.86438 -14033.50196 93 -89777.21321 137631.86438 94 90481.59227 -89777.21321 95 10678.39226 90481.59227 96 19473.69877 10678.39226 97 21519.12241 19473.69877 98 -48344.28494 21519.12241 99 42245.68748 -48344.28494 100 -21816.82929 42245.68748 101 -133542.53336 -21816.82929 102 -108868.71359 -133542.53336 103 -98741.69172 -108868.71359 104 39063.39389 -98741.69172 105 -59967.40542 39063.39389 106 -289570.08055 -59967.40542 107 198107.49533 -289570.08055 108 -86218.57294 198107.49533 109 -66492.21629 -86218.57294 110 83864.87293 -66492.21629 111 -222070.80804 83864.87293 112 -84451.81816 -222070.80804 113 -19997.95631 -84451.81816 114 112369.79821 -19997.95631 115 39988.90315 112369.79821 116 -24563.04824 39988.90315 117 71140.54934 -24563.04824 118 -86911.82797 71140.54934 119 3663.56051 -86911.82797 120 -65922.30668 3663.56051 121 9489.98057 -65922.30668 122 -3418.89374 9489.98057 123 152158.54615 -3418.89374 124 29286.34957 152158.54615 125 -138918.74718 29286.34957 126 4695.80258 -138918.74718 127 5100.36378 4695.80258 128 14817.85613 5100.36378 129 4695.80258 14817.85613 130 4973.42744 4695.80258 131 65208.58145 4973.42744 132 4695.80258 65208.58145 133 5357.41047 4695.80258 134 -1120.95879 5357.41047 135 3882.97829 -1120.95879 136 9075.68325 3882.97829 137 4364.04313 9075.68325 138 444.25026 4364.04313 139 4695.80258 444.25026 140 17633.90268 4695.80258 141 4695.80258 17633.90268 142 -7587.46063 4695.80258 143 5383.56487 -7587.46063 144 14587.67173 5383.56487 145 10257.96358 14587.67173 146 4695.80258 10257.96358 147 5344.44639 4695.80258 148 13569.03592 5344.44639 149 4560.05297 13569.03592 150 -16424.65558 4560.05297 151 -4530.63985 -16424.65558 152 5628.27428 -4530.63985 153 -9949.03909 5628.27428 154 4418.84197 -9949.03909 155 4695.80258 4418.84197 156 4695.80258 4695.80258 157 2967.75607 4695.80258 158 24390.36985 2967.75607 159 4695.80258 24390.36985 160 4695.80258 4695.80258 161 4695.80258 4695.80258 162 4695.80258 4695.80258 163 -5759.11887 4695.80258 164 -11139.17833 -5759.11887 165 4146.91131 -11139.17833 166 4695.80258 4146.91131 167 100101.84568 4695.80258 168 17631.71406 100101.84568 169 4695.80258 17631.71406 170 4695.80258 4695.80258 171 502.88764 4695.80258 172 4695.80258 502.88764 173 -2747.76674 4695.80258 174 4695.80258 -2747.76674 175 5279.74355 4695.80258 176 21636.35992 5279.74355 177 95888.64101 21636.35992 178 4695.80258 95888.64101 179 4695.80258 4695.80258 180 26546.37872 4695.80258 181 69353.44777 26546.37872 182 5194.43625 69353.44777 183 27469.56632 5194.43625 184 82604.48615 27469.56632 185 71777.73429 82604.48615 186 20793.09613 71777.73429 187 505.91678 20793.09613 188 -28647.56861 505.91678 189 4695.80258 -28647.56861 190 10293.81670 4695.80258 191 4695.80258 10293.81670 192 3547.21593 4695.80258 193 4695.80258 3547.21593 194 5099.07171 4695.80258 195 96745.06614 5099.07171 196 -159756.58648 96745.06614 197 96745.06614 -159756.58648 198 96745.06614 96745.06614 199 -11466.83054 96745.06614 200 96745.06614 -11466.83054 201 96745.06614 96745.06614 202 96745.06614 96745.06614 203 48580.16067 96745.06614 204 60559.55258 48580.16067 205 84396.23774 60559.55258 206 92599.63230 84396.23774 207 96745.06614 92599.63230 208 34240.26905 96745.06614 209 74513.23008 34240.26905 210 4695.80258 74513.23008 211 4695.80258 4695.80258 212 25925.93172 4695.80258 213 -24615.08018 25925.93172 214 -11571.51316 -24615.08018 215 5068.79126 -11571.51316 216 5836.24263 5068.79126 217 6492.62189 5836.24263 218 5068.79126 6492.62189 219 15978.94819 5068.79126 220 -14524.58874 15978.94819 221 4695.80258 -14524.58874 222 4882.91782 4695.80258 223 2842.93476 4882.91782 224 4815.88543 2842.93476 225 5637.80160 4815.88543 226 5237.15542 5637.80160 227 4695.80258 5237.15542 228 -7592.79784 4695.80258 229 4695.80258 -7592.79784 230 -18346.43070 4695.80258 231 8642.79244 -18346.43070 232 4695.80258 8642.79244 233 -8722.67173 4695.80258 234 4695.80258 -8722.67173 235 4749.02911 4695.80258 236 4695.80258 4749.02911 237 6664.52600 4695.80258 238 4695.80258 6664.52600 239 4695.80258 4695.80258 240 4695.80258 4695.80258 241 1072.29490 4695.80258 242 4695.80258 1072.29490 243 7849.74472 4695.80258 244 -20413.98966 7849.74472 245 96745.06614 -20413.98966 246 96745.06614 96745.06614 247 -1625.30469 96745.06614 248 96745.06614 -1625.30469 249 96745.06614 96745.06614 250 1105.51459 96745.06614 251 96745.06614 1105.51459 252 -10712.38062 96745.06614 253 -25693.29890 -10712.38062 254 -100020.13876 -25693.29890 255 4695.80258 -100020.13876 256 9540.54090 4695.80258 257 15504.80685 9540.54090 258 96745.06614 15504.80685 259 4695.80258 96745.06614 260 96745.06614 4695.80258 261 96745.06614 96745.06614 262 11524.03289 96745.06614 263 25124.49007 11524.03289 264 96745.06614 25124.49007 265 4695.80258 96745.06614 266 96745.06614 4695.80258 267 39068.86267 96745.06614 268 96745.06614 39068.86267 269 7614.44532 96745.06614 270 96745.06614 7614.44532 271 96745.06614 96745.06614 272 96745.06614 96745.06614 273 2280.74545 96745.06614 274 4695.80258 2280.74545 275 6856.73328 4695.80258 276 5750.88219 6856.73328 277 4695.80258 5750.88219 278 -32509.02284 4695.80258 279 -13492.53170 -32509.02284 280 96745.06614 -13492.53170 281 4695.80258 96745.06614 282 13981.74577 4695.80258 283 4695.80258 13981.74577 284 1474.75604 4695.80258 285 4695.80258 1474.75604 286 73875.87396 4695.80258 287 5305.87880 73875.87396 288 96745.06614 5305.87880 289 -4616.91129 96745.06614 290 96745.06614 -4616.91129 291 11869.04635 96745.06614 292 17785.58720 11869.04635 293 4695.80258 17785.58720 294 -4372.82474 4695.80258 295 90735.14542 -4372.82474 296 4695.80258 90735.14542 297 96745.06614 4695.80258 298 4695.80258 96745.06614 299 6374.25164 4695.80258 300 83624.34209 6374.25164 301 -7133.77660 83624.34209 302 4882.29692 -7133.77660 303 96745.06614 4882.29692 304 32437.27291 96745.06614 305 6909.15144 32437.27291 306 96745.06614 6909.15144 307 44212.05015 96745.06614 308 11845.65282 44212.05015 309 96745.06614 11845.65282 310 96745.06614 96745.06614 311 -10213.42054 96745.06614 312 96745.06614 -10213.42054 313 66041.66334 96745.06614 314 6370.75503 66041.66334 315 18828.45632 6370.75503 316 5508.08560 18828.45632 317 -12520.74431 5508.08560 318 5617.12230 -12520.74431 319 4695.80258 5617.12230 320 4695.80258 4695.80258 321 4695.80258 4695.80258 322 4695.80258 4695.80258 323 96745.06614 4695.80258 324 4695.80258 96745.06614 325 1889.38269 4695.80258 326 96745.06614 1889.38269 327 96745.06614 96745.06614 328 4695.80258 96745.06614 329 4695.80258 4695.80258 330 4695.80258 4695.80258 331 4695.80258 4695.80258 332 4695.80258 4695.80258 333 3953.39132 4695.80258 334 4871.95470 3953.39132 335 3846.84201 4871.95470 336 -15467.40013 3846.84201 337 4695.80258 -15467.40013 338 2916.79186 4695.80258 339 10119.59876 2916.79186 340 -18672.05287 10119.59876 341 4695.80258 -18672.05287 342 4695.80258 4695.80258 343 4695.80258 4695.80258 344 -3140.29155 4695.80258 345 4695.80258 -3140.29155 346 -4006.90833 4695.80258 347 4695.80258 -4006.90833 348 4695.80258 4695.80258 349 2493.78740 4695.80258 350 11018.66125 2493.78740 351 3600.11020 11018.66125 352 16051.07917 3600.11020 353 476.58495 16051.07917 354 4695.80258 476.58495 355 4695.80258 4695.80258 356 -22396.84628 4695.80258 357 4695.80258 -22396.84628 358 4695.80258 4695.80258 359 4695.80258 4695.80258 360 5450.93162 4695.80258 361 4695.80258 5450.93162 362 4695.80258 4695.80258 363 -22742.53307 4695.80258 364 8531.29201 -22742.53307 365 4695.80258 8531.29201 366 4695.80258 4695.80258 367 4695.80258 4695.80258 368 4695.80258 4695.80258 369 5073.19778 4695.80258 370 4695.80258 5073.19778 371 13502.47832 4695.80258 372 -4122.43855 13502.47832 373 -11364.54599 -4122.43855 374 4695.80258 -11364.54599 375 4695.80258 4695.80258 376 -5249.64067 4695.80258 377 10434.70666 -5249.64067 378 -127.53011 10434.70666 379 4695.80258 -127.53011 380 4695.80258 4695.80258 381 -16852.39645 4695.80258 382 5978.97241 -16852.39645 383 4695.80258 5978.97241 384 4730.19031 4695.80258 385 -93173.15808 4730.19031 386 18295.44481 -93173.15808 387 -1201.64187 18295.44481 388 4333.03834 -1201.64187 389 -87.88601 4333.03834 390 7724.02211 -87.88601 391 -17053.27366 7724.02211 392 -49494.53350 -17053.27366 393 40647.49443 -49494.53350 394 -40778.26494 40647.49443 395 -24144.92966 -40778.26494 396 -72002.46594 -24144.92966 397 -432306.25565 -72002.46594 398 -21770.74622 -432306.25565 399 -135659.00978 -21770.74622 400 -29640.24458 -135659.00978 401 13728.65001 -29640.24458 402 -12827.11030 13728.65001 403 -29957.78483 -12827.11030 404 -42211.45613 -29957.78483 405 -6534.19128 -42211.45613 406 -8555.89563 -6534.19128 407 -5487.64777 -8555.89563 408 -8564.90722 -5487.64777 409 -99513.27125 -8564.90722 410 -14573.70679 -99513.27125 411 -157855.32051 -14573.70679 412 -13610.02368 -157855.32051 413 68913.23945 -13610.02368 414 -188433.40848 68913.23945 415 -135107.61599 -188433.40848 416 132680.84269 -135107.61599 417 -5736.85380 132680.84269 418 -53424.84627 -5736.85380 419 -77168.91215 -53424.84627 420 -79330.08351 -77168.91215 421 -93066.08147 -79330.08351 422 -731095.45915 -93066.08147 423 -372773.51927 -731095.45915 424 -35145.58987 -372773.51927 425 -331902.94857 -35145.58987 426 -266234.14336 -331902.94857 427 -322175.21212 -266234.14336 428 -133587.07324 -322175.21212 429 -30715.63764 -133587.07324 430 -332780.81006 -30715.63764 431 NA -332780.81006 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 2275662.88080 208308.99237 [2,] 795828.19898 2275662.88080 [3,] -270987.19826 795828.19898 [4,] -809715.40781 -270987.19826 [5,] -22266.82622 -809715.40781 [6,] 800116.23761 -22266.82622 [7,] -49155.83317 800116.23761 [8,] -691527.85593 -49155.83317 [9,] 520274.61695 -691527.85593 [10,] 300395.15083 520274.61695 [11,] 657029.06624 300395.15083 [12,] -223632.43021 657029.06624 [13,] -214806.51231 -223632.43021 [14,] -584732.95539 -214806.51231 [15,] 567057.94638 -584732.95539 [16,] -120685.92990 567057.94638 [17,] -555208.39746 -120685.92990 [18,] -692034.64645 -555208.39746 [19,] -22517.65388 -692034.64645 [20,] 147556.65102 -22517.65388 [21,] -198905.77113 147556.65102 [22,] 144552.96763 -198905.77113 [23,] -85472.26929 144552.96763 [24,] -172515.67738 -85472.26929 [25,] 199630.89796 -172515.67738 [26,] -806793.29430 199630.89796 [27,] -91537.29089 -806793.29430 [28,] 493285.69665 -91537.29089 [29,] 102851.24420 493285.69665 [30,] 321503.16649 102851.24420 [31,] -104258.33967 321503.16649 [32,] 255545.92790 -104258.33967 [33,] 160686.08391 255545.92790 [34,] -132980.09112 160686.08391 [35,] -66128.03087 -132980.09112 [36,] 238335.88236 -66128.03087 [37,] 142898.18103 238335.88236 [38,] 44953.18110 142898.18103 [39,] -135648.29772 44953.18110 [40,] -143145.37863 -135648.29772 [41,] -91583.43595 -143145.37863 [42,] 121975.57705 -91583.43595 [43,] -56834.07852 121975.57705 [44,] 44665.36706 -56834.07852 [45,] 2854.95073 44665.36706 [46,] -109747.17134 2854.95073 [47,] -353406.30019 -109747.17134 [48,] 28280.34861 -353406.30019 [49,] -149859.43355 28280.34861 [50,] 1830.56568 -149859.43355 [51,] -154893.59582 1830.56568 [52,] -163816.82868 -154893.59582 [53,] 108069.30970 -163816.82868 [54,] -398627.52899 108069.30970 [55,] 96414.42623 -398627.52899 [56,] 344119.88167 96414.42623 [57,] 65897.34530 344119.88167 [58,] -39937.02809 65897.34530 [59,] -268342.71677 -39937.02809 [60,] 30684.17243 -268342.71677 [61,] -165191.81686 30684.17243 [62,] 54293.11167 -165191.81686 [63,] -229559.34975 54293.11167 [64,] 44287.34309 -229559.34975 [65,] -30222.19008 44287.34309 [66,] -89225.87105 -30222.19008 [67,] 12463.09493 -89225.87105 [68,] -62985.93646 12463.09493 [69,] -110309.30452 -62985.93646 [70,] -206946.64271 -110309.30452 [71,] 90240.21711 -206946.64271 [72,] -88273.72819 90240.21711 [73,] -21316.72675 -88273.72819 [74,] -30243.60472 -21316.72675 [75,] -66184.57835 -30243.60472 [76,] -44934.51284 -66184.57835 [77,] -101207.15920 -44934.51284 [78,] -11678.29645 -101207.15920 [79,] 63643.05378 -11678.29645 [80,] 34641.61231 63643.05378 [81,] -73790.49159 34641.61231 [82,] -104685.81508 -73790.49159 [83,] 29937.67397 -104685.81508 [84,] -12725.14358 29937.67397 [85,] -167364.74107 -12725.14358 [86,] 91819.60560 -167364.74107 [87,] -38862.44579 91819.60560 [88,] -91394.42511 -38862.44579 [89,] -5660.57195 -91394.42511 [90,] 13129.61159 -5660.57195 [91,] -14033.50196 13129.61159 [92,] 137631.86438 -14033.50196 [93,] -89777.21321 137631.86438 [94,] 90481.59227 -89777.21321 [95,] 10678.39226 90481.59227 [96,] 19473.69877 10678.39226 [97,] 21519.12241 19473.69877 [98,] -48344.28494 21519.12241 [99,] 42245.68748 -48344.28494 [100,] -21816.82929 42245.68748 [101,] -133542.53336 -21816.82929 [102,] -108868.71359 -133542.53336 [103,] -98741.69172 -108868.71359 [104,] 39063.39389 -98741.69172 [105,] -59967.40542 39063.39389 [106,] -289570.08055 -59967.40542 [107,] 198107.49533 -289570.08055 [108,] -86218.57294 198107.49533 [109,] -66492.21629 -86218.57294 [110,] 83864.87293 -66492.21629 [111,] -222070.80804 83864.87293 [112,] -84451.81816 -222070.80804 [113,] -19997.95631 -84451.81816 [114,] 112369.79821 -19997.95631 [115,] 39988.90315 112369.79821 [116,] -24563.04824 39988.90315 [117,] 71140.54934 -24563.04824 [118,] -86911.82797 71140.54934 [119,] 3663.56051 -86911.82797 [120,] -65922.30668 3663.56051 [121,] 9489.98057 -65922.30668 [122,] -3418.89374 9489.98057 [123,] 152158.54615 -3418.89374 [124,] 29286.34957 152158.54615 [125,] -138918.74718 29286.34957 [126,] 4695.80258 -138918.74718 [127,] 5100.36378 4695.80258 [128,] 14817.85613 5100.36378 [129,] 4695.80258 14817.85613 [130,] 4973.42744 4695.80258 [131,] 65208.58145 4973.42744 [132,] 4695.80258 65208.58145 [133,] 5357.41047 4695.80258 [134,] -1120.95879 5357.41047 [135,] 3882.97829 -1120.95879 [136,] 9075.68325 3882.97829 [137,] 4364.04313 9075.68325 [138,] 444.25026 4364.04313 [139,] 4695.80258 444.25026 [140,] 17633.90268 4695.80258 [141,] 4695.80258 17633.90268 [142,] -7587.46063 4695.80258 [143,] 5383.56487 -7587.46063 [144,] 14587.67173 5383.56487 [145,] 10257.96358 14587.67173 [146,] 4695.80258 10257.96358 [147,] 5344.44639 4695.80258 [148,] 13569.03592 5344.44639 [149,] 4560.05297 13569.03592 [150,] -16424.65558 4560.05297 [151,] -4530.63985 -16424.65558 [152,] 5628.27428 -4530.63985 [153,] -9949.03909 5628.27428 [154,] 4418.84197 -9949.03909 [155,] 4695.80258 4418.84197 [156,] 4695.80258 4695.80258 [157,] 2967.75607 4695.80258 [158,] 24390.36985 2967.75607 [159,] 4695.80258 24390.36985 [160,] 4695.80258 4695.80258 [161,] 4695.80258 4695.80258 [162,] 4695.80258 4695.80258 [163,] -5759.11887 4695.80258 [164,] -11139.17833 -5759.11887 [165,] 4146.91131 -11139.17833 [166,] 4695.80258 4146.91131 [167,] 100101.84568 4695.80258 [168,] 17631.71406 100101.84568 [169,] 4695.80258 17631.71406 [170,] 4695.80258 4695.80258 [171,] 502.88764 4695.80258 [172,] 4695.80258 502.88764 [173,] -2747.76674 4695.80258 [174,] 4695.80258 -2747.76674 [175,] 5279.74355 4695.80258 [176,] 21636.35992 5279.74355 [177,] 95888.64101 21636.35992 [178,] 4695.80258 95888.64101 [179,] 4695.80258 4695.80258 [180,] 26546.37872 4695.80258 [181,] 69353.44777 26546.37872 [182,] 5194.43625 69353.44777 [183,] 27469.56632 5194.43625 [184,] 82604.48615 27469.56632 [185,] 71777.73429 82604.48615 [186,] 20793.09613 71777.73429 [187,] 505.91678 20793.09613 [188,] -28647.56861 505.91678 [189,] 4695.80258 -28647.56861 [190,] 10293.81670 4695.80258 [191,] 4695.80258 10293.81670 [192,] 3547.21593 4695.80258 [193,] 4695.80258 3547.21593 [194,] 5099.07171 4695.80258 [195,] 96745.06614 5099.07171 [196,] -159756.58648 96745.06614 [197,] 96745.06614 -159756.58648 [198,] 96745.06614 96745.06614 [199,] -11466.83054 96745.06614 [200,] 96745.06614 -11466.83054 [201,] 96745.06614 96745.06614 [202,] 96745.06614 96745.06614 [203,] 48580.16067 96745.06614 [204,] 60559.55258 48580.16067 [205,] 84396.23774 60559.55258 [206,] 92599.63230 84396.23774 [207,] 96745.06614 92599.63230 [208,] 34240.26905 96745.06614 [209,] 74513.23008 34240.26905 [210,] 4695.80258 74513.23008 [211,] 4695.80258 4695.80258 [212,] 25925.93172 4695.80258 [213,] -24615.08018 25925.93172 [214,] -11571.51316 -24615.08018 [215,] 5068.79126 -11571.51316 [216,] 5836.24263 5068.79126 [217,] 6492.62189 5836.24263 [218,] 5068.79126 6492.62189 [219,] 15978.94819 5068.79126 [220,] -14524.58874 15978.94819 [221,] 4695.80258 -14524.58874 [222,] 4882.91782 4695.80258 [223,] 2842.93476 4882.91782 [224,] 4815.88543 2842.93476 [225,] 5637.80160 4815.88543 [226,] 5237.15542 5637.80160 [227,] 4695.80258 5237.15542 [228,] -7592.79784 4695.80258 [229,] 4695.80258 -7592.79784 [230,] -18346.43070 4695.80258 [231,] 8642.79244 -18346.43070 [232,] 4695.80258 8642.79244 [233,] -8722.67173 4695.80258 [234,] 4695.80258 -8722.67173 [235,] 4749.02911 4695.80258 [236,] 4695.80258 4749.02911 [237,] 6664.52600 4695.80258 [238,] 4695.80258 6664.52600 [239,] 4695.80258 4695.80258 [240,] 4695.80258 4695.80258 [241,] 1072.29490 4695.80258 [242,] 4695.80258 1072.29490 [243,] 7849.74472 4695.80258 [244,] -20413.98966 7849.74472 [245,] 96745.06614 -20413.98966 [246,] 96745.06614 96745.06614 [247,] -1625.30469 96745.06614 [248,] 96745.06614 -1625.30469 [249,] 96745.06614 96745.06614 [250,] 1105.51459 96745.06614 [251,] 96745.06614 1105.51459 [252,] -10712.38062 96745.06614 [253,] -25693.29890 -10712.38062 [254,] -100020.13876 -25693.29890 [255,] 4695.80258 -100020.13876 [256,] 9540.54090 4695.80258 [257,] 15504.80685 9540.54090 [258,] 96745.06614 15504.80685 [259,] 4695.80258 96745.06614 [260,] 96745.06614 4695.80258 [261,] 96745.06614 96745.06614 [262,] 11524.03289 96745.06614 [263,] 25124.49007 11524.03289 [264,] 96745.06614 25124.49007 [265,] 4695.80258 96745.06614 [266,] 96745.06614 4695.80258 [267,] 39068.86267 96745.06614 [268,] 96745.06614 39068.86267 [269,] 7614.44532 96745.06614 [270,] 96745.06614 7614.44532 [271,] 96745.06614 96745.06614 [272,] 96745.06614 96745.06614 [273,] 2280.74545 96745.06614 [274,] 4695.80258 2280.74545 [275,] 6856.73328 4695.80258 [276,] 5750.88219 6856.73328 [277,] 4695.80258 5750.88219 [278,] -32509.02284 4695.80258 [279,] -13492.53170 -32509.02284 [280,] 96745.06614 -13492.53170 [281,] 4695.80258 96745.06614 [282,] 13981.74577 4695.80258 [283,] 4695.80258 13981.74577 [284,] 1474.75604 4695.80258 [285,] 4695.80258 1474.75604 [286,] 73875.87396 4695.80258 [287,] 5305.87880 73875.87396 [288,] 96745.06614 5305.87880 [289,] -4616.91129 96745.06614 [290,] 96745.06614 -4616.91129 [291,] 11869.04635 96745.06614 [292,] 17785.58720 11869.04635 [293,] 4695.80258 17785.58720 [294,] -4372.82474 4695.80258 [295,] 90735.14542 -4372.82474 [296,] 4695.80258 90735.14542 [297,] 96745.06614 4695.80258 [298,] 4695.80258 96745.06614 [299,] 6374.25164 4695.80258 [300,] 83624.34209 6374.25164 [301,] -7133.77660 83624.34209 [302,] 4882.29692 -7133.77660 [303,] 96745.06614 4882.29692 [304,] 32437.27291 96745.06614 [305,] 6909.15144 32437.27291 [306,] 96745.06614 6909.15144 [307,] 44212.05015 96745.06614 [308,] 11845.65282 44212.05015 [309,] 96745.06614 11845.65282 [310,] 96745.06614 96745.06614 [311,] -10213.42054 96745.06614 [312,] 96745.06614 -10213.42054 [313,] 66041.66334 96745.06614 [314,] 6370.75503 66041.66334 [315,] 18828.45632 6370.75503 [316,] 5508.08560 18828.45632 [317,] -12520.74431 5508.08560 [318,] 5617.12230 -12520.74431 [319,] 4695.80258 5617.12230 [320,] 4695.80258 4695.80258 [321,] 4695.80258 4695.80258 [322,] 4695.80258 4695.80258 [323,] 96745.06614 4695.80258 [324,] 4695.80258 96745.06614 [325,] 1889.38269 4695.80258 [326,] 96745.06614 1889.38269 [327,] 96745.06614 96745.06614 [328,] 4695.80258 96745.06614 [329,] 4695.80258 4695.80258 [330,] 4695.80258 4695.80258 [331,] 4695.80258 4695.80258 [332,] 4695.80258 4695.80258 [333,] 3953.39132 4695.80258 [334,] 4871.95470 3953.39132 [335,] 3846.84201 4871.95470 [336,] -15467.40013 3846.84201 [337,] 4695.80258 -15467.40013 [338,] 2916.79186 4695.80258 [339,] 10119.59876 2916.79186 [340,] -18672.05287 10119.59876 [341,] 4695.80258 -18672.05287 [342,] 4695.80258 4695.80258 [343,] 4695.80258 4695.80258 [344,] -3140.29155 4695.80258 [345,] 4695.80258 -3140.29155 [346,] -4006.90833 4695.80258 [347,] 4695.80258 -4006.90833 [348,] 4695.80258 4695.80258 [349,] 2493.78740 4695.80258 [350,] 11018.66125 2493.78740 [351,] 3600.11020 11018.66125 [352,] 16051.07917 3600.11020 [353,] 476.58495 16051.07917 [354,] 4695.80258 476.58495 [355,] 4695.80258 4695.80258 [356,] -22396.84628 4695.80258 [357,] 4695.80258 -22396.84628 [358,] 4695.80258 4695.80258 [359,] 4695.80258 4695.80258 [360,] 5450.93162 4695.80258 [361,] 4695.80258 5450.93162 [362,] 4695.80258 4695.80258 [363,] -22742.53307 4695.80258 [364,] 8531.29201 -22742.53307 [365,] 4695.80258 8531.29201 [366,] 4695.80258 4695.80258 [367,] 4695.80258 4695.80258 [368,] 4695.80258 4695.80258 [369,] 5073.19778 4695.80258 [370,] 4695.80258 5073.19778 [371,] 13502.47832 4695.80258 [372,] -4122.43855 13502.47832 [373,] -11364.54599 -4122.43855 [374,] 4695.80258 -11364.54599 [375,] 4695.80258 4695.80258 [376,] -5249.64067 4695.80258 [377,] 10434.70666 -5249.64067 [378,] -127.53011 10434.70666 [379,] 4695.80258 -127.53011 [380,] 4695.80258 4695.80258 [381,] -16852.39645 4695.80258 [382,] 5978.97241 -16852.39645 [383,] 4695.80258 5978.97241 [384,] 4730.19031 4695.80258 [385,] -93173.15808 4730.19031 [386,] 18295.44481 -93173.15808 [387,] -1201.64187 18295.44481 [388,] 4333.03834 -1201.64187 [389,] -87.88601 4333.03834 [390,] 7724.02211 -87.88601 [391,] -17053.27366 7724.02211 [392,] -49494.53350 -17053.27366 [393,] 40647.49443 -49494.53350 [394,] -40778.26494 40647.49443 [395,] -24144.92966 -40778.26494 [396,] -72002.46594 -24144.92966 [397,] -432306.25565 -72002.46594 [398,] -21770.74622 -432306.25565 [399,] -135659.00978 -21770.74622 [400,] -29640.24458 -135659.00978 [401,] 13728.65001 -29640.24458 [402,] -12827.11030 13728.65001 [403,] -29957.78483 -12827.11030 [404,] -42211.45613 -29957.78483 [405,] -6534.19128 -42211.45613 [406,] -8555.89563 -6534.19128 [407,] -5487.64777 -8555.89563 [408,] -8564.90722 -5487.64777 [409,] -99513.27125 -8564.90722 [410,] -14573.70679 -99513.27125 [411,] -157855.32051 -14573.70679 [412,] -13610.02368 -157855.32051 [413,] 68913.23945 -13610.02368 [414,] -188433.40848 68913.23945 [415,] -135107.61599 -188433.40848 [416,] 132680.84269 -135107.61599 [417,] -5736.85380 132680.84269 [418,] -53424.84627 -5736.85380 [419,] -77168.91215 -53424.84627 [420,] -79330.08351 -77168.91215 [421,] -93066.08147 -79330.08351 [422,] -731095.45915 -93066.08147 [423,] -372773.51927 -731095.45915 [424,] -35145.58987 -372773.51927 [425,] -331902.94857 -35145.58987 [426,] -266234.14336 -331902.94857 [427,] -322175.21212 -266234.14336 [428,] -133587.07324 -322175.21212 [429,] -30715.63764 -133587.07324 [430,] -332780.81006 -30715.63764 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 2275662.88080 208308.99237 2 795828.19898 2275662.88080 3 -270987.19826 795828.19898 4 -809715.40781 -270987.19826 5 -22266.82622 -809715.40781 6 800116.23761 -22266.82622 7 -49155.83317 800116.23761 8 -691527.85593 -49155.83317 9 520274.61695 -691527.85593 10 300395.15083 520274.61695 11 657029.06624 300395.15083 12 -223632.43021 657029.06624 13 -214806.51231 -223632.43021 14 -584732.95539 -214806.51231 15 567057.94638 -584732.95539 16 -120685.92990 567057.94638 17 -555208.39746 -120685.92990 18 -692034.64645 -555208.39746 19 -22517.65388 -692034.64645 20 147556.65102 -22517.65388 21 -198905.77113 147556.65102 22 144552.96763 -198905.77113 23 -85472.26929 144552.96763 24 -172515.67738 -85472.26929 25 199630.89796 -172515.67738 26 -806793.29430 199630.89796 27 -91537.29089 -806793.29430 28 493285.69665 -91537.29089 29 102851.24420 493285.69665 30 321503.16649 102851.24420 31 -104258.33967 321503.16649 32 255545.92790 -104258.33967 33 160686.08391 255545.92790 34 -132980.09112 160686.08391 35 -66128.03087 -132980.09112 36 238335.88236 -66128.03087 37 142898.18103 238335.88236 38 44953.18110 142898.18103 39 -135648.29772 44953.18110 40 -143145.37863 -135648.29772 41 -91583.43595 -143145.37863 42 121975.57705 -91583.43595 43 -56834.07852 121975.57705 44 44665.36706 -56834.07852 45 2854.95073 44665.36706 46 -109747.17134 2854.95073 47 -353406.30019 -109747.17134 48 28280.34861 -353406.30019 49 -149859.43355 28280.34861 50 1830.56568 -149859.43355 51 -154893.59582 1830.56568 52 -163816.82868 -154893.59582 53 108069.30970 -163816.82868 54 -398627.52899 108069.30970 55 96414.42623 -398627.52899 56 344119.88167 96414.42623 57 65897.34530 344119.88167 58 -39937.02809 65897.34530 59 -268342.71677 -39937.02809 60 30684.17243 -268342.71677 61 -165191.81686 30684.17243 62 54293.11167 -165191.81686 63 -229559.34975 54293.11167 64 44287.34309 -229559.34975 65 -30222.19008 44287.34309 66 -89225.87105 -30222.19008 67 12463.09493 -89225.87105 68 -62985.93646 12463.09493 69 -110309.30452 -62985.93646 70 -206946.64271 -110309.30452 71 90240.21711 -206946.64271 72 -88273.72819 90240.21711 73 -21316.72675 -88273.72819 74 -30243.60472 -21316.72675 75 -66184.57835 -30243.60472 76 -44934.51284 -66184.57835 77 -101207.15920 -44934.51284 78 -11678.29645 -101207.15920 79 63643.05378 -11678.29645 80 34641.61231 63643.05378 81 -73790.49159 34641.61231 82 -104685.81508 -73790.49159 83 29937.67397 -104685.81508 84 -12725.14358 29937.67397 85 -167364.74107 -12725.14358 86 91819.60560 -167364.74107 87 -38862.44579 91819.60560 88 -91394.42511 -38862.44579 89 -5660.57195 -91394.42511 90 13129.61159 -5660.57195 91 -14033.50196 13129.61159 92 137631.86438 -14033.50196 93 -89777.21321 137631.86438 94 90481.59227 -89777.21321 95 10678.39226 90481.59227 96 19473.69877 10678.39226 97 21519.12241 19473.69877 98 -48344.28494 21519.12241 99 42245.68748 -48344.28494 100 -21816.82929 42245.68748 101 -133542.53336 -21816.82929 102 -108868.71359 -133542.53336 103 -98741.69172 -108868.71359 104 39063.39389 -98741.69172 105 -59967.40542 39063.39389 106 -289570.08055 -59967.40542 107 198107.49533 -289570.08055 108 -86218.57294 198107.49533 109 -66492.21629 -86218.57294 110 83864.87293 -66492.21629 111 -222070.80804 83864.87293 112 -84451.81816 -222070.80804 113 -19997.95631 -84451.81816 114 112369.79821 -19997.95631 115 39988.90315 112369.79821 116 -24563.04824 39988.90315 117 71140.54934 -24563.04824 118 -86911.82797 71140.54934 119 3663.56051 -86911.82797 120 -65922.30668 3663.56051 121 9489.98057 -65922.30668 122 -3418.89374 9489.98057 123 152158.54615 -3418.89374 124 29286.34957 152158.54615 125 -138918.74718 29286.34957 126 4695.80258 -138918.74718 127 5100.36378 4695.80258 128 14817.85613 5100.36378 129 4695.80258 14817.85613 130 4973.42744 4695.80258 131 65208.58145 4973.42744 132 4695.80258 65208.58145 133 5357.41047 4695.80258 134 -1120.95879 5357.41047 135 3882.97829 -1120.95879 136 9075.68325 3882.97829 137 4364.04313 9075.68325 138 444.25026 4364.04313 139 4695.80258 444.25026 140 17633.90268 4695.80258 141 4695.80258 17633.90268 142 -7587.46063 4695.80258 143 5383.56487 -7587.46063 144 14587.67173 5383.56487 145 10257.96358 14587.67173 146 4695.80258 10257.96358 147 5344.44639 4695.80258 148 13569.03592 5344.44639 149 4560.05297 13569.03592 150 -16424.65558 4560.05297 151 -4530.63985 -16424.65558 152 5628.27428 -4530.63985 153 -9949.03909 5628.27428 154 4418.84197 -9949.03909 155 4695.80258 4418.84197 156 4695.80258 4695.80258 157 2967.75607 4695.80258 158 24390.36985 2967.75607 159 4695.80258 24390.36985 160 4695.80258 4695.80258 161 4695.80258 4695.80258 162 4695.80258 4695.80258 163 -5759.11887 4695.80258 164 -11139.17833 -5759.11887 165 4146.91131 -11139.17833 166 4695.80258 4146.91131 167 100101.84568 4695.80258 168 17631.71406 100101.84568 169 4695.80258 17631.71406 170 4695.80258 4695.80258 171 502.88764 4695.80258 172 4695.80258 502.88764 173 -2747.76674 4695.80258 174 4695.80258 -2747.76674 175 5279.74355 4695.80258 176 21636.35992 5279.74355 177 95888.64101 21636.35992 178 4695.80258 95888.64101 179 4695.80258 4695.80258 180 26546.37872 4695.80258 181 69353.44777 26546.37872 182 5194.43625 69353.44777 183 27469.56632 5194.43625 184 82604.48615 27469.56632 185 71777.73429 82604.48615 186 20793.09613 71777.73429 187 505.91678 20793.09613 188 -28647.56861 505.91678 189 4695.80258 -28647.56861 190 10293.81670 4695.80258 191 4695.80258 10293.81670 192 3547.21593 4695.80258 193 4695.80258 3547.21593 194 5099.07171 4695.80258 195 96745.06614 5099.07171 196 -159756.58648 96745.06614 197 96745.06614 -159756.58648 198 96745.06614 96745.06614 199 -11466.83054 96745.06614 200 96745.06614 -11466.83054 201 96745.06614 96745.06614 202 96745.06614 96745.06614 203 48580.16067 96745.06614 204 60559.55258 48580.16067 205 84396.23774 60559.55258 206 92599.63230 84396.23774 207 96745.06614 92599.63230 208 34240.26905 96745.06614 209 74513.23008 34240.26905 210 4695.80258 74513.23008 211 4695.80258 4695.80258 212 25925.93172 4695.80258 213 -24615.08018 25925.93172 214 -11571.51316 -24615.08018 215 5068.79126 -11571.51316 216 5836.24263 5068.79126 217 6492.62189 5836.24263 218 5068.79126 6492.62189 219 15978.94819 5068.79126 220 -14524.58874 15978.94819 221 4695.80258 -14524.58874 222 4882.91782 4695.80258 223 2842.93476 4882.91782 224 4815.88543 2842.93476 225 5637.80160 4815.88543 226 5237.15542 5637.80160 227 4695.80258 5237.15542 228 -7592.79784 4695.80258 229 4695.80258 -7592.79784 230 -18346.43070 4695.80258 231 8642.79244 -18346.43070 232 4695.80258 8642.79244 233 -8722.67173 4695.80258 234 4695.80258 -8722.67173 235 4749.02911 4695.80258 236 4695.80258 4749.02911 237 6664.52600 4695.80258 238 4695.80258 6664.52600 239 4695.80258 4695.80258 240 4695.80258 4695.80258 241 1072.29490 4695.80258 242 4695.80258 1072.29490 243 7849.74472 4695.80258 244 -20413.98966 7849.74472 245 96745.06614 -20413.98966 246 96745.06614 96745.06614 247 -1625.30469 96745.06614 248 96745.06614 -1625.30469 249 96745.06614 96745.06614 250 1105.51459 96745.06614 251 96745.06614 1105.51459 252 -10712.38062 96745.06614 253 -25693.29890 -10712.38062 254 -100020.13876 -25693.29890 255 4695.80258 -100020.13876 256 9540.54090 4695.80258 257 15504.80685 9540.54090 258 96745.06614 15504.80685 259 4695.80258 96745.06614 260 96745.06614 4695.80258 261 96745.06614 96745.06614 262 11524.03289 96745.06614 263 25124.49007 11524.03289 264 96745.06614 25124.49007 265 4695.80258 96745.06614 266 96745.06614 4695.80258 267 39068.86267 96745.06614 268 96745.06614 39068.86267 269 7614.44532 96745.06614 270 96745.06614 7614.44532 271 96745.06614 96745.06614 272 96745.06614 96745.06614 273 2280.74545 96745.06614 274 4695.80258 2280.74545 275 6856.73328 4695.80258 276 5750.88219 6856.73328 277 4695.80258 5750.88219 278 -32509.02284 4695.80258 279 -13492.53170 -32509.02284 280 96745.06614 -13492.53170 281 4695.80258 96745.06614 282 13981.74577 4695.80258 283 4695.80258 13981.74577 284 1474.75604 4695.80258 285 4695.80258 1474.75604 286 73875.87396 4695.80258 287 5305.87880 73875.87396 288 96745.06614 5305.87880 289 -4616.91129 96745.06614 290 96745.06614 -4616.91129 291 11869.04635 96745.06614 292 17785.58720 11869.04635 293 4695.80258 17785.58720 294 -4372.82474 4695.80258 295 90735.14542 -4372.82474 296 4695.80258 90735.14542 297 96745.06614 4695.80258 298 4695.80258 96745.06614 299 6374.25164 4695.80258 300 83624.34209 6374.25164 301 -7133.77660 83624.34209 302 4882.29692 -7133.77660 303 96745.06614 4882.29692 304 32437.27291 96745.06614 305 6909.15144 32437.27291 306 96745.06614 6909.15144 307 44212.05015 96745.06614 308 11845.65282 44212.05015 309 96745.06614 11845.65282 310 96745.06614 96745.06614 311 -10213.42054 96745.06614 312 96745.06614 -10213.42054 313 66041.66334 96745.06614 314 6370.75503 66041.66334 315 18828.45632 6370.75503 316 5508.08560 18828.45632 317 -12520.74431 5508.08560 318 5617.12230 -12520.74431 319 4695.80258 5617.12230 320 4695.80258 4695.80258 321 4695.80258 4695.80258 322 4695.80258 4695.80258 323 96745.06614 4695.80258 324 4695.80258 96745.06614 325 1889.38269 4695.80258 326 96745.06614 1889.38269 327 96745.06614 96745.06614 328 4695.80258 96745.06614 329 4695.80258 4695.80258 330 4695.80258 4695.80258 331 4695.80258 4695.80258 332 4695.80258 4695.80258 333 3953.39132 4695.80258 334 4871.95470 3953.39132 335 3846.84201 4871.95470 336 -15467.40013 3846.84201 337 4695.80258 -15467.40013 338 2916.79186 4695.80258 339 10119.59876 2916.79186 340 -18672.05287 10119.59876 341 4695.80258 -18672.05287 342 4695.80258 4695.80258 343 4695.80258 4695.80258 344 -3140.29155 4695.80258 345 4695.80258 -3140.29155 346 -4006.90833 4695.80258 347 4695.80258 -4006.90833 348 4695.80258 4695.80258 349 2493.78740 4695.80258 350 11018.66125 2493.78740 351 3600.11020 11018.66125 352 16051.07917 3600.11020 353 476.58495 16051.07917 354 4695.80258 476.58495 355 4695.80258 4695.80258 356 -22396.84628 4695.80258 357 4695.80258 -22396.84628 358 4695.80258 4695.80258 359 4695.80258 4695.80258 360 5450.93162 4695.80258 361 4695.80258 5450.93162 362 4695.80258 4695.80258 363 -22742.53307 4695.80258 364 8531.29201 -22742.53307 365 4695.80258 8531.29201 366 4695.80258 4695.80258 367 4695.80258 4695.80258 368 4695.80258 4695.80258 369 5073.19778 4695.80258 370 4695.80258 5073.19778 371 13502.47832 4695.80258 372 -4122.43855 13502.47832 373 -11364.54599 -4122.43855 374 4695.80258 -11364.54599 375 4695.80258 4695.80258 376 -5249.64067 4695.80258 377 10434.70666 -5249.64067 378 -127.53011 10434.70666 379 4695.80258 -127.53011 380 4695.80258 4695.80258 381 -16852.39645 4695.80258 382 5978.97241 -16852.39645 383 4695.80258 5978.97241 384 4730.19031 4695.80258 385 -93173.15808 4730.19031 386 18295.44481 -93173.15808 387 -1201.64187 18295.44481 388 4333.03834 -1201.64187 389 -87.88601 4333.03834 390 7724.02211 -87.88601 391 -17053.27366 7724.02211 392 -49494.53350 -17053.27366 393 40647.49443 -49494.53350 394 -40778.26494 40647.49443 395 -24144.92966 -40778.26494 396 -72002.46594 -24144.92966 397 -432306.25565 -72002.46594 398 -21770.74622 -432306.25565 399 -135659.00978 -21770.74622 400 -29640.24458 -135659.00978 401 13728.65001 -29640.24458 402 -12827.11030 13728.65001 403 -29957.78483 -12827.11030 404 -42211.45613 -29957.78483 405 -6534.19128 -42211.45613 406 -8555.89563 -6534.19128 407 -5487.64777 -8555.89563 408 -8564.90722 -5487.64777 409 -99513.27125 -8564.90722 410 -14573.70679 -99513.27125 411 -157855.32051 -14573.70679 412 -13610.02368 -157855.32051 413 68913.23945 -13610.02368 414 -188433.40848 68913.23945 415 -135107.61599 -188433.40848 416 132680.84269 -135107.61599 417 -5736.85380 132680.84269 418 -53424.84627 -5736.85380 419 -77168.91215 -53424.84627 420 -79330.08351 -77168.91215 421 -93066.08147 -79330.08351 422 -731095.45915 -93066.08147 423 -372773.51927 -731095.45915 424 -35145.58987 -372773.51927 425 -331902.94857 -35145.58987 426 -266234.14336 -331902.94857 427 -322175.21212 -266234.14336 428 -133587.07324 -322175.21212 429 -30715.63764 -133587.07324 430 -332780.81006 -30715.63764 > plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals') > lines(lowess(z)) > abline(lm(z)) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/73qfr1291225443.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/83qfr1291225443.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/93qfr1291225443.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0)) > plot(mylm, las = 1, sub='Residual Diagnostics') > par(opar) > dev.off() null device 1 > if (n > n25) { + postscript(file="/var/www/html/rcomp/tmp/10w0fd1291225443.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) + plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint') + grid() + 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, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE) > a<-table.row.end(a) > myeq <- colnames(x)[1] > myeq <- paste(myeq, '[t] = ', sep='') > for (i in 1:k){ + if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '') + myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ') + if (rownames(mysum$coefficients)[i] != '(Intercept)') { + myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='') + if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='') + } + } > myeq <- paste(myeq, ' + e[t]') > a<-table.row.start(a) > a<-table.element(a, myeq) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/11zivi1291225443.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Variable',header=TRUE) > a<-table.element(a,'Parameter',header=TRUE) > a<-table.element(a,'S.D.',header=TRUE) > a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE) > a<-table.element(a,'2-tail p-value',header=TRUE) > a<-table.element(a,'1-tail p-value',header=TRUE) > a<-table.row.end(a) > for (i in 1:k){ + a<-table.row.start(a) + a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE) + a<-table.element(a,mysum$coefficients[i,1]) + a<-table.element(a, round(mysum$coefficients[i,2],6)) + a<-table.element(a, round(mysum$coefficients[i,3],4)) + a<-table.element(a, round(mysum$coefficients[i,4],6)) + a<-table.element(a, round(mysum$coefficients[i,4]/2,6)) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/12k1u61291225443.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple R',1,TRUE) > a<-table.element(a, sqrt(mysum$r.squared)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'R-squared',1,TRUE) > a<-table.element(a, mysum$r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Adjusted R-squared',1,TRUE) > a<-table.element(a, mysum$adj.r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (value)',1,TRUE) > a<-table.element(a, mysum$fstatistic[1]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[2]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[3]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'p-value',1,TRUE) > a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3])) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Residual Standard Deviation',1,TRUE) > a<-table.element(a, mysum$sigma) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Sum Squared Residuals',1,TRUE) > a<-table.element(a, sum(myerror*myerror)) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/13yssx1291225443.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Time or Index', 1, TRUE) > a<-table.element(a, 'Actuals', 1, TRUE) > a<-table.element(a, 'Interpolation
Forecast', 1, TRUE) > a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,i, 1, TRUE) + a<-table.element(a,x[i]) + a<-table.element(a,x[i]-mysum$resid[i]) + a<-table.element(a,mysum$resid[i]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/142t831291225443.tab") > if (n > n25) { + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'p-values',header=TRUE) + a<-table.element(a,'Alternative Hypothesis',3,header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'breakpoint index',header=TRUE) + a<-table.element(a,'greater',header=TRUE) + a<-table.element(a,'2-sided',header=TRUE) + a<-table.element(a,'less',header=TRUE) + a<-table.row.end(a) + for (mypoint in kp3:nmkm3) { + a<-table.row.start(a) + a<-table.element(a,mypoint,header=TRUE) + a<-table.element(a,gqarr[mypoint-kp3+1,1]) + a<-table.element(a,gqarr[mypoint-kp3+1,2]) + a<-table.element(a,gqarr[mypoint-kp3+1,3]) + a<-table.row.end(a) + } + a<-table.end(a) + table.save(a,file="/var/www/html/rcomp/tmp/155tp91291225443.tab") + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'Description',header=TRUE) + a<-table.element(a,'# significant tests',header=TRUE) + a<-table.element(a,'% significant tests',header=TRUE) + a<-table.element(a,'OK/NOK',header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'1% type I error level',header=TRUE) + a<-table.element(a,numsignificant1) + a<-table.element(a,numsignificant1/numgqtests) + if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'5% type I error level',header=TRUE) + a<-table.element(a,numsignificant5) + a<-table.element(a,numsignificant5/numgqtests) + if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'10% type I error level',header=TRUE) + a<-table.element(a,numsignificant10) + a<-table.element(a,numsignificant10/numgqtests) + if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.end(a) + table.save(a,file="/var/www/html/rcomp/tmp/169c5f1291225443.tab") + } > > try(system("convert tmp/1ph0j1291225443.ps tmp/1ph0j1291225443.png",intern=TRUE)) character(0) > try(system("convert tmp/2ph0j1291225443.ps tmp/2ph0j1291225443.png",intern=TRUE)) character(0) > try(system("convert tmp/3ph0j1291225443.ps tmp/3ph0j1291225443.png",intern=TRUE)) character(0) > try(system("convert tmp/40qhm1291225443.ps tmp/40qhm1291225443.png",intern=TRUE)) character(0) > try(system("convert tmp/50qhm1291225443.ps tmp/50qhm1291225443.png",intern=TRUE)) character(0) > try(system("convert tmp/6shy71291225443.ps tmp/6shy71291225443.png",intern=TRUE)) character(0) > try(system("convert tmp/73qfr1291225443.ps tmp/73qfr1291225443.png",intern=TRUE)) character(0) > try(system("convert tmp/83qfr1291225443.ps tmp/83qfr1291225443.png",intern=TRUE)) character(0) > try(system("convert tmp/93qfr1291225443.ps tmp/93qfr1291225443.png",intern=TRUE)) character(0) > try(system("convert tmp/10w0fd1291225443.ps tmp/10w0fd1291225443.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 14.328 2.371 61.470