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(6282154 + ,807 + ,1 + ,213118 + ,162556 + ,4321023 + ,444 + ,1 + ,81767 + ,29790 + ,4111912 + ,412 + ,1 + ,153198 + ,87550 + ,223193 + ,428 + ,0 + ,-26007 + ,84738 + ,1491348 + ,315 + ,1 + ,126942 + ,54660 + ,1629616 + ,168 + ,1 + ,157214 + ,42634 + ,1398893 + ,263 + ,0 + ,129352 + ,40949 + ,1926517 + ,267 + ,1 + ,234817 + ,45187 + ,983660 + ,228 + ,1 + ,60448 + ,37704 + ,1443586 + ,129 + ,1 + ,47818 + ,16275 + ,1073089 + ,104 + ,0 + ,245546 + ,25830 + ,984885 + ,122 + ,0 + ,48020 + ,12679 + ,1405225 + ,393 + ,1 + ,-1710 + ,18014 + ,227132 + ,190 + ,0 + ,32648 + ,43556 + ,929118 + ,280 + ,1 + ,95350 + ,24811 + ,1071292 + ,63 + ,0 + ,151352 + ,6575 + ,638830 + ,102 + ,0 + ,288170 + ,7123 + ,856956 + ,265 + ,1 + ,114337 + ,21950 + ,992426 + ,234 + ,1 + ,37884 + ,37597 + ,444477 + ,277 + ,0 + ,122844 + ,17821 + ,857217 + ,73 + ,1 + ,82340 + ,12988 + ,711969 + ,67 + ,1 + ,79801 + ,22330 + ,702380 + ,103 + ,0 + ,165548 + ,13326 + ,358589 + ,290 + ,0 + ,116384 + ,16189 + ,297978 + ,83 + ,0 + ,134028 + ,7146 + ,585715 + ,56 + ,0 + ,63838 + ,15824 + ,657954 + ,236 + ,1 + ,74996 + ,27664 + ,209458 + ,73 + ,0 + ,31080 + ,11920 + ,786690 + ,34 + ,0 + ,32168 + ,8568 + ,439798 + ,139 + ,0 + ,49857 + ,14416 + ,688779 + ,26 + ,1 + ,87161 + ,3369 + ,574339 + ,70 + ,1 + ,106113 + ,11819 + ,741409 + ,40 + ,1 + ,80570 + ,6984 + ,597793 + ,42 + ,1 + ,102129 + ,4519 + ,644190 + ,12 + ,0 + ,301670 + ,2220 + ,377934 + ,211 + ,0 + ,102313 + ,18562 + ,640273 + ,74 + ,0 + ,88577 + ,10327 + ,697458 + ,80 + ,1 + ,112477 + ,5336 + ,550608 + ,83 + ,1 + ,191778 + ,2365 + ,207393 + ,131 + ,0 + ,79804 + ,4069 + ,301607 + ,203 + ,0 + ,128294 + ,8636 + ,345783 + ,56 + ,0 + ,96448 + ,13718 + ,501749 + ,89 + ,0 + ,93811 + ,4525 + ,379983 + ,88 + ,0 + ,117520 + ,6869 + ,387475 + ,39 + ,0 + ,69159 + ,4628 + ,377305 + ,25 + ,1 + ,101792 + ,3689 + ,370837 + ,49 + ,1 + ,210568 + ,4891 + ,430866 + ,149 + ,1 + ,136996 + ,7489 + ,469107 + ,58 + ,0 + ,121920 + ,4901 + ,194493 + ,41 + ,0 + ,76403 + ,2284 + ,530670 + ,90 + ,1 + ,108094 + ,3160 + ,518365 + ,136 + ,1 + ,134759 + ,4150 + ,491303 + ,97 + ,1 + ,188873 + ,7285 + ,527021 + ,63 + ,1 + ,146216 + ,1134 + ,233773 + ,114 + ,1 + ,156608 + ,4658 + ,405972 + ,77 + ,0 + ,61348 + ,2384 + ,652925 + ,6 + ,0 + ,50350 + ,3748 + ,446211 + ,47 + ,0 + ,87720 + ,5371 + ,341340 + ,51 + ,0 + ,99489 + ,1285 + ,387699 + ,85 + ,1 + ,87419 + ,9327 + ,493408 + ,43 + ,1 + ,94355 + ,5565 + ,146494 + ,32 + ,0 + ,60326 + ,1528 + ,414462 + ,25 + ,1 + ,94670 + ,3122 + ,364304 + ,77 + ,1 + ,82425 + ,7561 + ,355178 + ,54 + ,0 + ,59017 + ,2675 + ,357760 + ,251 + ,0 + ,90829 + ,13253 + ,261216 + ,15 + ,0 + ,80791 + ,880 + ,397144 + ,44 + ,1 + ,100423 + ,2053 + ,374943 + ,73 + ,0 + ,131116 + ,1424 + ,424898 + ,85 + ,1 + ,100269 + ,4036 + ,202055 + ,49 + ,1 + ,27330 + ,3045 + ,378525 + ,38 + ,0 + ,39039 + ,5119 + ,310768 + ,35 + ,0 + ,106885 + ,1431 + ,325738 + ,9 + ,0 + ,79285 + ,554 + ,394510 + ,34 + ,0 + ,118881 + ,1975 + ,247060 + ,20 + ,1 + ,77623 + ,1765 + ,368078 + ,29 + ,0 + ,114768 + ,1012 + ,236761 + ,11 + ,0 + ,74015 + ,810 + ,312378 + ,52 + ,0 + ,69465 + ,1280 + ,339836 + ,13 + ,1 + ,117869 + ,666 + ,347385 + ,29 + ,0 + ,60982 + ,1380 + ,426280 + ,66 + ,1 + ,90131 + ,4677 + ,352850 + ,33 + ,0 + ,138971 + ,876 + ,301881 + ,15 + ,0 + ,39625 + ,814 + ,377516 + ,15 + ,0 + ,102725 + ,514 + ,357312 + ,68 + ,1 + ,64239 + ,5692 + ,458343 + ,100 + ,0 + ,90262 + ,3642 + ,354228 + ,13 + ,0 + ,103960 + ,540 + ,308636 + ,45 + ,0 + ,106611 + ,2099 + ,386212 + ,14 + ,0 + ,103345 + ,567 + ,393343 + ,36 + ,0 + ,95551 + ,2001 + ,378509 + ,40 + ,1 + ,82903 + ,2949 + ,452469 + ,68 + ,0 + ,63593 + ,2253 + ,364839 + ,29 + ,1 + ,126910 + ,6533 + ,358649 + ,43 + ,0 + ,37527 + ,1889 + ,376641 + ,30 + ,1 + ,60247 + ,3055 + ,429112 + ,9 + ,0 + ,112995 + ,272 + ,330546 + ,22 + ,1 + ,70184 + ,1414 + ,403560 + ,19 + ,0 + ,130140 + ,2564 + ,317892 + ,9 + ,1 + ,73221 + ,1383 + ,307528 + ,31 + ,1 + ,76114 + ,1261 + ,235133 + ,19 + ,0 + ,90534 + ,975 + ,299243 + ,55 + ,0 + ,108479 + ,3366 + ,314073 + ,8 + ,0 + ,113761 + ,576 + ,368186 + ,28 + ,0 + ,68696 + ,1686 + ,269661 + ,29 + ,0 + ,71561 + ,746 + ,125390 + ,48 + ,1 + ,59831 + ,3192 + ,510834 + ,16 + ,1 + ,97890 + ,2045 + ,321896 + ,47 + ,0 + ,101481 + ,5702 + ,249898 + ,20 + ,1 + ,72954 + ,1932 + ,408881 + ,22 + ,0 + ,67939 + ,936 + ,158492 + ,33 + ,1 + ,48022 + ,3437 + ,292154 + ,44 + ,0 + ,86111 + ,5131 + ,289513 + ,13 + ,1 + ,74020 + ,2397 + ,378049 + ,6 + ,1 + ,57530 + ,1389 + ,343466 + ,35 + ,0 + ,56364 + ,1503 + ,332743 + ,8 + ,0 + ,84990 + ,402 + ,442882 + ,17 + ,0 + ,88590 + ,2239 + ,214215 + ,11 + ,1 + ,77200 + ,2234 + ,315688 + ,21 + ,0 + ,61262 + ,837 + ,375195 + ,92 + ,0 + ,110309 + ,10579 + ,334280 + ,12 + ,0 + ,67000 + ,875 + ,355864 + ,112 + ,0 + ,93099 + ,1585 + ,480382 + ,25 + ,1 + ,107577 + ,1659 + ,353058 + ,17 + ,1 + ,62920 + ,2647 + ,217193 + ,23 + ,1 + ,75832 + ,3294 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,314533 + ,10 + ,0 + ,60793 + ,94 + ,318056 + ,23 + ,0 + ,57935 + ,422 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,314353 + ,7 + ,0 + ,60630 + ,34 + ,369448 + ,25 + ,0 + ,55637 + ,1558 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,312846 + ,20 + ,0 + ,60887 + ,43 + ,312075 + ,4 + ,0 + ,60720 + ,645 + ,315009 + ,4 + ,0 + ,60505 + ,316 + ,318903 + ,10 + ,0 + ,60945 + ,115 + ,314887 + ,1 + ,0 + ,60720 + ,5 + ,314913 + ,4 + ,0 + ,60720 + ,897 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,325506 + ,8 + ,0 + ,58990 + ,389 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,298568 + ,11 + ,0 + ,56750 + ,1002 + ,315834 + ,4 + ,0 + ,60894 + ,36 + ,329784 + ,15 + ,0 + ,63346 + ,460 + ,312878 + ,9 + ,0 + ,56535 + ,309 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,314987 + ,7 + ,0 + ,60835 + ,9 + ,325249 + ,2 + ,0 + ,60720 + ,271 + ,315877 + ,0 + ,0 + ,61016 + ,14 + ,291650 + ,7 + ,0 + ,58650 + ,520 + ,305959 + ,46 + ,0 + ,60438 + ,1766 + ,315380 + ,5 + ,0 + ,60720 + ,0 + ,297765 + ,7 + ,0 + ,58625 + ,458 + ,315245 + ,2 + ,0 + ,60938 + ,20 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315236 + ,2 + ,0 + ,61490 + ,98 + ,336425 + ,5 + ,0 + ,60845 + ,405 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,306268 + ,7 + ,0 + ,60830 + ,483 + ,302187 + ,24 + ,0 + ,63261 + ,454 + ,314882 + ,1 + ,0 + ,60720 + ,47 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,382712 + ,18 + ,0 + ,45689 + ,757 + ,341570 + ,55 + ,0 + ,60720 + ,4655 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,312412 + ,3 + ,0 + ,61564 + ,36 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,309596 + ,9 + ,0 + ,61938 + ,203 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315547 + ,8 + ,0 + ,60951 + ,126 + ,313267 + ,113 + ,0 + ,60720 + ,400 + ,316176 + ,0 + ,1 + ,60745 + ,71 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,359335 + ,19 + ,0 + ,71642 + ,972 + ,330068 + ,11 + ,1 + ,71641 + ,531 + ,314289 + ,25 + ,0 + ,55792 + ,2461 + ,297413 + ,16 + ,1 + ,71873 + ,378 + ,314806 + ,5 + ,1 + ,62555 + ,23 + ,333210 + ,11 + ,1 + ,60370 + ,638 + ,352108 + ,23 + ,1 + ,64873 + ,2300 + ,313332 + ,6 + ,0 + ,62041 + ,149 + ,291787 + ,5 + ,0 + ,65745 + ,226 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,318745 + ,7 + ,0 + ,59500 + ,275 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315366 + ,7 + ,0 + ,61630 + ,141 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315688 + ,3 + ,0 + ,60890 + ,28 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,409642 + ,89 + ,1 + ,113521 + ,4980 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,269587 + ,19 + ,1 + ,80045 + ,472 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,300962 + ,12 + ,1 + ,50804 + ,203 + ,325479 + ,12 + ,1 + ,87390 + ,496 + ,316155 + ,5 + ,1 + ,61656 + ,10 + ,318574 + ,2 + ,1 + ,65688 + ,63 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,343613 + ,26 + ,1 + ,48522 + ,1136 + ,306948 + ,3 + ,1 + ,60720 + ,265 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,330059 + ,11 + ,0 + ,57640 + ,267 + ,288985 + ,10 + ,0 + ,61977 + ,474 + ,304485 + ,5 + ,0 + ,62620 + ,534 + ,315380 + ,2 + ,0 + ,60720 + ,0 + ,315688 + ,6 + ,0 + ,60831 + ,15 + ,317736 + ,7 + ,0 + ,60646 + ,397 + ,315380 + ,2 + ,0 + ,60720 + ,0 + ,322331 + ,28 + ,0 + ,56225 + ,1866 + ,296656 + ,3 + ,0 + ,60510 + ,288 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315354 + ,1 + ,0 + ,60698 + ,3 + ,312161 + ,20 + ,0 + ,60720 + ,468 + ,315576 + ,1 + ,0 + ,60805 + ,20 + ,314922 + ,22 + ,0 + ,61404 + ,278 + ,314551 + ,9 + ,0 + ,60720 + ,61 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,312339 + ,2 + ,0 + ,65276 + ,192 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,298700 + ,7 + ,0 + ,63915 + ,317 + ,321376 + ,9 + ,0 + ,60720 + ,738 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,303230 + ,13 + ,0 + ,61686 + ,368 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315487 + ,0 + ,0 + ,60743 + ,2 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315793 + ,6 + ,0 + ,60349 + ,53 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,312887 + ,3 + ,0 + ,61360 + ,94 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315637 + ,7 + ,0 + ,59818 + ,24 + ,324385 + ,2 + ,0 + ,72680 + ,2332 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,308989 + ,15 + ,0 + ,61808 + ,131 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,296702 + ,9 + ,0 + ,53110 + ,206 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,307322 + ,1 + ,0 + ,64245 + ,167 + ,304376 + ,38 + ,0 + ,73007 + ,622 + ,253588 + ,57 + ,0 + ,82732 + ,2328 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,309560 + ,7 + ,0 + ,54820 + ,365 + ,298466 + ,26 + ,0 + ,47705 + ,364 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,343929 + ,13 + ,0 + ,72835 + ,226 + ,331955 + ,10 + ,0 + ,58856 + ,307 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,381180 + ,9 + ,0 + ,77655 + ,188 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,331420 + ,26 + ,0 + ,69817 + ,138 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,310201 + ,19 + ,0 + ,60798 + ,125 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,320016 + ,12 + ,0 + ,62452 + ,282 + ,320398 + ,23 + ,0 + ,64175 + ,335 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,291841 + ,29 + ,1 + ,67440 + ,1324 + ,310670 + ,8 + ,0 + ,68136 + ,176 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,313491 + ,26 + ,0 + ,56726 + ,249 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,331323 + ,9 + ,0 + ,70811 + ,333 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,319210 + ,5 + ,1 + ,60720 + ,601 + ,318098 + ,3 + ,0 + ,62045 + ,30 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,292754 + ,13 + ,0 + ,54323 + ,249 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,325176 + ,12 + ,0 + ,62841 + ,165 + ,365959 + ,19 + ,0 + ,81125 + ,453 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,302409 + ,10 + ,0 + ,59506 + ,53 + ,340968 + ,9 + ,1 + ,59365 + ,382 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,9 + ,0 + ,60720 + ,0 + ,313164 + ,4 + ,1 + ,60798 + ,30 + ,301164 + ,1 + ,0 + ,58790 + ,290 + ,315380 + ,1 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,344425 + ,14 + ,0 + ,61808 + ,366 + ,315394 + ,12 + ,0 + ,60735 + ,2 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,316647 + ,19 + ,1 + ,64016 + ,209 + ,309836 + ,17 + ,0 + ,54683 + ,384 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,346611 + ,32 + ,0 + ,87192 + ,365 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,322031 + ,14 + ,1 + ,64107 + ,49 + ,315656 + ,8 + ,0 + ,60761 + ,3 + ,339445 + ,4 + ,0 + ,65990 + ,133 + ,314964 + ,0 + ,0 + ,59988 + ,32 + ,297141 + ,20 + ,0 + ,61167 + ,368 + ,315372 + ,5 + ,0 + ,60719 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,312502 + ,1 + ,0 + ,60722 + ,22 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,315380 + ,0 + ,1 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,313729 + ,4 + ,0 + ,60379 + ,96 + ,315388 + ,1 + ,0 + ,60727 + ,1 + ,315371 + ,4 + ,0 + ,60720 + ,314 + ,296139 + ,20 + ,0 + ,60925 + ,844 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,313880 + ,1 + ,0 + ,60896 + ,26 + ,317698 + ,10 + ,0 + ,59734 + ,125 + ,295580 + ,12 + ,0 + ,62969 + ,304 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,308256 + ,13 + ,0 + ,60720 + ,621 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,303677 + ,3 + ,0 + ,59118 + ,119 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,319369 + ,10 + ,0 + ,60720 + ,1595 + ,318690 + ,3 + ,0 + ,58598 + ,312 + ,314049 + ,7 + ,0 + ,61124 + ,60 + ,325699 + ,10 + ,0 + ,59595 + ,587 + ,314210 + ,1 + ,0 + ,62065 + ,135 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,322378 + ,15 + ,0 + ,78780 + ,514 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315398 + ,4 + ,0 + ,60722 + ,1 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,308336 + ,28 + ,1 + ,61600 + ,1763 + ,316386 + ,9 + ,0 + ,59635 + ,180 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315553 + ,7 + ,0 + ,60720 + ,218 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,323361 + ,7 + ,0 + ,59781 + ,448 + ,336639 + ,7 + ,0 + ,76644 + ,227 + ,307424 + ,3 + ,0 + ,64820 + ,174 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,295370 + ,11 + ,0 + ,56178 + ,121 + ,322340 + ,7 + ,0 + ,60436 + ,607 + ,319864 + ,10 + ,0 + ,60720 + ,2212 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,317291 + ,18 + ,0 + ,73433 + ,530 + ,280398 + ,14 + ,0 + ,41477 + ,571 + ,315380 + ,0 + ,0 + ,60720 + ,0 + ,317330 + ,12 + ,0 + ,62700 + ,78 + ,238125 + ,29 + ,0 + ,67804 + ,2489 + ,327071 + ,3 + ,0 + ,59661 + ,131 + ,309038 + ,6 + ,0 + ,58620 + ,923 + ,314210 + ,3 + ,0 + ,60398 + ,72 + ,307930 + ,8 + ,0 + ,58580 + ,572 + ,322327 + ,10 + ,0 + ,62710 + ,397 + ,292136 + ,6 + ,0 + ,59325 + ,450 + ,263276 + ,8 + ,0 + ,60950 + ,622 + ,367655 + ,6 + ,0 + ,68060 + ,694 + ,283910 + ,9 + ,1 + ,83620 + ,3425 + ,283587 + ,8 + ,0 + ,58456 + ,562 + ,243650 + ,26 + ,0 + ,52811 + ,4917 + ,438493 + ,239 + ,1 + ,121173 + ,1442 + ,296261 + ,7 + ,0 + ,63870 + ,529 + ,230621 + ,41 + ,1 + ,21001 + ,2126 + ,304252 + ,3 + ,0 + ,70415 + ,1061 + ,333505 + ,8 + ,0 + ,64230 + ,776 + ,296919 + ,6 + ,0 + ,59190 + ,611 + ,278990 + ,21 + ,1 + ,69351 + ,1526 + ,276898 + ,7 + ,0 + ,64270 + ,592 + ,327007 + ,11 + ,0 + ,70694 + ,1182 + ,317046 + ,11 + ,0 + ,68005 + ,621 + ,304555 + ,12 + ,0 + ,58930 + ,989 + ,298096 + ,9 + ,0 + ,58320 + ,438 + ,231861 + ,3 + ,0 + ,69980 + ,726 + ,309422 + ,57 + ,0 + ,69863 + ,1303 + ,286963 + ,21 + ,1 + ,63255 + ,7419 + ,269753 + ,15 + ,1 + ,57320 + ,1164 + ,448243 + ,32 + ,1 + ,75230 + ,3310 + ,165404 + ,11 + ,0 + ,79420 + ,1920 + ,204325 + ,2 + ,0 + ,73490 + ,965 + ,407159 + ,23 + ,0 + ,35250 + ,3256 + ,290476 + ,20 + ,1 + ,62285 + ,1135 + ,275311 + ,24 + ,0 + ,69206 + ,1270 + ,246541 + ,1 + ,0 + ,65920 + ,661 + ,253468 + ,1 + ,0 + ,69770 + ,1013 + ,240897 + ,74 + ,0 + ,72683 + ,2844 + ,-83265 + ,68 + ,1 + ,-14545 + ,11528 + ,-42143 + ,20 + ,0 + ,55830 + ,6526 + ,272713 + ,20 + ,0 + ,55174 + ,2264 + ,215362 + ,82 + ,1 + ,67038 + ,5109 + ,42754 + ,21 + ,0 + ,51252 + ,3999 + ,306275 + ,244 + ,0 + ,157278 + ,35624 + ,253537 + ,32 + ,0 + ,79510 + ,9252 + ,372631 + ,86 + ,0 + ,77440 + ,15236 + ,-7170 + ,69 + ,0 + ,27284 + ,18073) + ,dim=c(5 + ,431) + ,dimnames=list(c('wealth' + ,'Orders' + ,'Group' + ,'Dividends' + ,'Costs ') + ,1:431)) > y <- array(NA,dim=c(5,431),dimnames=list(c('wealth','Orders','Group','Dividends','Costs '),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 = '1' > #'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 Orders Group Dividends Costs\r 1 6282154 807 1 213118 162556 2 4321023 444 1 81767 29790 3 4111912 412 1 153198 87550 4 223193 428 0 -26007 84738 5 1491348 315 1 126942 54660 6 1629616 168 1 157214 42634 7 1398893 263 0 129352 40949 8 1926517 267 1 234817 45187 9 983660 228 1 60448 37704 10 1443586 129 1 47818 16275 11 1073089 104 0 245546 25830 12 984885 122 0 48020 12679 13 1405225 393 1 -1710 18014 14 227132 190 0 32648 43556 15 929118 280 1 95350 24811 16 1071292 63 0 151352 6575 17 638830 102 0 288170 7123 18 856956 265 1 114337 21950 19 992426 234 1 37884 37597 20 444477 277 0 122844 17821 21 857217 73 1 82340 12988 22 711969 67 1 79801 22330 23 702380 103 0 165548 13326 24 358589 290 0 116384 16189 25 297978 83 0 134028 7146 26 585715 56 0 63838 15824 27 657954 236 1 74996 27664 28 209458 73 0 31080 11920 29 786690 34 0 32168 8568 30 439798 139 0 49857 14416 31 688779 26 1 87161 3369 32 574339 70 1 106113 11819 33 741409 40 1 80570 6984 34 597793 42 1 102129 4519 35 644190 12 0 301670 2220 36 377934 211 0 102313 18562 37 640273 74 0 88577 10327 38 697458 80 1 112477 5336 39 550608 83 1 191778 2365 40 207393 131 0 79804 4069 41 301607 203 0 128294 8636 42 345783 56 0 96448 13718 43 501749 89 0 93811 4525 44 379983 88 0 117520 6869 45 387475 39 0 69159 4628 46 377305 25 1 101792 3689 47 370837 49 1 210568 4891 48 430866 149 1 136996 7489 49 469107 58 0 121920 4901 50 194493 41 0 76403 2284 51 530670 90 1 108094 3160 52 518365 136 1 134759 4150 53 491303 97 1 188873 7285 54 527021 63 1 146216 1134 55 233773 114 1 156608 4658 56 405972 77 0 61348 2384 57 652925 6 0 50350 3748 58 446211 47 0 87720 5371 59 341340 51 0 99489 1285 60 387699 85 1 87419 9327 61 493408 43 1 94355 5565 62 146494 32 0 60326 1528 63 414462 25 1 94670 3122 64 364304 77 1 82425 7561 65 355178 54 0 59017 2675 66 357760 251 0 90829 13253 67 261216 15 0 80791 880 68 397144 44 1 100423 2053 69 374943 73 0 131116 1424 70 424898 85 1 100269 4036 71 202055 49 1 27330 3045 72 378525 38 0 39039 5119 73 310768 35 0 106885 1431 74 325738 9 0 79285 554 75 394510 34 0 118881 1975 76 247060 20 1 77623 1765 77 368078 29 0 114768 1012 78 236761 11 0 74015 810 79 312378 52 0 69465 1280 80 339836 13 1 117869 666 81 347385 29 0 60982 1380 82 426280 66 1 90131 4677 83 352850 33 0 138971 876 84 301881 15 0 39625 814 85 377516 15 0 102725 514 86 357312 68 1 64239 5692 87 458343 100 0 90262 3642 88 354228 13 0 103960 540 89 308636 45 0 106611 2099 90 386212 14 0 103345 567 91 393343 36 0 95551 2001 92 378509 40 1 82903 2949 93 452469 68 0 63593 2253 94 364839 29 1 126910 6533 95 358649 43 0 37527 1889 96 376641 30 1 60247 3055 97 429112 9 0 112995 272 98 330546 22 1 70184 1414 99 403560 19 0 130140 2564 100 317892 9 1 73221 1383 101 307528 31 1 76114 1261 102 235133 19 0 90534 975 103 299243 55 0 108479 3366 104 314073 8 0 113761 576 105 368186 28 0 68696 1686 106 269661 29 0 71561 746 107 125390 48 1 59831 3192 108 510834 16 1 97890 2045 109 321896 47 0 101481 5702 110 249898 20 1 72954 1932 111 408881 22 0 67939 936 112 158492 33 1 48022 3437 113 292154 44 0 86111 5131 114 289513 13 1 74020 2397 115 378049 6 1 57530 1389 116 343466 35 0 56364 1503 117 332743 8 0 84990 402 118 442882 17 0 88590 2239 119 214215 11 1 77200 2234 120 315688 21 0 61262 837 121 375195 92 0 110309 10579 122 334280 12 0 67000 875 123 355864 112 0 93099 1585 124 480382 25 1 107577 1659 125 353058 17 1 62920 2647 126 217193 23 1 75832 3294 127 315380 0 0 60720 0 128 314533 10 0 60793 94 129 318056 23 0 57935 422 130 315380 0 0 60720 0 131 314353 7 0 60630 34 132 369448 25 0 55637 1558 133 315380 0 0 60720 0 134 312846 20 0 60887 43 135 312075 4 0 60720 645 136 315009 4 0 60505 316 137 318903 10 0 60945 115 138 314887 1 0 60720 5 139 314913 4 0 60720 897 140 315380 0 0 60720 0 141 325506 8 0 58990 389 142 315380 0 0 60720 0 143 298568 11 0 56750 1002 144 315834 4 0 60894 36 145 329784 15 0 63346 460 146 312878 9 0 56535 309 147 315380 0 0 60720 0 148 314987 7 0 60835 9 149 325249 2 0 60720 271 150 315877 0 0 61016 14 151 291650 7 0 58650 520 152 305959 46 0 60438 1766 153 315380 5 0 60720 0 154 297765 7 0 58625 458 155 315245 2 0 60938 20 156 315380 0 0 60720 0 157 315380 0 0 60720 0 158 315236 2 0 61490 98 159 336425 5 0 60845 405 160 315380 0 0 60720 0 161 315380 0 0 60720 0 162 315380 0 0 60720 0 163 315380 0 0 60720 0 164 306268 7 0 60830 483 165 302187 24 0 63261 454 166 314882 1 0 60720 47 167 315380 0 0 60720 0 168 382712 18 0 45689 757 169 341570 55 0 60720 4655 170 315380 0 0 60720 0 171 315380 0 0 60720 0 172 312412 3 0 61564 36 173 315380 0 0 60720 0 174 309596 9 0 61938 203 175 315380 0 0 60720 0 176 315547 8 0 60951 126 177 313267 113 0 60720 400 178 316176 0 1 60745 71 179 315380 0 0 60720 0 180 315380 0 0 60720 0 181 359335 19 0 71642 972 182 330068 11 1 71641 531 183 314289 25 0 55792 2461 184 297413 16 1 71873 378 185 314806 5 1 62555 23 186 333210 11 1 60370 638 187 352108 23 1 64873 2300 188 313332 6 0 62041 149 189 291787 5 0 65745 226 190 315380 0 0 60720 0 191 318745 7 0 59500 275 192 315380 0 0 60720 0 193 315366 7 0 61630 141 194 315380 0 0 60720 0 195 315688 3 0 60890 28 196 315380 0 1 60720 0 197 409642 89 1 113521 4980 198 315380 0 1 60720 0 199 315380 0 1 60720 0 200 269587 19 1 80045 472 201 315380 0 1 60720 0 202 315380 0 1 60720 0 203 315380 0 1 60720 0 204 300962 12 1 50804 203 205 325479 12 1 87390 496 206 316155 5 1 61656 10 207 318574 2 1 65688 63 208 315380 0 1 60720 0 209 343613 26 1 48522 1136 210 306948 3 1 60720 265 211 315380 0 0 60720 0 212 315380 0 0 60720 0 213 330059 11 0 57640 267 214 288985 10 0 61977 474 215 304485 5 0 62620 534 216 315380 2 0 60720 0 217 315688 6 0 60831 15 218 317736 7 0 60646 397 219 315380 2 0 60720 0 220 322331 28 0 56225 1866 221 296656 3 0 60510 288 222 315380 0 0 60720 0 223 315354 1 0 60698 3 224 312161 20 0 60720 468 225 315576 1 0 60805 20 226 314922 22 0 61404 278 227 314551 9 0 60720 61 228 315380 0 0 60720 0 229 312339 2 0 65276 192 230 315380 0 0 60720 0 231 298700 7 0 63915 317 232 321376 9 0 60720 738 233 315380 0 0 60720 0 234 303230 13 0 61686 368 235 315380 0 0 60720 0 236 315487 0 0 60743 2 237 315380 0 0 60720 0 238 315793 6 0 60349 53 239 315380 0 0 60720 0 240 315380 0 0 60720 0 241 315380 0 0 60720 0 242 312887 3 0 61360 94 243 315380 0 0 60720 0 244 315637 7 0 59818 24 245 324385 2 0 72680 2332 246 315380 0 1 60720 0 247 315380 0 1 60720 0 248 308989 15 0 61808 131 249 315380 0 1 60720 0 250 315380 0 1 60720 0 251 296702 9 0 53110 206 252 315380 0 1 60720 0 253 307322 1 0 64245 167 254 304376 38 0 73007 622 255 253588 57 0 82732 2328 256 315380 0 0 60720 0 257 309560 7 0 54820 365 258 298466 26 0 47705 364 259 315380 0 1 60720 0 260 315380 0 0 60720 0 261 315380 0 1 60720 0 262 315380 0 1 60720 0 263 343929 13 0 72835 226 264 331955 10 0 58856 307 265 315380 0 1 60720 0 266 315380 0 0 60720 0 267 315380 0 1 60720 0 268 381180 9 0 77655 188 269 315380 0 1 60720 0 270 331420 26 0 69817 138 271 315380 0 1 60720 0 272 315380 0 1 60720 0 273 315380 0 1 60720 0 274 310201 19 0 60798 125 275 315380 0 0 60720 0 276 320016 12 0 62452 282 277 320398 23 0 64175 335 278 315380 0 0 60720 0 279 291841 29 1 67440 1324 280 310670 8 0 68136 176 281 315380 0 1 60720 0 282 315380 0 0 60720 0 283 313491 26 0 56726 249 284 315380 0 0 60720 0 285 331323 9 0 70811 333 286 315380 0 0 60720 0 287 319210 5 1 60720 601 288 318098 3 0 62045 30 289 315380 0 1 60720 0 290 292754 13 0 54323 249 291 315380 0 1 60720 0 292 325176 12 0 62841 165 293 365959 19 0 81125 453 294 315380 0 0 60720 0 295 302409 10 0 59506 53 296 340968 9 1 59365 382 297 315380 0 0 60720 0 298 315380 0 1 60720 0 299 315380 0 0 60720 0 300 315380 9 0 60720 0 301 313164 4 1 60798 30 302 301164 1 0 58790 290 303 315380 1 0 60720 0 304 315380 0 1 60720 0 305 344425 14 0 61808 366 306 315394 12 0 60735 2 307 315380 0 1 60720 0 308 316647 19 1 64016 209 309 309836 17 0 54683 384 310 315380 0 1 60720 0 311 315380 0 1 60720 0 312 346611 32 0 87192 365 313 315380 0 1 60720 0 314 322031 14 1 64107 49 315 315656 8 0 60761 3 316 339445 4 0 65990 133 317 314964 0 0 59988 32 318 297141 20 0 61167 368 319 315372 5 0 60719 1 320 315380 0 0 60720 0 321 315380 0 0 60720 0 322 315380 0 0 60720 0 323 315380 0 0 60720 0 324 315380 0 1 60720 0 325 315380 0 0 60720 0 326 312502 1 0 60722 22 327 315380 0 1 60720 0 328 315380 0 1 60720 0 329 315380 0 0 60720 0 330 315380 0 0 60720 0 331 315380 0 0 60720 0 332 315380 0 0 60720 0 333 315380 0 0 60720 0 334 313729 4 0 60379 96 335 315388 1 0 60727 1 336 315371 4 0 60720 314 337 296139 20 0 60925 844 338 315380 0 0 60720 0 339 313880 1 0 60896 26 340 317698 10 0 59734 125 341 295580 12 0 62969 304 342 315380 0 0 60720 0 343 315380 0 0 60720 0 344 315380 0 0 60720 0 345 308256 13 0 60720 621 346 315380 0 0 60720 0 347 303677 3 0 59118 119 348 315380 0 0 60720 0 349 315380 0 0 60720 0 350 319369 10 0 60720 1595 351 318690 3 0 58598 312 352 314049 7 0 61124 60 353 325699 10 0 59595 587 354 314210 1 0 62065 135 355 315380 0 0 60720 0 356 315380 0 0 60720 0 357 322378 15 0 78780 514 358 315380 0 0 60720 0 359 315380 0 0 60720 0 360 315380 0 0 60720 0 361 315398 4 0 60722 1 362 315380 0 0 60720 0 363 315380 0 0 60720 0 364 308336 28 1 61600 1763 365 316386 9 0 59635 180 366 315380 0 0 60720 0 367 315380 0 0 60720 0 368 315380 0 0 60720 0 369 315380 0 0 60720 0 370 315553 7 0 60720 218 371 315380 0 0 60720 0 372 323361 7 0 59781 448 373 336639 7 0 76644 227 374 307424 3 0 64820 174 375 315380 0 0 60720 0 376 315380 0 0 60720 0 377 295370 11 0 56178 121 378 322340 7 0 60436 607 379 319864 10 0 60720 2212 380 315380 0 0 60720 0 381 315380 0 0 60720 0 382 317291 18 0 73433 530 383 280398 14 0 41477 571 384 315380 0 0 60720 0 385 317330 12 0 62700 78 386 238125 29 0 67804 2489 387 327071 3 0 59661 131 388 309038 6 0 58620 923 389 314210 3 0 60398 72 390 307930 8 0 58580 572 391 322327 10 0 62710 397 392 292136 6 0 59325 450 393 263276 8 0 60950 622 394 367655 6 0 68060 694 395 283910 9 1 83620 3425 396 283587 8 0 58456 562 397 243650 26 0 52811 4917 398 438493 239 1 121173 1442 399 296261 7 0 63870 529 400 230621 41 1 21001 2126 401 304252 3 0 70415 1061 402 333505 8 0 64230 776 403 296919 6 0 59190 611 404 278990 21 1 69351 1526 405 276898 7 0 64270 592 406 327007 11 0 70694 1182 407 317046 11 0 68005 621 408 304555 12 0 58930 989 409 298096 9 0 58320 438 410 231861 3 0 69980 726 411 309422 57 0 69863 1303 412 286963 21 1 63255 7419 413 269753 15 1 57320 1164 414 448243 32 1 75230 3310 415 165404 11 0 79420 1920 416 204325 2 0 73490 965 417 407159 23 0 35250 3256 418 290476 20 1 62285 1135 419 275311 24 0 69206 1270 420 246541 1 0 65920 661 421 253468 1 0 69770 1013 422 240897 74 0 72683 2844 423 -83265 68 1 -14545 11528 424 -42143 20 0 55830 6526 425 272713 20 0 55174 2264 426 215362 82 1 67038 5109 427 42754 21 0 51252 3999 428 306275 244 0 157278 35624 429 253537 32 0 79510 9252 430 372631 86 0 77440 15236 431 -7170 69 0 27284 18073 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Orders Group Dividends `Costs\r` 1.573e+05 9.510e+02 4.756e+04 1.518e+00 2.204e+01 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -2169550 -26507 35623 65929 2913139 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 1.573e+05 3.090e+04 5.089 5.4e-07 *** Orders 9.510e+02 3.506e+02 2.713 0.006947 ** Group 4.756e+04 2.793e+04 1.703 0.089345 . Dividends 1.518e+00 4.116e-01 3.689 0.000254 *** `Costs\r` 2.204e+01 2.150e+00 10.253 < 2e-16 *** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 253000 on 426 degrees of freedom Multiple R-squared: 0.658, Adjusted R-squared: 0.6547 F-statistic: 204.9 on 4 and 426 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 8.929724e-44 4.464862e-44 [2,] 1.0000000 1.710130e-43 8.550649e-44 [3,] 1.0000000 1.472796e-56 7.363978e-57 [4,] 1.0000000 4.388800e-60 2.194400e-60 [5,] 1.0000000 6.105175e-74 3.052588e-74 [6,] 1.0000000 8.510729e-98 4.255365e-98 [7,] 1.0000000 1.060595e-99 5.302973e-100 [8,] 1.0000000 8.189103e-111 4.094551e-111 [9,] 1.0000000 1.455076e-120 7.275382e-121 [10,] 1.0000000 3.847186e-126 1.923593e-126 [11,] 1.0000000 1.133439e-134 5.667195e-135 [12,] 1.0000000 1.932293e-139 9.661465e-140 [13,] 1.0000000 1.562321e-142 7.811607e-143 [14,] 1.0000000 6.183486e-148 3.091743e-148 [15,] 1.0000000 1.613759e-149 8.068793e-150 [16,] 1.0000000 2.276656e-150 1.138328e-150 [17,] 1.0000000 2.100393e-152 1.050196e-152 [18,] 1.0000000 4.716624e-152 2.358312e-152 [19,] 1.0000000 3.782349e-157 1.891175e-157 [20,] 1.0000000 2.669433e-162 1.334716e-162 [21,] 1.0000000 7.852758e-163 3.926379e-163 [22,] 1.0000000 2.042638e-179 1.021319e-179 [23,] 1.0000000 9.704061e-180 4.852030e-180 [24,] 1.0000000 8.316716e-184 4.158358e-184 [25,] 1.0000000 2.684244e-185 1.342122e-185 [26,] 1.0000000 2.821387e-193 1.410693e-193 [27,] 1.0000000 3.137464e-195 1.568732e-195 [28,] 1.0000000 3.129579e-194 1.564790e-194 [29,] 1.0000000 6.222566e-194 3.111283e-194 [30,] 1.0000000 1.561918e-200 7.809590e-201 [31,] 1.0000000 4.094542e-206 2.047271e-206 [32,] 1.0000000 1.694215e-206 8.471075e-207 [33,] 1.0000000 9.272081e-207 4.636041e-207 [34,] 1.0000000 2.923527e-207 1.461763e-207 [35,] 1.0000000 1.286388e-206 6.431941e-207 [36,] 1.0000000 1.569793e-207 7.848964e-208 [37,] 1.0000000 1.744535e-206 8.722676e-207 [38,] 1.0000000 7.524827e-207 3.762413e-207 [39,] 1.0000000 5.955238e-206 2.977619e-206 [40,] 1.0000000 5.505826e-208 2.752913e-208 [41,] 1.0000000 1.962118e-208 9.810589e-209 [42,] 1.0000000 3.513849e-208 1.756924e-208 [43,] 1.0000000 1.164878e-208 5.824391e-209 [44,] 1.0000000 4.067565e-209 2.033783e-209 [45,] 1.0000000 3.537442e-209 1.768721e-209 [46,] 1.0000000 3.931178e-209 1.965589e-209 [47,] 1.0000000 1.181033e-208 5.905163e-209 [48,] 1.0000000 9.862900e-213 4.931450e-213 [49,] 1.0000000 8.382714e-213 4.191357e-213 [50,] 1.0000000 4.190671e-227 2.095336e-227 [51,] 1.0000000 8.613017e-228 4.306509e-228 [52,] 1.0000000 5.429237e-227 2.714619e-227 [53,] 1.0000000 8.915173e-227 4.457587e-227 [54,] 1.0000000 2.834720e-228 1.417360e-228 [55,] 1.0000000 6.806395e-230 3.403198e-230 [56,] 1.0000000 2.652220e-229 1.326110e-229 [57,] 1.0000000 1.001598e-228 5.007990e-229 [58,] 1.0000000 3.221348e-228 1.610674e-228 [59,] 1.0000000 6.432763e-228 3.216382e-228 [60,] 1.0000000 1.089958e-227 5.449789e-228 [61,] 1.0000000 1.117477e-226 5.587384e-227 [62,] 1.0000000 5.704019e-226 2.852010e-226 [63,] 1.0000000 2.524624e-225 1.262312e-225 [64,] 1.0000000 1.877316e-224 9.386578e-225 [65,] 1.0000000 8.390189e-226 4.195095e-226 [66,] 1.0000000 2.237458e-225 1.118729e-225 [67,] 1.0000000 1.134072e-224 5.670358e-225 [68,] 1.0000000 1.010447e-223 5.052234e-224 [69,] 1.0000000 2.371207e-223 1.185603e-223 [70,] 1.0000000 1.581880e-222 7.909402e-223 [71,] 1.0000000 1.595748e-222 7.978741e-223 [72,] 1.0000000 1.272961e-221 6.364805e-222 [73,] 1.0000000 4.949137e-221 2.474569e-221 [74,] 1.0000000 1.947523e-220 9.737613e-221 [75,] 1.0000000 3.071400e-220 1.535700e-220 [76,] 1.0000000 2.706247e-220 1.353123e-220 [77,] 1.0000000 1.000773e-219 5.003864e-220 [78,] 1.0000000 7.101104e-219 3.550552e-219 [79,] 1.0000000 1.642677e-218 8.213386e-219 [80,] 1.0000000 1.183264e-218 5.916319e-219 [81,] 1.0000000 6.844389e-218 3.422195e-218 [82,] 1.0000000 1.390227e-217 6.951136e-218 [83,] 1.0000000 1.030231e-216 5.151155e-217 [84,] 1.0000000 7.673293e-216 3.836647e-216 [85,] 1.0000000 5.126159e-215 2.563079e-215 [86,] 1.0000000 1.372404e-216 6.862019e-217 [87,] 1.0000000 1.071483e-215 5.357414e-216 [88,] 1.0000000 2.983363e-216 1.491681e-216 [89,] 1.0000000 3.073925e-216 1.536962e-216 [90,] 1.0000000 2.120654e-215 1.060327e-215 [91,] 1.0000000 2.348028e-214 1.174014e-214 [92,] 1.0000000 1.824227e-213 9.121133e-214 [93,] 1.0000000 2.064612e-212 1.032306e-212 [94,] 1.0000000 1.946156e-211 9.730779e-212 [95,] 1.0000000 1.420118e-212 7.100592e-213 [96,] 1.0000000 1.331733e-212 6.658663e-213 [97,] 1.0000000 1.913082e-213 9.565411e-214 [98,] 1.0000000 6.987772e-213 3.493886e-213 [99,] 1.0000000 2.149489e-212 1.074745e-212 [100,] 1.0000000 1.128984e-214 5.644918e-215 [101,] 1.0000000 5.620141e-216 2.810071e-216 [102,] 1.0000000 4.368869e-215 2.184434e-215 [103,] 1.0000000 9.236168e-215 4.618084e-215 [104,] 1.0000000 5.210613e-215 2.605307e-215 [105,] 1.0000000 2.884043e-215 1.442021e-215 [106,] 1.0000000 2.258819e-214 1.129409e-214 [107,] 1.0000000 1.881515e-213 9.407576e-214 [108,] 1.0000000 1.568024e-213 7.840122e-214 [109,] 1.0000000 5.293107e-213 2.646554e-213 [110,] 1.0000000 3.201947e-212 1.600974e-212 [111,] 1.0000000 2.139583e-212 1.069792e-212 [112,] 1.0000000 2.138771e-213 1.069386e-213 [113,] 1.0000000 1.872131e-212 9.360656e-213 [114,] 1.0000000 1.081818e-211 5.409088e-212 [115,] 1.0000000 8.438831e-211 4.219416e-211 [116,] 1.0000000 7.663780e-210 3.831890e-210 [117,] 1.0000000 1.241985e-209 6.209925e-210 [118,] 1.0000000 3.889246e-209 1.944623e-209 [119,] 1.0000000 1.130349e-209 5.651744e-210 [120,] 1.0000000 8.681176e-209 4.340588e-209 [121,] 1.0000000 7.245568e-208 3.622784e-208 [122,] 1.0000000 5.985974e-207 2.992987e-207 [123,] 1.0000000 4.768247e-206 2.384124e-206 [124,] 1.0000000 4.015630e-205 2.007815e-205 [125,] 1.0000000 2.980651e-205 1.490326e-205 [126,] 1.0000000 2.466344e-204 1.233172e-204 [127,] 1.0000000 2.287487e-203 1.143744e-203 [128,] 1.0000000 2.017381e-202 1.008690e-202 [129,] 1.0000000 1.742127e-201 8.710637e-202 [130,] 1.0000000 1.520844e-200 7.604218e-201 [131,] 1.0000000 1.309392e-199 6.546959e-200 [132,] 1.0000000 1.140130e-198 5.700652e-199 [133,] 1.0000000 9.838625e-198 4.919313e-198 [134,] 1.0000000 7.395459e-197 3.697730e-197 [135,] 1.0000000 6.439118e-196 3.219559e-196 [136,] 1.0000000 6.165250e-195 3.082625e-195 [137,] 1.0000000 5.479407e-194 2.739703e-194 [138,] 1.0000000 4.652079e-193 2.326040e-193 [139,] 1.0000000 3.960683e-192 1.980342e-192 [140,] 1.0000000 3.485863e-191 1.742932e-191 [141,] 1.0000000 3.160173e-190 1.580086e-190 [142,] 1.0000000 2.531805e-189 1.265903e-189 [143,] 1.0000000 2.241069e-188 1.120534e-188 [144,] 1.0000000 2.085650e-187 1.042825e-187 [145,] 1.0000000 2.103201e-186 1.051600e-186 [146,] 1.0000000 1.884981e-185 9.424907e-186 [147,] 1.0000000 1.768562e-184 8.842812e-185 [148,] 1.0000000 1.574150e-183 7.870749e-184 [149,] 1.0000000 1.387979e-182 6.939896e-183 [150,] 1.0000000 1.223813e-181 6.119066e-182 [151,] 1.0000000 1.095283e-180 5.476414e-181 [152,] 1.0000000 7.241611e-180 3.620805e-180 [153,] 1.0000000 6.380519e-179 3.190260e-179 [154,] 1.0000000 5.617652e-178 2.808826e-178 [155,] 1.0000000 4.941539e-177 2.470769e-177 [156,] 1.0000000 4.342227e-176 2.171114e-176 [157,] 1.0000000 4.015465e-175 2.007733e-175 [158,] 1.0000000 3.670896e-174 1.835448e-174 [159,] 1.0000000 3.222541e-173 1.611271e-173 [160,] 1.0000000 2.810352e-172 1.405176e-172 [161,] 1.0000000 3.212993e-173 1.606497e-173 [162,] 1.0000000 5.545004e-173 2.772502e-173 [163,] 1.0000000 4.904097e-172 2.452049e-172 [164,] 1.0000000 4.329083e-171 2.164542e-171 [165,] 1.0000000 3.912148e-170 1.956074e-170 [166,] 1.0000000 3.438805e-169 1.719403e-169 [167,] 1.0000000 3.151383e-168 1.575692e-168 [168,] 1.0000000 2.756594e-167 1.378297e-167 [169,] 1.0000000 2.443510e-166 1.221755e-166 [170,] 1.0000000 2.316665e-165 1.158332e-165 [171,] 1.0000000 2.180605e-164 1.090303e-164 [172,] 1.0000000 1.880971e-163 9.404853e-164 [173,] 1.0000000 1.618505e-162 8.092524e-163 [174,] 1.0000000 1.035890e-161 5.179450e-162 [175,] 1.0000000 9.663808e-161 4.831904e-161 [176,] 1.0000000 5.433790e-160 2.716895e-160 [177,] 1.0000000 3.059999e-159 1.530000e-159 [178,] 1.0000000 2.826486e-158 1.413243e-158 [179,] 1.0000000 2.071980e-157 1.035990e-157 [180,] 1.0000000 8.547270e-157 4.273635e-157 [181,] 1.0000000 7.404772e-156 3.702386e-156 [182,] 1.0000000 5.203063e-155 2.601532e-155 [183,] 1.0000000 4.353463e-154 2.176731e-154 [184,] 1.0000000 3.454640e-153 1.727320e-153 [185,] 1.0000000 2.873879e-152 1.436940e-152 [186,] 1.0000000 2.429975e-151 1.214988e-151 [187,] 1.0000000 2.009169e-150 1.004585e-150 [188,] 1.0000000 1.663834e-149 8.319169e-150 [189,] 1.0000000 1.463920e-148 7.319598e-149 [190,] 1.0000000 1.055848e-147 5.279239e-148 [191,] 1.0000000 9.200750e-147 4.600375e-147 [192,] 1.0000000 7.980481e-146 3.990240e-146 [193,] 1.0000000 6.884047e-146 3.442024e-146 [194,] 1.0000000 5.994865e-145 2.997433e-145 [195,] 1.0000000 5.196791e-144 2.598396e-144 [196,] 1.0000000 4.484394e-143 2.242197e-143 [197,] 1.0000000 3.623067e-142 1.811533e-142 [198,] 1.0000000 1.398935e-141 6.994675e-142 [199,] 1.0000000 1.199565e-140 5.997823e-141 [200,] 1.0000000 1.006251e-139 5.031253e-140 [201,] 1.0000000 8.529467e-139 4.264734e-139 [202,] 1.0000000 1.461810e-138 7.309050e-139 [203,] 1.0000000 1.238013e-137 6.190064e-138 [204,] 1.0000000 9.757329e-137 4.878665e-137 [205,] 1.0000000 7.665771e-136 3.832885e-136 [206,] 1.0000000 4.442382e-135 2.221191e-135 [207,] 1.0000000 3.322306e-134 1.661153e-134 [208,] 1.0000000 2.682326e-133 1.341163e-133 [209,] 1.0000000 2.086211e-132 1.043106e-132 [210,] 1.0000000 1.622627e-131 8.113137e-132 [211,] 1.0000000 1.219075e-130 6.095377e-131 [212,] 1.0000000 9.381112e-130 4.690556e-130 [213,] 1.0000000 4.307616e-129 2.153808e-129 [214,] 1.0000000 3.352272e-128 1.676136e-128 [215,] 1.0000000 2.551118e-127 1.275559e-127 [216,] 1.0000000 1.935829e-126 9.679143e-127 [217,] 1.0000000 1.488273e-125 7.441367e-126 [218,] 1.0000000 1.120540e-124 5.602700e-125 [219,] 1.0000000 8.556595e-124 4.278297e-124 [220,] 1.0000000 6.450997e-123 3.225498e-123 [221,] 1.0000000 4.798742e-122 2.399371e-122 [222,] 1.0000000 3.627399e-121 1.813699e-121 [223,] 1.0000000 2.678659e-120 1.339330e-120 [224,] 1.0000000 1.937733e-119 9.688663e-120 [225,] 1.0000000 1.301685e-118 6.508424e-119 [226,] 1.0000000 9.507079e-118 4.753540e-118 [227,] 1.0000000 7.176897e-117 3.588449e-117 [228,] 1.0000000 5.200686e-116 2.600343e-116 [229,] 1.0000000 3.752854e-115 1.876427e-115 [230,] 1.0000000 2.698230e-114 1.349115e-114 [231,] 1.0000000 1.922043e-113 9.610217e-114 [232,] 1.0000000 1.370863e-112 6.854316e-113 [233,] 1.0000000 9.738268e-112 4.869134e-112 [234,] 1.0000000 6.889927e-111 3.444964e-111 [235,] 1.0000000 4.938277e-110 2.469138e-110 [236,] 1.0000000 3.465376e-109 1.732688e-109 [237,] 1.0000000 2.393388e-108 1.196694e-108 [238,] 1.0000000 1.723177e-107 8.615886e-108 [239,] 1.0000000 1.247045e-106 6.235225e-107 [240,] 1.0000000 8.981264e-106 4.490632e-106 [241,] 1.0000000 6.373885e-105 3.186943e-105 [242,] 1.0000000 4.546689e-104 2.273345e-104 [243,] 1.0000000 3.227538e-103 1.613769e-103 [244,] 1.0000000 2.168018e-102 1.084009e-102 [245,] 1.0000000 1.524638e-101 7.623188e-102 [246,] 1.0000000 1.038565e-100 5.192823e-101 [247,] 1.0000000 5.954458e-100 2.977229e-100 [248,] 1.0000000 5.159507e-100 2.579754e-100 [249,] 1.0000000 3.470232e-99 1.735116e-99 [250,] 1.0000000 2.118732e-98 1.059366e-98 [251,] 1.0000000 1.164181e-97 5.820905e-98 [252,] 1.0000000 8.038866e-97 4.019433e-97 [253,] 1.0000000 5.315619e-96 2.657809e-96 [254,] 1.0000000 3.635138e-95 1.817569e-95 [255,] 1.0000000 2.473386e-94 1.236693e-94 [256,] 1.0000000 1.636538e-93 8.182691e-94 [257,] 1.0000000 8.323230e-93 4.161615e-93 [258,] 1.0000000 5.590268e-92 2.795134e-92 [259,] 1.0000000 3.599310e-91 1.799655e-91 [260,] 1.0000000 2.393728e-90 1.196864e-90 [261,] 1.0000000 1.244427e-89 6.222137e-90 [262,] 1.0000000 8.203833e-89 4.101917e-89 [263,] 1.0000000 5.354740e-88 2.677370e-88 [264,] 1.0000000 3.494243e-87 1.747121e-87 [265,] 1.0000000 2.268154e-86 1.134077e-86 [266,] 1.0000000 1.464461e-85 7.322305e-86 [267,] 1.0000000 9.324322e-85 4.662161e-85 [268,] 1.0000000 5.758235e-84 2.879117e-84 [269,] 1.0000000 3.526792e-83 1.763396e-83 [270,] 1.0000000 2.197237e-82 1.098618e-82 [271,] 1.0000000 1.337620e-81 6.688100e-82 [272,] 1.0000000 7.175639e-81 3.587819e-81 [273,] 1.0000000 4.282514e-80 2.141257e-80 [274,] 1.0000000 2.664198e-79 1.332099e-79 [275,] 1.0000000 1.593667e-78 7.968334e-79 [276,] 1.0000000 8.914938e-78 4.457469e-78 [277,] 1.0000000 5.281743e-77 2.640872e-77 [278,] 1.0000000 3.196466e-76 1.598233e-76 [279,] 1.0000000 1.875505e-75 9.377524e-76 [280,] 1.0000000 1.096519e-74 5.482597e-75 [281,] 1.0000000 6.395781e-74 3.197890e-74 [282,] 1.0000000 3.819060e-73 1.909530e-73 [283,] 1.0000000 2.241909e-72 1.120954e-72 [284,] 1.0000000 1.324254e-71 6.621270e-72 [285,] 1.0000000 7.362502e-71 3.681251e-71 [286,] 1.0000000 4.208816e-70 2.104408e-70 [287,] 1.0000000 2.371638e-69 1.185819e-69 [288,] 1.0000000 1.367691e-68 6.838456e-69 [289,] 1.0000000 5.971010e-68 2.985505e-68 [290,] 1.0000000 3.319440e-67 1.659720e-67 [291,] 1.0000000 1.885851e-66 9.429255e-67 [292,] 1.0000000 1.037447e-65 5.187236e-66 [293,] 1.0000000 5.698661e-65 2.849331e-65 [294,] 1.0000000 3.199901e-64 1.599951e-64 [295,] 1.0000000 1.773839e-63 8.869194e-64 [296,] 1.0000000 9.556411e-63 4.778206e-63 [297,] 1.0000000 5.247472e-62 2.623736e-62 [298,] 1.0000000 2.133624e-61 1.066812e-61 [299,] 1.0000000 1.136810e-60 5.684050e-61 [300,] 1.0000000 6.149827e-60 3.074914e-60 [301,] 1.0000000 3.339058e-59 1.669529e-59 [302,] 1.0000000 1.596429e-58 7.982144e-59 [303,] 1.0000000 8.488333e-58 4.244166e-58 [304,] 1.0000000 4.479978e-57 2.239989e-57 [305,] 1.0000000 2.100607e-56 1.050303e-56 [306,] 1.0000000 1.094907e-55 5.474537e-56 [307,] 1.0000000 5.705074e-55 2.852537e-55 [308,] 1.0000000 2.890970e-54 1.445485e-54 [309,] 1.0000000 1.365145e-53 6.825723e-54 [310,] 1.0000000 6.768636e-53 3.384318e-53 [311,] 1.0000000 3.465018e-52 1.732509e-52 [312,] 1.0000000 1.715208e-51 8.576040e-52 [313,] 1.0000000 8.423252e-51 4.211626e-51 [314,] 1.0000000 4.111810e-50 2.055905e-50 [315,] 1.0000000 1.995060e-49 9.975302e-50 [316,] 1.0000000 9.621166e-49 4.810583e-49 [317,] 1.0000000 4.681033e-48 2.340517e-48 [318,] 1.0000000 2.230061e-47 1.115031e-47 [319,] 1.0000000 1.066973e-46 5.334867e-47 [320,] 1.0000000 5.075600e-46 2.537800e-46 [321,] 1.0000000 2.385956e-45 1.192978e-45 [322,] 1.0000000 1.108746e-44 5.543729e-45 [323,] 1.0000000 5.119272e-44 2.559636e-44 [324,] 1.0000000 2.348372e-43 1.174186e-43 [325,] 1.0000000 1.070245e-42 5.351224e-43 [326,] 1.0000000 4.845424e-42 2.422712e-42 [327,] 1.0000000 2.183570e-41 1.091785e-41 [328,] 1.0000000 9.759100e-41 4.879550e-41 [329,] 1.0000000 4.293055e-40 2.146527e-40 [330,] 1.0000000 1.960476e-39 9.802379e-40 [331,] 1.0000000 8.581514e-39 4.290757e-39 [332,] 1.0000000 3.757470e-38 1.878735e-38 [333,] 1.0000000 1.584330e-37 7.921651e-38 [334,] 1.0000000 6.971298e-37 3.485649e-37 [335,] 1.0000000 2.969279e-36 1.484640e-36 [336,] 1.0000000 1.255578e-35 6.277890e-36 [337,] 1.0000000 5.270589e-35 2.635294e-35 [338,] 1.0000000 2.225750e-34 1.112875e-34 [339,] 1.0000000 9.205314e-34 4.602657e-34 [340,] 1.0000000 3.864565e-33 1.932283e-33 [341,] 1.0000000 1.574401e-32 7.872007e-33 [342,] 1.0000000 6.364734e-32 3.182367e-32 [343,] 1.0000000 2.313304e-31 1.156652e-31 [344,] 1.0000000 8.664499e-31 4.332249e-31 [345,] 1.0000000 3.458423e-30 1.729211e-30 [346,] 1.0000000 1.218564e-29 6.092819e-30 [347,] 1.0000000 4.808086e-29 2.404043e-29 [348,] 1.0000000 1.853821e-28 9.269104e-29 [349,] 1.0000000 7.087858e-28 3.543929e-28 [350,] 1.0000000 2.703424e-27 1.351712e-27 [351,] 1.0000000 1.018846e-26 5.094228e-27 [352,] 1.0000000 3.806428e-26 1.903214e-26 [353,] 1.0000000 1.409578e-25 7.047889e-26 [354,] 1.0000000 5.181652e-25 2.590826e-25 [355,] 1.0000000 1.884555e-24 9.422777e-25 [356,] 1.0000000 6.791133e-24 3.395566e-24 [357,] 1.0000000 2.417417e-23 1.208709e-23 [358,] 1.0000000 8.387977e-23 4.193988e-23 [359,] 1.0000000 2.942338e-22 1.471169e-22 [360,] 1.0000000 1.022085e-21 5.110426e-22 [361,] 1.0000000 3.515363e-21 1.757682e-21 [362,] 1.0000000 1.196928e-20 5.984640e-21 [363,] 1.0000000 4.013921e-20 2.006960e-20 [364,] 1.0000000 1.338684e-19 6.693420e-20 [365,] 1.0000000 4.103891e-19 2.051945e-19 [366,] 1.0000000 1.417499e-18 7.087495e-19 [367,] 1.0000000 4.830696e-18 2.415348e-18 [368,] 1.0000000 1.545519e-17 7.727594e-18 [369,] 1.0000000 4.888215e-17 2.444108e-17 [370,] 1.0000000 1.575527e-16 7.877633e-17 [371,] 1.0000000 4.574695e-16 2.287347e-16 [372,] 1.0000000 1.225756e-15 6.128782e-16 [373,] 1.0000000 3.696777e-15 1.848389e-15 [374,] 1.0000000 1.100606e-14 5.503029e-15 [375,] 1.0000000 3.499313e-14 1.749657e-14 [376,] 1.0000000 8.946392e-14 4.473196e-14 [377,] 1.0000000 2.565646e-13 1.282823e-13 [378,] 1.0000000 7.391409e-13 3.695705e-13 [379,] 1.0000000 1.996102e-12 9.980509e-13 [380,] 1.0000000 5.092436e-12 2.546218e-12 [381,] 1.0000000 1.350930e-11 6.754652e-12 [382,] 1.0000000 3.638123e-11 1.819062e-11 [383,] 1.0000000 9.534922e-11 4.767461e-11 [384,] 1.0000000 2.439596e-10 1.219798e-10 [385,] 1.0000000 6.685122e-10 3.342561e-10 [386,] 1.0000000 1.890270e-09 9.451348e-10 [387,] 1.0000000 3.651164e-09 1.825582e-09 [388,] 1.0000000 8.801424e-09 4.400712e-09 [389,] 1.0000000 2.319614e-08 1.159807e-08 [390,] 1.0000000 5.791891e-08 2.895945e-08 [391,] 0.9999999 1.200941e-07 6.004706e-08 [392,] 0.9999998 3.061282e-07 1.530641e-07 [393,] 0.9999997 6.854829e-07 3.427415e-07 [394,] 0.9999991 1.714987e-06 8.574935e-07 [395,] 0.9999982 3.530216e-06 1.765108e-06 [396,] 0.9999960 7.901662e-06 3.950831e-06 [397,] 0.9999908 1.848380e-05 9.241898e-06 [398,] 0.9999780 4.393486e-05 2.196743e-05 [399,] 0.9999532 9.353515e-05 4.676758e-05 [400,] 0.9999007 1.985446e-04 9.927231e-05 [401,] 0.9998100 3.799579e-04 1.899790e-04 [402,] 0.9996361 7.278703e-04 3.639352e-04 [403,] 0.9992154 1.569141e-03 7.845706e-04 [404,] 0.9983970 3.206047e-03 1.603024e-03 [405,] 0.9968513 6.297453e-03 3.148727e-03 [406,] 0.9936497 1.270069e-02 6.350347e-03 [407,] 0.9937102 1.257957e-02 6.289785e-03 [408,] 0.9899248 2.015036e-02 1.007518e-02 [409,] 0.9814504 3.709920e-02 1.854960e-02 [410,] 0.9948193 1.036136e-02 5.180679e-03 [411,] 0.9885426 2.291473e-02 1.145737e-02 [412,] 0.9748603 5.027941e-02 2.513971e-02 [413,] 0.9465972 1.068056e-01 5.340278e-02 [414,] 0.8937335 2.125330e-01 1.062665e-01 [415,] 0.7986671 4.026659e-01 2.013329e-01 [416,] 0.6517514 6.964971e-01 3.482486e-01 > postscript(file="/var/www/html/rcomp/tmp/1xpby1291221210.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/2xpby1291221210.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/38ga11291221210.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/48ga11291221210.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/58ga11291221210.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 6 1402961.479 2913138.769 1352753.355 -2169549.620 -410679.272 86519.825 7 8 9 10 11 12 -107543.117 115170.675 -360898.474 684728.075 -125278.630 359203.174 13 14 15 16 17 18 432172.621 -1120524.455 -233672.311 479383.913 -209978.041 -257330.290 19 20 21 22 23 24 -321220.597 -555562.003 171657.086 -169963.511 -97939.098 -608028.961 25 26 27 28 29 30 -299233.585 -70546.516 -494980.590 -327174.386 359384.866 -243129.762 31 32 33 34 35 36 252635.029 -118693.439 222270.011 98356.917 -31446.011 -544502.394 37 38 39 40 41 42 50507.376 128163.400 -76450.561 -285308.135 -433861.615 -313566.191 43 44 45 46 47 48 17670.150 -190813.141 -13896.029 -87156.361 -308099.912 -288734.002 49 50 51 52 53 54 -36458.009 -168107.866 6487.705 -111871.450 -253114.256 15295.099 55 56 57 58 59 60 -419913.622 29789.549 330892.738 -7327.732 -43800.516 -236281.701 61 62 63 64 65 66 -18233.594 -166473.087 -26687.129 -205557.027 -2007.474 -468248.697 67 68 69 70 71 72 -52372.158 -47243.728 -82202.671 -101959.036 -157977.521 13012.828 73 74 75 76 77 78 -73604.619 27328.585 -19117.089 -133538.174 -13321.181 -61192.121 79 80 81 82 83 84 -28017.638 -70984.910 39537.333 -81245.903 -66102.720 52250.348 85 86 87 88 89 90 38693.345 -135178.356 -11343.467 14859.058 -99555.689 46230.652 91 92 93 94 95 96 12662.905 -55225.344 84324.520 -204255.765 61879.176 -15521.425 97 98 99 100 101 102 85736.827 -32922.067 -25880.808 -37141.056 -70129.747 -99146.107 103 104 105 106 107 108 -149224.227 -36215.540 42831.684 -40273.045 -286278.469 97095.955 109 110 111 112 113 114 -159832.658 -127292.519 106914.799 -226382.812 -150798.315 -92889.550 115 116 117 118 119 120 49560.374 34212.597 29973.288 85593.041 -167520.659 26992.853 121 122 123 124 125 126 -270239.165 44593.971 -84197.513 51886.210 -21806.226 -197243.845 127 128 129 130 131 132 65929.058 53389.293 41658.487 65929.058 57632.350 69595.804 133 134 135 136 137 138 65929.058 43174.000 44601.864 55114.720 57065.590 64374.858 139 140 141 142 143 144 41884.819 65929.058 62498.832 65929.058 22596.029 61521.371 145 146 147 148 149 150 51941.111 54410.801 65929.058 58506.192 67922.218 65668.025 151 152 153 154 155 156 27222.295 -25738.300 61174.154 34741.970 63120.229 65929.058 157 158 159 160 161 162 65929.058 60553.705 73101.614 65929.058 65929.058 65929.058 163 164 165 166 167 168 65929.058 39346.008 16046.603 63444.017 65929.058 122277.892 169 170 171 172 173 174 -62798.915 65929.058 65929.058 58033.087 65929.058 45262.036 175 176 177 178 179 180 65929.058 55359.960 -52462.303 17565.671 65929.058 65929.058 181 182 183 184 185 186 53805.874 -5686.749 -5704.113 -40076.194 9750.736 12209.395 187 188 189 190 191 192 -23778.155 52884.959 24969.757 65929.058 58427.473 65929.058 193 194 195 196 197 198 54768.353 65929.058 62508.776 18372.740 -161948.935 18372.740 199 200 201 202 203 204 18372.740 -85234.859 18372.740 18372.740 18372.740 3123.603 205 206 207 208 209 210 -34367.008 12751.263 10733.070 18372.740 15358.754 1246.183 211 212 213 214 215 216 65929.058 65929.058 68937.951 17666.963 35622.966 64027.096 217 218 219 220 221 222 60031.984 52989.146 64027.096 11943.597 38322.337 65929.058 223 224 225 226 227 228 64919.348 33373.929 64604.144 37382.774 55196.557 65929.058 229 230 231 232 233 234 49836.275 65929.058 30753.308 47097.884 65929.058 31837.494 235 236 237 238 239 240 65929.058 65957.049 65929.058 60031.142 65929.058 65929.058 241 242 243 244 245 246 65929.058 57539.279 65929.058 60369.653 3466.936 18372.740 247 248 249 250 251 252 18372.740 40733.687 18372.740 18372.740 45705.516 18372.740 253 254 255 256 257 258 47886.726 -13578.909 -134807.851 65929.058 54364.204 36026.365 259 260 261 262 263 264 18372.740 65929.058 18372.740 18372.740 58739.116 69056.918 265 266 267 268 269 270 18372.740 65929.058 18372.740 93313.467 18372.740 40389.473 271 272 273 274 275 276 18372.740 18372.740 18372.740 39807.516 65929.058 50307.221 277 278 279 280 281 282 36444.067 65929.058 -72133.754 38471.735 18372.740 65929.058 283 284 285 286 287 288 39889.761 65929.058 50651.402 65929.058 4199.495 63121.044 289 290 291 292 293 294 18372.740 35164.002 18372.740 57455.728 57472.520 65929.058 295 296 297 298 299 300 44123.151 29038.476 65929.058 18372.740 65929.058 57370.231 301 302 303 304 305 306 11573.074 47299.698 64978.077 18372.740 71940.367 54464.426 307 308 309 310 311 312 18372.740 -8040.394 44919.571 18372.740 18372.740 18490.060 313 314 315 316 317 318 18372.740 5487.357 58468.829 75256.829 65919.056 19879.629 319 320 321 322 323 324 61145.628 65929.058 65929.058 65929.058 65929.058 18372.740 325 326 327 328 329 330 65929.058 61612.076 18372.740 18372.740 65929.058 65929.058 331 332 333 334 335 336 65929.058 65929.058 65929.058 58875.669 64953.405 55194.372 337 338 339 340 341 342 8752.196 65929.058 62637.715 57478.821 24601.292 65929.058 343 344 345 346 347 348 65929.058 65929.058 32753.089 65929.058 51182.227 65929.058 349 350 351 352 353 354 65929.058 25248.341 62730.287 56005.167 55506.617 58790.037 355 356 357 358 359 360 65929.058 65929.058 19911.198 65929.058 65929.058 65929.058 361 362 363 364 365 366 62118.054 65929.058 65929.058 -55498.104 56055.704 65929.058 367 368 369 370 371 372 65929.058 65929.058 65929.058 54639.637 65929.058 58803.246 373 374 375 376 377 378 51349.727 45059.431 65929.058 65929.058 39687.115 53282.786 379 380 381 382 383 384 12142.297 65929.058 65929.058 19736.941 34263.075 65929.058 385 386 387 388 389 390 51741.622 -104527.163 73487.260 36723.162 60807.854 41511.317 391 392 393 394 395 396 51593.419 29177.487 -7843.262 86055.384 -131925.428 17577.025 397 398 399 400 401 402 -126907.668 -209372.064 23709.340 -91936.031 13839.628 54010.943 403 404 405 406 407 408 30616.402 -84731.246 2350.256 25895.871 32384.186 24608.709 409 410 411 412 413 414 34074.969 -50506.225 -36839.836 -197406.940 -62015.747 25808.647 415 416 417 418 419 420 -165224.220 -87688.976 102732.060 -52946.773 -37843.496 -26327.090 421 422 423 424 425 426 -33004.999 -159782.676 -584784.875 -447047.074 -37244.261 -281840.274 427 428 429 430 431 -300445.488 -1107109.319 -258823.827 -319850.175 -669870.658 > postscript(file="/var/www/html/rcomp/tmp/61p941291221210.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 1402961.479 NA 1 2913138.769 1402961.479 2 1352753.355 2913138.769 3 -2169549.620 1352753.355 4 -410679.272 -2169549.620 5 86519.825 -410679.272 6 -107543.117 86519.825 7 115170.675 -107543.117 8 -360898.474 115170.675 9 684728.075 -360898.474 10 -125278.630 684728.075 11 359203.174 -125278.630 12 432172.621 359203.174 13 -1120524.455 432172.621 14 -233672.311 -1120524.455 15 479383.913 -233672.311 16 -209978.041 479383.913 17 -257330.290 -209978.041 18 -321220.597 -257330.290 19 -555562.003 -321220.597 20 171657.086 -555562.003 21 -169963.511 171657.086 22 -97939.098 -169963.511 23 -608028.961 -97939.098 24 -299233.585 -608028.961 25 -70546.516 -299233.585 26 -494980.590 -70546.516 27 -327174.386 -494980.590 28 359384.866 -327174.386 29 -243129.762 359384.866 30 252635.029 -243129.762 31 -118693.439 252635.029 32 222270.011 -118693.439 33 98356.917 222270.011 34 -31446.011 98356.917 35 -544502.394 -31446.011 36 50507.376 -544502.394 37 128163.400 50507.376 38 -76450.561 128163.400 39 -285308.135 -76450.561 40 -433861.615 -285308.135 41 -313566.191 -433861.615 42 17670.150 -313566.191 43 -190813.141 17670.150 44 -13896.029 -190813.141 45 -87156.361 -13896.029 46 -308099.912 -87156.361 47 -288734.002 -308099.912 48 -36458.009 -288734.002 49 -168107.866 -36458.009 50 6487.705 -168107.866 51 -111871.450 6487.705 52 -253114.256 -111871.450 53 15295.099 -253114.256 54 -419913.622 15295.099 55 29789.549 -419913.622 56 330892.738 29789.549 57 -7327.732 330892.738 58 -43800.516 -7327.732 59 -236281.701 -43800.516 60 -18233.594 -236281.701 61 -166473.087 -18233.594 62 -26687.129 -166473.087 63 -205557.027 -26687.129 64 -2007.474 -205557.027 65 -468248.697 -2007.474 66 -52372.158 -468248.697 67 -47243.728 -52372.158 68 -82202.671 -47243.728 69 -101959.036 -82202.671 70 -157977.521 -101959.036 71 13012.828 -157977.521 72 -73604.619 13012.828 73 27328.585 -73604.619 74 -19117.089 27328.585 75 -133538.174 -19117.089 76 -13321.181 -133538.174 77 -61192.121 -13321.181 78 -28017.638 -61192.121 79 -70984.910 -28017.638 80 39537.333 -70984.910 81 -81245.903 39537.333 82 -66102.720 -81245.903 83 52250.348 -66102.720 84 38693.345 52250.348 85 -135178.356 38693.345 86 -11343.467 -135178.356 87 14859.058 -11343.467 88 -99555.689 14859.058 89 46230.652 -99555.689 90 12662.905 46230.652 91 -55225.344 12662.905 92 84324.520 -55225.344 93 -204255.765 84324.520 94 61879.176 -204255.765 95 -15521.425 61879.176 96 85736.827 -15521.425 97 -32922.067 85736.827 98 -25880.808 -32922.067 99 -37141.056 -25880.808 100 -70129.747 -37141.056 101 -99146.107 -70129.747 102 -149224.227 -99146.107 103 -36215.540 -149224.227 104 42831.684 -36215.540 105 -40273.045 42831.684 106 -286278.469 -40273.045 107 97095.955 -286278.469 108 -159832.658 97095.955 109 -127292.519 -159832.658 110 106914.799 -127292.519 111 -226382.812 106914.799 112 -150798.315 -226382.812 113 -92889.550 -150798.315 114 49560.374 -92889.550 115 34212.597 49560.374 116 29973.288 34212.597 117 85593.041 29973.288 118 -167520.659 85593.041 119 26992.853 -167520.659 120 -270239.165 26992.853 121 44593.971 -270239.165 122 -84197.513 44593.971 123 51886.210 -84197.513 124 -21806.226 51886.210 125 -197243.845 -21806.226 126 65929.058 -197243.845 127 53389.293 65929.058 128 41658.487 53389.293 129 65929.058 41658.487 130 57632.350 65929.058 131 69595.804 57632.350 132 65929.058 69595.804 133 43174.000 65929.058 134 44601.864 43174.000 135 55114.720 44601.864 136 57065.590 55114.720 137 64374.858 57065.590 138 41884.819 64374.858 139 65929.058 41884.819 140 62498.832 65929.058 141 65929.058 62498.832 142 22596.029 65929.058 143 61521.371 22596.029 144 51941.111 61521.371 145 54410.801 51941.111 146 65929.058 54410.801 147 58506.192 65929.058 148 67922.218 58506.192 149 65668.025 67922.218 150 27222.295 65668.025 151 -25738.300 27222.295 152 61174.154 -25738.300 153 34741.970 61174.154 154 63120.229 34741.970 155 65929.058 63120.229 156 65929.058 65929.058 157 60553.705 65929.058 158 73101.614 60553.705 159 65929.058 73101.614 160 65929.058 65929.058 161 65929.058 65929.058 162 65929.058 65929.058 163 39346.008 65929.058 164 16046.603 39346.008 165 63444.017 16046.603 166 65929.058 63444.017 167 122277.892 65929.058 168 -62798.915 122277.892 169 65929.058 -62798.915 170 65929.058 65929.058 171 58033.087 65929.058 172 65929.058 58033.087 173 45262.036 65929.058 174 65929.058 45262.036 175 55359.960 65929.058 176 -52462.303 55359.960 177 17565.671 -52462.303 178 65929.058 17565.671 179 65929.058 65929.058 180 53805.874 65929.058 181 -5686.749 53805.874 182 -5704.113 -5686.749 183 -40076.194 -5704.113 184 9750.736 -40076.194 185 12209.395 9750.736 186 -23778.155 12209.395 187 52884.959 -23778.155 188 24969.757 52884.959 189 65929.058 24969.757 190 58427.473 65929.058 191 65929.058 58427.473 192 54768.353 65929.058 193 65929.058 54768.353 194 62508.776 65929.058 195 18372.740 62508.776 196 -161948.935 18372.740 197 18372.740 -161948.935 198 18372.740 18372.740 199 -85234.859 18372.740 200 18372.740 -85234.859 201 18372.740 18372.740 202 18372.740 18372.740 203 3123.603 18372.740 204 -34367.008 3123.603 205 12751.263 -34367.008 206 10733.070 12751.263 207 18372.740 10733.070 208 15358.754 18372.740 209 1246.183 15358.754 210 65929.058 1246.183 211 65929.058 65929.058 212 68937.951 65929.058 213 17666.963 68937.951 214 35622.966 17666.963 215 64027.096 35622.966 216 60031.984 64027.096 217 52989.146 60031.984 218 64027.096 52989.146 219 11943.597 64027.096 220 38322.337 11943.597 221 65929.058 38322.337 222 64919.348 65929.058 223 33373.929 64919.348 224 64604.144 33373.929 225 37382.774 64604.144 226 55196.557 37382.774 227 65929.058 55196.557 228 49836.275 65929.058 229 65929.058 49836.275 230 30753.308 65929.058 231 47097.884 30753.308 232 65929.058 47097.884 233 31837.494 65929.058 234 65929.058 31837.494 235 65957.049 65929.058 236 65929.058 65957.049 237 60031.142 65929.058 238 65929.058 60031.142 239 65929.058 65929.058 240 65929.058 65929.058 241 57539.279 65929.058 242 65929.058 57539.279 243 60369.653 65929.058 244 3466.936 60369.653 245 18372.740 3466.936 246 18372.740 18372.740 247 40733.687 18372.740 248 18372.740 40733.687 249 18372.740 18372.740 250 45705.516 18372.740 251 18372.740 45705.516 252 47886.726 18372.740 253 -13578.909 47886.726 254 -134807.851 -13578.909 255 65929.058 -134807.851 256 54364.204 65929.058 257 36026.365 54364.204 258 18372.740 36026.365 259 65929.058 18372.740 260 18372.740 65929.058 261 18372.740 18372.740 262 58739.116 18372.740 263 69056.918 58739.116 264 18372.740 69056.918 265 65929.058 18372.740 266 18372.740 65929.058 267 93313.467 18372.740 268 18372.740 93313.467 269 40389.473 18372.740 270 18372.740 40389.473 271 18372.740 18372.740 272 18372.740 18372.740 273 39807.516 18372.740 274 65929.058 39807.516 275 50307.221 65929.058 276 36444.067 50307.221 277 65929.058 36444.067 278 -72133.754 65929.058 279 38471.735 -72133.754 280 18372.740 38471.735 281 65929.058 18372.740 282 39889.761 65929.058 283 65929.058 39889.761 284 50651.402 65929.058 285 65929.058 50651.402 286 4199.495 65929.058 287 63121.044 4199.495 288 18372.740 63121.044 289 35164.002 18372.740 290 18372.740 35164.002 291 57455.728 18372.740 292 57472.520 57455.728 293 65929.058 57472.520 294 44123.151 65929.058 295 29038.476 44123.151 296 65929.058 29038.476 297 18372.740 65929.058 298 65929.058 18372.740 299 57370.231 65929.058 300 11573.074 57370.231 301 47299.698 11573.074 302 64978.077 47299.698 303 18372.740 64978.077 304 71940.367 18372.740 305 54464.426 71940.367 306 18372.740 54464.426 307 -8040.394 18372.740 308 44919.571 -8040.394 309 18372.740 44919.571 310 18372.740 18372.740 311 18490.060 18372.740 312 18372.740 18490.060 313 5487.357 18372.740 314 58468.829 5487.357 315 75256.829 58468.829 316 65919.056 75256.829 317 19879.629 65919.056 318 61145.628 19879.629 319 65929.058 61145.628 320 65929.058 65929.058 321 65929.058 65929.058 322 65929.058 65929.058 323 18372.740 65929.058 324 65929.058 18372.740 325 61612.076 65929.058 326 18372.740 61612.076 327 18372.740 18372.740 328 65929.058 18372.740 329 65929.058 65929.058 330 65929.058 65929.058 331 65929.058 65929.058 332 65929.058 65929.058 333 58875.669 65929.058 334 64953.405 58875.669 335 55194.372 64953.405 336 8752.196 55194.372 337 65929.058 8752.196 338 62637.715 65929.058 339 57478.821 62637.715 340 24601.292 57478.821 341 65929.058 24601.292 342 65929.058 65929.058 343 65929.058 65929.058 344 32753.089 65929.058 345 65929.058 32753.089 346 51182.227 65929.058 347 65929.058 51182.227 348 65929.058 65929.058 349 25248.341 65929.058 350 62730.287 25248.341 351 56005.167 62730.287 352 55506.617 56005.167 353 58790.037 55506.617 354 65929.058 58790.037 355 65929.058 65929.058 356 19911.198 65929.058 357 65929.058 19911.198 358 65929.058 65929.058 359 65929.058 65929.058 360 62118.054 65929.058 361 65929.058 62118.054 362 65929.058 65929.058 363 -55498.104 65929.058 364 56055.704 -55498.104 365 65929.058 56055.704 366 65929.058 65929.058 367 65929.058 65929.058 368 65929.058 65929.058 369 54639.637 65929.058 370 65929.058 54639.637 371 58803.246 65929.058 372 51349.727 58803.246 373 45059.431 51349.727 374 65929.058 45059.431 375 65929.058 65929.058 376 39687.115 65929.058 377 53282.786 39687.115 378 12142.297 53282.786 379 65929.058 12142.297 380 65929.058 65929.058 381 19736.941 65929.058 382 34263.075 19736.941 383 65929.058 34263.075 384 51741.622 65929.058 385 -104527.163 51741.622 386 73487.260 -104527.163 387 36723.162 73487.260 388 60807.854 36723.162 389 41511.317 60807.854 390 51593.419 41511.317 391 29177.487 51593.419 392 -7843.262 29177.487 393 86055.384 -7843.262 394 -131925.428 86055.384 395 17577.025 -131925.428 396 -126907.668 17577.025 397 -209372.064 -126907.668 398 23709.340 -209372.064 399 -91936.031 23709.340 400 13839.628 -91936.031 401 54010.943 13839.628 402 30616.402 54010.943 403 -84731.246 30616.402 404 2350.256 -84731.246 405 25895.871 2350.256 406 32384.186 25895.871 407 24608.709 32384.186 408 34074.969 24608.709 409 -50506.225 34074.969 410 -36839.836 -50506.225 411 -197406.940 -36839.836 412 -62015.747 -197406.940 413 25808.647 -62015.747 414 -165224.220 25808.647 415 -87688.976 -165224.220 416 102732.060 -87688.976 417 -52946.773 102732.060 418 -37843.496 -52946.773 419 -26327.090 -37843.496 420 -33004.999 -26327.090 421 -159782.676 -33004.999 422 -584784.875 -159782.676 423 -447047.074 -584784.875 424 -37244.261 -447047.074 425 -281840.274 -37244.261 426 -300445.488 -281840.274 427 -1107109.319 -300445.488 428 -258823.827 -1107109.319 429 -319850.175 -258823.827 430 -669870.658 -319850.175 431 NA -669870.658 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 2913138.769 1402961.479 [2,] 1352753.355 2913138.769 [3,] -2169549.620 1352753.355 [4,] -410679.272 -2169549.620 [5,] 86519.825 -410679.272 [6,] -107543.117 86519.825 [7,] 115170.675 -107543.117 [8,] -360898.474 115170.675 [9,] 684728.075 -360898.474 [10,] -125278.630 684728.075 [11,] 359203.174 -125278.630 [12,] 432172.621 359203.174 [13,] -1120524.455 432172.621 [14,] -233672.311 -1120524.455 [15,] 479383.913 -233672.311 [16,] -209978.041 479383.913 [17,] -257330.290 -209978.041 [18,] -321220.597 -257330.290 [19,] -555562.003 -321220.597 [20,] 171657.086 -555562.003 [21,] -169963.511 171657.086 [22,] -97939.098 -169963.511 [23,] -608028.961 -97939.098 [24,] -299233.585 -608028.961 [25,] -70546.516 -299233.585 [26,] -494980.590 -70546.516 [27,] -327174.386 -494980.590 [28,] 359384.866 -327174.386 [29,] -243129.762 359384.866 [30,] 252635.029 -243129.762 [31,] -118693.439 252635.029 [32,] 222270.011 -118693.439 [33,] 98356.917 222270.011 [34,] -31446.011 98356.917 [35,] -544502.394 -31446.011 [36,] 50507.376 -544502.394 [37,] 128163.400 50507.376 [38,] -76450.561 128163.400 [39,] -285308.135 -76450.561 [40,] -433861.615 -285308.135 [41,] -313566.191 -433861.615 [42,] 17670.150 -313566.191 [43,] -190813.141 17670.150 [44,] -13896.029 -190813.141 [45,] -87156.361 -13896.029 [46,] -308099.912 -87156.361 [47,] -288734.002 -308099.912 [48,] -36458.009 -288734.002 [49,] -168107.866 -36458.009 [50,] 6487.705 -168107.866 [51,] -111871.450 6487.705 [52,] -253114.256 -111871.450 [53,] 15295.099 -253114.256 [54,] -419913.622 15295.099 [55,] 29789.549 -419913.622 [56,] 330892.738 29789.549 [57,] -7327.732 330892.738 [58,] -43800.516 -7327.732 [59,] -236281.701 -43800.516 [60,] -18233.594 -236281.701 [61,] -166473.087 -18233.594 [62,] -26687.129 -166473.087 [63,] -205557.027 -26687.129 [64,] -2007.474 -205557.027 [65,] -468248.697 -2007.474 [66,] -52372.158 -468248.697 [67,] -47243.728 -52372.158 [68,] -82202.671 -47243.728 [69,] -101959.036 -82202.671 [70,] -157977.521 -101959.036 [71,] 13012.828 -157977.521 [72,] -73604.619 13012.828 [73,] 27328.585 -73604.619 [74,] -19117.089 27328.585 [75,] -133538.174 -19117.089 [76,] -13321.181 -133538.174 [77,] -61192.121 -13321.181 [78,] -28017.638 -61192.121 [79,] -70984.910 -28017.638 [80,] 39537.333 -70984.910 [81,] -81245.903 39537.333 [82,] -66102.720 -81245.903 [83,] 52250.348 -66102.720 [84,] 38693.345 52250.348 [85,] -135178.356 38693.345 [86,] -11343.467 -135178.356 [87,] 14859.058 -11343.467 [88,] -99555.689 14859.058 [89,] 46230.652 -99555.689 [90,] 12662.905 46230.652 [91,] -55225.344 12662.905 [92,] 84324.520 -55225.344 [93,] -204255.765 84324.520 [94,] 61879.176 -204255.765 [95,] -15521.425 61879.176 [96,] 85736.827 -15521.425 [97,] -32922.067 85736.827 [98,] -25880.808 -32922.067 [99,] -37141.056 -25880.808 [100,] -70129.747 -37141.056 [101,] -99146.107 -70129.747 [102,] -149224.227 -99146.107 [103,] -36215.540 -149224.227 [104,] 42831.684 -36215.540 [105,] -40273.045 42831.684 [106,] -286278.469 -40273.045 [107,] 97095.955 -286278.469 [108,] -159832.658 97095.955 [109,] -127292.519 -159832.658 [110,] 106914.799 -127292.519 [111,] -226382.812 106914.799 [112,] -150798.315 -226382.812 [113,] -92889.550 -150798.315 [114,] 49560.374 -92889.550 [115,] 34212.597 49560.374 [116,] 29973.288 34212.597 [117,] 85593.041 29973.288 [118,] -167520.659 85593.041 [119,] 26992.853 -167520.659 [120,] -270239.165 26992.853 [121,] 44593.971 -270239.165 [122,] -84197.513 44593.971 [123,] 51886.210 -84197.513 [124,] -21806.226 51886.210 [125,] -197243.845 -21806.226 [126,] 65929.058 -197243.845 [127,] 53389.293 65929.058 [128,] 41658.487 53389.293 [129,] 65929.058 41658.487 [130,] 57632.350 65929.058 [131,] 69595.804 57632.350 [132,] 65929.058 69595.804 [133,] 43174.000 65929.058 [134,] 44601.864 43174.000 [135,] 55114.720 44601.864 [136,] 57065.590 55114.720 [137,] 64374.858 57065.590 [138,] 41884.819 64374.858 [139,] 65929.058 41884.819 [140,] 62498.832 65929.058 [141,] 65929.058 62498.832 [142,] 22596.029 65929.058 [143,] 61521.371 22596.029 [144,] 51941.111 61521.371 [145,] 54410.801 51941.111 [146,] 65929.058 54410.801 [147,] 58506.192 65929.058 [148,] 67922.218 58506.192 [149,] 65668.025 67922.218 [150,] 27222.295 65668.025 [151,] -25738.300 27222.295 [152,] 61174.154 -25738.300 [153,] 34741.970 61174.154 [154,] 63120.229 34741.970 [155,] 65929.058 63120.229 [156,] 65929.058 65929.058 [157,] 60553.705 65929.058 [158,] 73101.614 60553.705 [159,] 65929.058 73101.614 [160,] 65929.058 65929.058 [161,] 65929.058 65929.058 [162,] 65929.058 65929.058 [163,] 39346.008 65929.058 [164,] 16046.603 39346.008 [165,] 63444.017 16046.603 [166,] 65929.058 63444.017 [167,] 122277.892 65929.058 [168,] -62798.915 122277.892 [169,] 65929.058 -62798.915 [170,] 65929.058 65929.058 [171,] 58033.087 65929.058 [172,] 65929.058 58033.087 [173,] 45262.036 65929.058 [174,] 65929.058 45262.036 [175,] 55359.960 65929.058 [176,] -52462.303 55359.960 [177,] 17565.671 -52462.303 [178,] 65929.058 17565.671 [179,] 65929.058 65929.058 [180,] 53805.874 65929.058 [181,] -5686.749 53805.874 [182,] -5704.113 -5686.749 [183,] -40076.194 -5704.113 [184,] 9750.736 -40076.194 [185,] 12209.395 9750.736 [186,] -23778.155 12209.395 [187,] 52884.959 -23778.155 [188,] 24969.757 52884.959 [189,] 65929.058 24969.757 [190,] 58427.473 65929.058 [191,] 65929.058 58427.473 [192,] 54768.353 65929.058 [193,] 65929.058 54768.353 [194,] 62508.776 65929.058 [195,] 18372.740 62508.776 [196,] -161948.935 18372.740 [197,] 18372.740 -161948.935 [198,] 18372.740 18372.740 [199,] -85234.859 18372.740 [200,] 18372.740 -85234.859 [201,] 18372.740 18372.740 [202,] 18372.740 18372.740 [203,] 3123.603 18372.740 [204,] -34367.008 3123.603 [205,] 12751.263 -34367.008 [206,] 10733.070 12751.263 [207,] 18372.740 10733.070 [208,] 15358.754 18372.740 [209,] 1246.183 15358.754 [210,] 65929.058 1246.183 [211,] 65929.058 65929.058 [212,] 68937.951 65929.058 [213,] 17666.963 68937.951 [214,] 35622.966 17666.963 [215,] 64027.096 35622.966 [216,] 60031.984 64027.096 [217,] 52989.146 60031.984 [218,] 64027.096 52989.146 [219,] 11943.597 64027.096 [220,] 38322.337 11943.597 [221,] 65929.058 38322.337 [222,] 64919.348 65929.058 [223,] 33373.929 64919.348 [224,] 64604.144 33373.929 [225,] 37382.774 64604.144 [226,] 55196.557 37382.774 [227,] 65929.058 55196.557 [228,] 49836.275 65929.058 [229,] 65929.058 49836.275 [230,] 30753.308 65929.058 [231,] 47097.884 30753.308 [232,] 65929.058 47097.884 [233,] 31837.494 65929.058 [234,] 65929.058 31837.494 [235,] 65957.049 65929.058 [236,] 65929.058 65957.049 [237,] 60031.142 65929.058 [238,] 65929.058 60031.142 [239,] 65929.058 65929.058 [240,] 65929.058 65929.058 [241,] 57539.279 65929.058 [242,] 65929.058 57539.279 [243,] 60369.653 65929.058 [244,] 3466.936 60369.653 [245,] 18372.740 3466.936 [246,] 18372.740 18372.740 [247,] 40733.687 18372.740 [248,] 18372.740 40733.687 [249,] 18372.740 18372.740 [250,] 45705.516 18372.740 [251,] 18372.740 45705.516 [252,] 47886.726 18372.740 [253,] -13578.909 47886.726 [254,] -134807.851 -13578.909 [255,] 65929.058 -134807.851 [256,] 54364.204 65929.058 [257,] 36026.365 54364.204 [258,] 18372.740 36026.365 [259,] 65929.058 18372.740 [260,] 18372.740 65929.058 [261,] 18372.740 18372.740 [262,] 58739.116 18372.740 [263,] 69056.918 58739.116 [264,] 18372.740 69056.918 [265,] 65929.058 18372.740 [266,] 18372.740 65929.058 [267,] 93313.467 18372.740 [268,] 18372.740 93313.467 [269,] 40389.473 18372.740 [270,] 18372.740 40389.473 [271,] 18372.740 18372.740 [272,] 18372.740 18372.740 [273,] 39807.516 18372.740 [274,] 65929.058 39807.516 [275,] 50307.221 65929.058 [276,] 36444.067 50307.221 [277,] 65929.058 36444.067 [278,] -72133.754 65929.058 [279,] 38471.735 -72133.754 [280,] 18372.740 38471.735 [281,] 65929.058 18372.740 [282,] 39889.761 65929.058 [283,] 65929.058 39889.761 [284,] 50651.402 65929.058 [285,] 65929.058 50651.402 [286,] 4199.495 65929.058 [287,] 63121.044 4199.495 [288,] 18372.740 63121.044 [289,] 35164.002 18372.740 [290,] 18372.740 35164.002 [291,] 57455.728 18372.740 [292,] 57472.520 57455.728 [293,] 65929.058 57472.520 [294,] 44123.151 65929.058 [295,] 29038.476 44123.151 [296,] 65929.058 29038.476 [297,] 18372.740 65929.058 [298,] 65929.058 18372.740 [299,] 57370.231 65929.058 [300,] 11573.074 57370.231 [301,] 47299.698 11573.074 [302,] 64978.077 47299.698 [303,] 18372.740 64978.077 [304,] 71940.367 18372.740 [305,] 54464.426 71940.367 [306,] 18372.740 54464.426 [307,] -8040.394 18372.740 [308,] 44919.571 -8040.394 [309,] 18372.740 44919.571 [310,] 18372.740 18372.740 [311,] 18490.060 18372.740 [312,] 18372.740 18490.060 [313,] 5487.357 18372.740 [314,] 58468.829 5487.357 [315,] 75256.829 58468.829 [316,] 65919.056 75256.829 [317,] 19879.629 65919.056 [318,] 61145.628 19879.629 [319,] 65929.058 61145.628 [320,] 65929.058 65929.058 [321,] 65929.058 65929.058 [322,] 65929.058 65929.058 [323,] 18372.740 65929.058 [324,] 65929.058 18372.740 [325,] 61612.076 65929.058 [326,] 18372.740 61612.076 [327,] 18372.740 18372.740 [328,] 65929.058 18372.740 [329,] 65929.058 65929.058 [330,] 65929.058 65929.058 [331,] 65929.058 65929.058 [332,] 65929.058 65929.058 [333,] 58875.669 65929.058 [334,] 64953.405 58875.669 [335,] 55194.372 64953.405 [336,] 8752.196 55194.372 [337,] 65929.058 8752.196 [338,] 62637.715 65929.058 [339,] 57478.821 62637.715 [340,] 24601.292 57478.821 [341,] 65929.058 24601.292 [342,] 65929.058 65929.058 [343,] 65929.058 65929.058 [344,] 32753.089 65929.058 [345,] 65929.058 32753.089 [346,] 51182.227 65929.058 [347,] 65929.058 51182.227 [348,] 65929.058 65929.058 [349,] 25248.341 65929.058 [350,] 62730.287 25248.341 [351,] 56005.167 62730.287 [352,] 55506.617 56005.167 [353,] 58790.037 55506.617 [354,] 65929.058 58790.037 [355,] 65929.058 65929.058 [356,] 19911.198 65929.058 [357,] 65929.058 19911.198 [358,] 65929.058 65929.058 [359,] 65929.058 65929.058 [360,] 62118.054 65929.058 [361,] 65929.058 62118.054 [362,] 65929.058 65929.058 [363,] -55498.104 65929.058 [364,] 56055.704 -55498.104 [365,] 65929.058 56055.704 [366,] 65929.058 65929.058 [367,] 65929.058 65929.058 [368,] 65929.058 65929.058 [369,] 54639.637 65929.058 [370,] 65929.058 54639.637 [371,] 58803.246 65929.058 [372,] 51349.727 58803.246 [373,] 45059.431 51349.727 [374,] 65929.058 45059.431 [375,] 65929.058 65929.058 [376,] 39687.115 65929.058 [377,] 53282.786 39687.115 [378,] 12142.297 53282.786 [379,] 65929.058 12142.297 [380,] 65929.058 65929.058 [381,] 19736.941 65929.058 [382,] 34263.075 19736.941 [383,] 65929.058 34263.075 [384,] 51741.622 65929.058 [385,] -104527.163 51741.622 [386,] 73487.260 -104527.163 [387,] 36723.162 73487.260 [388,] 60807.854 36723.162 [389,] 41511.317 60807.854 [390,] 51593.419 41511.317 [391,] 29177.487 51593.419 [392,] -7843.262 29177.487 [393,] 86055.384 -7843.262 [394,] -131925.428 86055.384 [395,] 17577.025 -131925.428 [396,] -126907.668 17577.025 [397,] -209372.064 -126907.668 [398,] 23709.340 -209372.064 [399,] -91936.031 23709.340 [400,] 13839.628 -91936.031 [401,] 54010.943 13839.628 [402,] 30616.402 54010.943 [403,] -84731.246 30616.402 [404,] 2350.256 -84731.246 [405,] 25895.871 2350.256 [406,] 32384.186 25895.871 [407,] 24608.709 32384.186 [408,] 34074.969 24608.709 [409,] -50506.225 34074.969 [410,] -36839.836 -50506.225 [411,] -197406.940 -36839.836 [412,] -62015.747 -197406.940 [413,] 25808.647 -62015.747 [414,] -165224.220 25808.647 [415,] -87688.976 -165224.220 [416,] 102732.060 -87688.976 [417,] -52946.773 102732.060 [418,] -37843.496 -52946.773 [419,] -26327.090 -37843.496 [420,] -33004.999 -26327.090 [421,] -159782.676 -33004.999 [422,] -584784.875 -159782.676 [423,] -447047.074 -584784.875 [424,] -37244.261 -447047.074 [425,] -281840.274 -37244.261 [426,] -300445.488 -281840.274 [427,] -1107109.319 -300445.488 [428,] -258823.827 -1107109.319 [429,] -319850.175 -258823.827 [430,] -669870.658 -319850.175 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 2913138.769 1402961.479 2 1352753.355 2913138.769 3 -2169549.620 1352753.355 4 -410679.272 -2169549.620 5 86519.825 -410679.272 6 -107543.117 86519.825 7 115170.675 -107543.117 8 -360898.474 115170.675 9 684728.075 -360898.474 10 -125278.630 684728.075 11 359203.174 -125278.630 12 432172.621 359203.174 13 -1120524.455 432172.621 14 -233672.311 -1120524.455 15 479383.913 -233672.311 16 -209978.041 479383.913 17 -257330.290 -209978.041 18 -321220.597 -257330.290 19 -555562.003 -321220.597 20 171657.086 -555562.003 21 -169963.511 171657.086 22 -97939.098 -169963.511 23 -608028.961 -97939.098 24 -299233.585 -608028.961 25 -70546.516 -299233.585 26 -494980.590 -70546.516 27 -327174.386 -494980.590 28 359384.866 -327174.386 29 -243129.762 359384.866 30 252635.029 -243129.762 31 -118693.439 252635.029 32 222270.011 -118693.439 33 98356.917 222270.011 34 -31446.011 98356.917 35 -544502.394 -31446.011 36 50507.376 -544502.394 37 128163.400 50507.376 38 -76450.561 128163.400 39 -285308.135 -76450.561 40 -433861.615 -285308.135 41 -313566.191 -433861.615 42 17670.150 -313566.191 43 -190813.141 17670.150 44 -13896.029 -190813.141 45 -87156.361 -13896.029 46 -308099.912 -87156.361 47 -288734.002 -308099.912 48 -36458.009 -288734.002 49 -168107.866 -36458.009 50 6487.705 -168107.866 51 -111871.450 6487.705 52 -253114.256 -111871.450 53 15295.099 -253114.256 54 -419913.622 15295.099 55 29789.549 -419913.622 56 330892.738 29789.549 57 -7327.732 330892.738 58 -43800.516 -7327.732 59 -236281.701 -43800.516 60 -18233.594 -236281.701 61 -166473.087 -18233.594 62 -26687.129 -166473.087 63 -205557.027 -26687.129 64 -2007.474 -205557.027 65 -468248.697 -2007.474 66 -52372.158 -468248.697 67 -47243.728 -52372.158 68 -82202.671 -47243.728 69 -101959.036 -82202.671 70 -157977.521 -101959.036 71 13012.828 -157977.521 72 -73604.619 13012.828 73 27328.585 -73604.619 74 -19117.089 27328.585 75 -133538.174 -19117.089 76 -13321.181 -133538.174 77 -61192.121 -13321.181 78 -28017.638 -61192.121 79 -70984.910 -28017.638 80 39537.333 -70984.910 81 -81245.903 39537.333 82 -66102.720 -81245.903 83 52250.348 -66102.720 84 38693.345 52250.348 85 -135178.356 38693.345 86 -11343.467 -135178.356 87 14859.058 -11343.467 88 -99555.689 14859.058 89 46230.652 -99555.689 90 12662.905 46230.652 91 -55225.344 12662.905 92 84324.520 -55225.344 93 -204255.765 84324.520 94 61879.176 -204255.765 95 -15521.425 61879.176 96 85736.827 -15521.425 97 -32922.067 85736.827 98 -25880.808 -32922.067 99 -37141.056 -25880.808 100 -70129.747 -37141.056 101 -99146.107 -70129.747 102 -149224.227 -99146.107 103 -36215.540 -149224.227 104 42831.684 -36215.540 105 -40273.045 42831.684 106 -286278.469 -40273.045 107 97095.955 -286278.469 108 -159832.658 97095.955 109 -127292.519 -159832.658 110 106914.799 -127292.519 111 -226382.812 106914.799 112 -150798.315 -226382.812 113 -92889.550 -150798.315 114 49560.374 -92889.550 115 34212.597 49560.374 116 29973.288 34212.597 117 85593.041 29973.288 118 -167520.659 85593.041 119 26992.853 -167520.659 120 -270239.165 26992.853 121 44593.971 -270239.165 122 -84197.513 44593.971 123 51886.210 -84197.513 124 -21806.226 51886.210 125 -197243.845 -21806.226 126 65929.058 -197243.845 127 53389.293 65929.058 128 41658.487 53389.293 129 65929.058 41658.487 130 57632.350 65929.058 131 69595.804 57632.350 132 65929.058 69595.804 133 43174.000 65929.058 134 44601.864 43174.000 135 55114.720 44601.864 136 57065.590 55114.720 137 64374.858 57065.590 138 41884.819 64374.858 139 65929.058 41884.819 140 62498.832 65929.058 141 65929.058 62498.832 142 22596.029 65929.058 143 61521.371 22596.029 144 51941.111 61521.371 145 54410.801 51941.111 146 65929.058 54410.801 147 58506.192 65929.058 148 67922.218 58506.192 149 65668.025 67922.218 150 27222.295 65668.025 151 -25738.300 27222.295 152 61174.154 -25738.300 153 34741.970 61174.154 154 63120.229 34741.970 155 65929.058 63120.229 156 65929.058 65929.058 157 60553.705 65929.058 158 73101.614 60553.705 159 65929.058 73101.614 160 65929.058 65929.058 161 65929.058 65929.058 162 65929.058 65929.058 163 39346.008 65929.058 164 16046.603 39346.008 165 63444.017 16046.603 166 65929.058 63444.017 167 122277.892 65929.058 168 -62798.915 122277.892 169 65929.058 -62798.915 170 65929.058 65929.058 171 58033.087 65929.058 172 65929.058 58033.087 173 45262.036 65929.058 174 65929.058 45262.036 175 55359.960 65929.058 176 -52462.303 55359.960 177 17565.671 -52462.303 178 65929.058 17565.671 179 65929.058 65929.058 180 53805.874 65929.058 181 -5686.749 53805.874 182 -5704.113 -5686.749 183 -40076.194 -5704.113 184 9750.736 -40076.194 185 12209.395 9750.736 186 -23778.155 12209.395 187 52884.959 -23778.155 188 24969.757 52884.959 189 65929.058 24969.757 190 58427.473 65929.058 191 65929.058 58427.473 192 54768.353 65929.058 193 65929.058 54768.353 194 62508.776 65929.058 195 18372.740 62508.776 196 -161948.935 18372.740 197 18372.740 -161948.935 198 18372.740 18372.740 199 -85234.859 18372.740 200 18372.740 -85234.859 201 18372.740 18372.740 202 18372.740 18372.740 203 3123.603 18372.740 204 -34367.008 3123.603 205 12751.263 -34367.008 206 10733.070 12751.263 207 18372.740 10733.070 208 15358.754 18372.740 209 1246.183 15358.754 210 65929.058 1246.183 211 65929.058 65929.058 212 68937.951 65929.058 213 17666.963 68937.951 214 35622.966 17666.963 215 64027.096 35622.966 216 60031.984 64027.096 217 52989.146 60031.984 218 64027.096 52989.146 219 11943.597 64027.096 220 38322.337 11943.597 221 65929.058 38322.337 222 64919.348 65929.058 223 33373.929 64919.348 224 64604.144 33373.929 225 37382.774 64604.144 226 55196.557 37382.774 227 65929.058 55196.557 228 49836.275 65929.058 229 65929.058 49836.275 230 30753.308 65929.058 231 47097.884 30753.308 232 65929.058 47097.884 233 31837.494 65929.058 234 65929.058 31837.494 235 65957.049 65929.058 236 65929.058 65957.049 237 60031.142 65929.058 238 65929.058 60031.142 239 65929.058 65929.058 240 65929.058 65929.058 241 57539.279 65929.058 242 65929.058 57539.279 243 60369.653 65929.058 244 3466.936 60369.653 245 18372.740 3466.936 246 18372.740 18372.740 247 40733.687 18372.740 248 18372.740 40733.687 249 18372.740 18372.740 250 45705.516 18372.740 251 18372.740 45705.516 252 47886.726 18372.740 253 -13578.909 47886.726 254 -134807.851 -13578.909 255 65929.058 -134807.851 256 54364.204 65929.058 257 36026.365 54364.204 258 18372.740 36026.365 259 65929.058 18372.740 260 18372.740 65929.058 261 18372.740 18372.740 262 58739.116 18372.740 263 69056.918 58739.116 264 18372.740 69056.918 265 65929.058 18372.740 266 18372.740 65929.058 267 93313.467 18372.740 268 18372.740 93313.467 269 40389.473 18372.740 270 18372.740 40389.473 271 18372.740 18372.740 272 18372.740 18372.740 273 39807.516 18372.740 274 65929.058 39807.516 275 50307.221 65929.058 276 36444.067 50307.221 277 65929.058 36444.067 278 -72133.754 65929.058 279 38471.735 -72133.754 280 18372.740 38471.735 281 65929.058 18372.740 282 39889.761 65929.058 283 65929.058 39889.761 284 50651.402 65929.058 285 65929.058 50651.402 286 4199.495 65929.058 287 63121.044 4199.495 288 18372.740 63121.044 289 35164.002 18372.740 290 18372.740 35164.002 291 57455.728 18372.740 292 57472.520 57455.728 293 65929.058 57472.520 294 44123.151 65929.058 295 29038.476 44123.151 296 65929.058 29038.476 297 18372.740 65929.058 298 65929.058 18372.740 299 57370.231 65929.058 300 11573.074 57370.231 301 47299.698 11573.074 302 64978.077 47299.698 303 18372.740 64978.077 304 71940.367 18372.740 305 54464.426 71940.367 306 18372.740 54464.426 307 -8040.394 18372.740 308 44919.571 -8040.394 309 18372.740 44919.571 310 18372.740 18372.740 311 18490.060 18372.740 312 18372.740 18490.060 313 5487.357 18372.740 314 58468.829 5487.357 315 75256.829 58468.829 316 65919.056 75256.829 317 19879.629 65919.056 318 61145.628 19879.629 319 65929.058 61145.628 320 65929.058 65929.058 321 65929.058 65929.058 322 65929.058 65929.058 323 18372.740 65929.058 324 65929.058 18372.740 325 61612.076 65929.058 326 18372.740 61612.076 327 18372.740 18372.740 328 65929.058 18372.740 329 65929.058 65929.058 330 65929.058 65929.058 331 65929.058 65929.058 332 65929.058 65929.058 333 58875.669 65929.058 334 64953.405 58875.669 335 55194.372 64953.405 336 8752.196 55194.372 337 65929.058 8752.196 338 62637.715 65929.058 339 57478.821 62637.715 340 24601.292 57478.821 341 65929.058 24601.292 342 65929.058 65929.058 343 65929.058 65929.058 344 32753.089 65929.058 345 65929.058 32753.089 346 51182.227 65929.058 347 65929.058 51182.227 348 65929.058 65929.058 349 25248.341 65929.058 350 62730.287 25248.341 351 56005.167 62730.287 352 55506.617 56005.167 353 58790.037 55506.617 354 65929.058 58790.037 355 65929.058 65929.058 356 19911.198 65929.058 357 65929.058 19911.198 358 65929.058 65929.058 359 65929.058 65929.058 360 62118.054 65929.058 361 65929.058 62118.054 362 65929.058 65929.058 363 -55498.104 65929.058 364 56055.704 -55498.104 365 65929.058 56055.704 366 65929.058 65929.058 367 65929.058 65929.058 368 65929.058 65929.058 369 54639.637 65929.058 370 65929.058 54639.637 371 58803.246 65929.058 372 51349.727 58803.246 373 45059.431 51349.727 374 65929.058 45059.431 375 65929.058 65929.058 376 39687.115 65929.058 377 53282.786 39687.115 378 12142.297 53282.786 379 65929.058 12142.297 380 65929.058 65929.058 381 19736.941 65929.058 382 34263.075 19736.941 383 65929.058 34263.075 384 51741.622 65929.058 385 -104527.163 51741.622 386 73487.260 -104527.163 387 36723.162 73487.260 388 60807.854 36723.162 389 41511.317 60807.854 390 51593.419 41511.317 391 29177.487 51593.419 392 -7843.262 29177.487 393 86055.384 -7843.262 394 -131925.428 86055.384 395 17577.025 -131925.428 396 -126907.668 17577.025 397 -209372.064 -126907.668 398 23709.340 -209372.064 399 -91936.031 23709.340 400 13839.628 -91936.031 401 54010.943 13839.628 402 30616.402 54010.943 403 -84731.246 30616.402 404 2350.256 -84731.246 405 25895.871 2350.256 406 32384.186 25895.871 407 24608.709 32384.186 408 34074.969 24608.709 409 -50506.225 34074.969 410 -36839.836 -50506.225 411 -197406.940 -36839.836 412 -62015.747 -197406.940 413 25808.647 -62015.747 414 -165224.220 25808.647 415 -87688.976 -165224.220 416 102732.060 -87688.976 417 -52946.773 102732.060 418 -37843.496 -52946.773 419 -26327.090 -37843.496 420 -33004.999 -26327.090 421 -159782.676 -33004.999 422 -584784.875 -159782.676 423 -447047.074 -584784.875 424 -37244.261 -447047.074 425 -281840.274 -37244.261 426 -300445.488 -281840.274 427 -1107109.319 -300445.488 428 -258823.827 -1107109.319 429 -319850.175 -258823.827 430 -669870.658 -319850.175 > 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/7by871291221210.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/8by871291221210.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/9by871291221210.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/10mq8s1291221210.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/11pqoy1291221210.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/12brn41291221210.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/13hskx1291221210.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/14sjji1291221210.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/15e2zo1291221210.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/16abxx1291221210.tab") + } > > try(system("convert tmp/1xpby1291221210.ps tmp/1xpby1291221210.png",intern=TRUE)) character(0) > try(system("convert tmp/2xpby1291221210.ps tmp/2xpby1291221210.png",intern=TRUE)) character(0) > try(system("convert tmp/38ga11291221210.ps tmp/38ga11291221210.png",intern=TRUE)) character(0) > try(system("convert tmp/48ga11291221210.ps tmp/48ga11291221210.png",intern=TRUE)) character(0) > try(system("convert tmp/58ga11291221210.ps tmp/58ga11291221210.png",intern=TRUE)) character(0) > try(system("convert tmp/61p941291221210.ps tmp/61p941291221210.png",intern=TRUE)) character(0) > try(system("convert tmp/7by871291221210.ps tmp/7by871291221210.png",intern=TRUE)) character(0) > try(system("convert tmp/8by871291221210.ps tmp/8by871291221210.png",intern=TRUE)) character(0) > try(system("convert tmp/9by871291221210.ps tmp/9by871291221210.png",intern=TRUE)) character(0) > try(system("convert tmp/10mq8s1291221210.ps tmp/10mq8s1291221210.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 12.309 2.246 29.216