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(0 + ,162556 + ,0 + ,807 + ,0 + ,213118 + ,0 + ,6282154 + ,0 + ,29790 + ,0 + ,444 + ,0 + ,81767 + ,0 + ,4321023 + ,0 + ,87550 + ,0 + ,412 + ,0 + ,153198 + ,0 + ,4111912 + ,1 + ,84738 + ,254214 + ,428 + ,2140 + ,-26007 + ,-52014 + ,223193 + ,0 + ,54660 + ,0 + ,315 + ,0 + ,126942 + ,0 + ,1491348 + ,0 + ,42634 + ,0 + ,168 + ,0 + ,157214 + ,0 + ,1629616 + ,1 + ,40949 + ,122847 + ,263 + ,1315 + ,129352 + ,258704 + ,1398893 + ,0 + ,45187 + ,0 + ,267 + ,0 + ,234817 + ,0 + ,1926517 + ,0 + ,37704 + ,0 + ,228 + ,0 + ,60448 + ,0 + ,983660 + ,0 + ,16275 + ,0 + ,129 + ,0 + ,47818 + ,0 + ,1443586 + ,1 + ,25830 + ,77490 + ,104 + ,520 + ,245546 + ,491092 + ,1073089 + ,1 + ,12679 + ,38037 + ,122 + ,610 + ,48020 + ,96040 + ,984885 + ,0 + ,18014 + ,0 + ,393 + ,0 + ,-1710 + ,0 + ,1405225 + ,1 + ,43556 + ,130668 + ,190 + ,950 + ,32648 + ,65296 + ,227132 + ,0 + ,24811 + ,0 + ,280 + ,0 + ,95350 + ,0 + ,929118 + ,1 + ,6575 + ,19725 + ,63 + ,315 + ,151352 + ,302704 + ,1071292 + ,1 + ,7123 + ,21369 + ,102 + ,510 + ,288170 + ,576340 + ,638830 + ,0 + ,21950 + ,0 + ,265 + ,0 + ,114337 + ,0 + ,856956 + ,0 + ,37597 + ,0 + ,234 + ,0 + ,37884 + ,0 + ,992426 + ,1 + ,17821 + ,53463 + ,277 + ,1385 + ,122844 + ,245688 + ,444477 + ,0 + ,12988 + ,0 + ,73 + ,0 + ,82340 + ,0 + ,857217 + ,0 + ,22330 + ,0 + ,67 + ,0 + ,79801 + ,0 + ,711969 + ,1 + ,13326 + ,39978 + ,103 + ,515 + ,165548 + ,331096 + ,702380 + ,1 + ,16189 + ,48567 + ,290 + ,1450 + ,116384 + ,232768 + ,358589 + ,1 + ,7146 + ,21438 + ,83 + ,415 + ,134028 + ,268056 + ,297978 + ,1 + ,15824 + ,47472 + ,56 + ,280 + ,63838 + ,127676 + ,585715 + ,0 + ,27664 + ,0 + ,236 + ,0 + ,74996 + ,0 + ,657954 + ,1 + ,11920 + ,35760 + ,73 + ,365 + ,31080 + ,62160 + ,209458 + ,1 + ,8568 + ,25704 + ,34 + ,170 + ,32168 + ,64336 + ,786690 + ,1 + ,14416 + ,43248 + ,139 + ,695 + ,49857 + ,99714 + ,439798 + ,0 + ,3369 + ,0 + ,26 + ,0 + ,87161 + ,0 + ,688779 + ,0 + ,11819 + ,0 + ,70 + ,0 + ,106113 + ,0 + ,574339 + ,0 + ,6984 + ,0 + ,40 + ,0 + ,80570 + ,0 + ,741409 + ,0 + ,4519 + ,0 + ,42 + ,0 + ,102129 + ,0 + ,597793 + ,1 + ,2220 + ,6660 + ,12 + ,60 + ,301670 + ,603340 + ,644190 + ,1 + ,18562 + ,55686 + ,211 + ,1055 + ,102313 + ,204626 + ,377934 + ,1 + ,10327 + ,30981 + ,74 + ,370 + ,88577 + ,177154 + ,640273 + ,0 + ,5336 + ,0 + ,80 + ,0 + ,112477 + ,0 + ,697458 + ,0 + ,2365 + ,0 + ,83 + ,0 + ,191778 + ,0 + ,550608 + ,1 + ,4069 + ,12207 + ,131 + ,655 + ,79804 + ,159608 + ,207393 + ,1 + ,8636 + ,25908 + ,203 + ,1015 + ,128294 + ,256588 + ,301607 + ,1 + ,13718 + ,41154 + ,56 + ,280 + ,96448 + ,192896 + ,345783 + ,1 + ,4525 + ,13575 + ,89 + ,445 + ,93811 + ,187622 + ,501749 + ,1 + ,6869 + ,20607 + ,88 + ,440 + ,117520 + ,235040 + ,379983 + ,1 + ,4628 + ,13884 + ,39 + ,195 + ,69159 + ,138318 + ,387475 + ,0 + ,3689 + ,0 + ,25 + ,0 + ,101792 + ,0 + ,377305 + ,0 + ,4891 + ,0 + ,49 + ,0 + ,210568 + ,0 + ,370837 + ,0 + ,7489 + ,0 + ,149 + ,0 + ,136996 + ,0 + ,430866 + ,1 + ,4901 + ,14703 + ,58 + ,290 + ,121920 + ,243840 + ,469107 + ,1 + ,2284 + ,6852 + ,41 + ,205 + ,76403 + ,152806 + ,194493 + ,0 + ,3160 + ,0 + ,90 + ,0 + ,108094 + ,0 + ,530670 + ,0 + ,4150 + ,0 + ,136 + ,0 + ,134759 + ,0 + ,518365 + ,0 + ,7285 + ,0 + ,97 + ,0 + ,188873 + ,0 + ,491303 + ,0 + ,1134 + ,0 + ,63 + ,0 + ,146216 + ,0 + ,527021 + ,0 + ,4658 + ,0 + ,114 + ,0 + ,156608 + ,0 + ,233773 + ,1 + ,2384 + ,7152 + ,77 + ,385 + ,61348 + ,122696 + ,405972 + ,1 + ,3748 + ,11244 + ,6 + ,30 + ,50350 + ,100700 + ,652925 + ,1 + ,5371 + ,16113 + ,47 + ,235 + ,87720 + ,175440 + ,446211 + ,1 + ,1285 + ,3855 + ,51 + ,255 + ,99489 + ,198978 + ,341340 + ,0 + ,9327 + ,0 + ,85 + ,0 + ,87419 + ,0 + ,387699 + ,0 + ,5565 + ,0 + ,43 + ,0 + ,94355 + ,0 + ,493408 + ,1 + ,1528 + ,4584 + ,32 + ,160 + ,60326 + ,120652 + ,146494 + ,0 + ,3122 + ,0 + ,25 + ,0 + ,94670 + ,0 + ,414462 + ,0 + ,7561 + ,0 + ,77 + ,0 + ,82425 + ,0 + ,364304 + ,1 + ,2675 + ,8025 + ,54 + ,270 + ,59017 + ,118034 + ,355178 + ,1 + ,13253 + ,39759 + ,251 + ,1255 + ,90829 + ,181658 + ,357760 + ,1 + ,880 + ,2640 + ,15 + ,75 + ,80791 + ,161582 + ,261216 + ,0 + ,2053 + ,0 + ,44 + ,0 + ,100423 + ,0 + ,397144 + ,1 + ,1424 + ,4272 + ,73 + ,365 + ,131116 + ,262232 + ,374943 + ,0 + ,4036 + ,0 + ,85 + ,0 + ,100269 + ,0 + ,424898 + ,0 + ,3045 + ,0 + ,49 + ,0 + ,27330 + ,0 + ,202055 + ,1 + ,5119 + ,15357 + ,38 + ,190 + ,39039 + ,78078 + ,378525 + ,1 + ,1431 + ,4293 + ,35 + ,175 + ,106885 + ,213770 + ,310768 + ,1 + ,554 + ,1662 + ,9 + ,45 + ,79285 + ,158570 + ,325738 + ,1 + ,1975 + ,5925 + ,34 + ,170 + ,118881 + ,237762 + ,394510 + ,0 + ,1765 + ,0 + ,20 + ,0 + ,77623 + ,0 + ,247060 + ,1 + ,1012 + ,3036 + ,29 + ,145 + ,114768 + ,229536 + ,368078 + ,1 + ,810 + ,2430 + ,11 + ,55 + ,74015 + ,148030 + ,236761 + ,1 + ,1280 + ,3840 + ,52 + ,260 + ,69465 + ,138930 + ,312378 + ,0 + ,666 + ,0 + ,13 + ,0 + ,117869 + ,0 + ,339836 + ,1 + ,1380 + ,4140 + ,29 + ,145 + ,60982 + ,121964 + ,347385 + ,0 + ,4677 + ,0 + ,66 + ,0 + ,90131 + ,0 + ,426280 + ,1 + ,876 + ,2628 + ,33 + ,165 + ,138971 + ,277942 + ,352850 + ,1 + ,814 + ,2442 + ,15 + ,75 + ,39625 + ,79250 + ,301881 + ,1 + ,514 + ,1542 + ,15 + ,75 + ,102725 + ,205450 + ,377516 + ,0 + ,5692 + ,0 + ,68 + ,0 + ,64239 + ,0 + ,357312 + ,1 + ,3642 + ,10926 + ,100 + ,500 + ,90262 + ,180524 + ,458343 + ,1 + ,540 + ,1620 + ,13 + ,65 + ,103960 + ,207920 + ,354228 + ,1 + ,2099 + ,6297 + ,45 + ,225 + ,106611 + ,213222 + ,308636 + ,1 + ,567 + ,1701 + ,14 + ,70 + ,103345 + ,206690 + ,386212 + ,1 + ,2001 + ,6003 + ,36 + ,180 + ,95551 + ,191102 + ,393343 + ,0 + ,2949 + ,0 + ,40 + ,0 + ,82903 + ,0 + ,378509 + ,1 + ,2253 + ,6759 + ,68 + ,340 + ,63593 + ,127186 + ,452469 + ,0 + ,6533 + ,0 + ,29 + ,0 + ,126910 + ,0 + ,364839 + ,1 + ,1889 + ,5667 + ,43 + ,215 + ,37527 + ,75054 + ,358649 + ,0 + ,3055 + ,0 + ,30 + ,0 + ,60247 + ,0 + ,376641 + ,1 + ,272 + ,816 + ,9 + ,45 + ,112995 + ,225990 + ,429112 + ,0 + ,1414 + ,0 + ,22 + ,0 + ,70184 + ,0 + ,330546 + ,1 + ,2564 + ,7692 + ,19 + ,95 + ,130140 + ,260280 + ,403560 + ,0 + ,1383 + ,0 + ,9 + ,0 + ,73221 + ,0 + ,317892 + ,0 + ,1261 + ,0 + ,31 + ,0 + ,76114 + ,0 + ,307528 + ,1 + ,975 + ,2925 + ,19 + ,95 + ,90534 + ,181068 + ,235133 + ,1 + ,3366 + ,10098 + ,55 + ,275 + ,108479 + ,216958 + ,299243 + ,1 + ,576 + ,1728 + ,8 + ,40 + ,113761 + ,227522 + ,314073 + ,1 + ,1686 + ,5058 + ,28 + ,140 + ,68696 + ,137392 + ,368186 + ,1 + ,746 + ,2238 + ,29 + ,145 + ,71561 + ,143122 + ,269661 + ,0 + ,3192 + ,0 + ,48 + ,0 + ,59831 + ,0 + ,125390 + ,0 + ,2045 + ,0 + ,16 + ,0 + ,97890 + ,0 + ,510834 + ,1 + ,5702 + ,17106 + ,47 + ,235 + ,101481 + ,202962 + ,321896 + ,0 + ,1932 + ,0 + ,20 + ,0 + ,72954 + ,0 + ,249898 + ,1 + ,936 + ,2808 + ,22 + ,110 + ,67939 + ,135878 + ,408881 + ,0 + ,3437 + ,0 + ,33 + ,0 + ,48022 + ,0 + ,158492 + ,1 + ,5131 + ,15393 + ,44 + ,220 + ,86111 + ,172222 + ,292154 + ,0 + ,2397 + ,0 + ,13 + ,0 + ,74020 + ,0 + ,289513 + ,0 + ,1389 + ,0 + ,6 + ,0 + ,57530 + ,0 + ,378049 + ,1 + ,1503 + ,4509 + ,35 + ,175 + ,56364 + ,112728 + ,343466 + ,1 + ,402 + ,1206 + ,8 + ,40 + ,84990 + ,169980 + ,332743 + ,1 + ,2239 + ,6717 + ,17 + ,85 + ,88590 + ,177180 + ,442882 + ,0 + ,2234 + ,0 + ,11 + ,0 + ,77200 + ,0 + ,214215 + ,1 + ,837 + ,2511 + ,21 + ,105 + ,61262 + ,122524 + ,315688 + ,1 + ,10579 + ,31737 + ,92 + ,460 + ,110309 + ,220618 + ,375195 + ,1 + ,875 + ,2625 + ,12 + ,60 + ,67000 + ,134000 + ,334280 + ,1 + ,1585 + ,4755 + ,112 + ,560 + ,93099 + ,186198 + ,355864 + ,0 + ,1659 + ,0 + ,25 + ,0 + ,107577 + ,0 + ,480382 + ,0 + ,2647 + ,0 + ,17 + ,0 + ,62920 + ,0 + ,353058 + ,0 + ,3294 + ,0 + ,23 + ,0 + ,75832 + ,0 + ,217193 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,94 + ,282 + ,10 + ,50 + ,60793 + ,121586 + ,314533 + ,1 + ,422 + ,1266 + ,23 + ,115 + ,57935 + ,115870 + ,318056 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,34 + ,102 + ,7 + ,35 + ,60630 + ,121260 + ,314353 + ,1 + ,1558 + ,4674 + ,25 + ,125 + ,55637 + ,111274 + ,369448 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,43 + ,129 + ,20 + ,100 + ,60887 + ,121774 + ,312846 + ,1 + ,645 + ,1935 + ,4 + ,20 + ,60720 + ,121440 + ,312075 + ,1 + ,316 + ,948 + ,4 + ,20 + ,60505 + ,121010 + ,315009 + ,1 + ,115 + ,345 + ,10 + ,50 + ,60945 + ,121890 + ,318903 + ,1 + ,5 + ,15 + ,1 + ,5 + ,60720 + ,121440 + ,314887 + ,1 + ,897 + ,2691 + ,4 + ,20 + ,60720 + ,121440 + ,314913 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,389 + ,1167 + ,8 + ,40 + ,58990 + ,117980 + ,325506 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,1002 + ,3006 + ,11 + ,55 + ,56750 + ,113500 + ,298568 + ,1 + ,36 + ,108 + ,4 + ,20 + ,60894 + ,121788 + ,315834 + ,1 + ,460 + ,1380 + ,15 + ,75 + ,63346 + ,126692 + ,329784 + ,1 + ,309 + ,927 + ,9 + ,45 + ,56535 + ,113070 + ,312878 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,9 + ,27 + ,7 + ,35 + ,60835 + ,121670 + ,314987 + ,1 + ,271 + ,813 + ,2 + ,10 + ,60720 + ,121440 + ,325249 + ,1 + ,14 + ,42 + ,0 + ,0 + ,61016 + ,122032 + ,315877 + ,1 + ,520 + ,1560 + ,7 + ,35 + ,58650 + ,117300 + ,291650 + ,1 + ,1766 + ,5298 + ,46 + ,230 + ,60438 + ,120876 + ,305959 + ,1 + ,0 + ,0 + ,5 + ,25 + ,60720 + ,121440 + ,315380 + ,1 + ,458 + ,1374 + ,7 + ,35 + ,58625 + ,117250 + ,297765 + ,1 + ,20 + ,60 + ,2 + ,10 + ,60938 + ,121876 + ,315245 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,98 + ,294 + ,2 + ,10 + ,61490 + ,122980 + ,315236 + ,1 + ,405 + ,1215 + ,5 + ,25 + ,60845 + ,121690 + ,336425 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,483 + ,1449 + ,7 + ,35 + ,60830 + ,121660 + ,306268 + ,1 + ,454 + ,1362 + ,24 + ,120 + ,63261 + ,126522 + ,302187 + ,1 + ,47 + ,141 + ,1 + ,5 + ,60720 + ,121440 + ,314882 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,757 + ,2271 + ,18 + ,90 + ,45689 + ,91378 + ,382712 + ,1 + ,4655 + ,13965 + ,55 + ,275 + ,60720 + ,121440 + ,341570 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,36 + ,108 + ,3 + ,15 + ,61564 + ,123128 + ,312412 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,203 + ,609 + ,9 + ,45 + ,61938 + ,123876 + ,309596 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,126 + ,378 + ,8 + ,40 + ,60951 + ,121902 + ,315547 + ,1 + ,400 + ,1200 + ,113 + ,565 + ,60720 + ,121440 + ,313267 + ,0 + ,71 + ,0 + ,0 + ,0 + ,60745 + ,0 + ,316176 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,972 + ,2916 + ,19 + ,95 + ,71642 + ,143284 + ,359335 + ,0 + ,531 + ,0 + ,11 + ,0 + ,71641 + ,0 + ,330068 + ,1 + ,2461 + ,7383 + ,25 + ,125 + ,55792 + ,111584 + ,314289 + ,0 + ,378 + ,0 + ,16 + ,0 + ,71873 + ,0 + ,297413 + ,0 + ,23 + ,0 + ,5 + ,0 + ,62555 + ,0 + ,314806 + ,0 + ,638 + ,0 + ,11 + ,0 + ,60370 + ,0 + ,333210 + ,0 + ,2300 + ,0 + ,23 + ,0 + ,64873 + ,0 + ,352108 + ,1 + ,149 + ,447 + ,6 + ,30 + ,62041 + ,124082 + ,313332 + ,1 + ,226 + ,678 + ,5 + ,25 + ,65745 + ,131490 + ,291787 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,275 + ,825 + ,7 + ,35 + ,59500 + ,119000 + ,318745 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,141 + ,423 + ,7 + ,35 + ,61630 + ,123260 + ,315366 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,28 + ,84 + ,3 + ,15 + ,60890 + ,121780 + ,315688 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,4980 + ,0 + ,89 + ,0 + ,113521 + ,0 + ,409642 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,472 + ,0 + ,19 + ,0 + ,80045 + ,0 + ,269587 + ,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 + ,203 + ,0 + ,12 + ,0 + ,50804 + ,0 + ,300962 + ,0 + ,496 + ,0 + ,12 + ,0 + ,87390 + ,0 + ,325479 + ,0 + ,10 + ,0 + ,5 + ,0 + ,61656 + ,0 + ,316155 + ,0 + ,63 + ,0 + ,2 + ,0 + ,65688 + ,0 + ,318574 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,1136 + ,0 + ,26 + ,0 + ,48522 + ,0 + ,343613 + ,0 + ,265 + ,0 + ,3 + ,0 + ,60720 + ,0 + ,306948 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,267 + ,801 + ,11 + ,55 + ,57640 + ,115280 + ,330059 + ,1 + ,474 + ,1422 + ,10 + ,50 + ,61977 + ,123954 + ,288985 + ,1 + ,534 + ,1602 + ,5 + ,25 + ,62620 + ,125240 + ,304485 + ,1 + ,0 + ,0 + ,2 + ,10 + ,60720 + ,121440 + ,315380 + ,1 + ,15 + ,45 + ,6 + ,30 + ,60831 + ,121662 + ,315688 + ,1 + ,397 + ,1191 + ,7 + ,35 + ,60646 + ,121292 + ,317736 + ,1 + ,0 + ,0 + ,2 + ,10 + ,60720 + ,121440 + ,315380 + ,1 + ,1866 + ,5598 + ,28 + ,140 + ,56225 + ,112450 + ,322331 + ,1 + ,288 + ,864 + ,3 + ,15 + ,60510 + ,121020 + ,296656 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,3 + ,9 + ,1 + ,5 + ,60698 + ,121396 + ,315354 + ,1 + ,468 + ,1404 + ,20 + ,100 + ,60720 + ,121440 + ,312161 + ,1 + ,20 + ,60 + ,1 + ,5 + ,60805 + ,121610 + ,315576 + ,1 + ,278 + ,834 + ,22 + ,110 + ,61404 + ,122808 + ,314922 + ,1 + ,61 + ,183 + ,9 + ,45 + ,60720 + ,121440 + ,314551 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,192 + ,576 + ,2 + ,10 + ,65276 + ,130552 + ,312339 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,317 + ,951 + ,7 + ,35 + ,63915 + ,127830 + ,298700 + ,1 + ,738 + ,2214 + ,9 + ,45 + ,60720 + ,121440 + ,321376 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,368 + ,1104 + ,13 + ,65 + ,61686 + ,123372 + ,303230 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,2 + ,6 + ,0 + ,0 + ,60743 + ,121486 + ,315487 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,53 + ,159 + ,6 + ,30 + ,60349 + ,120698 + ,315793 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,94 + ,282 + ,3 + ,15 + ,61360 + ,122720 + ,312887 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,24 + ,72 + ,7 + ,35 + ,59818 + ,119636 + ,315637 + ,1 + ,2332 + ,6996 + ,2 + ,10 + ,72680 + ,145360 + ,324385 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,131 + ,393 + ,15 + ,75 + ,61808 + ,123616 + ,308989 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,206 + ,618 + ,9 + ,45 + ,53110 + ,106220 + ,296702 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,167 + ,501 + ,1 + ,5 + ,64245 + ,128490 + ,307322 + ,1 + ,622 + ,1866 + ,38 + ,190 + ,73007 + ,146014 + ,304376 + ,1 + ,2328 + ,6984 + ,57 + ,285 + ,82732 + ,165464 + ,253588 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,365 + ,1095 + ,7 + ,35 + ,54820 + ,109640 + ,309560 + ,1 + ,364 + ,1092 + ,26 + ,130 + ,47705 + ,95410 + ,298466 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,226 + ,678 + ,13 + ,65 + ,72835 + ,145670 + ,343929 + ,1 + ,307 + ,921 + ,10 + ,50 + ,58856 + ,117712 + ,331955 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,188 + ,564 + ,9 + ,45 + ,77655 + ,155310 + ,381180 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,138 + ,414 + ,26 + ,130 + ,69817 + ,139634 + ,331420 + ,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 + ,125 + ,375 + ,19 + ,95 + ,60798 + ,121596 + ,310201 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,282 + ,846 + ,12 + ,60 + ,62452 + ,124904 + ,320016 + ,1 + ,335 + ,1005 + ,23 + ,115 + ,64175 + ,128350 + ,320398 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,0 + ,1324 + ,0 + ,29 + ,0 + ,67440 + ,0 + ,291841 + ,1 + ,176 + ,528 + ,8 + ,40 + ,68136 + ,136272 + ,310670 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,249 + ,747 + ,26 + ,130 + ,56726 + ,113452 + ,313491 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,333 + ,999 + ,9 + ,45 + ,70811 + ,141622 + ,331323 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,0 + ,601 + ,0 + ,5 + ,0 + ,60720 + ,0 + ,319210 + ,1 + ,30 + ,90 + ,3 + ,15 + ,62045 + ,124090 + ,318098 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,249 + ,747 + ,13 + ,65 + ,54323 + ,108646 + ,292754 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,165 + ,495 + ,12 + ,60 + ,62841 + ,125682 + ,325176 + ,1 + ,453 + ,1359 + ,19 + ,95 + ,81125 + ,162250 + ,365959 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,53 + ,159 + ,10 + ,50 + ,59506 + ,119012 + ,302409 + ,0 + ,382 + ,0 + ,9 + ,0 + ,59365 + ,0 + ,340968 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,9 + ,45 + ,60720 + ,121440 + ,315380 + ,0 + ,30 + ,0 + ,4 + ,0 + ,60798 + ,0 + ,313164 + ,1 + ,290 + ,870 + ,1 + ,5 + ,58790 + ,117580 + ,301164 + ,1 + ,0 + ,0 + ,1 + ,5 + ,60720 + ,121440 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,366 + ,1098 + ,14 + ,70 + ,61808 + ,123616 + ,344425 + ,1 + ,2 + ,6 + ,12 + ,60 + ,60735 + ,121470 + ,315394 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,209 + ,0 + ,19 + ,0 + ,64016 + ,0 + ,316647 + ,1 + ,384 + ,1152 + ,17 + ,85 + ,54683 + ,109366 + ,309836 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,365 + ,1095 + ,32 + ,160 + ,87192 + ,174384 + ,346611 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,49 + ,0 + ,14 + ,0 + ,64107 + ,0 + ,322031 + ,1 + ,3 + ,9 + ,8 + ,40 + ,60761 + ,121522 + ,315656 + ,1 + ,133 + ,399 + ,4 + ,20 + ,65990 + ,131980 + ,339445 + ,1 + ,32 + ,96 + ,0 + ,0 + ,59988 + ,119976 + ,314964 + ,1 + ,368 + ,1104 + ,20 + ,100 + ,61167 + ,122334 + ,297141 + ,1 + ,1 + ,3 + ,5 + ,25 + ,60719 + ,121438 + ,315372 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,22 + ,66 + ,1 + ,5 + ,60722 + ,121444 + ,312502 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,96 + ,288 + ,4 + ,20 + ,60379 + ,120758 + ,313729 + ,1 + ,1 + ,3 + ,1 + ,5 + ,60727 + ,121454 + ,315388 + ,1 + ,314 + ,942 + ,4 + ,20 + ,60720 + ,121440 + ,315371 + ,1 + ,844 + ,2532 + ,20 + ,100 + ,60925 + ,121850 + ,296139 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,26 + ,78 + ,1 + ,5 + ,60896 + ,121792 + ,313880 + ,1 + ,125 + ,375 + ,10 + ,50 + ,59734 + ,119468 + ,317698 + ,1 + ,304 + ,912 + ,12 + ,60 + ,62969 + ,125938 + ,295580 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,621 + ,1863 + ,13 + ,65 + ,60720 + ,121440 + ,308256 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,119 + ,357 + ,3 + ,15 + ,59118 + ,118236 + ,303677 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,1595 + ,4785 + ,10 + ,50 + ,60720 + ,121440 + ,319369 + ,1 + ,312 + ,936 + ,3 + ,15 + ,58598 + ,117196 + ,318690 + ,1 + ,60 + ,180 + ,7 + ,35 + ,61124 + ,122248 + ,314049 + ,1 + ,587 + ,1761 + ,10 + ,50 + ,59595 + ,119190 + ,325699 + ,1 + ,135 + ,405 + ,1 + ,5 + ,62065 + ,124130 + ,314210 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,514 + ,1542 + ,15 + ,75 + ,78780 + ,157560 + ,322378 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,1 + ,3 + ,4 + ,20 + ,60722 + ,121444 + ,315398 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,0 + ,1763 + ,0 + ,28 + ,0 + ,61600 + ,0 + ,308336 + ,1 + ,180 + ,540 + ,9 + ,45 + ,59635 + ,119270 + ,316386 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,218 + ,654 + ,7 + ,35 + ,60720 + ,121440 + ,315553 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,448 + ,1344 + ,7 + ,35 + ,59781 + ,119562 + ,323361 + ,1 + ,227 + ,681 + ,7 + ,35 + ,76644 + ,153288 + ,336639 + ,1 + ,174 + ,522 + ,3 + ,15 + ,64820 + ,129640 + ,307424 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,121 + ,363 + ,11 + ,55 + ,56178 + ,112356 + ,295370 + ,1 + ,607 + ,1821 + ,7 + ,35 + ,60436 + ,120872 + ,322340 + ,1 + ,2212 + ,6636 + ,10 + ,50 + ,60720 + ,121440 + ,319864 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,530 + ,1590 + ,18 + ,90 + ,73433 + ,146866 + ,317291 + ,1 + ,571 + ,1713 + ,14 + ,70 + ,41477 + ,82954 + ,280398 + ,1 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,121440 + ,315380 + ,1 + ,78 + ,234 + ,12 + ,60 + ,62700 + ,125400 + ,317330 + ,1 + ,2489 + ,7467 + ,29 + ,145 + ,67804 + ,135608 + ,238125 + ,1 + ,131 + ,393 + ,3 + ,15 + ,59661 + ,119322 + ,327071 + ,1 + ,923 + ,2769 + ,6 + ,30 + ,58620 + ,117240 + ,309038 + ,1 + ,72 + ,216 + ,3 + ,15 + ,60398 + ,120796 + ,314210 + ,1 + ,572 + ,1716 + ,8 + ,40 + ,58580 + ,117160 + ,307930 + ,1 + ,397 + ,1191 + ,10 + ,50 + ,62710 + ,125420 + ,322327 + ,1 + ,450 + ,1350 + ,6 + ,30 + ,59325 + ,118650 + ,292136 + ,1 + ,622 + ,1866 + ,8 + ,40 + ,60950 + ,121900 + ,263276 + ,1 + ,694 + ,2082 + ,6 + ,30 + ,68060 + ,136120 + ,367655 + ,0 + ,3425 + ,0 + ,9 + ,0 + ,83620 + ,0 + ,283910 + ,1 + ,562 + ,1686 + ,8 + ,40 + ,58456 + ,116912 + ,283587 + ,1 + ,4917 + ,14751 + ,26 + ,130 + ,52811 + ,105622 + ,243650 + ,0 + ,1442 + ,0 + ,239 + ,0 + ,121173 + ,0 + ,438493 + ,1 + ,529 + ,1587 + ,7 + ,35 + ,63870 + ,127740 + ,296261 + ,0 + ,2126 + ,0 + ,41 + ,0 + ,21001 + ,0 + ,230621 + ,1 + ,1061 + ,3183 + ,3 + ,15 + ,70415 + ,140830 + ,304252 + ,1 + ,776 + ,2328 + ,8 + ,40 + ,64230 + ,128460 + ,333505 + ,1 + ,611 + ,1833 + ,6 + ,30 + ,59190 + ,118380 + ,296919 + ,0 + ,1526 + ,0 + ,21 + ,0 + ,69351 + ,0 + ,278990 + ,1 + ,592 + ,1776 + ,7 + ,35 + ,64270 + ,128540 + ,276898 + ,1 + ,1182 + ,3546 + ,11 + ,55 + ,70694 + ,141388 + ,327007 + ,1 + ,621 + ,1863 + ,11 + ,55 + ,68005 + ,136010 + ,317046 + ,1 + ,989 + ,2967 + ,12 + ,60 + ,58930 + ,117860 + ,304555 + ,1 + ,438 + ,1314 + ,9 + ,45 + ,58320 + ,116640 + ,298096 + ,1 + ,726 + ,2178 + ,3 + ,15 + ,69980 + ,139960 + ,231861 + ,1 + ,1303 + ,3909 + ,57 + ,285 + ,69863 + ,139726 + ,309422 + ,0 + ,7419 + ,0 + ,21 + ,0 + ,63255 + ,0 + ,286963 + ,0 + ,1164 + ,0 + ,15 + ,0 + ,57320 + ,0 + ,269753 + ,0 + ,3310 + ,0 + ,32 + ,0 + ,75230 + ,0 + ,448243 + ,1 + ,1920 + ,5760 + ,11 + ,55 + ,79420 + ,158840 + ,165404 + ,1 + ,965 + ,2895 + ,2 + ,10 + ,73490 + ,146980 + ,204325 + ,1 + ,3256 + ,9768 + ,23 + ,115 + ,35250 + ,70500 + ,407159 + ,0 + ,1135 + ,0 + ,20 + ,0 + ,62285 + ,0 + ,290476 + ,1 + ,1270 + ,3810 + ,24 + ,120 + ,69206 + ,138412 + ,275311 + ,1 + ,661 + ,1983 + ,1 + ,5 + ,65920 + ,131840 + ,246541 + ,1 + ,1013 + ,3039 + ,1 + ,5 + ,69770 + ,139540 + ,253468 + ,1 + ,2844 + ,8532 + ,74 + ,370 + ,72683 + ,145366 + ,240897 + ,0 + ,11528 + ,0 + ,68 + ,0 + ,-14545 + ,0 + ,-83265 + ,1 + ,6526 + ,19578 + ,20 + ,100 + ,55830 + ,111660 + ,-42143 + ,1 + ,2264 + ,6792 + ,20 + ,100 + ,55174 + ,110348 + ,272713 + ,0 + ,5109 + ,0 + ,82 + ,0 + ,67038 + ,0 + ,215362 + ,1 + ,3999 + ,11997 + ,21 + ,105 + ,51252 + ,102504 + ,42754 + ,1 + ,35624 + ,106872 + ,244 + ,1220 + ,157278 + ,314556 + ,306275 + ,1 + ,9252 + ,27756 + ,32 + ,160 + ,79510 + ,159020 + ,253537 + ,1 + ,15236 + ,45708 + ,86 + ,430 + ,77440 + ,154880 + ,372631 + ,1 + ,18073 + ,54219 + ,69 + ,345 + ,27284 + ,54568 + ,-7170) + ,dim=c(8 + ,431) + ,dimnames=list(c('Group' + ,'Costs' + ,'Costs_p' + ,'Orders' + ,'Orders_p' + ,'Dividends' + ,'Dividends_p' + ,'Wealth ') + ,1:431)) > y <- array(NA,dim=c(8,431),dimnames=list(c('Group','Costs','Costs_p','Orders','Orders_p','Dividends','Dividends_p','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\r Group Costs Costs_p Orders Orders_p Dividends Dividends_p 1 6282154 0 162556 0 807 0 213118 0 2 4321023 0 29790 0 444 0 81767 0 3 4111912 0 87550 0 412 0 153198 0 4 223193 1 84738 254214 428 2140 -26007 -52014 5 1491348 0 54660 0 315 0 126942 0 6 1629616 0 42634 0 168 0 157214 0 7 1398893 1 40949 122847 263 1315 129352 258704 8 1926517 0 45187 0 267 0 234817 0 9 983660 0 37704 0 228 0 60448 0 10 1443586 0 16275 0 129 0 47818 0 11 1073089 1 25830 77490 104 520 245546 491092 12 984885 1 12679 38037 122 610 48020 96040 13 1405225 0 18014 0 393 0 -1710 0 14 227132 1 43556 130668 190 950 32648 65296 15 929118 0 24811 0 280 0 95350 0 16 1071292 1 6575 19725 63 315 151352 302704 17 638830 1 7123 21369 102 510 288170 576340 18 856956 0 21950 0 265 0 114337 0 19 992426 0 37597 0 234 0 37884 0 20 444477 1 17821 53463 277 1385 122844 245688 21 857217 0 12988 0 73 0 82340 0 22 711969 0 22330 0 67 0 79801 0 23 702380 1 13326 39978 103 515 165548 331096 24 358589 1 16189 48567 290 1450 116384 232768 25 297978 1 7146 21438 83 415 134028 268056 26 585715 1 15824 47472 56 280 63838 127676 27 657954 0 27664 0 236 0 74996 0 28 209458 1 11920 35760 73 365 31080 62160 29 786690 1 8568 25704 34 170 32168 64336 30 439798 1 14416 43248 139 695 49857 99714 31 688779 0 3369 0 26 0 87161 0 32 574339 0 11819 0 70 0 106113 0 33 741409 0 6984 0 40 0 80570 0 34 597793 0 4519 0 42 0 102129 0 35 644190 1 2220 6660 12 60 301670 603340 36 377934 1 18562 55686 211 1055 102313 204626 37 640273 1 10327 30981 74 370 88577 177154 38 697458 0 5336 0 80 0 112477 0 39 550608 0 2365 0 83 0 191778 0 40 207393 1 4069 12207 131 655 79804 159608 41 301607 1 8636 25908 203 1015 128294 256588 42 345783 1 13718 41154 56 280 96448 192896 43 501749 1 4525 13575 89 445 93811 187622 44 379983 1 6869 20607 88 440 117520 235040 45 387475 1 4628 13884 39 195 69159 138318 46 377305 0 3689 0 25 0 101792 0 47 370837 0 4891 0 49 0 210568 0 48 430866 0 7489 0 149 0 136996 0 49 469107 1 4901 14703 58 290 121920 243840 50 194493 1 2284 6852 41 205 76403 152806 51 530670 0 3160 0 90 0 108094 0 52 518365 0 4150 0 136 0 134759 0 53 491303 0 7285 0 97 0 188873 0 54 527021 0 1134 0 63 0 146216 0 55 233773 0 4658 0 114 0 156608 0 56 405972 1 2384 7152 77 385 61348 122696 57 652925 1 3748 11244 6 30 50350 100700 58 446211 1 5371 16113 47 235 87720 175440 59 341340 1 1285 3855 51 255 99489 198978 60 387699 0 9327 0 85 0 87419 0 61 493408 0 5565 0 43 0 94355 0 62 146494 1 1528 4584 32 160 60326 120652 63 414462 0 3122 0 25 0 94670 0 64 364304 0 7561 0 77 0 82425 0 65 355178 1 2675 8025 54 270 59017 118034 66 357760 1 13253 39759 251 1255 90829 181658 67 261216 1 880 2640 15 75 80791 161582 68 397144 0 2053 0 44 0 100423 0 69 374943 1 1424 4272 73 365 131116 262232 70 424898 0 4036 0 85 0 100269 0 71 202055 0 3045 0 49 0 27330 0 72 378525 1 5119 15357 38 190 39039 78078 73 310768 1 1431 4293 35 175 106885 213770 74 325738 1 554 1662 9 45 79285 158570 75 394510 1 1975 5925 34 170 118881 237762 76 247060 0 1765 0 20 0 77623 0 77 368078 1 1012 3036 29 145 114768 229536 78 236761 1 810 2430 11 55 74015 148030 79 312378 1 1280 3840 52 260 69465 138930 80 339836 0 666 0 13 0 117869 0 81 347385 1 1380 4140 29 145 60982 121964 82 426280 0 4677 0 66 0 90131 0 83 352850 1 876 2628 33 165 138971 277942 84 301881 1 814 2442 15 75 39625 79250 85 377516 1 514 1542 15 75 102725 205450 86 357312 0 5692 0 68 0 64239 0 87 458343 1 3642 10926 100 500 90262 180524 88 354228 1 540 1620 13 65 103960 207920 89 308636 1 2099 6297 45 225 106611 213222 90 386212 1 567 1701 14 70 103345 206690 91 393343 1 2001 6003 36 180 95551 191102 92 378509 0 2949 0 40 0 82903 0 93 452469 1 2253 6759 68 340 63593 127186 94 364839 0 6533 0 29 0 126910 0 95 358649 1 1889 5667 43 215 37527 75054 96 376641 0 3055 0 30 0 60247 0 97 429112 1 272 816 9 45 112995 225990 98 330546 0 1414 0 22 0 70184 0 99 403560 1 2564 7692 19 95 130140 260280 100 317892 0 1383 0 9 0 73221 0 101 307528 0 1261 0 31 0 76114 0 102 235133 1 975 2925 19 95 90534 181068 103 299243 1 3366 10098 55 275 108479 216958 104 314073 1 576 1728 8 40 113761 227522 105 368186 1 1686 5058 28 140 68696 137392 106 269661 1 746 2238 29 145 71561 143122 107 125390 0 3192 0 48 0 59831 0 108 510834 0 2045 0 16 0 97890 0 109 321896 1 5702 17106 47 235 101481 202962 110 249898 0 1932 0 20 0 72954 0 111 408881 1 936 2808 22 110 67939 135878 112 158492 0 3437 0 33 0 48022 0 113 292154 1 5131 15393 44 220 86111 172222 114 289513 0 2397 0 13 0 74020 0 115 378049 0 1389 0 6 0 57530 0 116 343466 1 1503 4509 35 175 56364 112728 117 332743 1 402 1206 8 40 84990 169980 118 442882 1 2239 6717 17 85 88590 177180 119 214215 0 2234 0 11 0 77200 0 120 315688 1 837 2511 21 105 61262 122524 121 375195 1 10579 31737 92 460 110309 220618 122 334280 1 875 2625 12 60 67000 134000 123 355864 1 1585 4755 112 560 93099 186198 124 480382 0 1659 0 25 0 107577 0 125 353058 0 2647 0 17 0 62920 0 126 217193 0 3294 0 23 0 75832 0 127 315380 1 0 0 0 0 60720 121440 128 314533 1 94 282 10 50 60793 121586 129 318056 1 422 1266 23 115 57935 115870 130 315380 1 0 0 0 0 60720 121440 131 314353 1 34 102 7 35 60630 121260 132 369448 1 1558 4674 25 125 55637 111274 133 315380 1 0 0 0 0 60720 121440 134 312846 1 43 129 20 100 60887 121774 135 312075 1 645 1935 4 20 60720 121440 136 315009 1 316 948 4 20 60505 121010 137 318903 1 115 345 10 50 60945 121890 138 314887 1 5 15 1 5 60720 121440 139 314913 1 897 2691 4 20 60720 121440 140 315380 1 0 0 0 0 60720 121440 141 325506 1 389 1167 8 40 58990 117980 142 315380 1 0 0 0 0 60720 121440 143 298568 1 1002 3006 11 55 56750 113500 144 315834 1 36 108 4 20 60894 121788 145 329784 1 460 1380 15 75 63346 126692 146 312878 1 309 927 9 45 56535 113070 147 315380 1 0 0 0 0 60720 121440 148 314987 1 9 27 7 35 60835 121670 149 325249 1 271 813 2 10 60720 121440 150 315877 1 14 42 0 0 61016 122032 151 291650 1 520 1560 7 35 58650 117300 152 305959 1 1766 5298 46 230 60438 120876 153 315380 1 0 0 5 25 60720 121440 154 297765 1 458 1374 7 35 58625 117250 155 315245 1 20 60 2 10 60938 121876 156 315380 1 0 0 0 0 60720 121440 157 315380 1 0 0 0 0 60720 121440 158 315236 1 98 294 2 10 61490 122980 159 336425 1 405 1215 5 25 60845 121690 160 315380 1 0 0 0 0 60720 121440 161 315380 1 0 0 0 0 60720 121440 162 315380 1 0 0 0 0 60720 121440 163 315380 1 0 0 0 0 60720 121440 164 306268 1 483 1449 7 35 60830 121660 165 302187 1 454 1362 24 120 63261 126522 166 314882 1 47 141 1 5 60720 121440 167 315380 1 0 0 0 0 60720 121440 168 382712 1 757 2271 18 90 45689 91378 169 341570 1 4655 13965 55 275 60720 121440 170 315380 1 0 0 0 0 60720 121440 171 315380 1 0 0 0 0 60720 121440 172 312412 1 36 108 3 15 61564 123128 173 315380 1 0 0 0 0 60720 121440 174 309596 1 203 609 9 45 61938 123876 175 315380 1 0 0 0 0 60720 121440 176 315547 1 126 378 8 40 60951 121902 177 313267 1 400 1200 113 565 60720 121440 178 316176 0 71 0 0 0 60745 0 179 315380 1 0 0 0 0 60720 121440 180 315380 1 0 0 0 0 60720 121440 181 359335 1 972 2916 19 95 71642 143284 182 330068 0 531 0 11 0 71641 0 183 314289 1 2461 7383 25 125 55792 111584 184 297413 0 378 0 16 0 71873 0 185 314806 0 23 0 5 0 62555 0 186 333210 0 638 0 11 0 60370 0 187 352108 0 2300 0 23 0 64873 0 188 313332 1 149 447 6 30 62041 124082 189 291787 1 226 678 5 25 65745 131490 190 315380 1 0 0 0 0 60720 121440 191 318745 1 275 825 7 35 59500 119000 192 315380 1 0 0 0 0 60720 121440 193 315366 1 141 423 7 35 61630 123260 194 315380 1 0 0 0 0 60720 121440 195 315688 1 28 84 3 15 60890 121780 196 315380 0 0 0 0 0 60720 0 197 409642 0 4980 0 89 0 113521 0 198 315380 0 0 0 0 0 60720 0 199 315380 0 0 0 0 0 60720 0 200 269587 0 472 0 19 0 80045 0 201 315380 0 0 0 0 0 60720 0 202 315380 0 0 0 0 0 60720 0 203 315380 0 0 0 0 0 60720 0 204 300962 0 203 0 12 0 50804 0 205 325479 0 496 0 12 0 87390 0 206 316155 0 10 0 5 0 61656 0 207 318574 0 63 0 2 0 65688 0 208 315380 0 0 0 0 0 60720 0 209 343613 0 1136 0 26 0 48522 0 210 306948 0 265 0 3 0 60720 0 211 315380 1 0 0 0 0 60720 121440 212 315380 1 0 0 0 0 60720 121440 213 330059 1 267 801 11 55 57640 115280 214 288985 1 474 1422 10 50 61977 123954 215 304485 1 534 1602 5 25 62620 125240 216 315380 1 0 0 2 10 60720 121440 217 315688 1 15 45 6 30 60831 121662 218 317736 1 397 1191 7 35 60646 121292 219 315380 1 0 0 2 10 60720 121440 220 322331 1 1866 5598 28 140 56225 112450 221 296656 1 288 864 3 15 60510 121020 222 315380 1 0 0 0 0 60720 121440 223 315354 1 3 9 1 5 60698 121396 224 312161 1 468 1404 20 100 60720 121440 225 315576 1 20 60 1 5 60805 121610 226 314922 1 278 834 22 110 61404 122808 227 314551 1 61 183 9 45 60720 121440 228 315380 1 0 0 0 0 60720 121440 229 312339 1 192 576 2 10 65276 130552 230 315380 1 0 0 0 0 60720 121440 231 298700 1 317 951 7 35 63915 127830 232 321376 1 738 2214 9 45 60720 121440 233 315380 1 0 0 0 0 60720 121440 234 303230 1 368 1104 13 65 61686 123372 235 315380 1 0 0 0 0 60720 121440 236 315487 1 2 6 0 0 60743 121486 237 315380 1 0 0 0 0 60720 121440 238 315793 1 53 159 6 30 60349 120698 239 315380 1 0 0 0 0 60720 121440 240 315380 1 0 0 0 0 60720 121440 241 315380 1 0 0 0 0 60720 121440 242 312887 1 94 282 3 15 61360 122720 243 315380 1 0 0 0 0 60720 121440 244 315637 1 24 72 7 35 59818 119636 245 324385 1 2332 6996 2 10 72680 145360 246 315380 0 0 0 0 0 60720 0 247 315380 0 0 0 0 0 60720 0 248 308989 1 131 393 15 75 61808 123616 249 315380 0 0 0 0 0 60720 0 250 315380 0 0 0 0 0 60720 0 251 296702 1 206 618 9 45 53110 106220 252 315380 0 0 0 0 0 60720 0 253 307322 1 167 501 1 5 64245 128490 254 304376 1 622 1866 38 190 73007 146014 255 253588 1 2328 6984 57 285 82732 165464 256 315380 1 0 0 0 0 60720 121440 257 309560 1 365 1095 7 35 54820 109640 258 298466 1 364 1092 26 130 47705 95410 259 315380 0 0 0 0 0 60720 0 260 315380 1 0 0 0 0 60720 121440 261 315380 0 0 0 0 0 60720 0 262 315380 0 0 0 0 0 60720 0 263 343929 1 226 678 13 65 72835 145670 264 331955 1 307 921 10 50 58856 117712 265 315380 0 0 0 0 0 60720 0 266 315380 1 0 0 0 0 60720 121440 267 315380 0 0 0 0 0 60720 0 268 381180 1 188 564 9 45 77655 155310 269 315380 0 0 0 0 0 60720 0 270 331420 1 138 414 26 130 69817 139634 271 315380 0 0 0 0 0 60720 0 272 315380 0 0 0 0 0 60720 0 273 315380 0 0 0 0 0 60720 0 274 310201 1 125 375 19 95 60798 121596 275 315380 1 0 0 0 0 60720 121440 276 320016 1 282 846 12 60 62452 124904 277 320398 1 335 1005 23 115 64175 128350 278 315380 1 0 0 0 0 60720 121440 279 291841 0 1324 0 29 0 67440 0 280 310670 1 176 528 8 40 68136 136272 281 315380 0 0 0 0 0 60720 0 282 315380 1 0 0 0 0 60720 121440 283 313491 1 249 747 26 130 56726 113452 284 315380 1 0 0 0 0 60720 121440 285 331323 1 333 999 9 45 70811 141622 286 315380 1 0 0 0 0 60720 121440 287 319210 0 601 0 5 0 60720 0 288 318098 1 30 90 3 15 62045 124090 289 315380 0 0 0 0 0 60720 0 290 292754 1 249 747 13 65 54323 108646 291 315380 0 0 0 0 0 60720 0 292 325176 1 165 495 12 60 62841 125682 293 365959 1 453 1359 19 95 81125 162250 294 315380 1 0 0 0 0 60720 121440 295 302409 1 53 159 10 50 59506 119012 296 340968 0 382 0 9 0 59365 0 297 315380 1 0 0 0 0 60720 121440 298 315380 0 0 0 0 0 60720 0 299 315380 1 0 0 0 0 60720 121440 300 315380 1 0 0 9 45 60720 121440 301 313164 0 30 0 4 0 60798 0 302 301164 1 290 870 1 5 58790 117580 303 315380 1 0 0 1 5 60720 121440 304 315380 0 0 0 0 0 60720 0 305 344425 1 366 1098 14 70 61808 123616 306 315394 1 2 6 12 60 60735 121470 307 315380 0 0 0 0 0 60720 0 308 316647 0 209 0 19 0 64016 0 309 309836 1 384 1152 17 85 54683 109366 310 315380 0 0 0 0 0 60720 0 311 315380 0 0 0 0 0 60720 0 312 346611 1 365 1095 32 160 87192 174384 313 315380 0 0 0 0 0 60720 0 314 322031 0 49 0 14 0 64107 0 315 315656 1 3 9 8 40 60761 121522 316 339445 1 133 399 4 20 65990 131980 317 314964 1 32 96 0 0 59988 119976 318 297141 1 368 1104 20 100 61167 122334 319 315372 1 1 3 5 25 60719 121438 320 315380 1 0 0 0 0 60720 121440 321 315380 1 0 0 0 0 60720 121440 322 315380 1 0 0 0 0 60720 121440 323 315380 1 0 0 0 0 60720 121440 324 315380 0 0 0 0 0 60720 0 325 315380 1 0 0 0 0 60720 121440 326 312502 1 22 66 1 5 60722 121444 327 315380 0 0 0 0 0 60720 0 328 315380 0 0 0 0 0 60720 0 329 315380 1 0 0 0 0 60720 121440 330 315380 1 0 0 0 0 60720 121440 331 315380 1 0 0 0 0 60720 121440 332 315380 1 0 0 0 0 60720 121440 333 315380 1 0 0 0 0 60720 121440 334 313729 1 96 288 4 20 60379 120758 335 315388 1 1 3 1 5 60727 121454 336 315371 1 314 942 4 20 60720 121440 337 296139 1 844 2532 20 100 60925 121850 338 315380 1 0 0 0 0 60720 121440 339 313880 1 26 78 1 5 60896 121792 340 317698 1 125 375 10 50 59734 119468 341 295580 1 304 912 12 60 62969 125938 342 315380 1 0 0 0 0 60720 121440 343 315380 1 0 0 0 0 60720 121440 344 315380 1 0 0 0 0 60720 121440 345 308256 1 621 1863 13 65 60720 121440 346 315380 1 0 0 0 0 60720 121440 347 303677 1 119 357 3 15 59118 118236 348 315380 1 0 0 0 0 60720 121440 349 315380 1 0 0 0 0 60720 121440 350 319369 1 1595 4785 10 50 60720 121440 351 318690 1 312 936 3 15 58598 117196 352 314049 1 60 180 7 35 61124 122248 353 325699 1 587 1761 10 50 59595 119190 354 314210 1 135 405 1 5 62065 124130 355 315380 1 0 0 0 0 60720 121440 356 315380 1 0 0 0 0 60720 121440 357 322378 1 514 1542 15 75 78780 157560 358 315380 1 0 0 0 0 60720 121440 359 315380 1 0 0 0 0 60720 121440 360 315380 1 0 0 0 0 60720 121440 361 315398 1 1 3 4 20 60722 121444 362 315380 1 0 0 0 0 60720 121440 363 315380 1 0 0 0 0 60720 121440 364 308336 0 1763 0 28 0 61600 0 365 316386 1 180 540 9 45 59635 119270 366 315380 1 0 0 0 0 60720 121440 367 315380 1 0 0 0 0 60720 121440 368 315380 1 0 0 0 0 60720 121440 369 315380 1 0 0 0 0 60720 121440 370 315553 1 218 654 7 35 60720 121440 371 315380 1 0 0 0 0 60720 121440 372 323361 1 448 1344 7 35 59781 119562 373 336639 1 227 681 7 35 76644 153288 374 307424 1 174 522 3 15 64820 129640 375 315380 1 0 0 0 0 60720 121440 376 315380 1 0 0 0 0 60720 121440 377 295370 1 121 363 11 55 56178 112356 378 322340 1 607 1821 7 35 60436 120872 379 319864 1 2212 6636 10 50 60720 121440 380 315380 1 0 0 0 0 60720 121440 381 315380 1 0 0 0 0 60720 121440 382 317291 1 530 1590 18 90 73433 146866 383 280398 1 571 1713 14 70 41477 82954 384 315380 1 0 0 0 0 60720 121440 385 317330 1 78 234 12 60 62700 125400 386 238125 1 2489 7467 29 145 67804 135608 387 327071 1 131 393 3 15 59661 119322 388 309038 1 923 2769 6 30 58620 117240 389 314210 1 72 216 3 15 60398 120796 390 307930 1 572 1716 8 40 58580 117160 391 322327 1 397 1191 10 50 62710 125420 392 292136 1 450 1350 6 30 59325 118650 393 263276 1 622 1866 8 40 60950 121900 394 367655 1 694 2082 6 30 68060 136120 395 283910 0 3425 0 9 0 83620 0 396 283587 1 562 1686 8 40 58456 116912 397 243650 1 4917 14751 26 130 52811 105622 398 438493 0 1442 0 239 0 121173 0 399 296261 1 529 1587 7 35 63870 127740 400 230621 0 2126 0 41 0 21001 0 401 304252 1 1061 3183 3 15 70415 140830 402 333505 1 776 2328 8 40 64230 128460 403 296919 1 611 1833 6 30 59190 118380 404 278990 0 1526 0 21 0 69351 0 405 276898 1 592 1776 7 35 64270 128540 406 327007 1 1182 3546 11 55 70694 141388 407 317046 1 621 1863 11 55 68005 136010 408 304555 1 989 2967 12 60 58930 117860 409 298096 1 438 1314 9 45 58320 116640 410 231861 1 726 2178 3 15 69980 139960 411 309422 1 1303 3909 57 285 69863 139726 412 286963 0 7419 0 21 0 63255 0 413 269753 0 1164 0 15 0 57320 0 414 448243 0 3310 0 32 0 75230 0 415 165404 1 1920 5760 11 55 79420 158840 416 204325 1 965 2895 2 10 73490 146980 417 407159 1 3256 9768 23 115 35250 70500 418 290476 0 1135 0 20 0 62285 0 419 275311 1 1270 3810 24 120 69206 138412 420 246541 1 661 1983 1 5 65920 131840 421 253468 1 1013 3039 1 5 69770 139540 422 240897 1 2844 8532 74 370 72683 145366 423 -83265 0 11528 0 68 0 -14545 0 424 -42143 1 6526 19578 20 100 55830 111660 425 272713 1 2264 6792 20 100 55174 110348 426 215362 0 5109 0 82 0 67038 0 427 42754 1 3999 11997 21 105 51252 102504 428 306275 1 35624 106872 244 1220 157278 314556 429 253537 1 9252 27756 32 160 79510 159020 430 372631 1 15236 45708 86 430 77440 154880 431 -7170 1 18073 54219 69 345 27284 54568 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Group Costs Costs_p Orders Orders_p 2.091e+05 -1.373e+04 2.322e+01 -6.056e+00 2.549e+03 -5.471e+02 Dividends Dividends_p 1.567e-01 8.710e-01 > (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) 2.091e+05 4.111e+04 5.087 5.49e-07 *** Group -1.373e+04 4.962e+04 -0.277 0.78218 Costs 2.322e+01 1.992e+00 11.654 < 2e-16 *** Costs_p -6.056e+00 1.160e+00 -5.219 2.83e-07 *** Orders 2.549e+03 3.336e+02 7.641 1.46e-13 *** Orders_p -5.471e+02 1.065e+02 -5.137 4.26e-07 *** Dividends 1.567e-01 5.031e-01 0.312 0.75555 Dividends_p 8.710e-01 3.181e-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/1p6jf1291284860.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/2p6jf1291284860.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/3p6jf1291284860.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/40x0h1291284860.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/50x0h1291284860.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/6bphk1291284860.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/73gg51291284860.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/83gg51291284860.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/93gg51291284860.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/10wpx81291284860.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/11z7ew1291284860.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/12d0fx1291284861.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/13rsdo1291284861.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/14cabt1291284861.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/15gbah1291284861.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/16ul881291284861.tab") + } > > try(system("convert tmp/1p6jf1291284860.ps tmp/1p6jf1291284860.png",intern=TRUE)) character(0) > try(system("convert tmp/2p6jf1291284860.ps tmp/2p6jf1291284860.png",intern=TRUE)) character(0) > try(system("convert tmp/3p6jf1291284860.ps tmp/3p6jf1291284860.png",intern=TRUE)) character(0) > try(system("convert tmp/40x0h1291284860.ps tmp/40x0h1291284860.png",intern=TRUE)) character(0) > try(system("convert tmp/50x0h1291284860.ps tmp/50x0h1291284860.png",intern=TRUE)) character(0) > try(system("convert tmp/6bphk1291284860.ps tmp/6bphk1291284860.png",intern=TRUE)) character(0) > try(system("convert tmp/73gg51291284860.ps tmp/73gg51291284860.png",intern=TRUE)) character(0) > try(system("convert tmp/83gg51291284860.ps tmp/83gg51291284860.png",intern=TRUE)) character(0) > try(system("convert tmp/93gg51291284860.ps tmp/93gg51291284860.png",intern=TRUE)) character(0) > try(system("convert tmp/10wpx81291284860.ps tmp/10wpx81291284860.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 14.077 2.290 35.644