R version 2.9.0 (2009-04-17) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(0 + ,162556 + ,807 + ,213118 + ,6282154 + ,0 + ,29790 + ,444 + ,81767 + ,4321023 + ,0 + ,87550 + ,412 + ,153198 + ,4111912 + ,1 + ,84738 + ,428 + ,-26007 + ,223193 + ,0 + ,54660 + ,315 + ,126942 + ,1491348 + ,0 + ,42634 + ,168 + ,157214 + ,1629616 + ,1 + ,40949 + ,263 + ,129352 + ,1398893 + ,0 + ,45187 + ,267 + ,234817 + ,1926517 + ,0 + ,37704 + ,228 + ,60448 + ,983660 + ,0 + ,16275 + ,129 + ,47818 + ,1443586 + ,1 + ,25830 + ,104 + ,245546 + ,1073089 + ,1 + ,12679 + ,122 + ,48020 + ,984885 + ,0 + ,18014 + ,393 + ,-1710 + ,1405225 + ,1 + ,43556 + ,190 + ,32648 + ,227132 + ,0 + ,24811 + ,280 + ,95350 + ,929118 + ,1 + ,6575 + ,63 + ,151352 + ,1071292 + ,1 + ,7123 + ,102 + ,288170 + ,638830 + ,0 + ,21950 + ,265 + ,114337 + ,856956 + ,0 + ,37597 + ,234 + ,37884 + ,992426 + ,1 + ,17821 + ,277 + ,122844 + ,444477 + ,0 + ,12988 + ,73 + ,82340 + ,857217 + ,0 + ,22330 + ,67 + ,79801 + ,711969 + ,1 + ,13326 + ,103 + ,165548 + ,702380 + ,1 + ,16189 + ,290 + ,116384 + ,358589 + ,1 + ,7146 + ,83 + ,134028 + ,297978 + ,1 + ,15824 + ,56 + ,63838 + ,585715 + ,0 + ,27664 + ,236 + ,74996 + ,657954 + ,1 + ,11920 + ,73 + ,31080 + ,209458 + ,1 + ,8568 + ,34 + ,32168 + ,786690 + ,1 + ,14416 + ,139 + ,49857 + ,439798 + ,0 + ,3369 + ,26 + ,87161 + ,688779 + ,0 + ,11819 + ,70 + ,106113 + ,574339 + ,0 + ,6984 + ,40 + ,80570 + ,741409 + ,0 + ,4519 + ,42 + ,102129 + ,597793 + ,1 + ,2220 + ,12 + ,301670 + ,644190 + ,1 + ,18562 + ,211 + ,102313 + ,377934 + ,1 + ,10327 + ,74 + ,88577 + ,640273 + ,0 + ,5336 + ,80 + ,112477 + ,697458 + ,0 + ,2365 + ,83 + ,191778 + ,550608 + ,1 + ,4069 + ,131 + ,79804 + ,207393 + ,1 + ,8636 + ,203 + ,128294 + ,301607 + ,1 + ,13718 + ,56 + ,96448 + ,345783 + ,1 + ,4525 + ,89 + ,93811 + ,501749 + ,1 + ,6869 + ,88 + ,117520 + ,379983 + ,1 + ,4628 + ,39 + ,69159 + ,387475 + ,0 + ,3689 + ,25 + ,101792 + ,377305 + ,0 + ,4891 + ,49 + ,210568 + ,370837 + ,0 + ,7489 + ,149 + ,136996 + ,430866 + ,1 + ,4901 + ,58 + ,121920 + ,469107 + ,1 + ,2284 + ,41 + ,76403 + ,194493 + ,0 + ,3160 + ,90 + ,108094 + ,530670 + ,0 + ,4150 + ,136 + ,134759 + ,518365 + ,0 + ,7285 + ,97 + ,188873 + ,491303 + ,0 + ,1134 + ,63 + ,146216 + ,527021 + ,0 + ,4658 + ,114 + ,156608 + ,233773 + ,1 + ,2384 + ,77 + ,61348 + ,405972 + ,1 + ,3748 + ,6 + ,50350 + ,652925 + ,1 + ,5371 + ,47 + ,87720 + ,446211 + ,1 + ,1285 + ,51 + ,99489 + ,341340 + ,0 + ,9327 + ,85 + ,87419 + ,387699 + ,0 + ,5565 + ,43 + ,94355 + ,493408 + ,1 + ,1528 + ,32 + ,60326 + ,146494 + ,0 + ,3122 + ,25 + ,94670 + ,414462 + ,0 + ,7561 + ,77 + ,82425 + ,364304 + ,1 + ,2675 + ,54 + ,59017 + ,355178 + ,1 + ,13253 + ,251 + ,90829 + ,357760 + ,1 + ,880 + ,15 + ,80791 + ,261216 + ,0 + ,2053 + ,44 + ,100423 + ,397144 + ,1 + ,1424 + ,73 + ,131116 + ,374943 + ,0 + ,4036 + ,85 + ,100269 + ,424898 + ,0 + ,3045 + ,49 + ,27330 + ,202055 + ,1 + ,5119 + ,38 + ,39039 + ,378525 + ,1 + ,1431 + ,35 + ,106885 + ,310768 + ,1 + ,554 + ,9 + ,79285 + ,325738 + ,1 + ,1975 + ,34 + ,118881 + ,394510 + ,0 + ,1765 + ,20 + ,77623 + ,247060 + ,1 + ,1012 + ,29 + ,114768 + ,368078 + ,1 + ,810 + ,11 + ,74015 + ,236761 + ,1 + ,1280 + ,52 + ,69465 + ,312378 + ,0 + ,666 + ,13 + ,117869 + ,339836 + ,1 + ,1380 + ,29 + ,60982 + ,347385 + ,0 + ,4677 + ,66 + ,90131 + ,426280 + ,1 + ,876 + ,33 + ,138971 + ,352850 + ,1 + ,814 + ,15 + ,39625 + ,301881 + ,1 + ,514 + ,15 + ,102725 + ,377516 + ,0 + ,5692 + ,68 + ,64239 + ,357312 + ,1 + ,3642 + ,100 + ,90262 + ,458343 + ,1 + ,540 + ,13 + ,103960 + ,354228 + ,1 + ,2099 + ,45 + ,106611 + ,308636 + ,1 + ,567 + ,14 + ,103345 + ,386212 + ,1 + ,2001 + ,36 + ,95551 + ,393343 + ,0 + ,2949 + ,40 + ,82903 + ,378509 + ,1 + ,2253 + ,68 + ,63593 + ,452469 + ,0 + ,6533 + ,29 + ,126910 + ,364839 + ,1 + ,1889 + ,43 + ,37527 + ,358649 + ,0 + ,3055 + ,30 + ,60247 + ,376641 + ,1 + ,272 + ,9 + ,112995 + ,429112 + ,0 + ,1414 + ,22 + ,70184 + ,330546 + ,1 + ,2564 + ,19 + ,130140 + ,403560 + ,0 + ,1383 + ,9 + ,73221 + ,317892 + ,0 + ,1261 + ,31 + ,76114 + ,307528 + ,1 + ,975 + ,19 + ,90534 + ,235133 + ,1 + ,3366 + ,55 + ,108479 + ,299243 + ,1 + ,576 + ,8 + ,113761 + ,314073 + ,1 + ,1686 + ,28 + ,68696 + ,368186 + ,1 + ,746 + ,29 + ,71561 + ,269661 + ,0 + ,3192 + ,48 + ,59831 + ,125390 + ,0 + ,2045 + ,16 + ,97890 + ,510834 + ,1 + ,5702 + ,47 + ,101481 + ,321896 + ,0 + ,1932 + ,20 + ,72954 + ,249898 + ,1 + ,936 + ,22 + ,67939 + ,408881 + ,0 + ,3437 + ,33 + ,48022 + ,158492 + ,1 + ,5131 + ,44 + ,86111 + ,292154 + ,0 + ,2397 + ,13 + ,74020 + ,289513 + ,0 + ,1389 + ,6 + ,57530 + ,378049 + ,1 + ,1503 + ,35 + ,56364 + ,343466 + ,1 + ,402 + ,8 + ,84990 + ,332743 + ,1 + ,2239 + ,17 + ,88590 + ,442882 + ,0 + ,2234 + ,11 + ,77200 + ,214215 + ,1 + ,837 + ,21 + ,61262 + ,315688 + ,1 + ,10579 + ,92 + ,110309 + ,375195 + ,1 + ,875 + ,12 + ,67000 + ,334280 + ,1 + ,1585 + ,112 + ,93099 + ,355864 + ,0 + ,1659 + ,25 + ,107577 + ,480382 + ,0 + ,2647 + ,17 + ,62920 + ,353058 + ,0 + ,3294 + ,23 + ,75832 + ,217193 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,94 + ,10 + ,60793 + ,314533 + ,1 + ,422 + ,23 + ,57935 + ,318056 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,34 + ,7 + ,60630 + ,314353 + ,1 + ,1558 + ,25 + ,55637 + ,369448 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,43 + ,20 + ,60887 + ,312846 + ,1 + ,645 + ,4 + ,60720 + ,312075 + ,1 + ,316 + ,4 + ,60505 + ,315009 + ,1 + ,115 + ,10 + ,60945 + ,318903 + ,1 + ,5 + ,1 + ,60720 + ,314887 + ,1 + ,897 + ,4 + ,60720 + ,314913 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,389 + ,8 + ,58990 + ,325506 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,1002 + ,11 + ,56750 + ,298568 + ,1 + ,36 + ,4 + ,60894 + ,315834 + ,1 + ,460 + ,15 + ,63346 + ,329784 + ,1 + ,309 + ,9 + ,56535 + ,312878 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,9 + ,7 + ,60835 + ,314987 + ,1 + ,271 + ,2 + ,60720 + ,325249 + ,1 + ,14 + ,0 + ,61016 + ,315877 + ,1 + ,520 + ,7 + ,58650 + ,291650 + ,1 + ,1766 + ,46 + ,60438 + ,305959 + ,1 + ,0 + ,5 + ,60720 + ,315380 + ,1 + ,458 + ,7 + ,58625 + ,297765 + ,1 + ,20 + ,2 + ,60938 + ,315245 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,98 + ,2 + ,61490 + ,315236 + ,1 + ,405 + ,5 + ,60845 + ,336425 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,483 + ,7 + ,60830 + ,306268 + ,1 + ,454 + ,24 + ,63261 + ,302187 + ,1 + ,47 + ,1 + ,60720 + ,314882 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,757 + ,18 + ,45689 + ,382712 + ,1 + ,4655 + ,55 + ,60720 + ,341570 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,36 + ,3 + ,61564 + ,312412 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,203 + ,9 + ,61938 + ,309596 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,126 + ,8 + ,60951 + ,315547 + ,1 + ,400 + ,113 + ,60720 + ,313267 + ,0 + ,71 + ,0 + ,60745 + ,316176 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,972 + ,19 + ,71642 + ,359335 + ,0 + ,531 + ,11 + ,71641 + ,330068 + ,1 + ,2461 + ,25 + ,55792 + ,314289 + ,0 + ,378 + ,16 + ,71873 + ,297413 + ,0 + ,23 + ,5 + ,62555 + ,314806 + ,0 + ,638 + ,11 + ,60370 + ,333210 + ,0 + ,2300 + ,23 + ,64873 + ,352108 + ,1 + ,149 + ,6 + ,62041 + ,313332 + ,1 + ,226 + ,5 + ,65745 + ,291787 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,275 + ,7 + ,59500 + ,318745 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,141 + ,7 + ,61630 + ,315366 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,28 + ,3 + ,60890 + ,315688 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,4980 + ,89 + ,113521 + ,409642 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,472 + ,19 + ,80045 + ,269587 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,203 + ,12 + ,50804 + ,300962 + ,0 + ,496 + ,12 + ,87390 + ,325479 + ,0 + ,10 + ,5 + ,61656 + ,316155 + ,0 + ,63 + ,2 + ,65688 + ,318574 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1136 + ,26 + ,48522 + ,343613 + ,0 + ,265 + ,3 + ,60720 + ,306948 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,267 + ,11 + ,57640 + ,330059 + ,1 + ,474 + ,10 + ,61977 + ,288985 + ,1 + ,534 + ,5 + ,62620 + ,304485 + ,1 + ,0 + ,2 + ,60720 + ,315380 + ,1 + ,15 + ,6 + ,60831 + ,315688 + ,1 + ,397 + ,7 + ,60646 + ,317736 + ,1 + ,0 + ,2 + ,60720 + ,315380 + ,1 + ,1866 + ,28 + ,56225 + ,322331 + ,1 + ,288 + ,3 + ,60510 + ,296656 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,3 + ,1 + ,60698 + ,315354 + ,1 + ,468 + ,20 + ,60720 + ,312161 + ,1 + ,20 + ,1 + ,60805 + ,315576 + ,1 + ,278 + ,22 + ,61404 + ,314922 + ,1 + ,61 + ,9 + ,60720 + ,314551 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,192 + ,2 + ,65276 + ,312339 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,317 + ,7 + ,63915 + ,298700 + ,1 + ,738 + ,9 + ,60720 + ,321376 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,368 + ,13 + ,61686 + ,303230 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,2 + ,0 + ,60743 + ,315487 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,53 + ,6 + ,60349 + ,315793 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,94 + ,3 + ,61360 + ,312887 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,24 + ,7 + ,59818 + ,315637 + ,1 + ,2332 + ,2 + ,72680 + ,324385 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,131 + ,15 + ,61808 + ,308989 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,206 + ,9 + ,53110 + ,296702 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,167 + ,1 + ,64245 + ,307322 + ,1 + ,622 + ,38 + ,73007 + ,304376 + ,1 + ,2328 + ,57 + ,82732 + ,253588 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,365 + ,7 + ,54820 + ,309560 + ,1 + ,364 + ,26 + ,47705 + ,298466 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,226 + ,13 + ,72835 + ,343929 + ,1 + ,307 + ,10 + ,58856 + ,331955 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,188 + ,9 + ,77655 + ,381180 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,138 + ,26 + ,69817 + ,331420 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,125 + ,19 + ,60798 + ,310201 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,282 + ,12 + ,62452 + ,320016 + ,1 + ,335 + ,23 + ,64175 + ,320398 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1324 + ,29 + ,67440 + ,291841 + ,1 + ,176 + ,8 + ,68136 + ,310670 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,249 + ,26 + ,56726 + ,313491 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,333 + ,9 + ,70811 + ,331323 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,601 + ,5 + ,60720 + ,319210 + ,1 + ,30 + ,3 + ,62045 + ,318098 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,249 + ,13 + ,54323 + ,292754 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,165 + ,12 + ,62841 + ,325176 + ,1 + ,453 + ,19 + ,81125 + ,365959 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,53 + ,10 + ,59506 + ,302409 + ,0 + ,382 + ,9 + ,59365 + ,340968 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,9 + ,60720 + ,315380 + ,0 + ,30 + ,4 + ,60798 + ,313164 + ,1 + ,290 + ,1 + ,58790 + ,301164 + ,1 + ,0 + ,1 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,366 + ,14 + ,61808 + ,344425 + ,1 + ,2 + ,12 + ,60735 + ,315394 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,209 + ,19 + ,64016 + ,316647 + ,1 + ,384 + ,17 + ,54683 + ,309836 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,365 + ,32 + ,87192 + ,346611 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,49 + ,14 + ,64107 + ,322031 + ,1 + ,3 + ,8 + ,60761 + ,315656 + ,1 + ,133 + ,4 + ,65990 + ,339445 + ,1 + ,32 + ,0 + ,59988 + ,314964 + ,1 + ,368 + ,20 + ,61167 + ,297141 + ,1 + ,1 + ,5 + ,60719 + ,315372 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,22 + ,1 + ,60722 + ,312502 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,96 + ,4 + ,60379 + ,313729 + ,1 + ,1 + ,1 + ,60727 + ,315388 + ,1 + ,314 + ,4 + ,60720 + ,315371 + ,1 + ,844 + ,20 + ,60925 + ,296139 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,26 + ,1 + ,60896 + ,313880 + ,1 + ,125 + ,10 + ,59734 + ,317698 + ,1 + ,304 + ,12 + ,62969 + ,295580 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,621 + ,13 + ,60720 + ,308256 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,119 + ,3 + ,59118 + ,303677 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,1595 + ,10 + ,60720 + ,319369 + ,1 + ,312 + ,3 + ,58598 + ,318690 + ,1 + ,60 + ,7 + ,61124 + ,314049 + ,1 + ,587 + ,10 + ,59595 + ,325699 + ,1 + ,135 + ,1 + ,62065 + ,314210 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,514 + ,15 + ,78780 + ,322378 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,1 + ,4 + ,60722 + ,315398 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1763 + ,28 + ,61600 + ,308336 + ,1 + ,180 + ,9 + ,59635 + ,316386 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,218 + ,7 + ,60720 + ,315553 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,448 + ,7 + ,59781 + ,323361 + ,1 + ,227 + ,7 + ,76644 + ,336639 + ,1 + ,174 + ,3 + ,64820 + ,307424 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,121 + ,11 + ,56178 + ,295370 + ,1 + ,607 + ,7 + ,60436 + ,322340 + ,1 + ,2212 + ,10 + ,60720 + ,319864 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,530 + ,18 + ,73433 + ,317291 + ,1 + ,571 + ,14 + ,41477 + ,280398 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,78 + ,12 + ,62700 + ,317330 + ,1 + ,2489 + ,29 + ,67804 + ,238125 + ,1 + ,131 + ,3 + ,59661 + ,327071 + ,1 + ,923 + ,6 + ,58620 + ,309038 + ,1 + ,72 + ,3 + ,60398 + ,314210 + ,1 + ,572 + ,8 + ,58580 + ,307930 + ,1 + ,397 + ,10 + ,62710 + ,322327 + ,1 + ,450 + ,6 + ,59325 + ,292136 + ,1 + ,622 + ,8 + ,60950 + ,263276 + ,1 + ,694 + ,6 + ,68060 + ,367655 + ,0 + ,3425 + ,9 + ,83620 + ,283910 + ,1 + ,562 + ,8 + ,58456 + ,283587 + ,1 + ,4917 + ,26 + ,52811 + ,243650 + ,0 + ,1442 + ,239 + ,121173 + ,438493 + ,1 + ,529 + ,7 + ,63870 + ,296261 + ,0 + ,2126 + ,41 + ,21001 + ,230621 + ,1 + ,1061 + ,3 + ,70415 + ,304252 + ,1 + ,776 + ,8 + ,64230 + ,333505 + ,1 + ,611 + ,6 + ,59190 + ,296919 + ,0 + ,1526 + ,21 + ,69351 + ,278990 + ,1 + ,592 + ,7 + ,64270 + ,276898 + ,1 + ,1182 + ,11 + ,70694 + ,327007 + ,1 + ,621 + ,11 + ,68005 + ,317046 + ,1 + ,989 + ,12 + ,58930 + ,304555 + ,1 + ,438 + ,9 + ,58320 + ,298096 + ,1 + ,726 + ,3 + ,69980 + ,231861 + ,1 + ,1303 + ,57 + ,69863 + ,309422 + ,0 + ,7419 + ,21 + ,63255 + ,286963 + ,0 + ,1164 + ,15 + ,57320 + ,269753 + ,0 + ,3310 + ,32 + ,75230 + ,448243 + ,1 + ,1920 + ,11 + ,79420 + ,165404 + ,1 + ,965 + ,2 + ,73490 + ,204325 + ,1 + ,3256 + ,23 + ,35250 + ,407159 + ,0 + ,1135 + ,20 + ,62285 + ,290476 + ,1 + ,1270 + ,24 + ,69206 + ,275311 + ,1 + ,661 + ,1 + ,65920 + ,246541 + ,1 + ,1013 + ,1 + ,69770 + ,253468 + ,1 + ,2844 + ,74 + ,72683 + ,240897 + ,0 + ,11528 + ,68 + ,-14545 + ,-83265 + ,1 + ,6526 + ,20 + ,55830 + ,-42143 + ,1 + ,2264 + ,20 + ,55174 + ,272713 + ,0 + ,5109 + ,82 + ,67038 + ,215362 + ,1 + ,3999 + ,21 + ,51252 + ,42754 + ,1 + ,35624 + ,244 + ,157278 + ,306275 + ,1 + ,9252 + ,32 + ,79510 + ,253537 + ,1 + ,15236 + ,86 + ,77440 + ,372631 + ,1 + ,18073 + ,69 + ,27284 + ,-7170) + ,dim=c(5 + ,431) + ,dimnames=list(c('group' + ,'costs' + ,'orders' + ,'dividends' + ,'wealth ') + ,1:431)) > y <- array(NA,dim=c(5,431),dimnames=list(c('group','costs','orders','dividends','wealth '),1:431)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par3 = 'No Linear Trend' > par2 = 'Do not include Seasonal Dummies' > par1 = '5' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > library(lattice) > library(lmtest) Loading required package: zoo Attaching package: 'zoo' The following object(s) are masked from package:base : as.Date.numeric > n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test > par1 <- as.numeric(par1) > x <- t(y) > k <- length(x[1,]) > n <- length(x[,1]) > x1 <- cbind(x[,par1], x[,1:k!=par1]) > mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1]) > colnames(x1) <- mycolnames #colnames(x)[par1] > x <- x1 > if (par3 == 'First Differences'){ + x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep=''))) + for (i in 1:n-1) { + for (j in 1:k) { + x2[i,j] <- x[i+1,j] - x[i,j] + } + } + x <- x2 + } > if (par2 == 'Include Monthly Dummies'){ + x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep =''))) + for (i in 1:11){ + x2[seq(i,n,12),i] <- 1 + } + x <- cbind(x, x2) + } > if (par2 == 'Include Quarterly Dummies'){ + x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep =''))) + for (i in 1:3){ + x2[seq(i,n,4),i] <- 1 + } + x <- cbind(x, x2) + } > k <- length(x[1,]) > if (par3 == 'Linear Trend'){ + x <- cbind(x, c(1:n)) + colnames(x)[k+1] <- 't' + } > x wealth\r group costs orders dividends 1 6282154 0 162556 807 213118 2 4321023 0 29790 444 81767 3 4111912 0 87550 412 153198 4 223193 1 84738 428 -26007 5 1491348 0 54660 315 126942 6 1629616 0 42634 168 157214 7 1398893 1 40949 263 129352 8 1926517 0 45187 267 234817 9 983660 0 37704 228 60448 10 1443586 0 16275 129 47818 11 1073089 1 25830 104 245546 12 984885 1 12679 122 48020 13 1405225 0 18014 393 -1710 14 227132 1 43556 190 32648 15 929118 0 24811 280 95350 16 1071292 1 6575 63 151352 17 638830 1 7123 102 288170 18 856956 0 21950 265 114337 19 992426 0 37597 234 37884 20 444477 1 17821 277 122844 21 857217 0 12988 73 82340 22 711969 0 22330 67 79801 23 702380 1 13326 103 165548 24 358589 1 16189 290 116384 25 297978 1 7146 83 134028 26 585715 1 15824 56 63838 27 657954 0 27664 236 74996 28 209458 1 11920 73 31080 29 786690 1 8568 34 32168 30 439798 1 14416 139 49857 31 688779 0 3369 26 87161 32 574339 0 11819 70 106113 33 741409 0 6984 40 80570 34 597793 0 4519 42 102129 35 644190 1 2220 12 301670 36 377934 1 18562 211 102313 37 640273 1 10327 74 88577 38 697458 0 5336 80 112477 39 550608 0 2365 83 191778 40 207393 1 4069 131 79804 41 301607 1 8636 203 128294 42 345783 1 13718 56 96448 43 501749 1 4525 89 93811 44 379983 1 6869 88 117520 45 387475 1 4628 39 69159 46 377305 0 3689 25 101792 47 370837 0 4891 49 210568 48 430866 0 7489 149 136996 49 469107 1 4901 58 121920 50 194493 1 2284 41 76403 51 530670 0 3160 90 108094 52 518365 0 4150 136 134759 53 491303 0 7285 97 188873 54 527021 0 1134 63 146216 55 233773 0 4658 114 156608 56 405972 1 2384 77 61348 57 652925 1 3748 6 50350 58 446211 1 5371 47 87720 59 341340 1 1285 51 99489 60 387699 0 9327 85 87419 61 493408 0 5565 43 94355 62 146494 1 1528 32 60326 63 414462 0 3122 25 94670 64 364304 0 7561 77 82425 65 355178 1 2675 54 59017 66 357760 1 13253 251 90829 67 261216 1 880 15 80791 68 397144 0 2053 44 100423 69 374943 1 1424 73 131116 70 424898 0 4036 85 100269 71 202055 0 3045 49 27330 72 378525 1 5119 38 39039 73 310768 1 1431 35 106885 74 325738 1 554 9 79285 75 394510 1 1975 34 118881 76 247060 0 1765 20 77623 77 368078 1 1012 29 114768 78 236761 1 810 11 74015 79 312378 1 1280 52 69465 80 339836 0 666 13 117869 81 347385 1 1380 29 60982 82 426280 0 4677 66 90131 83 352850 1 876 33 138971 84 301881 1 814 15 39625 85 377516 1 514 15 102725 86 357312 0 5692 68 64239 87 458343 1 3642 100 90262 88 354228 1 540 13 103960 89 308636 1 2099 45 106611 90 386212 1 567 14 103345 91 393343 1 2001 36 95551 92 378509 0 2949 40 82903 93 452469 1 2253 68 63593 94 364839 0 6533 29 126910 95 358649 1 1889 43 37527 96 376641 0 3055 30 60247 97 429112 1 272 9 112995 98 330546 0 1414 22 70184 99 403560 1 2564 19 130140 100 317892 0 1383 9 73221 101 307528 0 1261 31 76114 102 235133 1 975 19 90534 103 299243 1 3366 55 108479 104 314073 1 576 8 113761 105 368186 1 1686 28 68696 106 269661 1 746 29 71561 107 125390 0 3192 48 59831 108 510834 0 2045 16 97890 109 321896 1 5702 47 101481 110 249898 0 1932 20 72954 111 408881 1 936 22 67939 112 158492 0 3437 33 48022 113 292154 1 5131 44 86111 114 289513 0 2397 13 74020 115 378049 0 1389 6 57530 116 343466 1 1503 35 56364 117 332743 1 402 8 84990 118 442882 1 2239 17 88590 119 214215 0 2234 11 77200 120 315688 1 837 21 61262 121 375195 1 10579 92 110309 122 334280 1 875 12 67000 123 355864 1 1585 112 93099 124 480382 0 1659 25 107577 125 353058 0 2647 17 62920 126 217193 0 3294 23 75832 127 315380 1 0 0 60720 128 314533 1 94 10 60793 129 318056 1 422 23 57935 130 315380 1 0 0 60720 131 314353 1 34 7 60630 132 369448 1 1558 25 55637 133 315380 1 0 0 60720 134 312846 1 43 20 60887 135 312075 1 645 4 60720 136 315009 1 316 4 60505 137 318903 1 115 10 60945 138 314887 1 5 1 60720 139 314913 1 897 4 60720 140 315380 1 0 0 60720 141 325506 1 389 8 58990 142 315380 1 0 0 60720 143 298568 1 1002 11 56750 144 315834 1 36 4 60894 145 329784 1 460 15 63346 146 312878 1 309 9 56535 147 315380 1 0 0 60720 148 314987 1 9 7 60835 149 325249 1 271 2 60720 150 315877 1 14 0 61016 151 291650 1 520 7 58650 152 305959 1 1766 46 60438 153 315380 1 0 5 60720 154 297765 1 458 7 58625 155 315245 1 20 2 60938 156 315380 1 0 0 60720 157 315380 1 0 0 60720 158 315236 1 98 2 61490 159 336425 1 405 5 60845 160 315380 1 0 0 60720 161 315380 1 0 0 60720 162 315380 1 0 0 60720 163 315380 1 0 0 60720 164 306268 1 483 7 60830 165 302187 1 454 24 63261 166 314882 1 47 1 60720 167 315380 1 0 0 60720 168 382712 1 757 18 45689 169 341570 1 4655 55 60720 170 315380 1 0 0 60720 171 315380 1 0 0 60720 172 312412 1 36 3 61564 173 315380 1 0 0 60720 174 309596 1 203 9 61938 175 315380 1 0 0 60720 176 315547 1 126 8 60951 177 313267 1 400 113 60720 178 316176 0 71 0 60745 179 315380 1 0 0 60720 180 315380 1 0 0 60720 181 359335 1 972 19 71642 182 330068 0 531 11 71641 183 314289 1 2461 25 55792 184 297413 0 378 16 71873 185 314806 0 23 5 62555 186 333210 0 638 11 60370 187 352108 0 2300 23 64873 188 313332 1 149 6 62041 189 291787 1 226 5 65745 190 315380 1 0 0 60720 191 318745 1 275 7 59500 192 315380 1 0 0 60720 193 315366 1 141 7 61630 194 315380 1 0 0 60720 195 315688 1 28 3 60890 196 315380 0 0 0 60720 197 409642 0 4980 89 113521 198 315380 0 0 0 60720 199 315380 0 0 0 60720 200 269587 0 472 19 80045 201 315380 0 0 0 60720 202 315380 0 0 0 60720 203 315380 0 0 0 60720 204 300962 0 203 12 50804 205 325479 0 496 12 87390 206 316155 0 10 5 61656 207 318574 0 63 2 65688 208 315380 0 0 0 60720 209 343613 0 1136 26 48522 210 306948 0 265 3 60720 211 315380 1 0 0 60720 212 315380 1 0 0 60720 213 330059 1 267 11 57640 214 288985 1 474 10 61977 215 304485 1 534 5 62620 216 315380 1 0 2 60720 217 315688 1 15 6 60831 218 317736 1 397 7 60646 219 315380 1 0 2 60720 220 322331 1 1866 28 56225 221 296656 1 288 3 60510 222 315380 1 0 0 60720 223 315354 1 3 1 60698 224 312161 1 468 20 60720 225 315576 1 20 1 60805 226 314922 1 278 22 61404 227 314551 1 61 9 60720 228 315380 1 0 0 60720 229 312339 1 192 2 65276 230 315380 1 0 0 60720 231 298700 1 317 7 63915 232 321376 1 738 9 60720 233 315380 1 0 0 60720 234 303230 1 368 13 61686 235 315380 1 0 0 60720 236 315487 1 2 0 60743 237 315380 1 0 0 60720 238 315793 1 53 6 60349 239 315380 1 0 0 60720 240 315380 1 0 0 60720 241 315380 1 0 0 60720 242 312887 1 94 3 61360 243 315380 1 0 0 60720 244 315637 1 24 7 59818 245 324385 1 2332 2 72680 246 315380 0 0 0 60720 247 315380 0 0 0 60720 248 308989 1 131 15 61808 249 315380 0 0 0 60720 250 315380 0 0 0 60720 251 296702 1 206 9 53110 252 315380 0 0 0 60720 253 307322 1 167 1 64245 254 304376 1 622 38 73007 255 253588 1 2328 57 82732 256 315380 1 0 0 60720 257 309560 1 365 7 54820 258 298466 1 364 26 47705 259 315380 0 0 0 60720 260 315380 1 0 0 60720 261 315380 0 0 0 60720 262 315380 0 0 0 60720 263 343929 1 226 13 72835 264 331955 1 307 10 58856 265 315380 0 0 0 60720 266 315380 1 0 0 60720 267 315380 0 0 0 60720 268 381180 1 188 9 77655 269 315380 0 0 0 60720 270 331420 1 138 26 69817 271 315380 0 0 0 60720 272 315380 0 0 0 60720 273 315380 0 0 0 60720 274 310201 1 125 19 60798 275 315380 1 0 0 60720 276 320016 1 282 12 62452 277 320398 1 335 23 64175 278 315380 1 0 0 60720 279 291841 0 1324 29 67440 280 310670 1 176 8 68136 281 315380 0 0 0 60720 282 315380 1 0 0 60720 283 313491 1 249 26 56726 284 315380 1 0 0 60720 285 331323 1 333 9 70811 286 315380 1 0 0 60720 287 319210 0 601 5 60720 288 318098 1 30 3 62045 289 315380 0 0 0 60720 290 292754 1 249 13 54323 291 315380 0 0 0 60720 292 325176 1 165 12 62841 293 365959 1 453 19 81125 294 315380 1 0 0 60720 295 302409 1 53 10 59506 296 340968 0 382 9 59365 297 315380 1 0 0 60720 298 315380 0 0 0 60720 299 315380 1 0 0 60720 300 315380 1 0 9 60720 301 313164 0 30 4 60798 302 301164 1 290 1 58790 303 315380 1 0 1 60720 304 315380 0 0 0 60720 305 344425 1 366 14 61808 306 315394 1 2 12 60735 307 315380 0 0 0 60720 308 316647 0 209 19 64016 309 309836 1 384 17 54683 310 315380 0 0 0 60720 311 315380 0 0 0 60720 312 346611 1 365 32 87192 313 315380 0 0 0 60720 314 322031 0 49 14 64107 315 315656 1 3 8 60761 316 339445 1 133 4 65990 317 314964 1 32 0 59988 318 297141 1 368 20 61167 319 315372 1 1 5 60719 320 315380 1 0 0 60720 321 315380 1 0 0 60720 322 315380 1 0 0 60720 323 315380 1 0 0 60720 324 315380 0 0 0 60720 325 315380 1 0 0 60720 326 312502 1 22 1 60722 327 315380 0 0 0 60720 328 315380 0 0 0 60720 329 315380 1 0 0 60720 330 315380 1 0 0 60720 331 315380 1 0 0 60720 332 315380 1 0 0 60720 333 315380 1 0 0 60720 334 313729 1 96 4 60379 335 315388 1 1 1 60727 336 315371 1 314 4 60720 337 296139 1 844 20 60925 338 315380 1 0 0 60720 339 313880 1 26 1 60896 340 317698 1 125 10 59734 341 295580 1 304 12 62969 342 315380 1 0 0 60720 343 315380 1 0 0 60720 344 315380 1 0 0 60720 345 308256 1 621 13 60720 346 315380 1 0 0 60720 347 303677 1 119 3 59118 348 315380 1 0 0 60720 349 315380 1 0 0 60720 350 319369 1 1595 10 60720 351 318690 1 312 3 58598 352 314049 1 60 7 61124 353 325699 1 587 10 59595 354 314210 1 135 1 62065 355 315380 1 0 0 60720 356 315380 1 0 0 60720 357 322378 1 514 15 78780 358 315380 1 0 0 60720 359 315380 1 0 0 60720 360 315380 1 0 0 60720 361 315398 1 1 4 60722 362 315380 1 0 0 60720 363 315380 1 0 0 60720 364 308336 0 1763 28 61600 365 316386 1 180 9 59635 366 315380 1 0 0 60720 367 315380 1 0 0 60720 368 315380 1 0 0 60720 369 315380 1 0 0 60720 370 315553 1 218 7 60720 371 315380 1 0 0 60720 372 323361 1 448 7 59781 373 336639 1 227 7 76644 374 307424 1 174 3 64820 375 315380 1 0 0 60720 376 315380 1 0 0 60720 377 295370 1 121 11 56178 378 322340 1 607 7 60436 379 319864 1 2212 10 60720 380 315380 1 0 0 60720 381 315380 1 0 0 60720 382 317291 1 530 18 73433 383 280398 1 571 14 41477 384 315380 1 0 0 60720 385 317330 1 78 12 62700 386 238125 1 2489 29 67804 387 327071 1 131 3 59661 388 309038 1 923 6 58620 389 314210 1 72 3 60398 390 307930 1 572 8 58580 391 322327 1 397 10 62710 392 292136 1 450 6 59325 393 263276 1 622 8 60950 394 367655 1 694 6 68060 395 283910 0 3425 9 83620 396 283587 1 562 8 58456 397 243650 1 4917 26 52811 398 438493 0 1442 239 121173 399 296261 1 529 7 63870 400 230621 0 2126 41 21001 401 304252 1 1061 3 70415 402 333505 1 776 8 64230 403 296919 1 611 6 59190 404 278990 0 1526 21 69351 405 276898 1 592 7 64270 406 327007 1 1182 11 70694 407 317046 1 621 11 68005 408 304555 1 989 12 58930 409 298096 1 438 9 58320 410 231861 1 726 3 69980 411 309422 1 1303 57 69863 412 286963 0 7419 21 63255 413 269753 0 1164 15 57320 414 448243 0 3310 32 75230 415 165404 1 1920 11 79420 416 204325 1 965 2 73490 417 407159 1 3256 23 35250 418 290476 0 1135 20 62285 419 275311 1 1270 24 69206 420 246541 1 661 1 65920 421 253468 1 1013 1 69770 422 240897 1 2844 74 72683 423 -83265 0 11528 68 -14545 424 -42143 1 6526 20 55830 425 272713 1 2264 20 55174 426 215362 0 5109 82 67038 427 42754 1 3999 21 51252 428 306275 1 35624 244 157278 429 253537 1 9252 32 79510 430 372631 1 15236 86 77440 431 -7170 1 18073 69 27284 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) group costs orders dividends 204815.678 -47556.317 22.044 950.981 1.518 > (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) 2.048e+05 3.849e+04 5.322 1.66e-07 *** group -4.756e+04 2.793e+04 -1.703 0.089345 . costs 2.204e+01 2.150e+00 10.253 < 2e-16 *** orders 9.510e+02 3.506e+02 2.713 0.006947 ** dividends 1.518e+00 4.116e-01 3.689 0.000254 *** --- 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/11m531291219010.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/21m531291219010.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/31m531291219010.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/4tv4o1291219010.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/5tv4o1291219010.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/64m3r1291219010.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/74m3r1291219010.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/8fw2c1291219010.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/9fw2c1291219010.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/1085kx1291219010.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/11t50l1291219010.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/12w6h81291219010.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/13byeh1291219010.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/14egd51291219010.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/15hhbb1291219010.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/163hsh1291219010.tab") + } > > try(system("convert tmp/11m531291219010.ps tmp/11m531291219010.png",intern=TRUE)) character(0) > try(system("convert tmp/21m531291219010.ps tmp/21m531291219010.png",intern=TRUE)) character(0) > try(system("convert tmp/31m531291219010.ps tmp/31m531291219010.png",intern=TRUE)) character(0) > try(system("convert tmp/4tv4o1291219010.ps tmp/4tv4o1291219010.png",intern=TRUE)) character(0) > try(system("convert tmp/5tv4o1291219010.ps tmp/5tv4o1291219010.png",intern=TRUE)) character(0) > try(system("convert tmp/64m3r1291219010.ps tmp/64m3r1291219010.png",intern=TRUE)) character(0) > try(system("convert tmp/74m3r1291219010.ps tmp/74m3r1291219010.png",intern=TRUE)) character(0) > try(system("convert tmp/8fw2c1291219010.ps tmp/8fw2c1291219010.png",intern=TRUE)) character(0) > try(system("convert tmp/9fw2c1291219010.ps tmp/9fw2c1291219010.png",intern=TRUE)) character(0) > try(system("convert tmp/1085kx1291219010.ps tmp/1085kx1291219010.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 12.354 2.383 27.160