R version 2.8.0 (2008-10-20) Copyright (C) 2008 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. Natural language support but running in an English locale 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 + ,1081 + ,213118 + ,6282154 + ,0 + ,29790 + ,309 + ,81767 + ,4321023 + ,0 + ,87550 + ,458 + ,153198 + ,4111912 + ,1 + ,84738 + ,588 + ,-26007 + ,223193 + ,0 + ,54660 + ,302 + ,126942 + ,1491348 + ,0 + ,42634 + ,156 + ,157214 + ,1629616 + ,1 + ,40949 + ,481 + ,129352 + ,1398893 + ,0 + ,45187 + ,353 + ,234817 + ,1926517 + ,0 + ,37704 + ,452 + ,60448 + ,983660 + ,0 + ,16275 + ,109 + ,47818 + ,1443586 + ,1 + ,25830 + ,115 + ,245546 + ,1073089 + ,1 + ,12679 + ,110 + ,48020 + ,984885 + ,0 + ,18014 + ,239 + ,-1710 + ,1405225 + ,1 + ,43556 + ,247 + ,32648 + ,227132 + ,0 + ,24811 + ,505 + ,95350 + ,929118 + ,1 + ,6575 + ,159 + ,151352 + ,1071292 + ,1 + ,7123 + ,109 + ,288170 + ,638830 + ,0 + ,21950 + ,519 + ,114337 + ,856956 + ,0 + ,37597 + ,248 + ,37884 + ,992426 + ,1 + ,17821 + ,373 + ,122844 + ,444477 + ,0 + ,12988 + ,119 + ,82340 + ,857217 + ,0 + ,22330 + ,84 + ,79801 + ,711969 + ,1 + ,13326 + ,102 + ,165548 + ,702380 + ,1 + ,16189 + ,295 + ,116384 + ,358589 + ,1 + ,7146 + ,105 + ,134028 + ,297978 + ,1 + ,15824 + ,64 + ,63838 + ,585715 + ,0 + ,27664 + ,282 + ,74996 + ,657954 + ,1 + ,11920 + ,182 + ,31080 + ,209458 + ,1 + ,8568 + ,37 + ,32168 + ,786690 + ,1 + ,14416 + ,361 + ,49857 + ,439798 + ,0 + ,3369 + ,28 + ,87161 + ,688779 + ,0 + ,11819 + ,85 + ,106113 + ,574339 + ,0 + ,6984 + ,45 + ,80570 + ,741409 + ,0 + ,4519 + ,49 + ,102129 + ,597793 + ,1 + ,2220 + ,22 + ,301670 + ,644190 + ,1 + ,18562 + ,155 + ,102313 + ,377934 + ,1 + ,10327 + ,91 + ,88577 + ,640273 + ,0 + ,5336 + ,81 + ,112477 + ,697458 + ,0 + ,2365 + ,79 + ,191778 + ,550608 + ,1 + ,4069 + ,145 + ,79804 + ,207393 + ,1 + ,8636 + ,855 + ,128294 + ,301607 + ,1 + ,13718 + ,61 + ,96448 + ,345783 + ,1 + ,4525 + ,226 + ,93811 + ,501749 + ,1 + ,6869 + ,105 + ,117520 + ,379983 + ,1 + ,4628 + ,62 + ,69159 + ,387475 + ,0 + ,3689 + ,25 + ,101792 + ,377305 + ,0 + ,4891 + ,217 + ,210568 + ,370837 + ,0 + ,7489 + ,322 + ,136996 + ,430866 + ,1 + ,4901 + ,84 + ,121920 + ,469107 + ,1 + ,2284 + ,33 + ,76403 + ,194493 + ,0 + ,3160 + ,108 + ,108094 + ,530670 + ,0 + ,4150 + ,150 + ,134759 + ,518365 + ,0 + ,7285 + ,115 + ,188873 + ,491303 + ,0 + ,1134 + ,162 + ,146216 + ,527021 + ,0 + ,4658 + ,158 + ,156608 + ,233773 + ,1 + ,2384 + ,97 + ,61348 + ,405972 + ,1 + ,3748 + ,9 + ,50350 + ,652925 + ,1 + ,5371 + ,66 + ,87720 + ,446211 + ,1 + ,1285 + ,107 + ,99489 + ,341340 + ,0 + ,9327 + ,101 + ,87419 + ,387699 + ,0 + ,5565 + ,47 + ,94355 + ,493408 + ,1 + ,1528 + ,38 + ,60326 + ,146494 + ,0 + ,3122 + ,34 + ,94670 + ,414462 + ,0 + ,7561 + ,87 + ,82425 + ,364304 + ,1 + ,2675 + ,79 + ,59017 + ,355178 + ,1 + ,13253 + ,947 + ,90829 + ,357760 + ,1 + ,880 + ,74 + ,80791 + ,261216 + ,0 + ,2053 + ,53 + ,100423 + ,397144 + ,1 + ,1424 + ,94 + ,131116 + ,374943 + ,0 + ,4036 + ,63 + ,100269 + ,424898 + ,0 + ,3045 + ,58 + ,27330 + ,202055 + ,1 + ,5119 + ,49 + ,39039 + ,378525 + ,1 + ,1431 + ,34 + ,106885 + ,310768 + ,1 + ,554 + ,11 + ,79285 + ,325738 + ,1 + ,1975 + ,35 + ,118881 + ,394510 + ,0 + ,1765 + ,20 + ,77623 + ,247060 + ,1 + ,1012 + ,47 + ,114768 + ,368078 + ,1 + ,810 + ,43 + ,74015 + ,236761 + ,1 + ,1280 + ,117 + ,69465 + ,312378 + ,0 + ,666 + ,171 + ,117869 + ,339836 + ,1 + ,1380 + ,26 + ,60982 + ,347385 + ,0 + ,4677 + ,75 + ,90131 + ,426280 + ,1 + ,876 + ,59 + ,138971 + ,352850 + ,1 + ,814 + ,18 + ,39625 + ,301881 + ,1 + ,514 + ,15 + ,102725 + ,377516 + ,0 + ,5692 + ,72 + ,64239 + ,357312 + ,1 + ,3642 + ,86 + ,90262 + ,458343 + ,1 + ,540 + ,14 + ,103960 + ,354228 + ,1 + ,2099 + ,64 + ,106611 + ,308636 + ,1 + ,567 + ,11 + ,103345 + ,386212 + ,1 + ,2001 + ,52 + ,95551 + ,393343 + ,0 + ,2949 + ,41 + ,82903 + ,378509 + ,1 + ,2253 + ,99 + ,63593 + ,452469 + ,0 + ,6533 + ,75 + ,126910 + ,364839 + ,1 + ,1889 + ,45 + ,37527 + ,358649 + ,0 + ,3055 + ,43 + ,60247 + ,376641 + ,1 + ,272 + ,8 + ,112995 + ,429112 + ,0 + ,1414 + ,198 + ,70184 + ,330546 + ,1 + ,2564 + ,22 + ,130140 + ,403560 + ,0 + ,1383 + ,11 + ,73221 + ,317892 + ,0 + ,1261 + ,33 + ,76114 + ,307528 + ,1 + ,975 + ,23 + ,90534 + ,235133 + ,1 + ,3366 + ,80 + ,108479 + ,299243 + ,1 + ,576 + ,18 + ,113761 + ,314073 + ,1 + ,1686 + ,40 + ,68696 + ,368186 + ,1 + ,746 + ,23 + ,71561 + ,269661 + ,0 + ,3192 + ,60 + ,59831 + ,125390 + ,0 + ,2045 + ,20 + ,97890 + ,510834 + ,1 + ,5702 + ,61 + ,101481 + ,321896 + ,0 + ,1932 + ,36 + ,72954 + ,249898 + ,1 + ,936 + ,30 + ,67939 + ,408881 + ,0 + ,3437 + ,47 + ,48022 + ,158492 + ,1 + ,5131 + ,71 + ,86111 + ,292154 + ,0 + ,2397 + ,14 + ,74020 + ,289513 + ,0 + ,1389 + ,9 + ,57530 + ,378049 + ,1 + ,1503 + ,39 + ,56364 + ,343466 + ,1 + ,402 + ,26 + ,84990 + ,332743 + ,1 + ,2239 + ,21 + ,88590 + ,442882 + ,0 + ,2234 + ,16 + ,77200 + ,214215 + ,1 + ,837 + ,69 + ,61262 + ,315688 + ,1 + ,10579 + ,92 + ,110309 + ,375195 + ,1 + ,875 + ,14 + ,67000 + ,334280 + ,1 + ,1585 + ,107 + ,93099 + ,355864 + ,0 + ,1659 + ,29 + ,107577 + ,480382 + ,0 + ,2647 + ,37 + ,62920 + ,353058 + ,0 + ,3294 + ,23 + ,75832 + ,217193 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,94 + ,7 + ,60793 + ,314533 + ,1 + ,422 + ,28 + ,57935 + ,318056 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,34 + ,8 + ,60630 + ,314353 + ,1 + ,1558 + ,63 + ,55637 + ,369448 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,43 + ,3 + ,60887 + ,312846 + ,1 + ,645 + ,5 + ,60720 + ,312075 + ,1 + ,316 + ,9 + ,60505 + ,315009 + ,1 + ,115 + ,13 + ,60945 + ,318903 + ,1 + ,5 + ,2 + ,60720 + ,314887 + ,1 + ,897 + ,5 + ,60720 + ,314913 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,389 + ,14 + ,58990 + ,325506 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,1002 + ,15 + ,56750 + ,298568 + ,1 + ,36 + ,3 + ,60894 + ,315834 + ,1 + ,460 + ,15 + ,63346 + ,329784 + ,1 + ,309 + ,11 + ,56535 + ,312878 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,9 + ,6 + ,60835 + ,314987 + ,1 + ,271 + ,2 + ,60720 + ,325249 + ,1 + ,14 + ,1 + ,61016 + ,315877 + ,1 + ,520 + ,10 + ,58650 + ,291650 + ,1 + ,1766 + ,73 + ,60438 + ,305959 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,458 + ,11 + ,58625 + ,297765 + ,1 + ,20 + ,3 + ,60938 + ,315245 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,98 + ,2 + ,61490 + ,315236 + ,1 + ,405 + ,7 + ,60845 + ,336425 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,483 + ,27 + ,60830 + ,306268 + ,1 + ,454 + ,51 + ,63261 + ,302187 + ,1 + ,47 + ,3 + ,60720 + ,314882 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,757 + ,19 + ,45689 + ,382712 + ,1 + ,4655 + ,393 + ,60720 + ,341570 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,36 + ,4 + ,61564 + ,312412 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,203 + ,9 + ,61938 + ,309596 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,126 + ,10 + ,60951 + ,315547 + ,1 + ,400 + ,152 + ,60720 + ,313267 + ,0 + ,71 + ,1 + ,60745 + ,316176 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,972 + ,34 + ,71642 + ,359335 + ,0 + ,531 + ,10 + ,71641 + ,330068 + ,1 + ,2461 + ,57 + ,55792 + ,314289 + ,0 + ,378 + ,52 + ,71873 + ,297413 + ,0 + ,23 + ,5 + ,62555 + ,314806 + ,0 + ,638 + ,14 + ,60370 + ,333210 + ,0 + ,2300 + ,29 + ,64873 + ,352108 + ,1 + ,149 + ,5 + ,62041 + ,313332 + ,1 + ,226 + ,5 + ,65745 + ,291787 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,275 + ,4 + ,59500 + ,318745 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,141 + ,6 + ,61630 + ,315366 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,28 + ,2 + ,60890 + ,315688 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,4980 + ,91 + ,113521 + ,409642 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,472 + ,20 + ,80045 + ,269587 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,203 + ,27 + ,50804 + ,300962 + ,0 + ,496 + ,17 + ,87390 + ,325479 + ,0 + ,10 + ,2 + ,61656 + ,316155 + ,0 + ,63 + ,4 + ,65688 + ,318574 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1136 + ,32 + ,48522 + ,343613 + ,0 + ,265 + ,31 + ,60720 + ,306948 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,267 + ,32 + ,57640 + ,330059 + ,1 + ,474 + ,20 + ,61977 + ,288985 + ,1 + ,534 + ,7 + ,62620 + ,304485 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,15 + ,8 + ,60831 + ,315688 + ,1 + ,397 + ,28 + ,60646 + ,317736 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,1866 + ,29 + ,56225 + ,322331 + ,1 + ,288 + ,4 + ,60510 + ,296656 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,3 + ,2 + ,60698 + ,315354 + ,1 + ,468 + ,21 + ,60720 + ,312161 + ,1 + ,20 + ,2 + ,60805 + ,315576 + ,1 + ,278 + ,26 + ,61404 + ,314922 + ,1 + ,61 + ,14 + ,60720 + ,314551 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,192 + ,4 + ,65276 + ,312339 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,317 + ,9 + ,63915 + ,298700 + ,1 + ,738 + ,10 + ,60720 + ,321376 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,368 + ,17 + ,61686 + ,303230 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,2 + ,1 + ,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 + ,8 + ,59818 + ,315637 + ,1 + ,2332 + ,4 + ,72680 + ,324385 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,131 + ,11 + ,61808 + ,308989 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,206 + ,9 + ,53110 + ,296702 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,167 + ,2 + ,64245 + ,307322 + ,1 + ,622 + ,73 + ,73007 + ,304376 + ,1 + ,2328 + ,94 + ,82732 + ,253588 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,365 + ,8 + ,54820 + ,309560 + ,1 + ,364 + ,35 + ,47705 + ,298466 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,226 + ,12 + ,72835 + ,343929 + ,1 + ,307 + ,15 + ,58856 + ,331955 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,188 + ,11 + ,77655 + ,381180 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,138 + ,6 + ,69817 + ,331420 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,125 + ,12 + ,60798 + ,310201 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,282 + ,30 + ,62452 + ,320016 + ,1 + ,335 + ,33 + ,64175 + ,320398 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1324 + ,117 + ,67440 + ,291841 + ,1 + ,176 + ,28 + ,68136 + ,310670 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,249 + ,72 + ,56726 + ,313491 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,333 + ,13 + ,70811 + ,331323 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,601 + ,6 + ,60720 + ,319210 + ,1 + ,30 + ,4 + ,62045 + ,318098 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,249 + ,62 + ,54323 + ,292754 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,165 + ,24 + ,62841 + ,325176 + ,1 + ,453 + ,21 + ,81125 + ,365959 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,53 + ,14 + ,59506 + ,302409 + ,0 + ,382 + ,21 + ,59365 + ,340968 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,30 + ,4 + ,60798 + ,313164 + ,1 + ,290 + ,2 + ,58790 + ,301164 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,366 + ,53 + ,61808 + ,344425 + ,1 + ,2 + ,9 + ,60735 + ,315394 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,209 + ,13 + ,64016 + ,316647 + ,1 + ,384 + ,22 + ,54683 + ,309836 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,365 + ,83 + ,87192 + ,346611 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,49 + ,8 + ,64107 + ,322031 + ,1 + ,3 + ,4 + ,60761 + ,315656 + ,1 + ,133 + ,14 + ,65990 + ,339445 + ,1 + ,32 + ,1 + ,59988 + ,314964 + ,1 + ,368 + ,17 + ,61167 + ,297141 + ,1 + ,1 + ,6 + ,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 + ,2 + ,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 + ,5 + ,60379 + ,313729 + ,1 + ,1 + ,2 + ,60727 + ,315388 + ,1 + ,314 + ,5 + ,60720 + ,315371 + ,1 + ,844 + ,78 + ,60925 + ,296139 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,26 + ,1 + ,60896 + ,313880 + ,1 + ,125 + ,13 + ,59734 + ,317698 + ,1 + ,304 + ,15 + ,62969 + ,295580 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,621 + ,48 + ,60720 + ,308256 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,119 + ,6 + ,59118 + ,303677 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,1595 + ,17 + ,60720 + ,319369 + ,1 + ,312 + ,14 + ,58598 + ,318690 + ,1 + ,60 + ,10 + ,61124 + ,314049 + ,1 + ,587 + ,12 + ,59595 + ,325699 + ,1 + ,135 + ,2 + ,62065 + ,314210 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,514 + ,52 + ,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 + ,24 + ,61600 + ,308336 + ,1 + ,180 + ,11 + ,59635 + ,316386 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,218 + ,21 + ,60720 + ,315553 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,448 + ,40 + ,59781 + ,323361 + ,1 + ,227 + ,9 + ,76644 + ,336639 + ,1 + ,174 + ,1 + ,64820 + ,307424 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,121 + ,24 + ,56178 + ,295370 + ,1 + ,607 + ,11 + ,60436 + ,322340 + ,1 + ,2212 + ,14 + ,60720 + ,319864 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,530 + ,60 + ,73433 + ,317291 + ,1 + ,571 + ,80 + ,41477 + ,280398 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,78 + ,16 + ,62700 + ,317330 + ,1 + ,2489 + ,40 + ,67804 + ,238125 + ,1 + ,131 + ,6 + ,59661 + ,327071 + ,1 + ,923 + ,8 + ,58620 + ,309038 + ,1 + ,72 + ,3 + ,60398 + ,314210 + ,1 + ,572 + ,16 + ,58580 + ,307930 + ,1 + ,397 + ,10 + ,62710 + ,322327 + ,1 + ,450 + ,8 + ,59325 + ,292136 + ,1 + ,622 + ,7 + ,60950 + ,263276 + ,1 + ,694 + ,8 + ,68060 + ,367655 + ,0 + ,3425 + ,12 + ,83620 + ,283910 + ,1 + ,562 + ,13 + ,58456 + ,283587 + ,1 + ,4917 + ,42 + ,52811 + ,243650 + ,0 + ,1442 + ,118 + ,121173 + ,438493 + ,1 + ,529 + ,9 + ,63870 + ,296261 + ,0 + ,2126 + ,138 + ,21001 + ,230621 + ,1 + ,1061 + ,5 + ,70415 + ,304252 + ,1 + ,776 + ,9 + ,64230 + ,333505 + ,1 + ,611 + ,8 + ,59190 + ,296919 + ,0 + ,1526 + ,25 + ,69351 + ,278990 + ,1 + ,592 + ,7 + ,64270 + ,276898 + ,1 + ,1182 + ,13 + ,70694 + ,327007 + ,1 + ,621 + ,16 + ,68005 + ,317046 + ,1 + ,989 + ,11 + ,58930 + ,304555 + ,1 + ,438 + ,11 + ,58320 + ,298096 + ,1 + ,726 + ,3 + ,69980 + ,231861 + ,1 + ,1303 + ,61 + ,69863 + ,309422 + ,0 + ,7419 + ,29 + ,63255 + ,286963 + ,0 + ,1164 + ,17 + ,57320 + ,269753 + ,0 + ,3310 + ,33 + ,75230 + ,448243 + ,1 + ,1920 + ,15 + ,79420 + ,165404 + ,1 + ,965 + ,3 + ,73490 + ,204325 + ,1 + ,3256 + ,66 + ,35250 + ,407159 + ,0 + ,1135 + ,17 + ,62285 + ,290476 + ,1 + ,1270 + ,26 + ,69206 + ,275311 + ,1 + ,661 + ,3 + ,65920 + ,246541 + ,1 + ,1013 + ,2 + ,69770 + ,253468 + ,1 + ,2844 + ,67 + ,72683 + ,240897 + ,0 + ,11528 + ,70 + ,-14545 + ,-83265 + ,1 + ,6526 + ,26 + ,55830 + ,-42143 + ,1 + ,2264 + ,24 + ,55174 + ,272713 + ,0 + ,5109 + ,97 + ,67038 + ,215362 + ,1 + ,3999 + ,30 + ,51252 + ,42754 + ,1 + ,35624 + ,223 + ,157278 + ,306275 + ,1 + ,9252 + ,48 + ,79510 + ,253537 + ,1 + ,15236 + ,90 + ,77440 + ,372631 + ,1 + ,18073 + ,180 + ,27284 + ,-7170) + ,dim=c(5 + ,431) + ,dimnames=list(c('Group' + ,'Costs' + ,'Trades' + ,'Dividends' + ,'Wealth') + ,1:431)) > y <- array(NA,dim=c(5,431),dimnames=list(c('Group','Costs','Trades','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 Group Costs Trades Dividends 1 6282154 0 162556 1081 213118 2 4321023 0 29790 309 81767 3 4111912 0 87550 458 153198 4 223193 1 84738 588 -26007 5 1491348 0 54660 302 126942 6 1629616 0 42634 156 157214 7 1398893 1 40949 481 129352 8 1926517 0 45187 353 234817 9 983660 0 37704 452 60448 10 1443586 0 16275 109 47818 11 1073089 1 25830 115 245546 12 984885 1 12679 110 48020 13 1405225 0 18014 239 -1710 14 227132 1 43556 247 32648 15 929118 0 24811 505 95350 16 1071292 1 6575 159 151352 17 638830 1 7123 109 288170 18 856956 0 21950 519 114337 19 992426 0 37597 248 37884 20 444477 1 17821 373 122844 21 857217 0 12988 119 82340 22 711969 0 22330 84 79801 23 702380 1 13326 102 165548 24 358589 1 16189 295 116384 25 297978 1 7146 105 134028 26 585715 1 15824 64 63838 27 657954 0 27664 282 74996 28 209458 1 11920 182 31080 29 786690 1 8568 37 32168 30 439798 1 14416 361 49857 31 688779 0 3369 28 87161 32 574339 0 11819 85 106113 33 741409 0 6984 45 80570 34 597793 0 4519 49 102129 35 644190 1 2220 22 301670 36 377934 1 18562 155 102313 37 640273 1 10327 91 88577 38 697458 0 5336 81 112477 39 550608 0 2365 79 191778 40 207393 1 4069 145 79804 41 301607 1 8636 855 128294 42 345783 1 13718 61 96448 43 501749 1 4525 226 93811 44 379983 1 6869 105 117520 45 387475 1 4628 62 69159 46 377305 0 3689 25 101792 47 370837 0 4891 217 210568 48 430866 0 7489 322 136996 49 469107 1 4901 84 121920 50 194493 1 2284 33 76403 51 530670 0 3160 108 108094 52 518365 0 4150 150 134759 53 491303 0 7285 115 188873 54 527021 0 1134 162 146216 55 233773 0 4658 158 156608 56 405972 1 2384 97 61348 57 652925 1 3748 9 50350 58 446211 1 5371 66 87720 59 341340 1 1285 107 99489 60 387699 0 9327 101 87419 61 493408 0 5565 47 94355 62 146494 1 1528 38 60326 63 414462 0 3122 34 94670 64 364304 0 7561 87 82425 65 355178 1 2675 79 59017 66 357760 1 13253 947 90829 67 261216 1 880 74 80791 68 397144 0 2053 53 100423 69 374943 1 1424 94 131116 70 424898 0 4036 63 100269 71 202055 0 3045 58 27330 72 378525 1 5119 49 39039 73 310768 1 1431 34 106885 74 325738 1 554 11 79285 75 394510 1 1975 35 118881 76 247060 0 1765 20 77623 77 368078 1 1012 47 114768 78 236761 1 810 43 74015 79 312378 1 1280 117 69465 80 339836 0 666 171 117869 81 347385 1 1380 26 60982 82 426280 0 4677 75 90131 83 352850 1 876 59 138971 84 301881 1 814 18 39625 85 377516 1 514 15 102725 86 357312 0 5692 72 64239 87 458343 1 3642 86 90262 88 354228 1 540 14 103960 89 308636 1 2099 64 106611 90 386212 1 567 11 103345 91 393343 1 2001 52 95551 92 378509 0 2949 41 82903 93 452469 1 2253 99 63593 94 364839 0 6533 75 126910 95 358649 1 1889 45 37527 96 376641 0 3055 43 60247 97 429112 1 272 8 112995 98 330546 0 1414 198 70184 99 403560 1 2564 22 130140 100 317892 0 1383 11 73221 101 307528 0 1261 33 76114 102 235133 1 975 23 90534 103 299243 1 3366 80 108479 104 314073 1 576 18 113761 105 368186 1 1686 40 68696 106 269661 1 746 23 71561 107 125390 0 3192 60 59831 108 510834 0 2045 20 97890 109 321896 1 5702 61 101481 110 249898 0 1932 36 72954 111 408881 1 936 30 67939 112 158492 0 3437 47 48022 113 292154 1 5131 71 86111 114 289513 0 2397 14 74020 115 378049 0 1389 9 57530 116 343466 1 1503 39 56364 117 332743 1 402 26 84990 118 442882 1 2239 21 88590 119 214215 0 2234 16 77200 120 315688 1 837 69 61262 121 375195 1 10579 92 110309 122 334280 1 875 14 67000 123 355864 1 1585 107 93099 124 480382 0 1659 29 107577 125 353058 0 2647 37 62920 126 217193 0 3294 23 75832 127 315380 1 0 0 60720 128 314533 1 94 7 60793 129 318056 1 422 28 57935 130 315380 1 0 0 60720 131 314353 1 34 8 60630 132 369448 1 1558 63 55637 133 315380 1 0 0 60720 134 312846 1 43 3 60887 135 312075 1 645 5 60720 136 315009 1 316 9 60505 137 318903 1 115 13 60945 138 314887 1 5 2 60720 139 314913 1 897 5 60720 140 315380 1 0 0 60720 141 325506 1 389 14 58990 142 315380 1 0 0 60720 143 298568 1 1002 15 56750 144 315834 1 36 3 60894 145 329784 1 460 15 63346 146 312878 1 309 11 56535 147 315380 1 0 0 60720 148 314987 1 9 6 60835 149 325249 1 271 2 60720 150 315877 1 14 1 61016 151 291650 1 520 10 58650 152 305959 1 1766 73 60438 153 315380 1 0 0 60720 154 297765 1 458 11 58625 155 315245 1 20 3 60938 156 315380 1 0 0 60720 157 315380 1 0 0 60720 158 315236 1 98 2 61490 159 336425 1 405 7 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 27 60830 165 302187 1 454 51 63261 166 314882 1 47 3 60720 167 315380 1 0 0 60720 168 382712 1 757 19 45689 169 341570 1 4655 393 60720 170 315380 1 0 0 60720 171 315380 1 0 0 60720 172 312412 1 36 4 61564 173 315380 1 0 0 60720 174 309596 1 203 9 61938 175 315380 1 0 0 60720 176 315547 1 126 10 60951 177 313267 1 400 152 60720 178 316176 0 71 1 60745 179 315380 1 0 0 60720 180 315380 1 0 0 60720 181 359335 1 972 34 71642 182 330068 0 531 10 71641 183 314289 1 2461 57 55792 184 297413 0 378 52 71873 185 314806 0 23 5 62555 186 333210 0 638 14 60370 187 352108 0 2300 29 64873 188 313332 1 149 5 62041 189 291787 1 226 5 65745 190 315380 1 0 0 60720 191 318745 1 275 4 59500 192 315380 1 0 0 60720 193 315366 1 141 6 61630 194 315380 1 0 0 60720 195 315688 1 28 2 60890 196 315380 0 0 0 60720 197 409642 0 4980 91 113521 198 315380 0 0 0 60720 199 315380 0 0 0 60720 200 269587 0 472 20 80045 201 315380 0 0 0 60720 202 315380 0 0 0 60720 203 315380 0 0 0 60720 204 300962 0 203 27 50804 205 325479 0 496 17 87390 206 316155 0 10 2 61656 207 318574 0 63 4 65688 208 315380 0 0 0 60720 209 343613 0 1136 32 48522 210 306948 0 265 31 60720 211 315380 1 0 0 60720 212 315380 1 0 0 60720 213 330059 1 267 32 57640 214 288985 1 474 20 61977 215 304485 1 534 7 62620 216 315380 1 0 0 60720 217 315688 1 15 8 60831 218 317736 1 397 28 60646 219 315380 1 0 0 60720 220 322331 1 1866 29 56225 221 296656 1 288 4 60510 222 315380 1 0 0 60720 223 315354 1 3 2 60698 224 312161 1 468 21 60720 225 315576 1 20 2 60805 226 314922 1 278 26 61404 227 314551 1 61 14 60720 228 315380 1 0 0 60720 229 312339 1 192 4 65276 230 315380 1 0 0 60720 231 298700 1 317 9 63915 232 321376 1 738 10 60720 233 315380 1 0 0 60720 234 303230 1 368 17 61686 235 315380 1 0 0 60720 236 315487 1 2 1 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 8 59818 245 324385 1 2332 4 72680 246 315380 0 0 0 60720 247 315380 0 0 0 60720 248 308989 1 131 11 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 2 64245 254 304376 1 622 73 73007 255 253588 1 2328 94 82732 256 315380 1 0 0 60720 257 309560 1 365 8 54820 258 298466 1 364 35 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 12 72835 264 331955 1 307 15 58856 265 315380 0 0 0 60720 266 315380 1 0 0 60720 267 315380 0 0 0 60720 268 381180 1 188 11 77655 269 315380 0 0 0 60720 270 331420 1 138 6 69817 271 315380 0 0 0 60720 272 315380 0 0 0 60720 273 315380 0 0 0 60720 274 310201 1 125 12 60798 275 315380 1 0 0 60720 276 320016 1 282 30 62452 277 320398 1 335 33 64175 278 315380 1 0 0 60720 279 291841 0 1324 117 67440 280 310670 1 176 28 68136 281 315380 0 0 0 60720 282 315380 1 0 0 60720 283 313491 1 249 72 56726 284 315380 1 0 0 60720 285 331323 1 333 13 70811 286 315380 1 0 0 60720 287 319210 0 601 6 60720 288 318098 1 30 4 62045 289 315380 0 0 0 60720 290 292754 1 249 62 54323 291 315380 0 0 0 60720 292 325176 1 165 24 62841 293 365959 1 453 21 81125 294 315380 1 0 0 60720 295 302409 1 53 14 59506 296 340968 0 382 21 59365 297 315380 1 0 0 60720 298 315380 0 0 0 60720 299 315380 1 0 0 60720 300 315380 1 0 0 60720 301 313164 0 30 4 60798 302 301164 1 290 2 58790 303 315380 1 0 0 60720 304 315380 0 0 0 60720 305 344425 1 366 53 61808 306 315394 1 2 9 60735 307 315380 0 0 0 60720 308 316647 0 209 13 64016 309 309836 1 384 22 54683 310 315380 0 0 0 60720 311 315380 0 0 0 60720 312 346611 1 365 83 87192 313 315380 0 0 0 60720 314 322031 0 49 8 64107 315 315656 1 3 4 60761 316 339445 1 133 14 65990 317 314964 1 32 1 59988 318 297141 1 368 17 61167 319 315372 1 1 6 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 2 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 5 60379 335 315388 1 1 2 60727 336 315371 1 314 5 60720 337 296139 1 844 78 60925 338 315380 1 0 0 60720 339 313880 1 26 1 60896 340 317698 1 125 13 59734 341 295580 1 304 15 62969 342 315380 1 0 0 60720 343 315380 1 0 0 60720 344 315380 1 0 0 60720 345 308256 1 621 48 60720 346 315380 1 0 0 60720 347 303677 1 119 6 59118 348 315380 1 0 0 60720 349 315380 1 0 0 60720 350 319369 1 1595 17 60720 351 318690 1 312 14 58598 352 314049 1 60 10 61124 353 325699 1 587 12 59595 354 314210 1 135 2 62065 355 315380 1 0 0 60720 356 315380 1 0 0 60720 357 322378 1 514 52 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 24 61600 365 316386 1 180 11 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 21 60720 371 315380 1 0 0 60720 372 323361 1 448 40 59781 373 336639 1 227 9 76644 374 307424 1 174 1 64820 375 315380 1 0 0 60720 376 315380 1 0 0 60720 377 295370 1 121 24 56178 378 322340 1 607 11 60436 379 319864 1 2212 14 60720 380 315380 1 0 0 60720 381 315380 1 0 0 60720 382 317291 1 530 60 73433 383 280398 1 571 80 41477 384 315380 1 0 0 60720 385 317330 1 78 16 62700 386 238125 1 2489 40 67804 387 327071 1 131 6 59661 388 309038 1 923 8 58620 389 314210 1 72 3 60398 390 307930 1 572 16 58580 391 322327 1 397 10 62710 392 292136 1 450 8 59325 393 263276 1 622 7 60950 394 367655 1 694 8 68060 395 283910 0 3425 12 83620 396 283587 1 562 13 58456 397 243650 1 4917 42 52811 398 438493 0 1442 118 121173 399 296261 1 529 9 63870 400 230621 0 2126 138 21001 401 304252 1 1061 5 70415 402 333505 1 776 9 64230 403 296919 1 611 8 59190 404 278990 0 1526 25 69351 405 276898 1 592 7 64270 406 327007 1 1182 13 70694 407 317046 1 621 16 68005 408 304555 1 989 11 58930 409 298096 1 438 11 58320 410 231861 1 726 3 69980 411 309422 1 1303 61 69863 412 286963 0 7419 29 63255 413 269753 0 1164 17 57320 414 448243 0 3310 33 75230 415 165404 1 1920 15 79420 416 204325 1 965 3 73490 417 407159 1 3256 66 35250 418 290476 0 1135 17 62285 419 275311 1 1270 26 69206 420 246541 1 661 3 65920 421 253468 1 1013 2 69770 422 240897 1 2844 67 72683 423 -83265 0 11528 70 -14545 424 -42143 1 6526 26 55830 425 272713 1 2264 24 55174 426 215362 0 5109 97 67038 427 42754 1 3999 30 51252 428 306275 1 35624 223 157278 429 253537 1 9252 48 79510 430 372631 1 15236 90 77440 431 -7170 1 18073 180 27284 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Group Costs Trades Dividends 205978.971 -55485.819 29.811 -407.150 1.926 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -2163938 -14344 35201 47931 3195284 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 2.060e+05 3.854e+04 5.345 1.48e-07 *** Group -5.549e+04 2.784e+04 -1.993 0.0469 * Costs 2.981e+01 1.578e+00 18.889 < 2e-16 *** Trades -4.072e+02 1.663e+02 -2.448 0.0148 * Dividends 1.926e+00 4.105e-01 4.692 3.65e-06 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 253400 on 426 degrees of freedom Multiple R-squared: 0.6569, Adjusted R-squared: 0.6537 F-statistic: 203.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 6.580461e-56 3.290230e-56 [2,] 1.0000000 9.697449e-93 4.848724e-93 [3,] 1.0000000 3.784523e-102 1.892262e-102 [4,] 1.0000000 2.398605e-104 1.199303e-104 [5,] 1.0000000 2.524261e-116 1.262131e-116 [6,] 1.0000000 2.796176e-132 1.398088e-132 [7,] 1.0000000 1.296946e-134 6.484732e-135 [8,] 1.0000000 1.043990e-142 5.219952e-143 [9,] 1.0000000 5.532265e-153 2.766133e-153 [10,] 1.0000000 1.188565e-152 5.942826e-153 [11,] 1.0000000 1.597306e-157 7.986528e-158 [12,] 1.0000000 1.622078e-164 8.110390e-165 [13,] 1.0000000 2.174270e-163 1.087135e-163 [14,] 1.0000000 4.837755e-169 2.418878e-169 [15,] 1.0000000 1.876296e-172 9.381478e-173 [16,] 1.0000000 5.172766e-173 2.586383e-173 [17,] 1.0000000 5.308774e-172 2.654387e-172 [18,] 1.0000000 1.870383e-171 9.351913e-172 [19,] 1.0000000 4.297350e-173 2.148675e-173 [20,] 1.0000000 1.219015e-177 6.095076e-178 [21,] 1.0000000 5.159968e-177 2.579984e-177 [22,] 1.0000000 6.552292e-190 3.276146e-190 [23,] 1.0000000 1.649405e-190 8.247026e-191 [24,] 1.0000000 1.876207e-194 9.381033e-195 [25,] 1.0000000 1.330203e-196 6.651013e-197 [26,] 1.0000000 3.554703e-204 1.777351e-204 [27,] 1.0000000 2.959186e-206 1.479593e-206 [28,] 1.0000000 1.866239e-205 9.331196e-206 [29,] 1.0000000 1.143810e-204 5.719049e-205 [30,] 1.0000000 1.635797e-210 8.178985e-211 [31,] 1.0000000 2.790585e-216 1.395292e-216 [32,] 1.0000000 2.423065e-216 1.211532e-216 [33,] 1.0000000 3.344570e-216 1.672285e-216 [34,] 1.0000000 6.743693e-216 3.371847e-216 [35,] 1.0000000 5.618652e-215 2.809326e-215 [36,] 1.0000000 6.314230e-216 3.157115e-216 [37,] 1.0000000 7.510434e-215 3.755217e-215 [38,] 1.0000000 1.254756e-214 6.273779e-215 [39,] 1.0000000 4.183044e-214 2.091522e-214 [40,] 1.0000000 6.006534e-217 3.003267e-217 [41,] 1.0000000 1.506022e-216 7.530108e-217 [42,] 1.0000000 4.509424e-216 2.254712e-216 [43,] 1.0000000 3.030670e-216 1.515335e-216 [44,] 1.0000000 6.898805e-217 3.449402e-217 [45,] 1.0000000 8.546720e-217 4.273360e-217 [46,] 1.0000000 1.213618e-216 6.068091e-217 [47,] 1.0000000 4.108265e-216 2.054132e-216 [48,] 1.0000000 1.821504e-219 9.107519e-220 [49,] 1.0000000 1.134672e-219 5.673360e-220 [50,] 1.0000000 2.818042e-232 1.409021e-232 [51,] 1.0000000 1.096219e-232 5.481096e-233 [52,] 1.0000000 8.301183e-232 4.150591e-232 [53,] 1.0000000 8.007338e-232 4.003669e-232 [54,] 1.0000000 2.876362e-233 1.438181e-233 [55,] 1.0000000 1.474903e-234 7.374514e-235 [56,] 1.0000000 5.924023e-234 2.962012e-234 [57,] 1.0000000 1.467914e-233 7.339571e-234 [58,] 1.0000000 4.940991e-233 2.470495e-233 [59,] 1.0000000 5.649985e-232 2.824993e-232 [60,] 1.0000000 1.276709e-231 6.383544e-232 [61,] 1.0000000 1.167663e-230 5.838316e-231 [62,] 1.0000000 7.212219e-230 3.606110e-230 [63,] 1.0000000 2.128623e-229 1.064312e-229 [64,] 1.0000000 1.713134e-228 8.565670e-229 [65,] 1.0000000 1.227108e-229 6.135539e-230 [66,] 1.0000000 4.201118e-229 2.100559e-229 [67,] 1.0000000 3.799274e-228 1.899637e-228 [68,] 1.0000000 4.250804e-227 2.125402e-227 [69,] 1.0000000 7.622571e-227 3.811286e-227 [70,] 1.0000000 6.129166e-226 3.064583e-226 [71,] 1.0000000 7.428973e-226 3.714486e-226 [72,] 1.0000000 6.530709e-225 3.265355e-225 [73,] 1.0000000 9.533083e-225 4.766542e-225 [74,] 1.0000000 4.964354e-224 2.482177e-224 [75,] 1.0000000 6.469603e-224 3.234801e-224 [76,] 1.0000000 4.289373e-224 2.144687e-224 [77,] 1.0000000 2.451482e-223 1.225741e-223 [78,] 1.0000000 2.370248e-222 1.185124e-222 [79,] 1.0000000 3.987638e-222 1.993819e-222 [80,] 1.0000000 1.416318e-222 7.081592e-223 [81,] 1.0000000 1.039626e-221 5.198128e-222 [82,] 1.0000000 2.049985e-221 1.024992e-221 [83,] 1.0000000 2.062041e-220 1.031020e-220 [84,] 1.0000000 1.847115e-219 9.235577e-220 [85,] 1.0000000 1.133793e-218 5.668964e-219 [86,] 1.0000000 2.454126e-220 1.227063e-220 [87,] 1.0000000 1.126815e-219 5.634075e-220 [88,] 1.0000000 3.228381e-220 1.614190e-220 [89,] 1.0000000 3.360509e-220 1.680254e-220 [90,] 1.0000000 3.273129e-219 1.636565e-219 [91,] 1.0000000 3.889547e-218 1.944773e-218 [92,] 1.0000000 3.315415e-217 1.657708e-217 [93,] 1.0000000 3.596537e-216 1.798269e-216 [94,] 1.0000000 3.249622e-215 1.624811e-215 [95,] 1.0000000 2.372608e-216 1.186304e-216 [96,] 1.0000000 1.928803e-216 9.644017e-217 [97,] 1.0000000 2.834180e-217 1.417090e-217 [98,] 1.0000000 1.308202e-216 6.541011e-217 [99,] 1.0000000 4.425533e-216 2.212766e-216 [100,] 1.0000000 1.770637e-218 8.853185e-219 [101,] 1.0000000 1.113593e-219 5.567967e-220 [102,] 1.0000000 8.177657e-219 4.088829e-219 [103,] 1.0000000 1.414773e-218 7.073863e-219 [104,] 1.0000000 1.055940e-218 5.279699e-219 [105,] 1.0000000 5.001559e-219 2.500779e-219 [106,] 1.0000000 3.706097e-218 1.853048e-218 [107,] 1.0000000 2.808432e-217 1.404216e-217 [108,] 1.0000000 2.568869e-217 1.284434e-217 [109,] 1.0000000 9.624234e-217 4.812117e-217 [110,] 1.0000000 6.720738e-216 3.360369e-216 [111,] 1.0000000 6.046519e-216 3.023259e-216 [112,] 1.0000000 4.683900e-217 2.341950e-217 [113,] 1.0000000 4.786938e-216 2.393469e-216 [114,] 1.0000000 2.798069e-215 1.399034e-215 [115,] 1.0000000 2.695036e-214 1.347518e-214 [116,] 1.0000000 2.225628e-213 1.112814e-213 [117,] 1.0000000 4.454487e-213 2.227244e-213 [118,] 1.0000000 1.515581e-212 7.577903e-213 [119,] 1.0000000 3.594020e-213 1.797010e-213 [120,] 1.0000000 3.632179e-212 1.816090e-212 [121,] 1.0000000 3.700289e-211 1.850145e-211 [122,] 1.0000000 3.446286e-210 1.723143e-210 [123,] 1.0000000 3.499927e-209 1.749963e-209 [124,] 1.0000000 3.566552e-208 1.783276e-208 [125,] 1.0000000 3.410846e-208 1.705423e-208 [126,] 1.0000000 3.497851e-207 1.748926e-207 [127,] 1.0000000 3.633494e-206 1.816747e-206 [128,] 1.0000000 3.805841e-205 1.902921e-205 [129,] 1.0000000 3.893032e-204 1.946516e-204 [130,] 1.0000000 3.898673e-203 1.949337e-203 [131,] 1.0000000 3.988873e-202 1.994436e-202 [132,] 1.0000000 4.012666e-201 2.006333e-201 [133,] 1.0000000 4.081896e-200 2.040948e-200 [134,] 1.0000000 3.488969e-199 1.744485e-199 [135,] 1.0000000 3.538448e-198 1.769224e-198 [136,] 1.0000000 3.719160e-197 1.859580e-197 [137,] 1.0000000 3.751240e-196 1.875620e-196 [138,] 1.0000000 3.462981e-195 1.731491e-195 [139,] 1.0000000 3.246167e-194 1.623083e-194 [140,] 1.0000000 3.253262e-193 1.626631e-193 [141,] 1.0000000 3.259631e-192 1.629815e-192 [142,] 1.0000000 2.928456e-191 1.464228e-191 [143,] 1.0000000 2.914539e-190 1.457269e-190 [144,] 1.0000000 2.941399e-189 1.470699e-189 [145,] 1.0000000 2.997652e-188 1.498826e-188 [146,] 1.0000000 2.951277e-187 1.475639e-187 [147,] 1.0000000 2.990784e-186 1.495392e-186 [148,] 1.0000000 2.933777e-185 1.466889e-185 [149,] 1.0000000 2.860569e-184 1.430285e-184 [150,] 1.0000000 2.780129e-183 1.390065e-183 [151,] 1.0000000 2.716643e-182 1.358322e-182 [152,] 1.0000000 1.942374e-181 9.711868e-182 [153,] 1.0000000 1.870041e-180 9.350207e-181 [154,] 1.0000000 1.794144e-179 8.970720e-180 [155,] 1.0000000 1.715259e-178 8.576296e-179 [156,] 1.0000000 1.633987e-177 8.169935e-178 [157,] 1.0000000 1.589681e-176 7.948406e-177 [158,] 1.0000000 1.446267e-175 7.231333e-176 [159,] 1.0000000 1.365745e-174 6.828725e-175 [160,] 1.0000000 1.282240e-173 6.411199e-174 [161,] 1.0000000 1.362101e-174 6.810505e-175 [162,] 1.0000000 4.685075e-174 2.342537e-174 [163,] 1.0000000 4.451372e-173 2.225686e-173 [164,] 1.0000000 4.213031e-172 2.106516e-172 [165,] 1.0000000 4.033156e-171 2.016578e-171 [166,] 1.0000000 3.787567e-170 1.893784e-170 [167,] 1.0000000 3.617722e-169 1.808861e-169 [168,] 1.0000000 3.370592e-168 1.685296e-168 [169,] 1.0000000 3.117500e-167 1.558750e-167 [170,] 1.0000000 2.696623e-166 1.348311e-166 [171,] 1.0000000 2.570902e-165 1.285451e-165 [172,] 1.0000000 2.357167e-164 1.178583e-164 [173,] 1.0000000 2.152478e-163 1.076239e-163 [174,] 1.0000000 1.404651e-162 7.023257e-163 [175,] 1.0000000 1.304927e-161 6.524635e-162 [176,] 1.0000000 7.202553e-161 3.601276e-161 [177,] 1.0000000 4.132589e-160 2.066295e-160 [178,] 1.0000000 3.843771e-159 1.921886e-159 [179,] 1.0000000 2.845878e-158 1.422939e-158 [180,] 1.0000000 1.165707e-157 5.828537e-158 [181,] 1.0000000 1.052981e-156 5.264906e-157 [182,] 1.0000000 7.604224e-156 3.802112e-156 [183,] 1.0000000 6.710470e-155 3.355235e-155 [184,] 1.0000000 5.549776e-154 2.774888e-154 [185,] 1.0000000 4.857509e-153 2.428754e-153 [186,] 1.0000000 4.259409e-152 2.129705e-152 [187,] 1.0000000 3.697245e-151 1.848623e-151 [188,] 1.0000000 3.193753e-150 1.596877e-150 [189,] 1.0000000 2.838005e-149 1.419002e-149 [190,] 1.0000000 2.157850e-148 1.078925e-148 [191,] 1.0000000 1.901158e-147 9.505790e-148 [192,] 1.0000000 1.667024e-146 8.335121e-147 [193,] 1.0000000 1.406492e-146 7.032462e-147 [194,] 1.0000000 1.238510e-145 6.192551e-146 [195,] 1.0000000 1.085448e-144 5.427242e-145 [196,] 1.0000000 9.468132e-144 4.734066e-144 [197,] 1.0000000 7.679675e-143 3.839838e-143 [198,] 1.0000000 2.923387e-142 1.461694e-142 [199,] 1.0000000 2.528231e-141 1.264115e-141 [200,] 1.0000000 2.134462e-140 1.067231e-140 [201,] 1.0000000 1.828339e-139 9.141693e-140 [202,] 1.0000000 3.151406e-139 1.575703e-139 [203,] 1.0000000 2.688550e-138 1.344275e-138 [204,] 1.0000000 2.208901e-137 1.104450e-137 [205,] 1.0000000 1.807269e-136 9.036346e-137 [206,] 1.0000000 1.060026e-135 5.300132e-136 [207,] 1.0000000 8.074529e-135 4.037265e-135 [208,] 1.0000000 6.686373e-134 3.343186e-134 [209,] 1.0000000 5.390174e-133 2.695087e-133 [210,] 1.0000000 4.311870e-132 2.155935e-132 [211,] 1.0000000 3.297804e-131 1.648902e-131 [212,] 1.0000000 2.625235e-130 1.312617e-130 [213,] 1.0000000 1.227528e-129 6.137642e-130 [214,] 1.0000000 9.775420e-129 4.887710e-129 [215,] 1.0000000 7.697290e-128 3.848645e-128 [216,] 1.0000000 6.028882e-127 3.014441e-127 [217,] 1.0000000 4.688149e-126 2.344075e-126 [218,] 1.0000000 3.639727e-125 1.819864e-125 [219,] 1.0000000 2.800355e-124 1.400178e-124 [220,] 1.0000000 2.151705e-123 1.075853e-123 [221,] 1.0000000 1.650697e-122 8.253484e-123 [222,] 1.0000000 1.278360e-121 6.391799e-122 [223,] 1.0000000 9.723369e-121 4.861685e-121 [224,] 1.0000000 7.140998e-120 3.570499e-120 [225,] 1.0000000 4.900579e-119 2.450290e-119 [226,] 1.0000000 3.680856e-118 1.840428e-118 [227,] 1.0000000 2.809766e-117 1.404883e-117 [228,] 1.0000000 2.091925e-116 1.045962e-116 [229,] 1.0000000 1.549537e-115 7.747684e-116 [230,] 1.0000000 1.143310e-114 5.716552e-115 [231,] 1.0000000 8.307959e-114 4.153980e-114 [232,] 1.0000000 6.074670e-113 3.037335e-113 [233,] 1.0000000 4.421512e-112 2.210756e-112 [234,] 1.0000000 3.203556e-111 1.601778e-111 [235,] 1.0000000 2.341111e-110 1.170555e-110 [236,] 1.0000000 1.680760e-109 8.403800e-110 [237,] 1.0000000 1.180110e-108 5.900549e-109 [238,] 1.0000000 8.518591e-108 4.259296e-108 [239,] 1.0000000 6.190722e-107 3.095361e-107 [240,] 1.0000000 4.477218e-106 2.238609e-106 [241,] 1.0000000 3.204745e-105 1.602372e-105 [242,] 1.0000000 2.295426e-104 1.147713e-104 [243,] 1.0000000 1.636100e-103 8.180498e-104 [244,] 1.0000000 1.113834e-102 5.569168e-103 [245,] 1.0000000 7.863750e-102 3.931875e-102 [246,] 1.0000000 5.443887e-101 2.721943e-101 [247,] 1.0000000 3.127684e-100 1.563842e-100 [248,] 1.0000000 2.535970e-100 1.267985e-100 [249,] 1.0000000 1.735667e-99 8.678335e-100 [250,] 1.0000000 1.069487e-98 5.347437e-99 [251,] 1.0000000 5.986435e-98 2.993218e-98 [252,] 1.0000000 4.144425e-97 2.072213e-97 [253,] 1.0000000 2.785412e-96 1.392706e-96 [254,] 1.0000000 1.909385e-95 9.546924e-96 [255,] 1.0000000 1.302166e-94 6.510832e-95 [256,] 1.0000000 8.718921e-94 4.359461e-94 [257,] 1.0000000 4.501405e-93 2.250703e-93 [258,] 1.0000000 3.029487e-92 1.514744e-92 [259,] 1.0000000 1.979688e-91 9.898441e-92 [260,] 1.0000000 1.318994e-90 6.594972e-91 [261,] 1.0000000 6.980669e-90 3.490335e-90 [262,] 1.0000000 4.609886e-89 2.304943e-89 [263,] 1.0000000 3.022311e-88 1.511155e-88 [264,] 1.0000000 1.975367e-87 9.876836e-88 [265,] 1.0000000 1.284132e-86 6.420662e-87 [266,] 1.0000000 8.302334e-86 4.151167e-86 [267,] 1.0000000 5.314259e-85 2.657129e-85 [268,] 1.0000000 3.325474e-84 1.662737e-84 [269,] 1.0000000 2.061703e-83 1.030852e-83 [270,] 1.0000000 1.292498e-82 6.462492e-83 [271,] 1.0000000 7.968747e-82 3.984373e-82 [272,] 1.0000000 3.903588e-81 1.951794e-81 [273,] 1.0000000 2.318534e-80 1.159267e-80 [274,] 1.0000000 1.444618e-79 7.223088e-80 [275,] 1.0000000 8.735332e-79 4.367666e-79 [276,] 1.0000000 5.099288e-78 2.549644e-78 [277,] 1.0000000 3.053613e-77 1.526806e-77 [278,] 1.0000000 1.866870e-76 9.334352e-77 [279,] 1.0000000 1.106537e-75 5.532685e-76 [280,] 1.0000000 6.429784e-75 3.214892e-75 [281,] 1.0000000 3.788418e-74 1.894209e-74 [282,] 1.0000000 2.263444e-73 1.131722e-73 [283,] 1.0000000 1.342261e-72 6.711306e-73 [284,] 1.0000000 7.930846e-72 3.965423e-72 [285,] 1.0000000 4.477237e-71 2.238618e-71 [286,] 1.0000000 2.586387e-70 1.293193e-70 [287,] 1.0000000 1.470859e-69 7.354296e-70 [288,] 1.0000000 8.539692e-69 4.269846e-69 [289,] 1.0000000 3.765639e-68 1.882819e-68 [290,] 1.0000000 2.111891e-67 1.055945e-67 [291,] 1.0000000 1.198781e-66 5.993906e-67 [292,] 1.0000000 6.650270e-66 3.325135e-66 [293,] 1.0000000 3.668631e-65 1.834315e-65 [294,] 1.0000000 2.060884e-64 1.030442e-64 [295,] 1.0000000 1.151110e-63 5.755551e-64 [296,] 1.0000000 6.245272e-63 3.122636e-63 [297,] 1.0000000 3.423197e-62 1.711598e-62 [298,] 1.0000000 1.477790e-61 7.388951e-62 [299,] 1.0000000 7.918983e-61 3.959492e-61 [300,] 1.0000000 4.276271e-60 2.138136e-60 [301,] 1.0000000 2.326807e-59 1.163404e-59 [302,] 1.0000000 1.125418e-58 5.627090e-59 [303,] 1.0000000 5.968006e-58 2.984003e-58 [304,] 1.0000000 3.139574e-57 1.569787e-57 [305,] 1.0000000 1.332911e-56 6.664553e-57 [306,] 1.0000000 6.903447e-56 3.451724e-56 [307,] 1.0000000 3.591003e-55 1.795501e-55 [308,] 1.0000000 1.829546e-54 9.147728e-55 [309,] 1.0000000 8.782550e-54 4.391275e-54 [310,] 1.0000000 4.378783e-53 2.189392e-53 [311,] 1.0000000 2.248769e-52 1.124385e-52 [312,] 1.0000000 1.121044e-51 5.605221e-52 [313,] 1.0000000 5.535517e-51 2.767759e-51 [314,] 1.0000000 2.716167e-50 1.358084e-50 [315,] 1.0000000 1.324332e-49 6.621658e-50 [316,] 1.0000000 6.415903e-49 3.207952e-49 [317,] 1.0000000 3.093586e-48 1.546793e-48 [318,] 1.0000000 1.479997e-47 7.399987e-48 [319,] 1.0000000 7.119930e-47 3.559965e-47 [320,] 1.0000000 3.349358e-46 1.674679e-46 [321,] 1.0000000 1.552927e-45 7.764635e-46 [322,] 1.0000000 7.241226e-45 3.620613e-45 [323,] 1.0000000 3.353885e-44 1.676943e-44 [324,] 1.0000000 1.542886e-43 7.714428e-44 [325,] 1.0000000 7.049247e-43 3.524623e-43 [326,] 1.0000000 3.198514e-42 1.599257e-42 [327,] 1.0000000 1.446390e-41 7.231949e-42 [328,] 1.0000000 6.483087e-41 3.241544e-41 [329,] 1.0000000 2.854850e-40 1.427425e-40 [330,] 1.0000000 1.283172e-39 6.415861e-40 [331,] 1.0000000 5.624803e-39 2.812401e-39 [332,] 1.0000000 2.468498e-38 1.234249e-38 [333,] 1.0000000 1.048704e-37 5.243518e-38 [334,] 1.0000000 4.633853e-37 2.316926e-37 [335,] 1.0000000 1.975253e-36 9.876264e-37 [336,] 1.0000000 8.356319e-36 4.178159e-36 [337,] 1.0000000 3.508191e-35 1.754095e-35 [338,] 1.0000000 1.516141e-34 7.580706e-35 [339,] 1.0000000 6.271530e-34 3.135765e-34 [340,] 1.0000000 2.645058e-33 1.322529e-33 [341,] 1.0000000 1.077222e-32 5.386108e-33 [342,] 1.0000000 4.351832e-32 2.175916e-32 [343,] 1.0000000 1.580431e-31 7.902153e-32 [344,] 1.0000000 5.981313e-31 2.990656e-31 [345,] 1.0000000 2.401539e-30 1.200770e-30 [346,] 1.0000000 8.459500e-30 4.229750e-30 [347,] 1.0000000 3.340945e-29 1.670472e-29 [348,] 1.0000000 1.285048e-28 6.425239e-29 [349,] 1.0000000 4.899308e-28 2.449654e-28 [350,] 1.0000000 1.763136e-27 8.815681e-28 [351,] 1.0000000 6.628002e-27 3.314001e-27 [352,] 1.0000000 2.468904e-26 1.234452e-26 [353,] 1.0000000 9.111468e-26 4.555734e-26 [354,] 1.0000000 3.349750e-25 1.674875e-25 [355,] 1.0000000 1.213088e-24 6.065441e-25 [356,] 1.0000000 4.350442e-24 2.175221e-24 [357,] 1.0000000 1.524975e-23 7.624877e-24 [358,] 1.0000000 5.313013e-23 2.656507e-23 [359,] 1.0000000 1.853148e-22 9.265742e-23 [360,] 1.0000000 6.396865e-22 3.198433e-22 [361,] 1.0000000 2.184853e-21 1.092427e-21 [362,] 1.0000000 7.382103e-21 3.691052e-21 [363,] 1.0000000 2.527319e-20 1.263659e-20 [364,] 1.0000000 8.359335e-20 4.179668e-20 [365,] 1.0000000 2.723883e-19 1.361941e-19 [366,] 1.0000000 9.501666e-19 4.750833e-19 [367,] 1.0000000 3.249154e-18 1.624577e-18 [368,] 1.0000000 1.030860e-17 5.154301e-18 [369,] 1.0000000 3.230398e-17 1.615199e-17 [370,] 1.0000000 1.059753e-16 5.298765e-17 [371,] 1.0000000 3.065219e-16 1.532609e-16 [372,] 1.0000000 7.980859e-16 3.990429e-16 [373,] 1.0000000 2.377767e-15 1.188883e-15 [374,] 1.0000000 6.983956e-15 3.491978e-15 [375,] 1.0000000 2.202736e-14 1.101368e-14 [376,] 1.0000000 6.534443e-14 3.267222e-14 [377,] 1.0000000 1.856724e-13 9.283621e-14 [378,] 1.0000000 5.428907e-13 2.714453e-13 [379,] 1.0000000 1.429329e-12 7.146643e-13 [380,] 1.0000000 3.629379e-12 1.814690e-12 [381,] 1.0000000 9.523941e-12 4.761970e-12 [382,] 1.0000000 2.547438e-11 1.273719e-11 [383,] 1.0000000 6.756486e-11 3.378243e-11 [384,] 1.0000000 1.730454e-10 8.652272e-11 [385,] 1.0000000 4.751507e-10 2.375753e-10 [386,] 1.0000000 1.358775e-09 6.793876e-10 [387,] 1.0000000 2.576872e-09 1.288436e-09 [388,] 1.0000000 6.635957e-09 3.317979e-09 [389,] 1.0000000 1.761673e-08 8.808366e-09 [390,] 1.0000000 4.432430e-08 2.216215e-08 [391,] 1.0000000 8.306709e-08 4.153354e-08 [392,] 0.9999999 2.131984e-07 1.065992e-07 [393,] 0.9999997 5.623451e-07 2.811726e-07 [394,] 0.9999993 1.408515e-06 7.042574e-07 [395,] 0.9999986 2.844751e-06 1.422376e-06 [396,] 0.9999969 6.263561e-06 3.131780e-06 [397,] 0.9999928 1.444125e-05 7.220627e-06 [398,] 0.9999828 3.446735e-05 1.723367e-05 [399,] 0.9999633 7.343965e-05 3.671983e-05 [400,] 0.9999209 1.581122e-04 7.905612e-05 [401,] 0.9998541 2.918114e-04 1.459057e-04 [402,] 0.9997261 5.477223e-04 2.738612e-04 [403,] 0.9993926 1.214886e-03 6.074429e-04 [404,] 0.9986910 2.618086e-03 1.309043e-03 [405,] 0.9975113 4.977359e-03 2.488680e-03 [406,] 0.9948779 1.024411e-02 5.122056e-03 [407,] 0.9943533 1.129332e-02 5.646660e-03 [408,] 0.9911964 1.760730e-02 8.803649e-03 [409,] 0.9835631 3.287384e-02 1.643692e-02 [410,] 0.9948437 1.031264e-02 5.156320e-03 [411,] 0.9884542 2.309153e-02 1.154577e-02 [412,] 0.9747599 5.048022e-02 2.524011e-02 [413,] 0.9463041 1.073918e-01 5.369590e-02 [414,] 0.8926270 2.147460e-01 1.073730e-01 [415,] 0.7962125 4.075750e-01 2.037875e-01 [416,] 0.6628701 6.742599e-01 3.371299e-01 > postscript(file="/var/www/html/freestat/rcomp/tmp/1qki21291414557.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/freestat/rcomp/tmp/2qki21291414557.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/freestat/rcomp/tmp/3jbi51291414557.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/freestat/rcomp/tmp/4jbi51291414557.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/freestat/rcomp/tmp/5jbi51291414557.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 1259827.443 3195283.878 1187360.002 -2163938.266 -465658.107 -86633.872 7 8 9 10 11 12 -25648.977 64892.532 -278718.282 704705.338 -273562.970 408709.063 13 14 15 16 17 18 762830.868 -1184133.917 5446.994 498000.209 -234688.381 -12297.139 19 20 21 22 23 24 -306359.379 -322029.866 153902.340 -279201.028 -122718.005 -378580.851 25 26 27 28 29 30 -280953.361 -133413.771 -402357.714 -282147.420 333879.073 -89503.531 31 32 33 34 35 36 225880.987 -153760.016 190360.572 80331.658 -144589.537 -459875.779 37 38 39 40 41 42 48357.958 148737.922 -63103.461 -159079.702 -5335.757 -374596.864 43 44 45 46 47 48 127681.941 -158893.721 -8951.678 -124535.591 -298183.103 -131141.638 49 50 51 52 53 54 -28126.641 -157816.795 66254.339 -9824.282 -248826.358 71559.952 55 56 57 58 59 60 -348387.682 105736.955 297382.198 -6488.100 4473.442 -223588.738 61 62 63 64 65 66 -41076.166 -150276.044 -53093.255 -190418.162 43429.078 22800.421 67 68 69 70 71 72 -40997.771 -41888.465 -32278.740 -68881.815 -123725.963 20183.871 73 74 75 76 77 78 -74418.783 10492.987 -29593.037 -152906.495 -14508.840 -62936.025 79 80 81 82 83 84 37562.856 -43408.390 48877.645 -62195.598 -67415.436 58126.423 85 86 87 88 89 90 19942.976 -112771.362 60434.071 -6906.064 -83722.349 24236.170 91 92 93 94 95 96 20323.606 -58373.820 152629.298 -249808.092 97881.580 -18948.548 97 98 99 100 101 102 56121.615 27844.604 -65081.319 -65872.062 -69214.165 -109444.008 103 104 105 106 107 108 -127969.710 -65383.903 51398.186 -31544.333 -266560.840 63483.018 109 110 111 112 113 114 -169212.146 -139539.338 111838.134 -223309.564 -148255.495 -124797.099 115 116 117 118 119 120 23515.058 55479.808 17147.822 63553.861 -200546.755 50336.225 121 122 123 124 125 126 -265684.483 34349.869 22362.228 29543.810 -37960.116 -223683.527 127 128 129 130 131 132 47930.758 46990.954 54791.004 47930.758 49320.735 90994.114 133 134 135 136 137 138 47930.758 45014.662 27433.327 42217.914 52885.045 48103.003 139 140 141 142 143 144 22758.922 47930.758 55492.576 47930.758 15002.087 48197.857 145 146 147 148 149 150 49670.808 48756.725 47930.758 49490.851 50535.241 47847.411 151 152 153 154 155 156 16757.609 16128.433 47930.758 25176.204 48001.084 47930.758 157 158 159 160 161 162 47930.758 44196.429 59511.531 47930.758 47930.758 47930.758 163 164 165 166 167 168 47930.758 35201.158 37073.804 47253.085 47930.758 129383.608 169 170 171 172 173 174 95359.943 47930.758 47930.758 43892.484 47930.758 37413.394 175 176 177 178 179 180 47930.758 47968.115 95780.130 -8516.655 47930.758 47930.758 181 182 183 184 185 186 55714.984 -25660.802 6174.204 -37101.259 -10313.459 -2370.308 187 188 189 190 191 192 -35584.633 40932.200 9957.267 47930.758 47076.205 47930.758 193 194 195 196 197 198 44403.489 47930.758 47890.901 -7555.060 -126404.718 -7555.060 199 200 201 202 203 204 -7555.060 -96498.844 -7555.060 -7555.060 -7555.060 2068.080 205 206 207 208 209 210 -56691.365 -8066.752 -14179.696 -7555.060 23336.527 -11265.356 211 212 213 214 215 216 47930.758 47930.758 73611.546 13127.108 20306.970 47930.758 217 218 219 220 221 222 50834.989 49994.478 47930.758 19719.600 22654.245 47930.758 223 224 225 226 227 228 48672.001 39310.301 48181.113 48453.678 50983.381 47930.758 229 230 231 232 233 234 32019.062 47930.758 19310.917 35997.643 47930.758 29871.149 235 236 237 238 239 240 47930.758 48340.985 47930.758 49921.272 47930.758 47930.758 241 242 243 244 245 246 47930.758 42624.223 47930.758 52466.884 -33992.012 -7555.060 247 248 249 250 251 252 -7555.060 40017.496 -7555.060 -7555.060 41434.050 -7555.060 253 254 255 256 257 258 28918.905 24439.537 -87387.956 47930.758 45851.206 59484.659 259 260 261 262 263 264 -7555.060 47930.758 -7555.060 -7555.060 51292.884 65051.344 265 266 267 268 269 270 -7555.060 47930.758 -7555.060 79985.494 -7555.060 44777.497 271 272 273 274 275 276 -7555.060 -7555.060 -7555.060 43760.928 47930.758 53038.424 277 278 279 280 281 282 49743.121 47930.758 -35871.190 35089.843 -7555.060 47930.758 283 284 285 286 287 288 75626.658 47930.758 39802.780 47930.758 -19198.651 48830.870 289 290 291 292 293 294 -7555.060 55446.705 -7555.060 58494.152 54252.289 47930.758 295 296 297 298 299 300 41418.222 17805.179 47930.758 -7555.060 47930.758 47930.758 301 302 303 304 305 306 -9187.034 29601.307 47930.758 -7555.060 85548.186 51520.595 307 308 309 310 311 312 -7555.060 -13574.240 51524.780 -7555.060 -7555.060 51085.001 313 314 315 316 317 318 -7555.060 -5631.489 49666.953 63580.147 48377.897 24781.823 319 320 321 322 323 324 50337.774 47930.758 47930.758 47930.758 47930.758 -7555.060 325 326 327 328 329 330 47930.758 45207.361 -7555.060 -7555.060 47930.758 47930.758 331 332 333 334 335 336 47930.758 47930.758 47930.758 46110.459 48709.764 40596.813 337 338 339 340 341 342 34892.012 47930.758 45723.816 53714.507 20843.505 47930.758 343 344 345 346 347 348 47930.758 47930.758 41837.252 47930.758 38208.833 47930.758 349 350 351 352 353 354 47930.758 11292.552 51727.086 48104.426 47803.348 40959.878 355 356 357 358 359 360 47930.758 47930.758 25991.297 47930.758 47930.758 47930.758 361 362 363 364 365 366 49543.695 47930.758 47930.758 -59079.502 50139.283 47930.758 367 368 369 370 371 372 47930.758 47930.758 47930.758 50155.084 47930.758 60651.036 373 374 375 376 377 378 35414.901 27297.538 47930.758 47930.758 42833.808 41821.079 379 380 381 382 383 384 -7827.369 47930.758 47930.758 33983.667 65563.601 47930.758 385 386 387 388 389 390 50256.107 -100883.026 60199.198 21375.207 46456.029 34065.162 391 392 393 394 395 396 43281.192 17215.934 -20308.732 68636.060 -180351.327 9037.666 397 398 399 400 401 402 -138046.328 4172.197 10638.634 -23000.886 -11465.171 39825.868 403 404 405 406 407 408 17459.373 -95882.735 -12187.230 10402.487 23566.412 15549.014 409 410 411 412 413 414 26690.904 -73845.864 10354.181 -250216.313 -74411.743 12120.541 415 416 417 418 419 420 -189194.453 -115267.526 118575.766 -62387.577 -35757.824 -49407.953 421 422 423 424 425 426 -60797.317 -107098.631 -576390.302 -484134.906 -41774.594 -232554.023 427 428 429 430 431 -313458.634 -1118357.214 -306374.089 -344582.470 -675705.948 > postscript(file="/var/www/html/freestat/rcomp/tmp/6c2hq1291414557.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 1259827.443 NA 1 3195283.878 1259827.443 2 1187360.002 3195283.878 3 -2163938.266 1187360.002 4 -465658.107 -2163938.266 5 -86633.872 -465658.107 6 -25648.977 -86633.872 7 64892.532 -25648.977 8 -278718.282 64892.532 9 704705.338 -278718.282 10 -273562.970 704705.338 11 408709.063 -273562.970 12 762830.868 408709.063 13 -1184133.917 762830.868 14 5446.994 -1184133.917 15 498000.209 5446.994 16 -234688.381 498000.209 17 -12297.139 -234688.381 18 -306359.379 -12297.139 19 -322029.866 -306359.379 20 153902.340 -322029.866 21 -279201.028 153902.340 22 -122718.005 -279201.028 23 -378580.851 -122718.005 24 -280953.361 -378580.851 25 -133413.771 -280953.361 26 -402357.714 -133413.771 27 -282147.420 -402357.714 28 333879.073 -282147.420 29 -89503.531 333879.073 30 225880.987 -89503.531 31 -153760.016 225880.987 32 190360.572 -153760.016 33 80331.658 190360.572 34 -144589.537 80331.658 35 -459875.779 -144589.537 36 48357.958 -459875.779 37 148737.922 48357.958 38 -63103.461 148737.922 39 -159079.702 -63103.461 40 -5335.757 -159079.702 41 -374596.864 -5335.757 42 127681.941 -374596.864 43 -158893.721 127681.941 44 -8951.678 -158893.721 45 -124535.591 -8951.678 46 -298183.103 -124535.591 47 -131141.638 -298183.103 48 -28126.641 -131141.638 49 -157816.795 -28126.641 50 66254.339 -157816.795 51 -9824.282 66254.339 52 -248826.358 -9824.282 53 71559.952 -248826.358 54 -348387.682 71559.952 55 105736.955 -348387.682 56 297382.198 105736.955 57 -6488.100 297382.198 58 4473.442 -6488.100 59 -223588.738 4473.442 60 -41076.166 -223588.738 61 -150276.044 -41076.166 62 -53093.255 -150276.044 63 -190418.162 -53093.255 64 43429.078 -190418.162 65 22800.421 43429.078 66 -40997.771 22800.421 67 -41888.465 -40997.771 68 -32278.740 -41888.465 69 -68881.815 -32278.740 70 -123725.963 -68881.815 71 20183.871 -123725.963 72 -74418.783 20183.871 73 10492.987 -74418.783 74 -29593.037 10492.987 75 -152906.495 -29593.037 76 -14508.840 -152906.495 77 -62936.025 -14508.840 78 37562.856 -62936.025 79 -43408.390 37562.856 80 48877.645 -43408.390 81 -62195.598 48877.645 82 -67415.436 -62195.598 83 58126.423 -67415.436 84 19942.976 58126.423 85 -112771.362 19942.976 86 60434.071 -112771.362 87 -6906.064 60434.071 88 -83722.349 -6906.064 89 24236.170 -83722.349 90 20323.606 24236.170 91 -58373.820 20323.606 92 152629.298 -58373.820 93 -249808.092 152629.298 94 97881.580 -249808.092 95 -18948.548 97881.580 96 56121.615 -18948.548 97 27844.604 56121.615 98 -65081.319 27844.604 99 -65872.062 -65081.319 100 -69214.165 -65872.062 101 -109444.008 -69214.165 102 -127969.710 -109444.008 103 -65383.903 -127969.710 104 51398.186 -65383.903 105 -31544.333 51398.186 106 -266560.840 -31544.333 107 63483.018 -266560.840 108 -169212.146 63483.018 109 -139539.338 -169212.146 110 111838.134 -139539.338 111 -223309.564 111838.134 112 -148255.495 -223309.564 113 -124797.099 -148255.495 114 23515.058 -124797.099 115 55479.808 23515.058 116 17147.822 55479.808 117 63553.861 17147.822 118 -200546.755 63553.861 119 50336.225 -200546.755 120 -265684.483 50336.225 121 34349.869 -265684.483 122 22362.228 34349.869 123 29543.810 22362.228 124 -37960.116 29543.810 125 -223683.527 -37960.116 126 47930.758 -223683.527 127 46990.954 47930.758 128 54791.004 46990.954 129 47930.758 54791.004 130 49320.735 47930.758 131 90994.114 49320.735 132 47930.758 90994.114 133 45014.662 47930.758 134 27433.327 45014.662 135 42217.914 27433.327 136 52885.045 42217.914 137 48103.003 52885.045 138 22758.922 48103.003 139 47930.758 22758.922 140 55492.576 47930.758 141 47930.758 55492.576 142 15002.087 47930.758 143 48197.857 15002.087 144 49670.808 48197.857 145 48756.725 49670.808 146 47930.758 48756.725 147 49490.851 47930.758 148 50535.241 49490.851 149 47847.411 50535.241 150 16757.609 47847.411 151 16128.433 16757.609 152 47930.758 16128.433 153 25176.204 47930.758 154 48001.084 25176.204 155 47930.758 48001.084 156 47930.758 47930.758 157 44196.429 47930.758 158 59511.531 44196.429 159 47930.758 59511.531 160 47930.758 47930.758 161 47930.758 47930.758 162 47930.758 47930.758 163 35201.158 47930.758 164 37073.804 35201.158 165 47253.085 37073.804 166 47930.758 47253.085 167 129383.608 47930.758 168 95359.943 129383.608 169 47930.758 95359.943 170 47930.758 47930.758 171 43892.484 47930.758 172 47930.758 43892.484 173 37413.394 47930.758 174 47930.758 37413.394 175 47968.115 47930.758 176 95780.130 47968.115 177 -8516.655 95780.130 178 47930.758 -8516.655 179 47930.758 47930.758 180 55714.984 47930.758 181 -25660.802 55714.984 182 6174.204 -25660.802 183 -37101.259 6174.204 184 -10313.459 -37101.259 185 -2370.308 -10313.459 186 -35584.633 -2370.308 187 40932.200 -35584.633 188 9957.267 40932.200 189 47930.758 9957.267 190 47076.205 47930.758 191 47930.758 47076.205 192 44403.489 47930.758 193 47930.758 44403.489 194 47890.901 47930.758 195 -7555.060 47890.901 196 -126404.718 -7555.060 197 -7555.060 -126404.718 198 -7555.060 -7555.060 199 -96498.844 -7555.060 200 -7555.060 -96498.844 201 -7555.060 -7555.060 202 -7555.060 -7555.060 203 2068.080 -7555.060 204 -56691.365 2068.080 205 -8066.752 -56691.365 206 -14179.696 -8066.752 207 -7555.060 -14179.696 208 23336.527 -7555.060 209 -11265.356 23336.527 210 47930.758 -11265.356 211 47930.758 47930.758 212 73611.546 47930.758 213 13127.108 73611.546 214 20306.970 13127.108 215 47930.758 20306.970 216 50834.989 47930.758 217 49994.478 50834.989 218 47930.758 49994.478 219 19719.600 47930.758 220 22654.245 19719.600 221 47930.758 22654.245 222 48672.001 47930.758 223 39310.301 48672.001 224 48181.113 39310.301 225 48453.678 48181.113 226 50983.381 48453.678 227 47930.758 50983.381 228 32019.062 47930.758 229 47930.758 32019.062 230 19310.917 47930.758 231 35997.643 19310.917 232 47930.758 35997.643 233 29871.149 47930.758 234 47930.758 29871.149 235 48340.985 47930.758 236 47930.758 48340.985 237 49921.272 47930.758 238 47930.758 49921.272 239 47930.758 47930.758 240 47930.758 47930.758 241 42624.223 47930.758 242 47930.758 42624.223 243 52466.884 47930.758 244 -33992.012 52466.884 245 -7555.060 -33992.012 246 -7555.060 -7555.060 247 40017.496 -7555.060 248 -7555.060 40017.496 249 -7555.060 -7555.060 250 41434.050 -7555.060 251 -7555.060 41434.050 252 28918.905 -7555.060 253 24439.537 28918.905 254 -87387.956 24439.537 255 47930.758 -87387.956 256 45851.206 47930.758 257 59484.659 45851.206 258 -7555.060 59484.659 259 47930.758 -7555.060 260 -7555.060 47930.758 261 -7555.060 -7555.060 262 51292.884 -7555.060 263 65051.344 51292.884 264 -7555.060 65051.344 265 47930.758 -7555.060 266 -7555.060 47930.758 267 79985.494 -7555.060 268 -7555.060 79985.494 269 44777.497 -7555.060 270 -7555.060 44777.497 271 -7555.060 -7555.060 272 -7555.060 -7555.060 273 43760.928 -7555.060 274 47930.758 43760.928 275 53038.424 47930.758 276 49743.121 53038.424 277 47930.758 49743.121 278 -35871.190 47930.758 279 35089.843 -35871.190 280 -7555.060 35089.843 281 47930.758 -7555.060 282 75626.658 47930.758 283 47930.758 75626.658 284 39802.780 47930.758 285 47930.758 39802.780 286 -19198.651 47930.758 287 48830.870 -19198.651 288 -7555.060 48830.870 289 55446.705 -7555.060 290 -7555.060 55446.705 291 58494.152 -7555.060 292 54252.289 58494.152 293 47930.758 54252.289 294 41418.222 47930.758 295 17805.179 41418.222 296 47930.758 17805.179 297 -7555.060 47930.758 298 47930.758 -7555.060 299 47930.758 47930.758 300 -9187.034 47930.758 301 29601.307 -9187.034 302 47930.758 29601.307 303 -7555.060 47930.758 304 85548.186 -7555.060 305 51520.595 85548.186 306 -7555.060 51520.595 307 -13574.240 -7555.060 308 51524.780 -13574.240 309 -7555.060 51524.780 310 -7555.060 -7555.060 311 51085.001 -7555.060 312 -7555.060 51085.001 313 -5631.489 -7555.060 314 49666.953 -5631.489 315 63580.147 49666.953 316 48377.897 63580.147 317 24781.823 48377.897 318 50337.774 24781.823 319 47930.758 50337.774 320 47930.758 47930.758 321 47930.758 47930.758 322 47930.758 47930.758 323 -7555.060 47930.758 324 47930.758 -7555.060 325 45207.361 47930.758 326 -7555.060 45207.361 327 -7555.060 -7555.060 328 47930.758 -7555.060 329 47930.758 47930.758 330 47930.758 47930.758 331 47930.758 47930.758 332 47930.758 47930.758 333 46110.459 47930.758 334 48709.764 46110.459 335 40596.813 48709.764 336 34892.012 40596.813 337 47930.758 34892.012 338 45723.816 47930.758 339 53714.507 45723.816 340 20843.505 53714.507 341 47930.758 20843.505 342 47930.758 47930.758 343 47930.758 47930.758 344 41837.252 47930.758 345 47930.758 41837.252 346 38208.833 47930.758 347 47930.758 38208.833 348 47930.758 47930.758 349 11292.552 47930.758 350 51727.086 11292.552 351 48104.426 51727.086 352 47803.348 48104.426 353 40959.878 47803.348 354 47930.758 40959.878 355 47930.758 47930.758 356 25991.297 47930.758 357 47930.758 25991.297 358 47930.758 47930.758 359 47930.758 47930.758 360 49543.695 47930.758 361 47930.758 49543.695 362 47930.758 47930.758 363 -59079.502 47930.758 364 50139.283 -59079.502 365 47930.758 50139.283 366 47930.758 47930.758 367 47930.758 47930.758 368 47930.758 47930.758 369 50155.084 47930.758 370 47930.758 50155.084 371 60651.036 47930.758 372 35414.901 60651.036 373 27297.538 35414.901 374 47930.758 27297.538 375 47930.758 47930.758 376 42833.808 47930.758 377 41821.079 42833.808 378 -7827.369 41821.079 379 47930.758 -7827.369 380 47930.758 47930.758 381 33983.667 47930.758 382 65563.601 33983.667 383 47930.758 65563.601 384 50256.107 47930.758 385 -100883.026 50256.107 386 60199.198 -100883.026 387 21375.207 60199.198 388 46456.029 21375.207 389 34065.162 46456.029 390 43281.192 34065.162 391 17215.934 43281.192 392 -20308.732 17215.934 393 68636.060 -20308.732 394 -180351.327 68636.060 395 9037.666 -180351.327 396 -138046.328 9037.666 397 4172.197 -138046.328 398 10638.634 4172.197 399 -23000.886 10638.634 400 -11465.171 -23000.886 401 39825.868 -11465.171 402 17459.373 39825.868 403 -95882.735 17459.373 404 -12187.230 -95882.735 405 10402.487 -12187.230 406 23566.412 10402.487 407 15549.014 23566.412 408 26690.904 15549.014 409 -73845.864 26690.904 410 10354.181 -73845.864 411 -250216.313 10354.181 412 -74411.743 -250216.313 413 12120.541 -74411.743 414 -189194.453 12120.541 415 -115267.526 -189194.453 416 118575.766 -115267.526 417 -62387.577 118575.766 418 -35757.824 -62387.577 419 -49407.953 -35757.824 420 -60797.317 -49407.953 421 -107098.631 -60797.317 422 -576390.302 -107098.631 423 -484134.906 -576390.302 424 -41774.594 -484134.906 425 -232554.023 -41774.594 426 -313458.634 -232554.023 427 -1118357.214 -313458.634 428 -306374.089 -1118357.214 429 -344582.470 -306374.089 430 -675705.948 -344582.470 431 NA -675705.948 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 3195283.878 1259827.443 [2,] 1187360.002 3195283.878 [3,] -2163938.266 1187360.002 [4,] -465658.107 -2163938.266 [5,] -86633.872 -465658.107 [6,] -25648.977 -86633.872 [7,] 64892.532 -25648.977 [8,] -278718.282 64892.532 [9,] 704705.338 -278718.282 [10,] -273562.970 704705.338 [11,] 408709.063 -273562.970 [12,] 762830.868 408709.063 [13,] -1184133.917 762830.868 [14,] 5446.994 -1184133.917 [15,] 498000.209 5446.994 [16,] -234688.381 498000.209 [17,] -12297.139 -234688.381 [18,] -306359.379 -12297.139 [19,] -322029.866 -306359.379 [20,] 153902.340 -322029.866 [21,] -279201.028 153902.340 [22,] -122718.005 -279201.028 [23,] -378580.851 -122718.005 [24,] -280953.361 -378580.851 [25,] -133413.771 -280953.361 [26,] -402357.714 -133413.771 [27,] -282147.420 -402357.714 [28,] 333879.073 -282147.420 [29,] -89503.531 333879.073 [30,] 225880.987 -89503.531 [31,] -153760.016 225880.987 [32,] 190360.572 -153760.016 [33,] 80331.658 190360.572 [34,] -144589.537 80331.658 [35,] -459875.779 -144589.537 [36,] 48357.958 -459875.779 [37,] 148737.922 48357.958 [38,] -63103.461 148737.922 [39,] -159079.702 -63103.461 [40,] -5335.757 -159079.702 [41,] -374596.864 -5335.757 [42,] 127681.941 -374596.864 [43,] -158893.721 127681.941 [44,] -8951.678 -158893.721 [45,] -124535.591 -8951.678 [46,] -298183.103 -124535.591 [47,] -131141.638 -298183.103 [48,] -28126.641 -131141.638 [49,] -157816.795 -28126.641 [50,] 66254.339 -157816.795 [51,] -9824.282 66254.339 [52,] -248826.358 -9824.282 [53,] 71559.952 -248826.358 [54,] -348387.682 71559.952 [55,] 105736.955 -348387.682 [56,] 297382.198 105736.955 [57,] -6488.100 297382.198 [58,] 4473.442 -6488.100 [59,] -223588.738 4473.442 [60,] -41076.166 -223588.738 [61,] -150276.044 -41076.166 [62,] -53093.255 -150276.044 [63,] -190418.162 -53093.255 [64,] 43429.078 -190418.162 [65,] 22800.421 43429.078 [66,] -40997.771 22800.421 [67,] -41888.465 -40997.771 [68,] -32278.740 -41888.465 [69,] -68881.815 -32278.740 [70,] -123725.963 -68881.815 [71,] 20183.871 -123725.963 [72,] -74418.783 20183.871 [73,] 10492.987 -74418.783 [74,] -29593.037 10492.987 [75,] -152906.495 -29593.037 [76,] -14508.840 -152906.495 [77,] -62936.025 -14508.840 [78,] 37562.856 -62936.025 [79,] -43408.390 37562.856 [80,] 48877.645 -43408.390 [81,] -62195.598 48877.645 [82,] -67415.436 -62195.598 [83,] 58126.423 -67415.436 [84,] 19942.976 58126.423 [85,] -112771.362 19942.976 [86,] 60434.071 -112771.362 [87,] -6906.064 60434.071 [88,] -83722.349 -6906.064 [89,] 24236.170 -83722.349 [90,] 20323.606 24236.170 [91,] -58373.820 20323.606 [92,] 152629.298 -58373.820 [93,] -249808.092 152629.298 [94,] 97881.580 -249808.092 [95,] -18948.548 97881.580 [96,] 56121.615 -18948.548 [97,] 27844.604 56121.615 [98,] -65081.319 27844.604 [99,] -65872.062 -65081.319 [100,] -69214.165 -65872.062 [101,] -109444.008 -69214.165 [102,] -127969.710 -109444.008 [103,] -65383.903 -127969.710 [104,] 51398.186 -65383.903 [105,] -31544.333 51398.186 [106,] -266560.840 -31544.333 [107,] 63483.018 -266560.840 [108,] -169212.146 63483.018 [109,] -139539.338 -169212.146 [110,] 111838.134 -139539.338 [111,] -223309.564 111838.134 [112,] -148255.495 -223309.564 [113,] -124797.099 -148255.495 [114,] 23515.058 -124797.099 [115,] 55479.808 23515.058 [116,] 17147.822 55479.808 [117,] 63553.861 17147.822 [118,] -200546.755 63553.861 [119,] 50336.225 -200546.755 [120,] -265684.483 50336.225 [121,] 34349.869 -265684.483 [122,] 22362.228 34349.869 [123,] 29543.810 22362.228 [124,] -37960.116 29543.810 [125,] -223683.527 -37960.116 [126,] 47930.758 -223683.527 [127,] 46990.954 47930.758 [128,] 54791.004 46990.954 [129,] 47930.758 54791.004 [130,] 49320.735 47930.758 [131,] 90994.114 49320.735 [132,] 47930.758 90994.114 [133,] 45014.662 47930.758 [134,] 27433.327 45014.662 [135,] 42217.914 27433.327 [136,] 52885.045 42217.914 [137,] 48103.003 52885.045 [138,] 22758.922 48103.003 [139,] 47930.758 22758.922 [140,] 55492.576 47930.758 [141,] 47930.758 55492.576 [142,] 15002.087 47930.758 [143,] 48197.857 15002.087 [144,] 49670.808 48197.857 [145,] 48756.725 49670.808 [146,] 47930.758 48756.725 [147,] 49490.851 47930.758 [148,] 50535.241 49490.851 [149,] 47847.411 50535.241 [150,] 16757.609 47847.411 [151,] 16128.433 16757.609 [152,] 47930.758 16128.433 [153,] 25176.204 47930.758 [154,] 48001.084 25176.204 [155,] 47930.758 48001.084 [156,] 47930.758 47930.758 [157,] 44196.429 47930.758 [158,] 59511.531 44196.429 [159,] 47930.758 59511.531 [160,] 47930.758 47930.758 [161,] 47930.758 47930.758 [162,] 47930.758 47930.758 [163,] 35201.158 47930.758 [164,] 37073.804 35201.158 [165,] 47253.085 37073.804 [166,] 47930.758 47253.085 [167,] 129383.608 47930.758 [168,] 95359.943 129383.608 [169,] 47930.758 95359.943 [170,] 47930.758 47930.758 [171,] 43892.484 47930.758 [172,] 47930.758 43892.484 [173,] 37413.394 47930.758 [174,] 47930.758 37413.394 [175,] 47968.115 47930.758 [176,] 95780.130 47968.115 [177,] -8516.655 95780.130 [178,] 47930.758 -8516.655 [179,] 47930.758 47930.758 [180,] 55714.984 47930.758 [181,] -25660.802 55714.984 [182,] 6174.204 -25660.802 [183,] -37101.259 6174.204 [184,] -10313.459 -37101.259 [185,] -2370.308 -10313.459 [186,] -35584.633 -2370.308 [187,] 40932.200 -35584.633 [188,] 9957.267 40932.200 [189,] 47930.758 9957.267 [190,] 47076.205 47930.758 [191,] 47930.758 47076.205 [192,] 44403.489 47930.758 [193,] 47930.758 44403.489 [194,] 47890.901 47930.758 [195,] -7555.060 47890.901 [196,] -126404.718 -7555.060 [197,] -7555.060 -126404.718 [198,] -7555.060 -7555.060 [199,] -96498.844 -7555.060 [200,] -7555.060 -96498.844 [201,] -7555.060 -7555.060 [202,] -7555.060 -7555.060 [203,] 2068.080 -7555.060 [204,] -56691.365 2068.080 [205,] -8066.752 -56691.365 [206,] -14179.696 -8066.752 [207,] -7555.060 -14179.696 [208,] 23336.527 -7555.060 [209,] -11265.356 23336.527 [210,] 47930.758 -11265.356 [211,] 47930.758 47930.758 [212,] 73611.546 47930.758 [213,] 13127.108 73611.546 [214,] 20306.970 13127.108 [215,] 47930.758 20306.970 [216,] 50834.989 47930.758 [217,] 49994.478 50834.989 [218,] 47930.758 49994.478 [219,] 19719.600 47930.758 [220,] 22654.245 19719.600 [221,] 47930.758 22654.245 [222,] 48672.001 47930.758 [223,] 39310.301 48672.001 [224,] 48181.113 39310.301 [225,] 48453.678 48181.113 [226,] 50983.381 48453.678 [227,] 47930.758 50983.381 [228,] 32019.062 47930.758 [229,] 47930.758 32019.062 [230,] 19310.917 47930.758 [231,] 35997.643 19310.917 [232,] 47930.758 35997.643 [233,] 29871.149 47930.758 [234,] 47930.758 29871.149 [235,] 48340.985 47930.758 [236,] 47930.758 48340.985 [237,] 49921.272 47930.758 [238,] 47930.758 49921.272 [239,] 47930.758 47930.758 [240,] 47930.758 47930.758 [241,] 42624.223 47930.758 [242,] 47930.758 42624.223 [243,] 52466.884 47930.758 [244,] -33992.012 52466.884 [245,] -7555.060 -33992.012 [246,] -7555.060 -7555.060 [247,] 40017.496 -7555.060 [248,] -7555.060 40017.496 [249,] -7555.060 -7555.060 [250,] 41434.050 -7555.060 [251,] -7555.060 41434.050 [252,] 28918.905 -7555.060 [253,] 24439.537 28918.905 [254,] -87387.956 24439.537 [255,] 47930.758 -87387.956 [256,] 45851.206 47930.758 [257,] 59484.659 45851.206 [258,] -7555.060 59484.659 [259,] 47930.758 -7555.060 [260,] -7555.060 47930.758 [261,] -7555.060 -7555.060 [262,] 51292.884 -7555.060 [263,] 65051.344 51292.884 [264,] -7555.060 65051.344 [265,] 47930.758 -7555.060 [266,] -7555.060 47930.758 [267,] 79985.494 -7555.060 [268,] -7555.060 79985.494 [269,] 44777.497 -7555.060 [270,] -7555.060 44777.497 [271,] -7555.060 -7555.060 [272,] -7555.060 -7555.060 [273,] 43760.928 -7555.060 [274,] 47930.758 43760.928 [275,] 53038.424 47930.758 [276,] 49743.121 53038.424 [277,] 47930.758 49743.121 [278,] -35871.190 47930.758 [279,] 35089.843 -35871.190 [280,] -7555.060 35089.843 [281,] 47930.758 -7555.060 [282,] 75626.658 47930.758 [283,] 47930.758 75626.658 [284,] 39802.780 47930.758 [285,] 47930.758 39802.780 [286,] -19198.651 47930.758 [287,] 48830.870 -19198.651 [288,] -7555.060 48830.870 [289,] 55446.705 -7555.060 [290,] -7555.060 55446.705 [291,] 58494.152 -7555.060 [292,] 54252.289 58494.152 [293,] 47930.758 54252.289 [294,] 41418.222 47930.758 [295,] 17805.179 41418.222 [296,] 47930.758 17805.179 [297,] -7555.060 47930.758 [298,] 47930.758 -7555.060 [299,] 47930.758 47930.758 [300,] -9187.034 47930.758 [301,] 29601.307 -9187.034 [302,] 47930.758 29601.307 [303,] -7555.060 47930.758 [304,] 85548.186 -7555.060 [305,] 51520.595 85548.186 [306,] -7555.060 51520.595 [307,] -13574.240 -7555.060 [308,] 51524.780 -13574.240 [309,] -7555.060 51524.780 [310,] -7555.060 -7555.060 [311,] 51085.001 -7555.060 [312,] -7555.060 51085.001 [313,] -5631.489 -7555.060 [314,] 49666.953 -5631.489 [315,] 63580.147 49666.953 [316,] 48377.897 63580.147 [317,] 24781.823 48377.897 [318,] 50337.774 24781.823 [319,] 47930.758 50337.774 [320,] 47930.758 47930.758 [321,] 47930.758 47930.758 [322,] 47930.758 47930.758 [323,] -7555.060 47930.758 [324,] 47930.758 -7555.060 [325,] 45207.361 47930.758 [326,] -7555.060 45207.361 [327,] -7555.060 -7555.060 [328,] 47930.758 -7555.060 [329,] 47930.758 47930.758 [330,] 47930.758 47930.758 [331,] 47930.758 47930.758 [332,] 47930.758 47930.758 [333,] 46110.459 47930.758 [334,] 48709.764 46110.459 [335,] 40596.813 48709.764 [336,] 34892.012 40596.813 [337,] 47930.758 34892.012 [338,] 45723.816 47930.758 [339,] 53714.507 45723.816 [340,] 20843.505 53714.507 [341,] 47930.758 20843.505 [342,] 47930.758 47930.758 [343,] 47930.758 47930.758 [344,] 41837.252 47930.758 [345,] 47930.758 41837.252 [346,] 38208.833 47930.758 [347,] 47930.758 38208.833 [348,] 47930.758 47930.758 [349,] 11292.552 47930.758 [350,] 51727.086 11292.552 [351,] 48104.426 51727.086 [352,] 47803.348 48104.426 [353,] 40959.878 47803.348 [354,] 47930.758 40959.878 [355,] 47930.758 47930.758 [356,] 25991.297 47930.758 [357,] 47930.758 25991.297 [358,] 47930.758 47930.758 [359,] 47930.758 47930.758 [360,] 49543.695 47930.758 [361,] 47930.758 49543.695 [362,] 47930.758 47930.758 [363,] -59079.502 47930.758 [364,] 50139.283 -59079.502 [365,] 47930.758 50139.283 [366,] 47930.758 47930.758 [367,] 47930.758 47930.758 [368,] 47930.758 47930.758 [369,] 50155.084 47930.758 [370,] 47930.758 50155.084 [371,] 60651.036 47930.758 [372,] 35414.901 60651.036 [373,] 27297.538 35414.901 [374,] 47930.758 27297.538 [375,] 47930.758 47930.758 [376,] 42833.808 47930.758 [377,] 41821.079 42833.808 [378,] -7827.369 41821.079 [379,] 47930.758 -7827.369 [380,] 47930.758 47930.758 [381,] 33983.667 47930.758 [382,] 65563.601 33983.667 [383,] 47930.758 65563.601 [384,] 50256.107 47930.758 [385,] -100883.026 50256.107 [386,] 60199.198 -100883.026 [387,] 21375.207 60199.198 [388,] 46456.029 21375.207 [389,] 34065.162 46456.029 [390,] 43281.192 34065.162 [391,] 17215.934 43281.192 [392,] -20308.732 17215.934 [393,] 68636.060 -20308.732 [394,] -180351.327 68636.060 [395,] 9037.666 -180351.327 [396,] -138046.328 9037.666 [397,] 4172.197 -138046.328 [398,] 10638.634 4172.197 [399,] -23000.886 10638.634 [400,] -11465.171 -23000.886 [401,] 39825.868 -11465.171 [402,] 17459.373 39825.868 [403,] -95882.735 17459.373 [404,] -12187.230 -95882.735 [405,] 10402.487 -12187.230 [406,] 23566.412 10402.487 [407,] 15549.014 23566.412 [408,] 26690.904 15549.014 [409,] -73845.864 26690.904 [410,] 10354.181 -73845.864 [411,] -250216.313 10354.181 [412,] -74411.743 -250216.313 [413,] 12120.541 -74411.743 [414,] -189194.453 12120.541 [415,] -115267.526 -189194.453 [416,] 118575.766 -115267.526 [417,] -62387.577 118575.766 [418,] -35757.824 -62387.577 [419,] -49407.953 -35757.824 [420,] -60797.317 -49407.953 [421,] -107098.631 -60797.317 [422,] -576390.302 -107098.631 [423,] -484134.906 -576390.302 [424,] -41774.594 -484134.906 [425,] -232554.023 -41774.594 [426,] -313458.634 -232554.023 [427,] -1118357.214 -313458.634 [428,] -306374.089 -1118357.214 [429,] -344582.470 -306374.089 [430,] -675705.948 -344582.470 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 3195283.878 1259827.443 2 1187360.002 3195283.878 3 -2163938.266 1187360.002 4 -465658.107 -2163938.266 5 -86633.872 -465658.107 6 -25648.977 -86633.872 7 64892.532 -25648.977 8 -278718.282 64892.532 9 704705.338 -278718.282 10 -273562.970 704705.338 11 408709.063 -273562.970 12 762830.868 408709.063 13 -1184133.917 762830.868 14 5446.994 -1184133.917 15 498000.209 5446.994 16 -234688.381 498000.209 17 -12297.139 -234688.381 18 -306359.379 -12297.139 19 -322029.866 -306359.379 20 153902.340 -322029.866 21 -279201.028 153902.340 22 -122718.005 -279201.028 23 -378580.851 -122718.005 24 -280953.361 -378580.851 25 -133413.771 -280953.361 26 -402357.714 -133413.771 27 -282147.420 -402357.714 28 333879.073 -282147.420 29 -89503.531 333879.073 30 225880.987 -89503.531 31 -153760.016 225880.987 32 190360.572 -153760.016 33 80331.658 190360.572 34 -144589.537 80331.658 35 -459875.779 -144589.537 36 48357.958 -459875.779 37 148737.922 48357.958 38 -63103.461 148737.922 39 -159079.702 -63103.461 40 -5335.757 -159079.702 41 -374596.864 -5335.757 42 127681.941 -374596.864 43 -158893.721 127681.941 44 -8951.678 -158893.721 45 -124535.591 -8951.678 46 -298183.103 -124535.591 47 -131141.638 -298183.103 48 -28126.641 -131141.638 49 -157816.795 -28126.641 50 66254.339 -157816.795 51 -9824.282 66254.339 52 -248826.358 -9824.282 53 71559.952 -248826.358 54 -348387.682 71559.952 55 105736.955 -348387.682 56 297382.198 105736.955 57 -6488.100 297382.198 58 4473.442 -6488.100 59 -223588.738 4473.442 60 -41076.166 -223588.738 61 -150276.044 -41076.166 62 -53093.255 -150276.044 63 -190418.162 -53093.255 64 43429.078 -190418.162 65 22800.421 43429.078 66 -40997.771 22800.421 67 -41888.465 -40997.771 68 -32278.740 -41888.465 69 -68881.815 -32278.740 70 -123725.963 -68881.815 71 20183.871 -123725.963 72 -74418.783 20183.871 73 10492.987 -74418.783 74 -29593.037 10492.987 75 -152906.495 -29593.037 76 -14508.840 -152906.495 77 -62936.025 -14508.840 78 37562.856 -62936.025 79 -43408.390 37562.856 80 48877.645 -43408.390 81 -62195.598 48877.645 82 -67415.436 -62195.598 83 58126.423 -67415.436 84 19942.976 58126.423 85 -112771.362 19942.976 86 60434.071 -112771.362 87 -6906.064 60434.071 88 -83722.349 -6906.064 89 24236.170 -83722.349 90 20323.606 24236.170 91 -58373.820 20323.606 92 152629.298 -58373.820 93 -249808.092 152629.298 94 97881.580 -249808.092 95 -18948.548 97881.580 96 56121.615 -18948.548 97 27844.604 56121.615 98 -65081.319 27844.604 99 -65872.062 -65081.319 100 -69214.165 -65872.062 101 -109444.008 -69214.165 102 -127969.710 -109444.008 103 -65383.903 -127969.710 104 51398.186 -65383.903 105 -31544.333 51398.186 106 -266560.840 -31544.333 107 63483.018 -266560.840 108 -169212.146 63483.018 109 -139539.338 -169212.146 110 111838.134 -139539.338 111 -223309.564 111838.134 112 -148255.495 -223309.564 113 -124797.099 -148255.495 114 23515.058 -124797.099 115 55479.808 23515.058 116 17147.822 55479.808 117 63553.861 17147.822 118 -200546.755 63553.861 119 50336.225 -200546.755 120 -265684.483 50336.225 121 34349.869 -265684.483 122 22362.228 34349.869 123 29543.810 22362.228 124 -37960.116 29543.810 125 -223683.527 -37960.116 126 47930.758 -223683.527 127 46990.954 47930.758 128 54791.004 46990.954 129 47930.758 54791.004 130 49320.735 47930.758 131 90994.114 49320.735 132 47930.758 90994.114 133 45014.662 47930.758 134 27433.327 45014.662 135 42217.914 27433.327 136 52885.045 42217.914 137 48103.003 52885.045 138 22758.922 48103.003 139 47930.758 22758.922 140 55492.576 47930.758 141 47930.758 55492.576 142 15002.087 47930.758 143 48197.857 15002.087 144 49670.808 48197.857 145 48756.725 49670.808 146 47930.758 48756.725 147 49490.851 47930.758 148 50535.241 49490.851 149 47847.411 50535.241 150 16757.609 47847.411 151 16128.433 16757.609 152 47930.758 16128.433 153 25176.204 47930.758 154 48001.084 25176.204 155 47930.758 48001.084 156 47930.758 47930.758 157 44196.429 47930.758 158 59511.531 44196.429 159 47930.758 59511.531 160 47930.758 47930.758 161 47930.758 47930.758 162 47930.758 47930.758 163 35201.158 47930.758 164 37073.804 35201.158 165 47253.085 37073.804 166 47930.758 47253.085 167 129383.608 47930.758 168 95359.943 129383.608 169 47930.758 95359.943 170 47930.758 47930.758 171 43892.484 47930.758 172 47930.758 43892.484 173 37413.394 47930.758 174 47930.758 37413.394 175 47968.115 47930.758 176 95780.130 47968.115 177 -8516.655 95780.130 178 47930.758 -8516.655 179 47930.758 47930.758 180 55714.984 47930.758 181 -25660.802 55714.984 182 6174.204 -25660.802 183 -37101.259 6174.204 184 -10313.459 -37101.259 185 -2370.308 -10313.459 186 -35584.633 -2370.308 187 40932.200 -35584.633 188 9957.267 40932.200 189 47930.758 9957.267 190 47076.205 47930.758 191 47930.758 47076.205 192 44403.489 47930.758 193 47930.758 44403.489 194 47890.901 47930.758 195 -7555.060 47890.901 196 -126404.718 -7555.060 197 -7555.060 -126404.718 198 -7555.060 -7555.060 199 -96498.844 -7555.060 200 -7555.060 -96498.844 201 -7555.060 -7555.060 202 -7555.060 -7555.060 203 2068.080 -7555.060 204 -56691.365 2068.080 205 -8066.752 -56691.365 206 -14179.696 -8066.752 207 -7555.060 -14179.696 208 23336.527 -7555.060 209 -11265.356 23336.527 210 47930.758 -11265.356 211 47930.758 47930.758 212 73611.546 47930.758 213 13127.108 73611.546 214 20306.970 13127.108 215 47930.758 20306.970 216 50834.989 47930.758 217 49994.478 50834.989 218 47930.758 49994.478 219 19719.600 47930.758 220 22654.245 19719.600 221 47930.758 22654.245 222 48672.001 47930.758 223 39310.301 48672.001 224 48181.113 39310.301 225 48453.678 48181.113 226 50983.381 48453.678 227 47930.758 50983.381 228 32019.062 47930.758 229 47930.758 32019.062 230 19310.917 47930.758 231 35997.643 19310.917 232 47930.758 35997.643 233 29871.149 47930.758 234 47930.758 29871.149 235 48340.985 47930.758 236 47930.758 48340.985 237 49921.272 47930.758 238 47930.758 49921.272 239 47930.758 47930.758 240 47930.758 47930.758 241 42624.223 47930.758 242 47930.758 42624.223 243 52466.884 47930.758 244 -33992.012 52466.884 245 -7555.060 -33992.012 246 -7555.060 -7555.060 247 40017.496 -7555.060 248 -7555.060 40017.496 249 -7555.060 -7555.060 250 41434.050 -7555.060 251 -7555.060 41434.050 252 28918.905 -7555.060 253 24439.537 28918.905 254 -87387.956 24439.537 255 47930.758 -87387.956 256 45851.206 47930.758 257 59484.659 45851.206 258 -7555.060 59484.659 259 47930.758 -7555.060 260 -7555.060 47930.758 261 -7555.060 -7555.060 262 51292.884 -7555.060 263 65051.344 51292.884 264 -7555.060 65051.344 265 47930.758 -7555.060 266 -7555.060 47930.758 267 79985.494 -7555.060 268 -7555.060 79985.494 269 44777.497 -7555.060 270 -7555.060 44777.497 271 -7555.060 -7555.060 272 -7555.060 -7555.060 273 43760.928 -7555.060 274 47930.758 43760.928 275 53038.424 47930.758 276 49743.121 53038.424 277 47930.758 49743.121 278 -35871.190 47930.758 279 35089.843 -35871.190 280 -7555.060 35089.843 281 47930.758 -7555.060 282 75626.658 47930.758 283 47930.758 75626.658 284 39802.780 47930.758 285 47930.758 39802.780 286 -19198.651 47930.758 287 48830.870 -19198.651 288 -7555.060 48830.870 289 55446.705 -7555.060 290 -7555.060 55446.705 291 58494.152 -7555.060 292 54252.289 58494.152 293 47930.758 54252.289 294 41418.222 47930.758 295 17805.179 41418.222 296 47930.758 17805.179 297 -7555.060 47930.758 298 47930.758 -7555.060 299 47930.758 47930.758 300 -9187.034 47930.758 301 29601.307 -9187.034 302 47930.758 29601.307 303 -7555.060 47930.758 304 85548.186 -7555.060 305 51520.595 85548.186 306 -7555.060 51520.595 307 -13574.240 -7555.060 308 51524.780 -13574.240 309 -7555.060 51524.780 310 -7555.060 -7555.060 311 51085.001 -7555.060 312 -7555.060 51085.001 313 -5631.489 -7555.060 314 49666.953 -5631.489 315 63580.147 49666.953 316 48377.897 63580.147 317 24781.823 48377.897 318 50337.774 24781.823 319 47930.758 50337.774 320 47930.758 47930.758 321 47930.758 47930.758 322 47930.758 47930.758 323 -7555.060 47930.758 324 47930.758 -7555.060 325 45207.361 47930.758 326 -7555.060 45207.361 327 -7555.060 -7555.060 328 47930.758 -7555.060 329 47930.758 47930.758 330 47930.758 47930.758 331 47930.758 47930.758 332 47930.758 47930.758 333 46110.459 47930.758 334 48709.764 46110.459 335 40596.813 48709.764 336 34892.012 40596.813 337 47930.758 34892.012 338 45723.816 47930.758 339 53714.507 45723.816 340 20843.505 53714.507 341 47930.758 20843.505 342 47930.758 47930.758 343 47930.758 47930.758 344 41837.252 47930.758 345 47930.758 41837.252 346 38208.833 47930.758 347 47930.758 38208.833 348 47930.758 47930.758 349 11292.552 47930.758 350 51727.086 11292.552 351 48104.426 51727.086 352 47803.348 48104.426 353 40959.878 47803.348 354 47930.758 40959.878 355 47930.758 47930.758 356 25991.297 47930.758 357 47930.758 25991.297 358 47930.758 47930.758 359 47930.758 47930.758 360 49543.695 47930.758 361 47930.758 49543.695 362 47930.758 47930.758 363 -59079.502 47930.758 364 50139.283 -59079.502 365 47930.758 50139.283 366 47930.758 47930.758 367 47930.758 47930.758 368 47930.758 47930.758 369 50155.084 47930.758 370 47930.758 50155.084 371 60651.036 47930.758 372 35414.901 60651.036 373 27297.538 35414.901 374 47930.758 27297.538 375 47930.758 47930.758 376 42833.808 47930.758 377 41821.079 42833.808 378 -7827.369 41821.079 379 47930.758 -7827.369 380 47930.758 47930.758 381 33983.667 47930.758 382 65563.601 33983.667 383 47930.758 65563.601 384 50256.107 47930.758 385 -100883.026 50256.107 386 60199.198 -100883.026 387 21375.207 60199.198 388 46456.029 21375.207 389 34065.162 46456.029 390 43281.192 34065.162 391 17215.934 43281.192 392 -20308.732 17215.934 393 68636.060 -20308.732 394 -180351.327 68636.060 395 9037.666 -180351.327 396 -138046.328 9037.666 397 4172.197 -138046.328 398 10638.634 4172.197 399 -23000.886 10638.634 400 -11465.171 -23000.886 401 39825.868 -11465.171 402 17459.373 39825.868 403 -95882.735 17459.373 404 -12187.230 -95882.735 405 10402.487 -12187.230 406 23566.412 10402.487 407 15549.014 23566.412 408 26690.904 15549.014 409 -73845.864 26690.904 410 10354.181 -73845.864 411 -250216.313 10354.181 412 -74411.743 -250216.313 413 12120.541 -74411.743 414 -189194.453 12120.541 415 -115267.526 -189194.453 416 118575.766 -115267.526 417 -62387.577 118575.766 418 -35757.824 -62387.577 419 -49407.953 -35757.824 420 -60797.317 -49407.953 421 -107098.631 -60797.317 422 -576390.302 -107098.631 423 -484134.906 -576390.302 424 -41774.594 -484134.906 425 -232554.023 -41774.594 426 -313458.634 -232554.023 427 -1118357.214 -313458.634 428 -306374.089 -1118357.214 429 -344582.470 -306374.089 430 -675705.948 -344582.470 > 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/freestat/rcomp/tmp/74cys1291414557.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/freestat/rcomp/tmp/84cys1291414557.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/freestat/rcomp/tmp/94cys1291414557.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/freestat/rcomp/tmp/10x3fd1291414557.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/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/freestat/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/freestat/rcomp/tmp/11i3e11291414557.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/freestat/rcomp/tmp/1234cp1291414557.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/freestat/rcomp/tmp/13zdag1291414557.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/freestat/rcomp/tmp/14le9m1291414557.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/freestat/rcomp/tmp/15oxpa1291414557.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/freestat/rcomp/tmp/16ax6y1291414557.tab") + } > > try(system("convert tmp/1qki21291414557.ps tmp/1qki21291414557.png",intern=TRUE)) character(0) > try(system("convert tmp/2qki21291414557.ps tmp/2qki21291414557.png",intern=TRUE)) character(0) > try(system("convert tmp/3jbi51291414557.ps tmp/3jbi51291414557.png",intern=TRUE)) character(0) > try(system("convert tmp/4jbi51291414557.ps tmp/4jbi51291414557.png",intern=TRUE)) character(0) > try(system("convert tmp/5jbi51291414557.ps tmp/5jbi51291414557.png",intern=TRUE)) character(0) > try(system("convert tmp/6c2hq1291414557.ps tmp/6c2hq1291414557.png",intern=TRUE)) character(0) > try(system("convert tmp/74cys1291414557.ps tmp/74cys1291414557.png",intern=TRUE)) character(0) > try(system("convert tmp/84cys1291414557.ps tmp/84cys1291414557.png",intern=TRUE)) character(0) > try(system("convert tmp/94cys1291414557.ps tmp/94cys1291414557.png",intern=TRUE)) character(0) > try(system("convert tmp/10x3fd1291414557.ps tmp/10x3fd1291414557.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 13.963 3.164 14.377