R version 2.9.0 (2009-04-17) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(162556 + ,807 + ,213118 + ,6282154 + ,56289 + ,29790 + ,444 + ,81767 + ,4321023 + ,28328 + ,87550 + ,412 + ,153198 + ,4111912 + ,17936 + ,84738 + ,428 + ,-26007 + ,223193 + ,4145 + ,54660 + ,315 + ,126942 + ,1491348 + ,3040 + ,42634 + ,168 + ,157214 + ,1629616 + ,2964 + ,40949 + ,263 + ,129352 + ,1398893 + ,2865 + ,45187 + ,267 + ,234817 + ,1926517 + ,2854 + ,37704 + ,228 + ,60448 + ,983660 + ,2167 + ,16275 + ,129 + ,47818 + ,1443586 + ,1974 + ,25830 + ,104 + ,245546 + ,1073089 + ,1910 + ,12679 + ,122 + ,48020 + ,984885 + ,1871 + ,18014 + ,393 + ,-1710 + ,1405225 + ,943 + ,43556 + ,190 + ,32648 + ,227132 + ,929 + ,24811 + ,280 + ,95350 + ,929118 + ,822 + ,6575 + ,63 + ,151352 + ,1071292 + ,819 + ,7123 + ,102 + ,288170 + ,638830 + ,769 + ,21950 + ,265 + ,114337 + ,856956 + ,745 + ,37597 + ,234 + ,37884 + ,992426 + ,652 + ,17821 + ,277 + ,122844 + ,444477 + ,643 + ,12988 + ,73 + ,82340 + ,857217 + ,601 + ,22330 + ,67 + ,79801 + ,711969 + ,446 + ,13326 + ,103 + ,165548 + ,702380 + ,436 + ,16189 + ,290 + ,116384 + ,358589 + ,379 + ,7146 + ,83 + ,134028 + ,297978 + ,305 + ,15824 + ,56 + ,63838 + ,585715 + ,284 + ,27664 + ,236 + ,74996 + ,657954 + ,247 + ,11920 + ,73 + ,31080 + ,209458 + ,238 + ,8568 + ,34 + ,32168 + ,786690 + ,223 + ,14416 + ,139 + ,49857 + ,439798 + ,220 + ,3369 + ,26 + ,87161 + ,688779 + ,217 + ,11819 + ,70 + ,106113 + ,574339 + ,199 + ,6984 + ,40 + ,80570 + ,741409 + ,195 + ,4519 + ,42 + ,102129 + ,597793 + ,163 + ,2220 + ,12 + ,301670 + ,644190 + ,154 + ,18562 + ,211 + ,102313 + ,377934 + ,143 + ,10327 + ,74 + ,88577 + ,640273 + ,135 + ,5336 + ,80 + ,112477 + ,697458 + ,132 + ,2365 + ,83 + ,191778 + ,550608 + ,120 + ,4069 + ,131 + ,79804 + ,207393 + ,119 + ,8636 + ,203 + ,128294 + ,301607 + ,108 + ,13718 + ,56 + ,96448 + ,345783 + ,104 + ,4525 + ,89 + ,93811 + ,501749 + ,101 + ,6869 + ,88 + ,117520 + ,379983 + ,90 + ,4628 + ,39 + ,69159 + ,387475 + ,85 + ,3689 + ,25 + ,101792 + ,377305 + ,74 + ,4891 + ,49 + ,210568 + ,370837 + ,73 + ,7489 + ,149 + ,136996 + ,430866 + ,70 + ,4901 + ,58 + ,121920 + ,469107 + ,67 + ,2284 + ,41 + ,76403 + ,194493 + ,66 + ,3160 + ,90 + ,108094 + ,530670 + ,66 + ,4150 + ,136 + ,134759 + ,518365 + ,66 + ,7285 + ,97 + ,188873 + ,491303 + ,59 + ,1134 + ,63 + ,146216 + ,527021 + ,58 + ,4658 + ,114 + ,156608 + ,233773 + ,58 + ,2384 + ,77 + ,61348 + ,405972 + ,54 + ,3748 + ,6 + ,50350 + ,652925 + ,53 + ,5371 + ,47 + ,87720 + ,446211 + ,49 + ,1285 + ,51 + ,99489 + ,341340 + ,49 + ,9327 + ,85 + ,87419 + ,387699 + ,44 + ,5565 + ,43 + ,94355 + ,493408 + ,39 + ,1528 + ,32 + ,60326 + ,146494 + ,38 + ,3122 + ,25 + ,94670 + ,414462 + ,38 + ,7561 + ,77 + ,82425 + ,364304 + ,37 + ,2675 + ,54 + ,59017 + ,355178 + ,36 + ,13253 + ,251 + ,90829 + ,357760 + ,35 + ,880 + ,15 + ,80791 + ,261216 + ,34 + ,2053 + ,44 + ,100423 + ,397144 + ,33 + ,1424 + ,73 + ,131116 + ,374943 + ,32 + ,4036 + ,85 + ,100269 + ,424898 + ,32 + ,3045 + ,49 + ,27330 + ,202055 + ,31 + ,5119 + ,38 + ,39039 + ,378525 + ,31 + ,1431 + ,35 + ,106885 + ,310768 + ,30 + ,554 + ,9 + ,79285 + ,325738 + ,30 + ,1975 + ,34 + ,118881 + ,394510 + ,30 + ,1765 + ,20 + ,77623 + ,247060 + ,28 + ,1012 + ,29 + ,114768 + ,368078 + ,26 + ,810 + ,11 + ,74015 + ,236761 + ,26 + ,1280 + ,52 + ,69465 + ,312378 + ,26 + ,666 + ,13 + ,117869 + ,339836 + ,25 + ,1380 + ,29 + ,60982 + ,347385 + ,25 + ,4677 + ,66 + ,90131 + ,426280 + ,24 + ,876 + ,33 + ,138971 + ,352850 + ,24 + ,814 + ,15 + ,39625 + ,301881 + ,23 + ,514 + ,15 + ,102725 + ,377516 + ,23 + ,5692 + ,68 + ,64239 + ,357312 + ,22 + ,3642 + ,100 + ,90262 + ,458343 + ,22 + ,540 + ,13 + ,103960 + ,354228 + ,22 + ,2099 + ,45 + ,106611 + ,308636 + ,20 + ,567 + ,14 + ,103345 + ,386212 + ,20 + ,2001 + ,36 + ,95551 + ,393343 + ,19 + ,2949 + ,40 + ,82903 + ,378509 + ,18 + ,2253 + ,68 + ,63593 + ,452469 + ,17 + ,6533 + ,29 + ,126910 + ,364839 + ,16 + ,1889 + ,43 + ,37527 + ,358649 + ,14 + ,3055 + ,30 + ,60247 + ,376641 + ,14 + ,272 + ,9 + ,112995 + ,429112 + ,14 + ,1414 + ,22 + ,70184 + ,330546 + ,13 + ,2564 + ,19 + ,130140 + ,403560 + ,13 + ,1383 + ,9 + ,73221 + ,317892 + ,13 + ,1261 + ,31 + ,76114 + ,307528 + ,13 + ,975 + ,19 + ,90534 + ,235133 + ,13 + ,3366 + ,55 + ,108479 + ,299243 + ,12 + ,576 + ,8 + ,113761 + ,314073 + ,12 + ,1686 + ,28 + ,68696 + ,368186 + ,11 + ,746 + ,29 + ,71561 + ,269661 + ,10 + ,3192 + ,48 + ,59831 + ,125390 + ,10 + ,2045 + ,16 + ,97890 + ,510834 + ,10 + ,5702 + ,47 + ,101481 + ,321896 + ,10 + ,1932 + ,20 + ,72954 + ,249898 + ,10 + ,936 + ,22 + ,67939 + ,408881 + ,9 + ,3437 + ,33 + ,48022 + ,158492 + ,9 + ,5131 + ,44 + ,86111 + ,292154 + ,9 + ,2397 + ,13 + ,74020 + ,289513 + ,6 + ,1389 + ,6 + ,57530 + ,378049 + ,6 + ,1503 + ,35 + ,56364 + ,343466 + ,6 + ,402 + ,8 + ,84990 + ,332743 + ,5 + ,2239 + ,17 + ,88590 + ,442882 + ,5 + ,2234 + ,11 + ,77200 + ,214215 + ,5 + ,837 + ,21 + ,61262 + ,315688 + ,4 + ,10579 + ,92 + ,110309 + ,375195 + ,4 + ,875 + ,12 + ,67000 + ,334280 + ,3 + ,1585 + ,112 + ,93099 + ,355864 + ,2 + ,1659 + ,25 + ,107577 + ,480382 + ,2 + ,2647 + ,17 + ,62920 + ,353058 + ,1 + ,3294 + ,23 + ,75832 + ,217193 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,94 + ,10 + ,60793 + ,314533 + ,0 + ,422 + ,23 + ,57935 + ,318056 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,34 + ,7 + ,60630 + ,314353 + ,0 + ,1558 + ,25 + ,55637 + ,369448 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,43 + ,20 + ,60887 + ,312846 + ,0 + ,645 + ,4 + ,60720 + ,312075 + ,0 + ,316 + ,4 + ,60505 + ,315009 + ,0 + ,115 + ,10 + ,60945 + ,318903 + ,0 + ,5 + ,1 + ,60720 + ,314887 + ,0 + ,897 + ,4 + ,60720 + ,314913 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,389 + ,8 + ,58990 + ,325506 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1002 + ,11 + ,56750 + ,298568 + ,0 + ,36 + ,4 + ,60894 + ,315834 + ,0 + ,460 + ,15 + ,63346 + ,329784 + ,0 + ,309 + ,9 + ,56535 + ,312878 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,9 + ,7 + ,60835 + ,314987 + ,0 + ,271 + ,2 + ,60720 + ,325249 + ,0 + ,14 + ,0 + ,61016 + ,315877 + ,0 + ,520 + ,7 + ,58650 + ,291650 + ,0 + ,1766 + ,46 + ,60438 + ,305959 + ,0 + ,0 + ,5 + ,60720 + ,315380 + ,0 + ,458 + ,7 + ,58625 + ,297765 + ,0 + ,20 + ,2 + ,60938 + ,315245 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,98 + ,2 + ,61490 + ,315236 + ,0 + ,405 + ,5 + ,60845 + ,336425 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,483 + ,7 + ,60830 + ,306268 + ,0 + ,454 + ,24 + ,63261 + ,302187 + ,0 + ,47 + ,1 + ,60720 + ,314882 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,757 + ,18 + ,45689 + ,382712 + ,0 + ,4655 + ,55 + ,60720 + ,341570 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,36 + ,3 + ,61564 + ,312412 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,203 + ,9 + ,61938 + ,309596 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,126 + ,8 + ,60951 + ,315547 + ,0 + ,400 + ,113 + ,60720 + ,313267 + ,0 + ,71 + ,0 + ,60745 + ,316176 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,972 + ,19 + ,71642 + ,359335 + ,0 + ,531 + ,11 + ,71641 + ,330068 + ,0 + ,2461 + ,25 + ,55792 + ,314289 + ,0 + ,378 + ,16 + ,71873 + ,297413 + ,0 + ,23 + ,5 + ,62555 + ,314806 + ,0 + ,638 + ,11 + ,60370 + ,333210 + ,0 + ,2300 + ,23 + ,64873 + ,352108 + ,0 + ,149 + ,6 + ,62041 + ,313332 + ,0 + ,226 + ,5 + ,65745 + ,291787 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,275 + ,7 + ,59500 + ,318745 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,141 + ,7 + ,61630 + ,315366 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,28 + ,3 + ,60890 + ,315688 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,4980 + ,89 + ,113521 + ,409642 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,472 + ,19 + ,80045 + ,269587 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,203 + ,12 + ,50804 + ,300962 + ,0 + ,496 + ,12 + ,87390 + ,325479 + ,0 + ,10 + ,5 + ,61656 + ,316155 + ,0 + ,63 + ,2 + ,65688 + ,318574 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1136 + ,26 + ,48522 + ,343613 + ,0 + ,265 + ,3 + ,60720 + ,306948 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,267 + ,11 + ,57640 + ,330059 + ,0 + ,474 + ,10 + ,61977 + ,288985 + ,0 + ,534 + ,5 + ,62620 + ,304485 + ,0 + ,0 + ,2 + ,60720 + ,315380 + ,0 + ,15 + ,6 + ,60831 + ,315688 + ,0 + ,397 + ,7 + ,60646 + ,317736 + ,0 + ,0 + ,2 + ,60720 + ,315380 + ,0 + ,1866 + ,28 + ,56225 + ,322331 + ,0 + ,288 + ,3 + ,60510 + ,296656 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,3 + ,1 + ,60698 + ,315354 + ,0 + ,468 + ,20 + ,60720 + ,312161 + ,0 + ,20 + ,1 + ,60805 + ,315576 + ,0 + ,278 + ,22 + ,61404 + ,314922 + ,0 + ,61 + ,9 + ,60720 + ,314551 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,192 + ,2 + ,65276 + ,312339 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,317 + ,7 + ,63915 + ,298700 + ,0 + ,738 + ,9 + ,60720 + ,321376 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,368 + ,13 + ,61686 + ,303230 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,2 + ,0 + ,60743 + ,315487 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,53 + ,6 + ,60349 + ,315793 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,94 + ,3 + ,61360 + ,312887 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,24 + ,7 + ,59818 + ,315637 + ,0 + ,2332 + ,2 + ,72680 + ,324385 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,131 + ,15 + ,61808 + ,308989 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,206 + ,9 + ,53110 + ,296702 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,167 + ,1 + ,64245 + ,307322 + ,0 + ,622 + ,38 + ,73007 + ,304376 + ,0 + ,2328 + ,57 + ,82732 + ,253588 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,365 + ,7 + ,54820 + ,309560 + ,0 + ,364 + ,26 + ,47705 + ,298466 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,226 + ,13 + ,72835 + ,343929 + ,0 + ,307 + ,10 + ,58856 + ,331955 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,188 + ,9 + ,77655 + ,381180 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,138 + ,26 + ,69817 + ,331420 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,125 + ,19 + ,60798 + ,310201 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,282 + ,12 + ,62452 + ,320016 + ,0 + ,335 + ,23 + ,64175 + ,320398 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1324 + ,29 + ,67440 + ,291841 + ,0 + ,176 + ,8 + ,68136 + ,310670 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,249 + ,26 + ,56726 + ,313491 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,333 + ,9 + ,70811 + ,331323 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,601 + ,5 + ,60720 + ,319210 + ,0 + ,30 + ,3 + ,62045 + ,318098 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,249 + ,13 + ,54323 + ,292754 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,165 + ,12 + ,62841 + ,325176 + ,0 + ,453 + ,19 + ,81125 + ,365959 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,53 + ,10 + ,59506 + ,302409 + ,0 + ,382 + ,9 + ,59365 + ,340968 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,9 + ,60720 + ,315380 + ,0 + ,30 + ,4 + ,60798 + ,313164 + ,0 + ,290 + ,1 + ,58790 + ,301164 + ,0 + ,0 + ,1 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,366 + ,14 + ,61808 + ,344425 + ,0 + ,2 + ,12 + ,60735 + ,315394 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,209 + ,19 + ,64016 + ,316647 + ,0 + ,384 + ,17 + ,54683 + ,309836 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,365 + ,32 + ,87192 + ,346611 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,49 + ,14 + ,64107 + ,322031 + ,0 + ,3 + ,8 + ,60761 + ,315656 + ,0 + ,133 + ,4 + ,65990 + ,339445 + ,0 + ,32 + ,0 + ,59988 + ,314964 + ,0 + ,368 + ,20 + ,61167 + ,297141 + ,0 + ,1 + ,5 + ,60719 + ,315372 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,22 + ,1 + ,60722 + ,312502 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,96 + ,4 + ,60379 + ,313729 + ,0 + ,1 + ,1 + ,60727 + ,315388 + ,0 + ,314 + ,4 + ,60720 + ,315371 + ,0 + ,844 + ,20 + ,60925 + ,296139 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,26 + ,1 + ,60896 + ,313880 + ,0 + ,125 + ,10 + ,59734 + ,317698 + ,0 + ,304 + ,12 + ,62969 + ,295580 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,621 + ,13 + ,60720 + ,308256 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,119 + ,3 + ,59118 + ,303677 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1595 + ,10 + ,60720 + ,319369 + ,0 + ,312 + ,3 + ,58598 + ,318690 + ,0 + ,60 + ,7 + ,61124 + ,314049 + ,0 + ,587 + ,10 + ,59595 + ,325699 + ,0 + ,135 + ,1 + ,62065 + ,314210 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,514 + ,15 + ,78780 + ,322378 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,4 + ,60722 + ,315398 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1763 + ,28 + ,61600 + ,308336 + ,0 + ,180 + ,9 + ,59635 + ,316386 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,218 + ,7 + ,60720 + ,315553 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,448 + ,7 + ,59781 + ,323361 + ,0 + ,227 + ,7 + ,76644 + ,336639 + ,0 + ,174 + ,3 + ,64820 + ,307424 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,121 + ,11 + ,56178 + ,295370 + ,0 + ,607 + ,7 + ,60436 + ,322340 + ,0 + ,2212 + ,10 + ,60720 + ,319864 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,530 + ,18 + ,73433 + ,317291 + ,0 + ,571 + ,14 + ,41477 + ,280398 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,78 + ,12 + ,62700 + ,317330 + ,0 + ,2489 + ,29 + ,67804 + ,238125 + ,0 + ,131 + ,3 + ,59661 + ,327071 + ,0 + ,923 + ,6 + ,58620 + ,309038 + ,0 + ,72 + ,3 + ,60398 + ,314210 + ,0 + ,572 + ,8 + ,58580 + ,307930 + ,-1 + ,397 + ,10 + ,62710 + ,322327 + ,-1 + ,450 + ,6 + ,59325 + ,292136 + ,-1 + ,622 + ,8 + ,60950 + ,263276 + ,-2 + ,694 + ,6 + ,68060 + ,367655 + ,-2 + ,3425 + ,9 + ,83620 + ,283910 + ,-3 + ,562 + ,8 + ,58456 + ,283587 + ,-3 + ,4917 + ,26 + ,52811 + ,243650 + ,-3 + ,1442 + ,239 + ,121173 + ,438493 + ,-3 + ,529 + ,7 + ,63870 + ,296261 + ,-3 + ,2126 + ,41 + ,21001 + ,230621 + ,-3 + ,1061 + ,3 + ,70415 + ,304252 + ,-4 + ,776 + ,8 + ,64230 + ,333505 + ,-4 + ,611 + ,6 + ,59190 + ,296919 + ,-4 + ,1526 + ,21 + ,69351 + ,278990 + ,-4 + ,592 + ,7 + ,64270 + ,276898 + ,-4 + ,1182 + ,11 + ,70694 + ,327007 + ,-5 + ,621 + ,11 + ,68005 + ,317046 + ,-5 + ,989 + ,12 + ,58930 + ,304555 + ,-5 + ,438 + ,9 + ,58320 + ,298096 + ,-5 + ,726 + ,3 + ,69980 + ,231861 + ,-5 + ,1303 + ,57 + ,69863 + ,309422 + ,-5 + ,7419 + ,21 + ,63255 + ,286963 + ,-5 + ,1164 + ,15 + ,57320 + ,269753 + ,-6 + ,3310 + ,32 + ,75230 + ,448243 + ,-6 + ,1920 + ,11 + ,79420 + ,165404 + ,-6 + ,965 + ,2 + ,73490 + ,204325 + ,-8 + ,3256 + ,23 + ,35250 + ,407159 + ,-8 + ,1135 + ,20 + ,62285 + ,290476 + ,-9 + ,1270 + ,24 + ,69206 + ,275311 + ,-10 + ,661 + ,1 + ,65920 + ,246541 + ,-12 + ,1013 + ,1 + ,69770 + ,253468 + ,-14 + ,2844 + ,74 + ,72683 + ,240897 + ,-15 + ,11528 + ,68 + ,-14545 + ,-83265 + ,-25 + ,6526 + ,20 + ,55830 + ,-42143 + ,-28 + ,2264 + ,20 + ,55174 + ,272713 + ,-31 + ,5109 + ,82 + ,67038 + ,215362 + ,-31 + ,3999 + ,21 + ,51252 + ,42754 + ,-58 + ,35624 + ,244 + ,157278 + ,306275 + ,-70 + ,9252 + ,32 + ,79510 + ,253537 + ,-129 + ,15236 + ,86 + ,77440 + ,372631 + ,-158 + ,18073 + ,69 + ,27284 + ,-7170 + ,-366) + ,dim=c(5 + ,431) + ,dimnames=list(c('kosten' + ,'orders' + ,'dividenden' + ,'rijkdom' + ,'CC-score') + ,1:431)) > y <- array(NA,dim=c(5,431),dimnames=list(c('kosten','orders','dividenden','rijkdom','CC-score'),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 CC-score kosten orders dividenden rijkdom 1 56289 162556 807 213118 6282154 2 28328 29790 444 81767 4321023 3 17936 87550 412 153198 4111912 4 4145 84738 428 -26007 223193 5 3040 54660 315 126942 1491348 6 2964 42634 168 157214 1629616 7 2865 40949 263 129352 1398893 8 2854 45187 267 234817 1926517 9 2167 37704 228 60448 983660 10 1974 16275 129 47818 1443586 11 1910 25830 104 245546 1073089 12 1871 12679 122 48020 984885 13 943 18014 393 -1710 1405225 14 929 43556 190 32648 227132 15 822 24811 280 95350 929118 16 819 6575 63 151352 1071292 17 769 7123 102 288170 638830 18 745 21950 265 114337 856956 19 652 37597 234 37884 992426 20 643 17821 277 122844 444477 21 601 12988 73 82340 857217 22 446 22330 67 79801 711969 23 436 13326 103 165548 702380 24 379 16189 290 116384 358589 25 305 7146 83 134028 297978 26 284 15824 56 63838 585715 27 247 27664 236 74996 657954 28 238 11920 73 31080 209458 29 223 8568 34 32168 786690 30 220 14416 139 49857 439798 31 217 3369 26 87161 688779 32 199 11819 70 106113 574339 33 195 6984 40 80570 741409 34 163 4519 42 102129 597793 35 154 2220 12 301670 644190 36 143 18562 211 102313 377934 37 135 10327 74 88577 640273 38 132 5336 80 112477 697458 39 120 2365 83 191778 550608 40 119 4069 131 79804 207393 41 108 8636 203 128294 301607 42 104 13718 56 96448 345783 43 101 4525 89 93811 501749 44 90 6869 88 117520 379983 45 85 4628 39 69159 387475 46 74 3689 25 101792 377305 47 73 4891 49 210568 370837 48 70 7489 149 136996 430866 49 67 4901 58 121920 469107 50 66 2284 41 76403 194493 51 66 3160 90 108094 530670 52 66 4150 136 134759 518365 53 59 7285 97 188873 491303 54 58 1134 63 146216 527021 55 58 4658 114 156608 233773 56 54 2384 77 61348 405972 57 53 3748 6 50350 652925 58 49 5371 47 87720 446211 59 49 1285 51 99489 341340 60 44 9327 85 87419 387699 61 39 5565 43 94355 493408 62 38 1528 32 60326 146494 63 38 3122 25 94670 414462 64 37 7561 77 82425 364304 65 36 2675 54 59017 355178 66 35 13253 251 90829 357760 67 34 880 15 80791 261216 68 33 2053 44 100423 397144 69 32 1424 73 131116 374943 70 32 4036 85 100269 424898 71 31 3045 49 27330 202055 72 31 5119 38 39039 378525 73 30 1431 35 106885 310768 74 30 554 9 79285 325738 75 30 1975 34 118881 394510 76 28 1765 20 77623 247060 77 26 1012 29 114768 368078 78 26 810 11 74015 236761 79 26 1280 52 69465 312378 80 25 666 13 117869 339836 81 25 1380 29 60982 347385 82 24 4677 66 90131 426280 83 24 876 33 138971 352850 84 23 814 15 39625 301881 85 23 514 15 102725 377516 86 22 5692 68 64239 357312 87 22 3642 100 90262 458343 88 22 540 13 103960 354228 89 20 2099 45 106611 308636 90 20 567 14 103345 386212 91 19 2001 36 95551 393343 92 18 2949 40 82903 378509 93 17 2253 68 63593 452469 94 16 6533 29 126910 364839 95 14 1889 43 37527 358649 96 14 3055 30 60247 376641 97 14 272 9 112995 429112 98 13 1414 22 70184 330546 99 13 2564 19 130140 403560 100 13 1383 9 73221 317892 101 13 1261 31 76114 307528 102 13 975 19 90534 235133 103 12 3366 55 108479 299243 104 12 576 8 113761 314073 105 11 1686 28 68696 368186 106 10 746 29 71561 269661 107 10 3192 48 59831 125390 108 10 2045 16 97890 510834 109 10 5702 47 101481 321896 110 10 1932 20 72954 249898 111 9 936 22 67939 408881 112 9 3437 33 48022 158492 113 9 5131 44 86111 292154 114 6 2397 13 74020 289513 115 6 1389 6 57530 378049 116 6 1503 35 56364 343466 117 5 402 8 84990 332743 118 5 2239 17 88590 442882 119 5 2234 11 77200 214215 120 4 837 21 61262 315688 121 4 10579 92 110309 375195 122 3 875 12 67000 334280 123 2 1585 112 93099 355864 124 2 1659 25 107577 480382 125 1 2647 17 62920 353058 126 1 3294 23 75832 217193 127 0 0 0 60720 315380 128 0 94 10 60793 314533 129 0 422 23 57935 318056 130 0 0 0 60720 315380 131 0 34 7 60630 314353 132 0 1558 25 55637 369448 133 0 0 0 60720 315380 134 0 43 20 60887 312846 135 0 645 4 60720 312075 136 0 316 4 60505 315009 137 0 115 10 60945 318903 138 0 5 1 60720 314887 139 0 897 4 60720 314913 140 0 0 0 60720 315380 141 0 389 8 58990 325506 142 0 0 0 60720 315380 143 0 1002 11 56750 298568 144 0 36 4 60894 315834 145 0 460 15 63346 329784 146 0 309 9 56535 312878 147 0 0 0 60720 315380 148 0 9 7 60835 314987 149 0 271 2 60720 325249 150 0 14 0 61016 315877 151 0 520 7 58650 291650 152 0 1766 46 60438 305959 153 0 0 5 60720 315380 154 0 458 7 58625 297765 155 0 20 2 60938 315245 156 0 0 0 60720 315380 157 0 0 0 60720 315380 158 0 98 2 61490 315236 159 0 405 5 60845 336425 160 0 0 0 60720 315380 161 0 0 0 60720 315380 162 0 0 0 60720 315380 163 0 0 0 60720 315380 164 0 483 7 60830 306268 165 0 454 24 63261 302187 166 0 47 1 60720 314882 167 0 0 0 60720 315380 168 0 757 18 45689 382712 169 0 4655 55 60720 341570 170 0 0 0 60720 315380 171 0 0 0 60720 315380 172 0 36 3 61564 312412 173 0 0 0 60720 315380 174 0 203 9 61938 309596 175 0 0 0 60720 315380 176 0 126 8 60951 315547 177 0 400 113 60720 313267 178 0 71 0 60745 316176 179 0 0 0 60720 315380 180 0 0 0 60720 315380 181 0 972 19 71642 359335 182 0 531 11 71641 330068 183 0 2461 25 55792 314289 184 0 378 16 71873 297413 185 0 23 5 62555 314806 186 0 638 11 60370 333210 187 0 2300 23 64873 352108 188 0 149 6 62041 313332 189 0 226 5 65745 291787 190 0 0 0 60720 315380 191 0 275 7 59500 318745 192 0 0 0 60720 315380 193 0 141 7 61630 315366 194 0 0 0 60720 315380 195 0 28 3 60890 315688 196 0 0 0 60720 315380 197 0 4980 89 113521 409642 198 0 0 0 60720 315380 199 0 0 0 60720 315380 200 0 472 19 80045 269587 201 0 0 0 60720 315380 202 0 0 0 60720 315380 203 0 0 0 60720 315380 204 0 203 12 50804 300962 205 0 496 12 87390 325479 206 0 10 5 61656 316155 207 0 63 2 65688 318574 208 0 0 0 60720 315380 209 0 1136 26 48522 343613 210 0 265 3 60720 306948 211 0 0 60720 315380 0 212 267 0 60720 315380 0 213 474 11 57640 330059 0 214 534 10 61977 288985 0 215 0 5 62620 304485 0 216 15 2 60720 315380 0 217 397 6 60831 315688 0 218 0 7 60646 317736 0 219 1866 2 60720 315380 0 220 288 28 56225 322331 0 221 0 3 60510 296656 0 222 3 0 60720 315380 0 223 468 1 60698 315354 0 224 20 20 60720 312161 0 225 278 1 60805 315576 0 226 61 22 61404 314922 0 227 0 9 60720 314551 0 228 192 0 60720 315380 0 229 0 2 65276 312339 0 230 317 0 60720 315380 0 231 738 7 63915 298700 0 232 0 9 60720 321376 0 233 368 0 60720 315380 0 234 0 13 61686 303230 0 235 2 0 60720 315380 0 236 0 0 60743 315487 0 237 53 0 60720 315380 0 238 0 6 60349 315793 0 239 0 0 60720 315380 0 240 0 0 60720 315380 0 241 94 0 60720 315380 0 242 0 3 61360 312887 0 243 24 0 60720 315380 0 244 2332 7 59818 315637 0 245 0 2 72680 324385 0 246 0 0 60720 315380 0 247 131 0 60720 315380 0 248 0 15 61808 308989 0 249 0 0 60720 315380 0 250 206 0 60720 315380 0 251 0 9 53110 296702 0 252 167 0 60720 315380 0 253 622 1 64245 307322 0 254 2328 38 73007 304376 0 255 0 57 82732 253588 0 256 365 0 60720 315380 0 257 364 7 54820 309560 0 258 0 26 47705 298466 0 259 0 0 60720 315380 0 260 0 0 60720 315380 0 261 0 0 60720 315380 0 262 226 0 60720 315380 0 263 307 13 72835 343929 0 264 0 10 58856 331955 0 265 0 0 60720 315380 0 266 0 0 60720 315380 0 267 188 0 60720 315380 0 268 0 9 77655 381180 0 269 138 0 60720 315380 0 270 0 26 69817 331420 0 271 0 0 60720 315380 0 272 0 0 60720 315380 0 273 125 0 60720 315380 0 274 0 19 60798 310201 0 275 282 0 60720 315380 0 276 335 12 62452 320016 0 277 0 23 64175 320398 0 278 1324 0 60720 315380 0 279 176 29 67440 291841 0 280 0 8 68136 310670 0 281 0 0 60720 315380 0 282 249 0 60720 315380 0 283 0 26 56726 313491 0 284 333 0 60720 315380 0 285 0 9 70811 331323 0 286 601 0 60720 315380 0 287 30 5 60720 319210 0 288 0 3 62045 318098 0 289 249 0 60720 315380 0 290 0 13 54323 292754 0 291 165 0 60720 315380 0 292 453 12 62841 325176 0 293 0 19 81125 365959 0 294 53 0 60720 315380 0 295 382 10 59506 302409 0 296 0 9 59365 340968 0 297 0 0 60720 315380 0 298 0 0 60720 315380 0 299 0 0 60720 315380 0 300 30 9 60720 315380 0 301 290 4 60798 313164 0 302 0 1 58790 301164 0 303 0 1 60720 315380 0 304 366 0 60720 315380 0 305 2 14 61808 344425 0 306 0 12 60735 315394 0 307 209 0 60720 315380 0 308 384 19 64016 316647 0 309 0 17 54683 309836 0 310 0 0 60720 315380 0 311 365 0 60720 315380 0 312 0 32 87192 346611 0 313 49 0 60720 315380 0 314 3 14 64107 322031 0 315 133 8 60761 315656 0 316 32 4 65990 339445 0 317 368 0 59988 314964 0 318 1 20 61167 297141 0 319 0 5 60719 315372 0 320 0 0 60720 315380 0 321 0 0 60720 315380 0 322 0 0 60720 315380 0 323 0 0 60720 315380 0 324 0 0 60720 315380 0 325 22 0 60720 315380 0 326 0 1 60722 312502 0 327 0 0 60720 315380 0 328 0 0 60720 315380 0 329 0 0 60720 315380 0 330 0 0 60720 315380 0 331 0 0 60720 315380 0 332 0 0 60720 315380 0 333 96 0 60720 315380 0 334 1 4 60379 313729 0 335 314 1 60727 315388 0 336 844 4 60720 315371 0 337 0 20 60925 296139 0 338 26 0 60720 315380 0 339 125 1 60896 313880 0 340 304 10 59734 317698 0 341 0 12 62969 295580 0 342 0 0 60720 315380 0 343 0 0 60720 315380 0 344 621 0 60720 315380 0 345 0 13 60720 308256 0 346 119 0 60720 315380 0 347 0 3 59118 303677 0 348 0 0 60720 315380 0 349 1595 0 60720 315380 0 350 312 10 60720 319369 0 351 60 3 58598 318690 0 352 587 7 61124 314049 0 353 135 10 59595 325699 0 354 0 1 62065 314210 0 355 0 0 60720 315380 0 356 514 0 60720 315380 0 357 0 15 78780 322378 0 358 0 0 60720 315380 0 359 0 0 60720 315380 0 360 1 0 60720 315380 0 361 0 4 60722 315398 0 362 0 0 60720 315380 0 363 1763 0 60720 315380 0 364 180 28 61600 308336 0 365 0 9 59635 316386 0 366 0 0 60720 315380 0 367 0 0 60720 315380 0 368 0 0 60720 315380 0 369 218 0 60720 315380 0 370 0 7 60720 315553 0 371 448 0 60720 315380 0 372 227 7 59781 323361 0 373 174 7 76644 336639 0 374 0 3 64820 307424 0 375 0 0 60720 315380 0 376 121 0 60720 315380 0 377 607 11 56178 295370 0 378 2212 7 60436 322340 0 379 0 10 60720 319864 0 380 0 0 60720 315380 0 381 530 0 60720 315380 0 382 571 18 73433 317291 0 383 0 14 41477 280398 0 384 78 0 60720 315380 0 385 2489 12 62700 317330 0 386 131 29 67804 238125 0 387 923 3 59661 327071 0 388 72 6 58620 309038 0 389 572 3 60398 314210 0 390 397 8 58580 307930 -1 391 450 10 62710 322327 -1 392 622 6 59325 292136 -1 393 694 8 60950 263276 -2 394 3425 6 68060 367655 -2 395 562 9 83620 283910 -3 396 4917 8 58456 283587 -3 397 1442 26 52811 243650 -3 398 529 239 121173 438493 -3 399 2126 7 63870 296261 -3 400 1061 41 21001 230621 -3 401 776 3 70415 304252 -4 402 611 8 64230 333505 -4 403 1526 6 59190 296919 -4 404 592 21 69351 278990 -4 405 1182 7 64270 276898 -4 406 621 11 70694 327007 -5 407 989 11 68005 317046 -5 408 438 12 58930 304555 -5 409 726 9 58320 298096 -5 410 1303 3 69980 231861 -5 411 7419 57 69863 309422 -5 412 1164 21 63255 286963 -5 413 3310 15 57320 269753 -6 414 1920 32 75230 448243 -6 415 965 11 79420 165404 -6 416 3256 2 73490 204325 -8 417 1135 23 35250 407159 -8 418 1270 20 62285 290476 -9 419 661 24 69206 275311 -10 420 1013 1 65920 246541 -12 421 2844 1 69770 253468 -14 422 11528 74 72683 240897 -15 423 6526 68 -14545 -83265 -25 424 2264 20 55830 -42143 -28 425 5109 20 55174 272713 -31 426 3999 82 67038 215362 -31 427 35624 21 51252 42754 -58 428 9252 244 157278 306275 -70 429 15236 32 79510 253537 -129 430 18073 86 77440 372631 -158 431 162556 69 27284 -7170 -366 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) kosten orders dividenden rijkdom 2.225e+03 6.075e-02 2.215e-01 -4.925e-02 5.793e-03 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -14401.8 -1061.7 -144.3 182.6 153933.4 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 2.225e+03 9.701e+02 2.294 0.022287 * kosten 6.075e-02 5.822e-02 1.043 0.297355 orders 2.215e-01 3.277e-02 6.758 4.62e-11 *** dividenden -4.925e-02 8.248e-03 -5.971 4.96e-09 *** rijkdom 5.793e-03 1.557e-03 3.721 0.000225 *** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 7992 on 426 degrees of freedom Multiple R-squared: 0.1639, Adjusted R-squared: 0.156 F-statistic: 20.88 on 4 and 426 DF, p-value: 9.929e-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,] 2.125167e-01 4.250334e-01 0.7874833 [2,] 1.298748e-01 2.597497e-01 0.8701252 [3,] 8.070540e-02 1.614108e-01 0.9192946 [4,] 1.132189e-01 2.264379e-01 0.8867811 [5,] 8.523411e-02 1.704682e-01 0.9147659 [6,] 7.627887e-02 1.525577e-01 0.9237211 [7,] 5.666892e-02 1.133378e-01 0.9433311 [8,] 3.196655e-02 6.393311e-02 0.9680334 [9,] 2.288814e-02 4.577629e-02 0.9771119 [10,] 2.338291e-02 4.676583e-02 0.9766171 [11,] 1.289508e-02 2.579017e-02 0.9871049 [12,] 7.005531e-03 1.401106e-02 0.9929945 [13,] 4.044243e-03 8.088486e-03 0.9959558 [14,] 2.606657e-03 5.213314e-03 0.9973933 [15,] 1.593433e-03 3.186867e-03 0.9984066 [16,] 9.500028e-04 1.900006e-03 0.9990500 [17,] 5.234383e-04 1.046877e-03 0.9994766 [18,] 4.890164e-04 9.780328e-04 0.9995110 [19,] 3.064900e-04 6.129800e-04 0.9996935 [20,] 1.488652e-04 2.977304e-04 0.9998511 [21,] 1.353729e-04 2.707459e-04 0.9998646 [22,] 6.953701e-05 1.390740e-04 0.9999305 [23,] 3.868887e-05 7.737774e-05 0.9999613 [24,] 2.109271e-05 4.218543e-05 0.9999789 [25,] 1.098999e-05 2.197998e-05 0.9999890 [26,] 5.324534e-06 1.064907e-05 0.9999947 [27,] 2.818646e-06 5.637291e-06 0.9999972 [28,] 1.829723e-06 3.659447e-06 0.9999982 [29,] 9.085901e-07 1.817180e-06 0.9999991 [30,] 4.203005e-07 8.406010e-07 0.9999996 [31,] 1.982056e-07 3.964112e-07 0.9999998 [32,] 1.082174e-07 2.164347e-07 0.9999999 [33,] 7.587511e-08 1.517502e-07 0.9999999 [34,] 3.899582e-08 7.799165e-08 1.0000000 [35,] 2.102492e-08 4.204985e-08 1.0000000 [36,] 1.008029e-08 2.016058e-08 1.0000000 [37,] 5.139158e-09 1.027832e-08 1.0000000 [38,] 2.736776e-09 5.473553e-09 1.0000000 [39,] 1.463929e-09 2.927857e-09 1.0000000 [40,] 7.785330e-10 1.557066e-09 1.0000000 [41,] 3.451712e-10 6.903424e-10 1.0000000 [42,] 1.583951e-10 3.167902e-10 1.0000000 [43,] 1.046328e-10 2.092655e-10 1.0000000 [44,] 4.589607e-11 9.179213e-11 1.0000000 [45,] 2.000633e-11 4.001267e-11 1.0000000 [46,] 9.019674e-12 1.803935e-11 1.0000000 [47,] 4.253908e-12 8.507817e-12 1.0000000 [48,] 2.162127e-12 4.324254e-12 1.0000000 [49,] 9.709263e-13 1.941853e-12 1.0000000 [50,] 4.053908e-13 8.107816e-13 1.0000000 [51,] 1.716025e-13 3.432049e-13 1.0000000 [52,] 8.134551e-14 1.626910e-13 1.0000000 [53,] 3.324879e-14 6.649758e-14 1.0000000 [54,] 1.348461e-14 2.696923e-14 1.0000000 [55,] 9.065586e-15 1.813117e-14 1.0000000 [56,] 3.821229e-15 7.642458e-15 1.0000000 [57,] 1.549156e-15 3.098312e-15 1.0000000 [58,] 6.562047e-16 1.312409e-15 1.0000000 [59,] 2.451638e-16 4.903276e-16 1.0000000 [60,] 1.197098e-16 2.394195e-16 1.0000000 [61,] 4.874148e-17 9.748295e-17 1.0000000 [62,] 2.054608e-17 4.109215e-17 1.0000000 [63,] 7.892697e-18 1.578539e-17 1.0000000 [64,] 4.283150e-18 8.566301e-18 1.0000000 [65,] 1.631164e-18 3.262328e-18 1.0000000 [66,] 6.865030e-19 1.373006e-18 1.0000000 [67,] 2.808978e-19 5.617956e-19 1.0000000 [68,] 1.110217e-19 2.220435e-19 1.0000000 [69,] 4.915851e-20 9.831702e-20 1.0000000 [70,] 1.940275e-20 3.880549e-20 1.0000000 [71,] 8.615550e-21 1.723110e-20 1.0000000 [72,] 3.329965e-21 6.659930e-21 1.0000000 [73,] 1.307677e-21 2.615353e-21 1.0000000 [74,] 4.759898e-22 9.519796e-22 1.0000000 [75,] 1.649766e-22 3.299531e-22 1.0000000 [76,] 6.549897e-23 1.309979e-22 1.0000000 [77,] 2.470359e-23 4.940717e-23 1.0000000 [78,] 9.032674e-24 1.806535e-23 1.0000000 [79,] 3.072709e-24 6.145417e-24 1.0000000 [80,] 1.063313e-24 2.126627e-24 1.0000000 [81,] 3.839309e-25 7.678617e-25 1.0000000 [82,] 1.391529e-25 2.783058e-25 1.0000000 [83,] 4.944555e-26 9.889109e-26 1.0000000 [84,] 1.688219e-26 3.376439e-26 1.0000000 [85,] 5.561028e-27 1.112206e-26 1.0000000 [86,] 1.816844e-27 3.633688e-27 1.0000000 [87,] 6.045227e-28 1.209045e-27 1.0000000 [88,] 1.940506e-28 3.881013e-28 1.0000000 [89,] 6.092513e-29 1.218503e-28 1.0000000 [90,] 2.245242e-29 4.490483e-29 1.0000000 [91,] 7.293385e-30 1.458677e-29 1.0000000 [92,] 2.642190e-30 5.284380e-30 1.0000000 [93,] 8.558370e-31 1.711674e-30 1.0000000 [94,] 2.805789e-31 5.611579e-31 1.0000000 [95,] 1.028237e-31 2.056474e-31 1.0000000 [96,] 3.398433e-32 6.796865e-32 1.0000000 [97,] 1.139655e-32 2.279310e-32 1.0000000 [98,] 3.453843e-33 6.907686e-33 1.0000000 [99,] 1.156283e-33 2.312565e-33 1.0000000 [100,] 6.609495e-34 1.321899e-33 1.0000000 [101,] 2.794897e-34 5.589794e-34 1.0000000 [102,] 8.499258e-35 1.699852e-34 1.0000000 [103,] 2.905298e-35 5.810596e-35 1.0000000 [104,] 8.901899e-36 1.780380e-35 1.0000000 [105,] 4.999616e-36 9.999232e-36 1.0000000 [106,] 1.575010e-36 3.150020e-36 1.0000000 [107,] 4.876073e-37 9.752145e-37 1.0000000 [108,] 1.393284e-37 2.786569e-37 1.0000000 [109,] 4.002322e-38 8.004643e-38 1.0000000 [110,] 1.191275e-38 2.382551e-38 1.0000000 [111,] 4.008940e-39 8.017881e-39 1.0000000 [112,] 1.508357e-39 3.016715e-39 1.0000000 [113,] 4.346884e-40 8.693768e-40 1.0000000 [114,] 1.210674e-40 2.421347e-40 1.0000000 [115,] 3.394383e-41 6.788767e-41 1.0000000 [116,] 1.024693e-41 2.049386e-41 1.0000000 [117,] 4.924799e-42 9.849597e-42 1.0000000 [118,] 1.335838e-42 2.671676e-42 1.0000000 [119,] 5.401411e-43 1.080282e-42 1.0000000 [120,] 1.496811e-43 2.993621e-43 1.0000000 [121,] 4.138401e-44 8.276802e-44 1.0000000 [122,] 1.135942e-44 2.271884e-44 1.0000000 [123,] 3.087154e-45 6.174308e-45 1.0000000 [124,] 8.365772e-46 1.673154e-45 1.0000000 [125,] 2.190263e-46 4.380527e-46 1.0000000 [126,] 5.843829e-47 1.168766e-46 1.0000000 [127,] 1.567106e-47 3.134211e-47 1.0000000 [128,] 4.141103e-48 8.282206e-48 1.0000000 [129,] 1.084649e-48 2.169299e-48 1.0000000 [130,] 2.828215e-49 5.656429e-49 1.0000000 [131,] 7.324404e-50 1.464881e-49 1.0000000 [132,] 1.884864e-50 3.769729e-50 1.0000000 [133,] 4.820098e-51 9.640196e-51 1.0000000 [134,] 1.217029e-51 2.434059e-51 1.0000000 [135,] 3.076106e-52 6.152212e-52 1.0000000 [136,] 8.021139e-53 1.604228e-52 1.0000000 [137,] 2.007301e-53 4.014603e-53 1.0000000 [138,] 4.979964e-54 9.959927e-54 1.0000000 [139,] 1.234700e-54 2.469401e-54 1.0000000 [140,] 3.033337e-55 6.066673e-55 1.0000000 [141,] 7.437712e-56 1.487542e-55 1.0000000 [142,] 1.795917e-56 3.591833e-56 1.0000000 [143,] 4.338172e-57 8.676345e-57 1.0000000 [144,] 1.085420e-57 2.170841e-57 1.0000000 [145,] 2.751911e-58 5.503821e-58 1.0000000 [146,] 6.558370e-59 1.311674e-58 1.0000000 [147,] 1.595341e-59 3.190682e-59 1.0000000 [148,] 3.755854e-60 7.511709e-60 1.0000000 [149,] 8.784980e-61 1.756996e-60 1.0000000 [150,] 2.043755e-61 4.087510e-61 1.0000000 [151,] 4.733893e-62 9.467786e-62 1.0000000 [152,] 1.082159e-62 2.164319e-62 1.0000000 [153,] 2.477840e-63 4.955681e-63 1.0000000 [154,] 5.643535e-64 1.128707e-63 1.0000000 [155,] 1.278599e-64 2.557198e-64 1.0000000 [156,] 2.881588e-65 5.763176e-65 1.0000000 [157,] 6.547938e-66 1.309588e-65 1.0000000 [158,] 1.506076e-66 3.012153e-66 1.0000000 [159,] 3.343327e-67 6.686653e-67 1.0000000 [160,] 7.379812e-68 1.475962e-67 1.0000000 [161,] 1.634993e-68 3.269986e-68 1.0000000 [162,] 4.199971e-69 8.399942e-69 1.0000000 [163,] 9.149807e-70 1.829961e-69 1.0000000 [164,] 1.983188e-70 3.966375e-70 1.0000000 [165,] 4.291285e-71 8.582570e-71 1.0000000 [166,] 9.207111e-72 1.841422e-71 1.0000000 [167,] 1.982459e-72 3.964918e-72 1.0000000 [168,] 4.210498e-73 8.420996e-73 1.0000000 [169,] 8.923368e-74 1.784674e-73 1.0000000 [170,] 2.064496e-74 4.128992e-74 1.0000000 [171,] 4.317101e-75 8.634201e-75 1.0000000 [172,] 8.992734e-76 1.798547e-75 1.0000000 [173,] 1.863978e-76 3.727956e-76 1.0000000 [174,] 3.923746e-77 7.847491e-77 1.0000000 [175,] 8.146689e-78 1.629338e-77 1.0000000 [176,] 1.962605e-78 3.925210e-78 1.0000000 [177,] 4.138857e-79 8.277715e-79 1.0000000 [178,] 8.424276e-80 1.684855e-79 1.0000000 [179,] 1.683556e-80 3.367113e-80 1.0000000 [180,] 3.436523e-81 6.873046e-81 1.0000000 [181,] 6.884869e-82 1.376974e-81 1.0000000 [182,] 1.415035e-82 2.830070e-82 1.0000000 [183,] 2.797711e-83 5.595422e-83 1.0000000 [184,] 5.493742e-84 1.098748e-83 1.0000000 [185,] 1.075712e-84 2.151424e-84 1.0000000 [186,] 2.101082e-85 4.202165e-85 1.0000000 [187,] 4.075116e-86 8.150232e-86 1.0000000 [188,] 7.875212e-87 1.575042e-86 1.0000000 [189,] 1.513430e-87 3.026860e-87 1.0000000 [190,] 3.504022e-88 7.008044e-88 1.0000000 [191,] 6.756597e-89 1.351319e-88 1.0000000 [192,] 1.301737e-89 2.603473e-89 1.0000000 [193,] 2.933585e-90 5.867171e-90 1.0000000 [194,] 5.619687e-91 1.123937e-90 1.0000000 [195,] 1.080354e-91 2.160709e-91 1.0000000 [196,] 2.096977e-92 4.193955e-92 1.0000000 [197,] 3.941754e-93 7.883507e-93 1.0000000 [198,] 7.454274e-94 1.490855e-93 1.0000000 [199,] 1.473384e-94 2.946768e-94 1.0000000 [200,] 3.152911e-95 6.305822e-95 1.0000000 [201,] 1.143386e-95 2.286772e-95 1.0000000 [202,] 1.572090e-94 3.144180e-94 1.0000000 [203,] 1.719520e-31 3.439040e-31 1.0000000 [204,] 3.179882e-31 6.359764e-31 1.0000000 [205,] 1.359689e-31 2.719378e-31 1.0000000 [206,] 6.022542e-32 1.204508e-31 1.0000000 [207,] 2.553133e-32 5.106265e-32 1.0000000 [208,] 1.101506e-32 2.203012e-32 1.0000000 [209,] 4.613961e-33 9.227921e-33 1.0000000 [210,] 1.933284e-33 3.866567e-33 1.0000000 [211,] 8.020659e-34 1.604132e-33 1.0000000 [212,] 4.720117e-34 9.440234e-34 1.0000000 [213,] 1.973781e-34 3.947562e-34 1.0000000 [214,] 8.337541e-35 1.667508e-34 1.0000000 [215,] 3.437695e-35 6.875391e-35 1.0000000 [216,] 1.399912e-35 2.799823e-35 1.0000000 [217,] 5.727277e-36 1.145455e-35 1.0000000 [218,] 2.300448e-36 4.600896e-36 1.0000000 [219,] 9.328884e-37 1.865777e-36 1.0000000 [220,] 3.745340e-37 7.490679e-37 1.0000000 [221,] 1.483174e-37 2.966349e-37 1.0000000 [222,] 6.147085e-38 1.229417e-37 1.0000000 [223,] 2.414678e-38 4.829357e-38 1.0000000 [224,] 9.593374e-39 1.918675e-38 1.0000000 [225,] 3.743117e-39 7.486234e-39 1.0000000 [226,] 1.452639e-39 2.905278e-39 1.0000000 [227,] 5.709466e-40 1.141893e-39 1.0000000 [228,] 2.199261e-40 4.398522e-40 1.0000000 [229,] 8.429253e-41 1.685851e-40 1.0000000 [230,] 3.206270e-41 6.412539e-41 1.0000000 [231,] 1.214409e-41 2.428818e-41 1.0000000 [232,] 4.584536e-42 9.169072e-42 1.0000000 [233,] 1.722015e-42 3.444031e-42 1.0000000 [234,] 6.414213e-43 1.282843e-42 1.0000000 [235,] 2.394097e-43 4.788194e-43 1.0000000 [236,] 8.849297e-44 1.769859e-43 1.0000000 [237,] 7.057370e-44 1.411474e-43 1.0000000 [238,] 3.093903e-44 6.187806e-44 1.0000000 [239,] 1.135660e-44 2.271320e-44 1.0000000 [240,] 4.130207e-45 8.260415e-45 1.0000000 [241,] 1.517328e-45 3.034656e-45 1.0000000 [242,] 5.486182e-46 1.097236e-45 1.0000000 [243,] 1.966769e-46 3.933537e-46 1.0000000 [244,] 7.167588e-47 1.433518e-46 1.0000000 [245,] 2.543980e-47 5.087961e-47 1.0000000 [246,] 9.127289e-48 1.825458e-47 1.0000000 [247,] 5.595379e-48 1.119076e-47 1.0000000 [248,] 4.996868e-48 9.993736e-48 1.0000000 [249,] 1.768589e-48 3.537179e-48 1.0000000 [250,] 6.498884e-49 1.299777e-48 1.0000000 [251,] 2.443329e-49 4.886658e-49 1.0000000 [252,] 8.486887e-50 1.697377e-49 1.0000000 [253,] 2.932990e-50 5.865980e-50 1.0000000 [254,] 1.008479e-50 2.016958e-50 1.0000000 [255,] 3.441039e-51 6.882079e-51 1.0000000 [256,] 1.330236e-51 2.660472e-51 1.0000000 [257,] 4.502382e-52 9.004763e-52 1.0000000 [258,] 1.517786e-52 3.035572e-52 1.0000000 [259,] 5.090485e-53 1.018097e-52 1.0000000 [260,] 1.694250e-53 3.388499e-53 1.0000000 [261,] 7.244327e-54 1.448865e-53 1.0000000 [262,] 2.388934e-54 4.777867e-54 1.0000000 [263,] 8.901083e-55 1.780217e-54 1.0000000 [264,] 2.912523e-55 5.825045e-55 1.0000000 [265,] 9.481096e-56 1.896219e-55 1.0000000 [266,] 3.065366e-56 6.130732e-56 1.0000000 [267,] 9.926118e-57 1.985224e-56 1.0000000 [268,] 3.195349e-57 6.390697e-57 1.0000000 [269,] 1.036771e-57 2.073542e-57 1.0000000 [270,] 3.404980e-58 6.809960e-58 1.0000000 [271,] 1.408534e-58 2.817067e-58 1.0000000 [272,] 4.872769e-59 9.745538e-59 1.0000000 [273,] 1.667587e-59 3.335174e-59 1.0000000 [274,] 5.193284e-60 1.038657e-59 1.0000000 [275,] 1.613981e-60 3.227962e-60 1.0000000 [276,] 4.991982e-61 9.983964e-61 1.0000000 [277,] 1.544566e-61 3.089132e-61 1.0000000 [278,] 5.353445e-62 1.070689e-61 1.0000000 [279,] 1.705631e-62 3.411263e-62 1.0000000 [280,] 5.149150e-63 1.029830e-62 1.0000000 [281,] 1.554649e-63 3.109299e-63 1.0000000 [282,] 4.661483e-64 9.322965e-64 1.0000000 [283,] 1.400195e-64 2.800391e-64 1.0000000 [284,] 4.136964e-65 8.273927e-65 1.0000000 [285,] 1.267296e-65 2.534592e-65 1.0000000 [286,] 5.972430e-66 1.194486e-65 1.0000000 [287,] 1.738281e-66 3.476561e-66 1.0000000 [288,] 5.118820e-67 1.023764e-66 1.0000000 [289,] 1.487909e-67 2.975818e-67 1.0000000 [290,] 4.264271e-68 8.528542e-68 1.0000000 [291,] 1.215266e-68 2.430531e-68 1.0000000 [292,] 3.443857e-69 6.887715e-69 1.0000000 [293,] 9.720773e-70 1.944155e-69 1.0000000 [294,] 2.748689e-70 5.497378e-70 1.0000000 [295,] 7.661794e-71 1.532359e-70 1.0000000 [296,] 2.122818e-71 4.245636e-71 1.0000000 [297,] 5.946448e-72 1.189290e-71 1.0000000 [298,] 1.658230e-72 3.316459e-72 1.0000000 [299,] 4.532656e-73 9.065311e-73 1.0000000 [300,] 1.231090e-73 2.462179e-73 1.0000000 [301,] 3.467139e-74 6.934278e-74 1.0000000 [302,] 9.405537e-75 1.881107e-74 1.0000000 [303,] 2.503292e-75 5.006585e-75 1.0000000 [304,] 6.737848e-76 1.347570e-75 1.0000000 [305,] 4.635612e-76 9.271223e-76 1.0000000 [306,] 1.215673e-76 2.431347e-76 1.0000000 [307,] 3.279483e-77 6.558966e-77 1.0000000 [308,] 8.548824e-78 1.709765e-77 1.0000000 [309,] 2.303353e-78 4.606706e-78 1.0000000 [310,] 6.052766e-79 1.210553e-78 1.0000000 [311,] 1.563370e-79 3.126740e-79 1.0000000 [312,] 3.962455e-80 7.924909e-80 1.0000000 [313,] 9.964648e-81 1.992930e-80 1.0000000 [314,] 2.490101e-81 4.980202e-81 1.0000000 [315,] 6.183183e-82 1.236637e-81 1.0000000 [316,] 1.525562e-82 3.051124e-82 1.0000000 [317,] 3.739836e-83 7.479671e-83 1.0000000 [318,] 9.106880e-84 1.821376e-83 1.0000000 [319,] 2.204781e-84 4.409562e-84 1.0000000 [320,] 5.300201e-85 1.060040e-84 1.0000000 [321,] 1.265751e-85 2.531503e-85 1.0000000 [322,] 3.002711e-86 6.005422e-86 1.0000000 [323,] 7.075667e-87 1.415133e-86 1.0000000 [324,] 1.656109e-87 3.312218e-87 1.0000000 [325,] 3.849965e-88 7.699931e-88 1.0000000 [326,] 8.897520e-89 1.779504e-88 1.0000000 [327,] 2.040650e-89 4.081300e-89 1.0000000 [328,] 4.739729e-90 9.479459e-90 1.0000000 [329,] 1.257923e-90 2.515846e-90 1.0000000 [330,] 2.867742e-91 5.735484e-91 1.0000000 [331,] 6.400568e-92 1.280114e-91 1.0000000 [332,] 1.422070e-92 2.844140e-92 1.0000000 [333,] 3.234830e-93 6.469660e-93 1.0000000 [334,] 7.302499e-94 1.460500e-93 1.0000000 [335,] 1.584240e-94 3.168479e-94 1.0000000 [336,] 3.411563e-95 6.823127e-95 1.0000000 [337,] 7.942871e-96 1.588574e-95 1.0000000 [338,] 1.698175e-96 3.396351e-96 1.0000000 [339,] 3.586179e-97 7.172358e-97 1.0000000 [340,] 7.468370e-98 1.493674e-97 1.0000000 [341,] 1.549369e-98 3.098739e-98 1.0000000 [342,] 5.799935e-99 1.159987e-98 1.0000000 [343,] 1.233566e-99 2.467132e-99 1.0000000 [344,] 2.511657e-100 5.023314e-100 1.0000000 [345,] 5.508980e-101 1.101796e-100 1.0000000 [346,] 1.125702e-101 2.251404e-101 1.0000000 [347,] 2.260961e-102 4.521922e-102 1.0000000 [348,] 4.456139e-103 8.912278e-103 1.0000000 [349,] 9.264020e-104 1.852804e-103 1.0000000 [350,] 3.225471e-104 6.450942e-104 1.0000000 [351,] 6.220965e-105 1.244193e-104 1.0000000 [352,] 1.189691e-105 2.379383e-105 1.0000000 [353,] 2.255730e-106 4.511461e-106 1.0000000 [354,] 4.253237e-107 8.506474e-107 1.0000000 [355,] 7.925552e-108 1.585110e-107 1.0000000 [356,] 3.224558e-108 6.449116e-108 1.0000000 [357,] 6.221687e-109 1.244337e-108 1.0000000 [358,] 1.136705e-109 2.273409e-109 1.0000000 [359,] 2.054750e-110 4.109500e-110 1.0000000 [360,] 3.680542e-111 7.361085e-111 1.0000000 [361,] 6.532314e-112 1.306463e-111 1.0000000 [362,] 1.166145e-112 2.332291e-112 1.0000000 [363,] 2.047116e-113 4.094233e-113 1.0000000 [364,] 3.749704e-114 7.499408e-114 1.0000000 [365,] 6.666230e-115 1.333246e-114 1.0000000 [366,] 1.813893e-115 3.627787e-115 1.0000000 [367,] 3.268060e-116 6.536119e-116 1.0000000 [368,] 5.469190e-117 1.093838e-116 1.0000000 [369,] 9.139755e-118 1.827951e-117 1.0000000 [370,] 1.678753e-118 3.357506e-118 1.0000000 [371,] 1.189169e-118 2.378339e-118 1.0000000 [372,] 1.974428e-119 3.948856e-119 1.0000000 [373,] 3.189461e-120 6.378921e-120 1.0000000 [374,] 5.640105e-121 1.128021e-120 1.0000000 [375,] 1.364488e-121 2.728977e-121 1.0000000 [376,] 2.611075e-122 5.222151e-122 1.0000000 [377,] 4.105721e-123 8.211441e-123 1.0000000 [378,] 3.859383e-123 7.718767e-123 1.0000000 [379,] 7.379655e-124 1.475931e-123 1.0000000 [380,] 1.655556e-124 3.311112e-124 1.0000000 [381,] 2.507606e-125 5.015212e-125 1.0000000 [382,] 4.419909e-126 8.839819e-126 1.0000000 [383,] 6.785805e-127 1.357161e-126 1.0000000 [384,] 1.138074e-127 2.276149e-127 1.0000000 [385,] 1.762340e-128 3.524680e-128 1.0000000 [386,] 2.444690e-129 4.889380e-129 1.0000000 [387,] 1.626040e-128 3.252080e-128 1.0000000 [388,] 5.083736e-129 1.016747e-128 1.0000000 [389,] 3.453745e-127 6.907490e-127 1.0000000 [390,] 6.878716e-128 1.375743e-127 1.0000000 [391,] 4.061757e-123 8.123514e-123 1.0000000 [392,] 1.649220e-123 3.298439e-123 1.0000000 [393,] 1.298665e-123 2.597331e-123 1.0000000 [394,] 2.095980e-124 4.191960e-124 1.0000000 [395,] 3.717499e-125 7.434998e-125 1.0000000 [396,] 8.612823e-126 1.722565e-125 1.0000000 [397,] 1.454992e-126 2.909985e-126 1.0000000 [398,] 2.328971e-127 4.657943e-127 1.0000000 [399,] 4.811506e-128 9.623012e-128 1.0000000 [400,] 1.057440e-128 2.114881e-128 1.0000000 [401,] 1.462247e-129 2.924493e-129 1.0000000 [402,] 2.153549e-130 4.307097e-130 1.0000000 [403,] 2.733865e-131 5.467730e-131 1.0000000 [404,] 1.267842e-124 2.535683e-124 1.0000000 [405,] 4.409326e-125 8.818653e-125 1.0000000 [406,] 1.514407e-124 3.028815e-124 1.0000000 [407,] 7.863393e-121 1.572679e-120 1.0000000 [408,] 1.156639e-121 2.313278e-121 1.0000000 [409,] 8.431701e-122 1.686340e-121 1.0000000 [410,] 5.847202e-118 1.169440e-117 1.0000000 [411,] 2.992952e-116 5.985904e-116 1.0000000 [412,] 9.868378e-115 1.973676e-114 1.0000000 [413,] 6.394888e-114 1.278978e-113 1.0000000 [414,] 6.690517e-110 1.338103e-109 1.0000000 [415,] 1.131871e-84 2.263741e-84 1.0000000 [416,] 5.888222e-51 1.177644e-50 1.0000000 > postscript(file="/var/www/html/rcomp/tmp/1menx1291297725.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index') > points(x[,1]-mysum$resid) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/2menx1291297725.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/3f5401291297725.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/4f5401291297725.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals') > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5f5401291297725.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 18112.030820 3189.149011 -5975.240229 -5896.762010 -4963.274386 6 7 8 9 10 -3586.218801 -3639.496594 -1771.217901 -5120.670902 -7276.418018 11 12 13 14 15 3969.352621 -4492.098869 -10688.605698 -3692.210291 -3658.988881 16 17 18 19 20 -571.565104 8580.303693 -2205.676265 -7792.585933 749.071535 21 22 23 24 25 -3340.109525 -3344.907924 1462.809853 760.744286 2502.076015 26 27 28 29 30 -3164.005247 -3829.109271 -2410.225134 -5503.404364 -3004.098847 31 32 33 34 35 -1916.099554 -860.802328 -2790.326038 -779.212044 8916.934765 36 37 38 39 40 -407.005211 -2080.687204 -935.992092 3988.236079 346.561017 41 42 43 44 45 1884.552701 -220.002535 -705.251195 1014.684702 -1268.559636 46 47 48 49 50 446.715966 5762.193305 1607.948969 818.275166 329.161773 51 52 53 54 55 -121.639573 1192.597994 3825.734858 1898.113007 3883.404235 56 57 58 59 60 -1663.508514 -3703.968891 -777.589508 656.930989 -707.185995 61 62 63 64 65 -745.130968 -164.659171 -120.864169 -715.556170 -1514.636577 66 67 68 69 70 -650.065442 217.786117 318.539775 1989.611071 19.621278 71 72 73 74 75 -2214.551765 -2783.762310 1173.976909 -213.043438 1246.830179 76 77 78 79 80 82.895705 1252.997086 22.872089 -676.919991 1592.900422 81 82 83 84 85 -1299.502738 -530.426761 2538.618480 -2052.259268 635.543149 86 87 88 89 90 -1470.195210 -656.375725 829.144472 1120.037178 609.702605 91 92 93 94 95 91.538895 -504.930538 -1849.356020 1524.364018 -2564.966168 96 97 98 99 100 -1618.168056 849.475477 -761.256688 1699.459918 -533.611160 101 102 103 104 105 -328.546992 821.087835 1179.268373 1533.402505 -1072.452550 106 107 108 109 110 -304.685718 -199.410741 -481.158937 561.232612 -191.644522 111 112 113 114 115 -1300.599370 -985.347261 10.897599 -399.340499 -1661.614585 116 117 118 119 120 -1532.042011 11.808421 -562.540244 202.842125 -1088.325075 121 122 123 124 125 375.001446 -914.745069 179.341066 185.809164 -1335.232702 126 127 128 129 130 47.161008 -1061.735897 -1061.158843 -1245.133403 -1061.735897 131 132 133 134 135 -1063.834583 -1725.494028 -1061.735897 -1045.872342 -1082.659339 136 137 138 139 140 -1090.258677 -1080.264849 -1059.405046 -1114.409716 -1061.735897 141 142 143 144 145 -1231.006356 -1061.735897 -1223.175503 -1058.869179 -1047.115252 146 147 148 149 150 -1274.122439 -1061.735897 -1055.892234 -1135.815656 -1050.887300 151 152 153 154 155 -1059.353040 -1138.519584 -1062.843185 -1092.243375 -1051.874989 156 157 158 159 160 -1061.735897 -1061.735897 -1029.374768 -1203.209376 -1061.735897 161 162 163 164 165 -1061.735897 -1061.735897 -1061.735897 -1034.423174 -893.054439 166 167 168 169 170 -1061.927612 -1061.735897 -2242.074409 -1508.435528 -1061.735897 171 172 173 174 175 -1061.735897 -1005.825011 -1061.735897 -982.565463 -1061.735897 176 177 178 179 180 -1060.752605 -1098.819809 -1069.429333 -1061.735897 -1061.735897 181 182 183 184 185 -841.713596 -643.649635 -1453.169244 -434.857855 -970.539261 186 187 188 189 190 -1223.462113 -1214.789853 -995.191188 -692.406099 -1061.735897 191 192 193 194 195 -1159.573192 -1061.735897 -1026.952296 -1061.735897 -1057.512909 196 197 198 199 200 -1061.735897 670.444619 -1061.735897 -1061.735897 122.450100 201 202 203 204 205 -1061.735897 -1061.735897 -1061.735897 -1481.573495 160.496827 206 207 208 209 210 -1021.841361 -839.829927 -1061.735897 -1900.832563 -1029.650703 211 212 213 214 215 -139.270215 127.729785 1739.108604 -1184.234238 -1096.934730 216 217 218 219 220 -124.391717 247.952847 -7.271879 1726.608283 1484.825340 221 222 223 224 225 -1015.125087 -136.270215 332.260570 -279.024713 129.498368 226 227 228 229 230 -253.640739 -180.646183 52.729785 -1298.125172 177.729785 231 232 233 234 235 -930.761739 155.492980 228.729785 -952.389572 -137.270215 236 237 238 239 240 -139.093865 -86.270215 -37.133233 -139.270215 -139.270215 241 242 243 244 245 -45.270215 -403.968451 -115.270215 2404.716804 -2344.517538 246 247 248 249 250 -139.270215 -8.270215 -695.891470 -139.270215 66.729785 251 252 253 254 255 625.561964 27.729785 -694.834735 -1076.587638 -8060.784634 256 257 258 259 260 225.729785 1244.262432 1908.386400 -139.270215 -139.270215 261 262 263 264 265 -139.270215 86.729785 -1109.947060 1089.257151 -139.270215 266 267 268 269 270 -139.270215 48.729785 -649.474595 -1.270215 -1365.460633 271 272 273 274 275 -139.270215 -139.270215 -14.270215 -412.769924 142.729785 276 277 278 279 280 39.764643 -658.661062 1184.729785 -2612.549884 -2014.058505 281 282 283 284 285 -139.270215 109.729785 650.616351 193.729785 -1589.334055 286 287 288 289 290 461.729785 79.057971 -299.019094 109.729785 162.247354 291 292 293 294 295 25.729785 325.753632 -2168.191760 -86.270215 -127.865024 296 297 298 299 300 1420.496699 -139.270215 -139.270215 -139.270215 -109.816972 301 302 303 304 305 24.072522 -412.072342 -139.330966 226.729785 1051.433365 306 307 308 309 310 -142.631572 69.729785 -423.947404 923.587931 -139.270215 311 312 313 314 315 225.729785 -4465.476016 -90.270215 -559.628047 -2.242662 316 317 318 319 320 -89.365450 370.348242 -1136.768660 -139.746521 -139.270215 321 322 323 324 325 -139.270215 -139.270215 -139.270215 -139.270215 -117.270215 326 327 328 329 330 -281.518718 -139.270215 -139.270215 -139.270215 -139.270215 331 332 333 334 335 -139.270215 -139.270215 -43.270215 -144.309852 173.512840 336 337 338 339 340 704.043521 -1133.525590 -113.270215 -127.184237 496.643627 341 342 343 344 345 -1613.230251 -139.270215 -139.270215 481.729785 -490.925235 346 347 348 349 350 -20.270215 -361.063762 -139.270215 1455.729785 368.585152 351 352 353 354 355 553.501829 292.282379 752.484755 -494.815251 -139.270215 356 357 358 359 360 374.729785 -3795.045538 -139.270215 -139.270215 -138.270215 361 362 363 364 365 -139.069613 -139.270215 1623.729785 -502.779063 150.011155 366 367 368 369 370 -139.270215 -139.270215 -139.270215 78.729785 -131.175020 371 372 373 374 375 308.729785 688.326689 -2445.155383 -1439.270720 -139.270215 376 377 378 379 380 -18.270215 487.406100 2477.986548 80.964476 -139.270215 381 382 383 384 385 390.729785 -2290.634823 2398.483194 -61.270215 2006.554552 386 387 388 389 390 -5383.735718 1593.866401 85.075320 446.233013 364.247630 391 392 393 394 395 211.575315 -353.489569 -2062.862285 4234.482677 -6199.111883 396 397 398 399 400 3712.798981 -480.110115 -6949.088627 347.097654 5540.850501 401 402 403 404 405 -2058.527317 586.628153 815.992878 -3252.172767 -1639.129780 406 407 408 409 410 -1155.084059 -682.175458 161.294988 266.453118 -5000.528162 411 412 413 414 415 4958.071007 -937.482146 1675.626363 5109.128002 -10702.652116 416 417 418 419 420 -5180.946020 11155.125275 -443.565059 -3332.403928 -3668.241362 421 422 423 424 425 -2348.678631 5066.650137 3417.018967 -14401.817777 4095.498813 426 427 428 429 430 -2470.243562 24153.403695 -12733.619516 7888.505936 17046.328143 431 153933.352077 > postscript(file="/var/www/html/rcomp/tmp/6pell1291297725.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 18112.030820 NA 1 3189.149011 18112.030820 2 -5975.240229 3189.149011 3 -5896.762010 -5975.240229 4 -4963.274386 -5896.762010 5 -3586.218801 -4963.274386 6 -3639.496594 -3586.218801 7 -1771.217901 -3639.496594 8 -5120.670902 -1771.217901 9 -7276.418018 -5120.670902 10 3969.352621 -7276.418018 11 -4492.098869 3969.352621 12 -10688.605698 -4492.098869 13 -3692.210291 -10688.605698 14 -3658.988881 -3692.210291 15 -571.565104 -3658.988881 16 8580.303693 -571.565104 17 -2205.676265 8580.303693 18 -7792.585933 -2205.676265 19 749.071535 -7792.585933 20 -3340.109525 749.071535 21 -3344.907924 -3340.109525 22 1462.809853 -3344.907924 23 760.744286 1462.809853 24 2502.076015 760.744286 25 -3164.005247 2502.076015 26 -3829.109271 -3164.005247 27 -2410.225134 -3829.109271 28 -5503.404364 -2410.225134 29 -3004.098847 -5503.404364 30 -1916.099554 -3004.098847 31 -860.802328 -1916.099554 32 -2790.326038 -860.802328 33 -779.212044 -2790.326038 34 8916.934765 -779.212044 35 -407.005211 8916.934765 36 -2080.687204 -407.005211 37 -935.992092 -2080.687204 38 3988.236079 -935.992092 39 346.561017 3988.236079 40 1884.552701 346.561017 41 -220.002535 1884.552701 42 -705.251195 -220.002535 43 1014.684702 -705.251195 44 -1268.559636 1014.684702 45 446.715966 -1268.559636 46 5762.193305 446.715966 47 1607.948969 5762.193305 48 818.275166 1607.948969 49 329.161773 818.275166 50 -121.639573 329.161773 51 1192.597994 -121.639573 52 3825.734858 1192.597994 53 1898.113007 3825.734858 54 3883.404235 1898.113007 55 -1663.508514 3883.404235 56 -3703.968891 -1663.508514 57 -777.589508 -3703.968891 58 656.930989 -777.589508 59 -707.185995 656.930989 60 -745.130968 -707.185995 61 -164.659171 -745.130968 62 -120.864169 -164.659171 63 -715.556170 -120.864169 64 -1514.636577 -715.556170 65 -650.065442 -1514.636577 66 217.786117 -650.065442 67 318.539775 217.786117 68 1989.611071 318.539775 69 19.621278 1989.611071 70 -2214.551765 19.621278 71 -2783.762310 -2214.551765 72 1173.976909 -2783.762310 73 -213.043438 1173.976909 74 1246.830179 -213.043438 75 82.895705 1246.830179 76 1252.997086 82.895705 77 22.872089 1252.997086 78 -676.919991 22.872089 79 1592.900422 -676.919991 80 -1299.502738 1592.900422 81 -530.426761 -1299.502738 82 2538.618480 -530.426761 83 -2052.259268 2538.618480 84 635.543149 -2052.259268 85 -1470.195210 635.543149 86 -656.375725 -1470.195210 87 829.144472 -656.375725 88 1120.037178 829.144472 89 609.702605 1120.037178 90 91.538895 609.702605 91 -504.930538 91.538895 92 -1849.356020 -504.930538 93 1524.364018 -1849.356020 94 -2564.966168 1524.364018 95 -1618.168056 -2564.966168 96 849.475477 -1618.168056 97 -761.256688 849.475477 98 1699.459918 -761.256688 99 -533.611160 1699.459918 100 -328.546992 -533.611160 101 821.087835 -328.546992 102 1179.268373 821.087835 103 1533.402505 1179.268373 104 -1072.452550 1533.402505 105 -304.685718 -1072.452550 106 -199.410741 -304.685718 107 -481.158937 -199.410741 108 561.232612 -481.158937 109 -191.644522 561.232612 110 -1300.599370 -191.644522 111 -985.347261 -1300.599370 112 10.897599 -985.347261 113 -399.340499 10.897599 114 -1661.614585 -399.340499 115 -1532.042011 -1661.614585 116 11.808421 -1532.042011 117 -562.540244 11.808421 118 202.842125 -562.540244 119 -1088.325075 202.842125 120 375.001446 -1088.325075 121 -914.745069 375.001446 122 179.341066 -914.745069 123 185.809164 179.341066 124 -1335.232702 185.809164 125 47.161008 -1335.232702 126 -1061.735897 47.161008 127 -1061.158843 -1061.735897 128 -1245.133403 -1061.158843 129 -1061.735897 -1245.133403 130 -1063.834583 -1061.735897 131 -1725.494028 -1063.834583 132 -1061.735897 -1725.494028 133 -1045.872342 -1061.735897 134 -1082.659339 -1045.872342 135 -1090.258677 -1082.659339 136 -1080.264849 -1090.258677 137 -1059.405046 -1080.264849 138 -1114.409716 -1059.405046 139 -1061.735897 -1114.409716 140 -1231.006356 -1061.735897 141 -1061.735897 -1231.006356 142 -1223.175503 -1061.735897 143 -1058.869179 -1223.175503 144 -1047.115252 -1058.869179 145 -1274.122439 -1047.115252 146 -1061.735897 -1274.122439 147 -1055.892234 -1061.735897 148 -1135.815656 -1055.892234 149 -1050.887300 -1135.815656 150 -1059.353040 -1050.887300 151 -1138.519584 -1059.353040 152 -1062.843185 -1138.519584 153 -1092.243375 -1062.843185 154 -1051.874989 -1092.243375 155 -1061.735897 -1051.874989 156 -1061.735897 -1061.735897 157 -1029.374768 -1061.735897 158 -1203.209376 -1029.374768 159 -1061.735897 -1203.209376 160 -1061.735897 -1061.735897 161 -1061.735897 -1061.735897 162 -1061.735897 -1061.735897 163 -1034.423174 -1061.735897 164 -893.054439 -1034.423174 165 -1061.927612 -893.054439 166 -1061.735897 -1061.927612 167 -2242.074409 -1061.735897 168 -1508.435528 -2242.074409 169 -1061.735897 -1508.435528 170 -1061.735897 -1061.735897 171 -1005.825011 -1061.735897 172 -1061.735897 -1005.825011 173 -982.565463 -1061.735897 174 -1061.735897 -982.565463 175 -1060.752605 -1061.735897 176 -1098.819809 -1060.752605 177 -1069.429333 -1098.819809 178 -1061.735897 -1069.429333 179 -1061.735897 -1061.735897 180 -841.713596 -1061.735897 181 -643.649635 -841.713596 182 -1453.169244 -643.649635 183 -434.857855 -1453.169244 184 -970.539261 -434.857855 185 -1223.462113 -970.539261 186 -1214.789853 -1223.462113 187 -995.191188 -1214.789853 188 -692.406099 -995.191188 189 -1061.735897 -692.406099 190 -1159.573192 -1061.735897 191 -1061.735897 -1159.573192 192 -1026.952296 -1061.735897 193 -1061.735897 -1026.952296 194 -1057.512909 -1061.735897 195 -1061.735897 -1057.512909 196 670.444619 -1061.735897 197 -1061.735897 670.444619 198 -1061.735897 -1061.735897 199 122.450100 -1061.735897 200 -1061.735897 122.450100 201 -1061.735897 -1061.735897 202 -1061.735897 -1061.735897 203 -1481.573495 -1061.735897 204 160.496827 -1481.573495 205 -1021.841361 160.496827 206 -839.829927 -1021.841361 207 -1061.735897 -839.829927 208 -1900.832563 -1061.735897 209 -1029.650703 -1900.832563 210 -139.270215 -1029.650703 211 127.729785 -139.270215 212 1739.108604 127.729785 213 -1184.234238 1739.108604 214 -1096.934730 -1184.234238 215 -124.391717 -1096.934730 216 247.952847 -124.391717 217 -7.271879 247.952847 218 1726.608283 -7.271879 219 1484.825340 1726.608283 220 -1015.125087 1484.825340 221 -136.270215 -1015.125087 222 332.260570 -136.270215 223 -279.024713 332.260570 224 129.498368 -279.024713 225 -253.640739 129.498368 226 -180.646183 -253.640739 227 52.729785 -180.646183 228 -1298.125172 52.729785 229 177.729785 -1298.125172 230 -930.761739 177.729785 231 155.492980 -930.761739 232 228.729785 155.492980 233 -952.389572 228.729785 234 -137.270215 -952.389572 235 -139.093865 -137.270215 236 -86.270215 -139.093865 237 -37.133233 -86.270215 238 -139.270215 -37.133233 239 -139.270215 -139.270215 240 -45.270215 -139.270215 241 -403.968451 -45.270215 242 -115.270215 -403.968451 243 2404.716804 -115.270215 244 -2344.517538 2404.716804 245 -139.270215 -2344.517538 246 -8.270215 -139.270215 247 -695.891470 -8.270215 248 -139.270215 -695.891470 249 66.729785 -139.270215 250 625.561964 66.729785 251 27.729785 625.561964 252 -694.834735 27.729785 253 -1076.587638 -694.834735 254 -8060.784634 -1076.587638 255 225.729785 -8060.784634 256 1244.262432 225.729785 257 1908.386400 1244.262432 258 -139.270215 1908.386400 259 -139.270215 -139.270215 260 -139.270215 -139.270215 261 86.729785 -139.270215 262 -1109.947060 86.729785 263 1089.257151 -1109.947060 264 -139.270215 1089.257151 265 -139.270215 -139.270215 266 48.729785 -139.270215 267 -649.474595 48.729785 268 -1.270215 -649.474595 269 -1365.460633 -1.270215 270 -139.270215 -1365.460633 271 -139.270215 -139.270215 272 -14.270215 -139.270215 273 -412.769924 -14.270215 274 142.729785 -412.769924 275 39.764643 142.729785 276 -658.661062 39.764643 277 1184.729785 -658.661062 278 -2612.549884 1184.729785 279 -2014.058505 -2612.549884 280 -139.270215 -2014.058505 281 109.729785 -139.270215 282 650.616351 109.729785 283 193.729785 650.616351 284 -1589.334055 193.729785 285 461.729785 -1589.334055 286 79.057971 461.729785 287 -299.019094 79.057971 288 109.729785 -299.019094 289 162.247354 109.729785 290 25.729785 162.247354 291 325.753632 25.729785 292 -2168.191760 325.753632 293 -86.270215 -2168.191760 294 -127.865024 -86.270215 295 1420.496699 -127.865024 296 -139.270215 1420.496699 297 -139.270215 -139.270215 298 -139.270215 -139.270215 299 -109.816972 -139.270215 300 24.072522 -109.816972 301 -412.072342 24.072522 302 -139.330966 -412.072342 303 226.729785 -139.330966 304 1051.433365 226.729785 305 -142.631572 1051.433365 306 69.729785 -142.631572 307 -423.947404 69.729785 308 923.587931 -423.947404 309 -139.270215 923.587931 310 225.729785 -139.270215 311 -4465.476016 225.729785 312 -90.270215 -4465.476016 313 -559.628047 -90.270215 314 -2.242662 -559.628047 315 -89.365450 -2.242662 316 370.348242 -89.365450 317 -1136.768660 370.348242 318 -139.746521 -1136.768660 319 -139.270215 -139.746521 320 -139.270215 -139.270215 321 -139.270215 -139.270215 322 -139.270215 -139.270215 323 -139.270215 -139.270215 324 -117.270215 -139.270215 325 -281.518718 -117.270215 326 -139.270215 -281.518718 327 -139.270215 -139.270215 328 -139.270215 -139.270215 329 -139.270215 -139.270215 330 -139.270215 -139.270215 331 -139.270215 -139.270215 332 -43.270215 -139.270215 333 -144.309852 -43.270215 334 173.512840 -144.309852 335 704.043521 173.512840 336 -1133.525590 704.043521 337 -113.270215 -1133.525590 338 -127.184237 -113.270215 339 496.643627 -127.184237 340 -1613.230251 496.643627 341 -139.270215 -1613.230251 342 -139.270215 -139.270215 343 481.729785 -139.270215 344 -490.925235 481.729785 345 -20.270215 -490.925235 346 -361.063762 -20.270215 347 -139.270215 -361.063762 348 1455.729785 -139.270215 349 368.585152 1455.729785 350 553.501829 368.585152 351 292.282379 553.501829 352 752.484755 292.282379 353 -494.815251 752.484755 354 -139.270215 -494.815251 355 374.729785 -139.270215 356 -3795.045538 374.729785 357 -139.270215 -3795.045538 358 -139.270215 -139.270215 359 -138.270215 -139.270215 360 -139.069613 -138.270215 361 -139.270215 -139.069613 362 1623.729785 -139.270215 363 -502.779063 1623.729785 364 150.011155 -502.779063 365 -139.270215 150.011155 366 -139.270215 -139.270215 367 -139.270215 -139.270215 368 78.729785 -139.270215 369 -131.175020 78.729785 370 308.729785 -131.175020 371 688.326689 308.729785 372 -2445.155383 688.326689 373 -1439.270720 -2445.155383 374 -139.270215 -1439.270720 375 -18.270215 -139.270215 376 487.406100 -18.270215 377 2477.986548 487.406100 378 80.964476 2477.986548 379 -139.270215 80.964476 380 390.729785 -139.270215 381 -2290.634823 390.729785 382 2398.483194 -2290.634823 383 -61.270215 2398.483194 384 2006.554552 -61.270215 385 -5383.735718 2006.554552 386 1593.866401 -5383.735718 387 85.075320 1593.866401 388 446.233013 85.075320 389 364.247630 446.233013 390 211.575315 364.247630 391 -353.489569 211.575315 392 -2062.862285 -353.489569 393 4234.482677 -2062.862285 394 -6199.111883 4234.482677 395 3712.798981 -6199.111883 396 -480.110115 3712.798981 397 -6949.088627 -480.110115 398 347.097654 -6949.088627 399 5540.850501 347.097654 400 -2058.527317 5540.850501 401 586.628153 -2058.527317 402 815.992878 586.628153 403 -3252.172767 815.992878 404 -1639.129780 -3252.172767 405 -1155.084059 -1639.129780 406 -682.175458 -1155.084059 407 161.294988 -682.175458 408 266.453118 161.294988 409 -5000.528162 266.453118 410 4958.071007 -5000.528162 411 -937.482146 4958.071007 412 1675.626363 -937.482146 413 5109.128002 1675.626363 414 -10702.652116 5109.128002 415 -5180.946020 -10702.652116 416 11155.125275 -5180.946020 417 -443.565059 11155.125275 418 -3332.403928 -443.565059 419 -3668.241362 -3332.403928 420 -2348.678631 -3668.241362 421 5066.650137 -2348.678631 422 3417.018967 5066.650137 423 -14401.817777 3417.018967 424 4095.498813 -14401.817777 425 -2470.243562 4095.498813 426 24153.403695 -2470.243562 427 -12733.619516 24153.403695 428 7888.505936 -12733.619516 429 17046.328143 7888.505936 430 153933.352077 17046.328143 431 NA 153933.352077 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 3189.149011 18112.030820 [2,] -5975.240229 3189.149011 [3,] -5896.762010 -5975.240229 [4,] -4963.274386 -5896.762010 [5,] -3586.218801 -4963.274386 [6,] -3639.496594 -3586.218801 [7,] -1771.217901 -3639.496594 [8,] -5120.670902 -1771.217901 [9,] -7276.418018 -5120.670902 [10,] 3969.352621 -7276.418018 [11,] -4492.098869 3969.352621 [12,] -10688.605698 -4492.098869 [13,] -3692.210291 -10688.605698 [14,] -3658.988881 -3692.210291 [15,] -571.565104 -3658.988881 [16,] 8580.303693 -571.565104 [17,] -2205.676265 8580.303693 [18,] -7792.585933 -2205.676265 [19,] 749.071535 -7792.585933 [20,] -3340.109525 749.071535 [21,] -3344.907924 -3340.109525 [22,] 1462.809853 -3344.907924 [23,] 760.744286 1462.809853 [24,] 2502.076015 760.744286 [25,] -3164.005247 2502.076015 [26,] -3829.109271 -3164.005247 [27,] -2410.225134 -3829.109271 [28,] -5503.404364 -2410.225134 [29,] -3004.098847 -5503.404364 [30,] -1916.099554 -3004.098847 [31,] -860.802328 -1916.099554 [32,] -2790.326038 -860.802328 [33,] -779.212044 -2790.326038 [34,] 8916.934765 -779.212044 [35,] -407.005211 8916.934765 [36,] -2080.687204 -407.005211 [37,] -935.992092 -2080.687204 [38,] 3988.236079 -935.992092 [39,] 346.561017 3988.236079 [40,] 1884.552701 346.561017 [41,] -220.002535 1884.552701 [42,] -705.251195 -220.002535 [43,] 1014.684702 -705.251195 [44,] -1268.559636 1014.684702 [45,] 446.715966 -1268.559636 [46,] 5762.193305 446.715966 [47,] 1607.948969 5762.193305 [48,] 818.275166 1607.948969 [49,] 329.161773 818.275166 [50,] -121.639573 329.161773 [51,] 1192.597994 -121.639573 [52,] 3825.734858 1192.597994 [53,] 1898.113007 3825.734858 [54,] 3883.404235 1898.113007 [55,] -1663.508514 3883.404235 [56,] -3703.968891 -1663.508514 [57,] -777.589508 -3703.968891 [58,] 656.930989 -777.589508 [59,] -707.185995 656.930989 [60,] -745.130968 -707.185995 [61,] -164.659171 -745.130968 [62,] -120.864169 -164.659171 [63,] -715.556170 -120.864169 [64,] -1514.636577 -715.556170 [65,] -650.065442 -1514.636577 [66,] 217.786117 -650.065442 [67,] 318.539775 217.786117 [68,] 1989.611071 318.539775 [69,] 19.621278 1989.611071 [70,] -2214.551765 19.621278 [71,] -2783.762310 -2214.551765 [72,] 1173.976909 -2783.762310 [73,] -213.043438 1173.976909 [74,] 1246.830179 -213.043438 [75,] 82.895705 1246.830179 [76,] 1252.997086 82.895705 [77,] 22.872089 1252.997086 [78,] -676.919991 22.872089 [79,] 1592.900422 -676.919991 [80,] -1299.502738 1592.900422 [81,] -530.426761 -1299.502738 [82,] 2538.618480 -530.426761 [83,] -2052.259268 2538.618480 [84,] 635.543149 -2052.259268 [85,] -1470.195210 635.543149 [86,] -656.375725 -1470.195210 [87,] 829.144472 -656.375725 [88,] 1120.037178 829.144472 [89,] 609.702605 1120.037178 [90,] 91.538895 609.702605 [91,] -504.930538 91.538895 [92,] -1849.356020 -504.930538 [93,] 1524.364018 -1849.356020 [94,] -2564.966168 1524.364018 [95,] -1618.168056 -2564.966168 [96,] 849.475477 -1618.168056 [97,] -761.256688 849.475477 [98,] 1699.459918 -761.256688 [99,] -533.611160 1699.459918 [100,] -328.546992 -533.611160 [101,] 821.087835 -328.546992 [102,] 1179.268373 821.087835 [103,] 1533.402505 1179.268373 [104,] -1072.452550 1533.402505 [105,] -304.685718 -1072.452550 [106,] -199.410741 -304.685718 [107,] -481.158937 -199.410741 [108,] 561.232612 -481.158937 [109,] -191.644522 561.232612 [110,] -1300.599370 -191.644522 [111,] -985.347261 -1300.599370 [112,] 10.897599 -985.347261 [113,] -399.340499 10.897599 [114,] -1661.614585 -399.340499 [115,] -1532.042011 -1661.614585 [116,] 11.808421 -1532.042011 [117,] -562.540244 11.808421 [118,] 202.842125 -562.540244 [119,] -1088.325075 202.842125 [120,] 375.001446 -1088.325075 [121,] -914.745069 375.001446 [122,] 179.341066 -914.745069 [123,] 185.809164 179.341066 [124,] -1335.232702 185.809164 [125,] 47.161008 -1335.232702 [126,] -1061.735897 47.161008 [127,] -1061.158843 -1061.735897 [128,] -1245.133403 -1061.158843 [129,] -1061.735897 -1245.133403 [130,] -1063.834583 -1061.735897 [131,] -1725.494028 -1063.834583 [132,] -1061.735897 -1725.494028 [133,] -1045.872342 -1061.735897 [134,] -1082.659339 -1045.872342 [135,] -1090.258677 -1082.659339 [136,] -1080.264849 -1090.258677 [137,] -1059.405046 -1080.264849 [138,] -1114.409716 -1059.405046 [139,] -1061.735897 -1114.409716 [140,] -1231.006356 -1061.735897 [141,] -1061.735897 -1231.006356 [142,] -1223.175503 -1061.735897 [143,] -1058.869179 -1223.175503 [144,] -1047.115252 -1058.869179 [145,] -1274.122439 -1047.115252 [146,] -1061.735897 -1274.122439 [147,] -1055.892234 -1061.735897 [148,] -1135.815656 -1055.892234 [149,] -1050.887300 -1135.815656 [150,] -1059.353040 -1050.887300 [151,] -1138.519584 -1059.353040 [152,] -1062.843185 -1138.519584 [153,] -1092.243375 -1062.843185 [154,] -1051.874989 -1092.243375 [155,] -1061.735897 -1051.874989 [156,] -1061.735897 -1061.735897 [157,] -1029.374768 -1061.735897 [158,] -1203.209376 -1029.374768 [159,] -1061.735897 -1203.209376 [160,] -1061.735897 -1061.735897 [161,] -1061.735897 -1061.735897 [162,] -1061.735897 -1061.735897 [163,] -1034.423174 -1061.735897 [164,] -893.054439 -1034.423174 [165,] -1061.927612 -893.054439 [166,] -1061.735897 -1061.927612 [167,] -2242.074409 -1061.735897 [168,] -1508.435528 -2242.074409 [169,] -1061.735897 -1508.435528 [170,] -1061.735897 -1061.735897 [171,] -1005.825011 -1061.735897 [172,] -1061.735897 -1005.825011 [173,] -982.565463 -1061.735897 [174,] -1061.735897 -982.565463 [175,] -1060.752605 -1061.735897 [176,] -1098.819809 -1060.752605 [177,] -1069.429333 -1098.819809 [178,] -1061.735897 -1069.429333 [179,] -1061.735897 -1061.735897 [180,] -841.713596 -1061.735897 [181,] -643.649635 -841.713596 [182,] -1453.169244 -643.649635 [183,] -434.857855 -1453.169244 [184,] -970.539261 -434.857855 [185,] -1223.462113 -970.539261 [186,] -1214.789853 -1223.462113 [187,] -995.191188 -1214.789853 [188,] -692.406099 -995.191188 [189,] -1061.735897 -692.406099 [190,] -1159.573192 -1061.735897 [191,] -1061.735897 -1159.573192 [192,] -1026.952296 -1061.735897 [193,] -1061.735897 -1026.952296 [194,] -1057.512909 -1061.735897 [195,] -1061.735897 -1057.512909 [196,] 670.444619 -1061.735897 [197,] -1061.735897 670.444619 [198,] -1061.735897 -1061.735897 [199,] 122.450100 -1061.735897 [200,] -1061.735897 122.450100 [201,] -1061.735897 -1061.735897 [202,] -1061.735897 -1061.735897 [203,] -1481.573495 -1061.735897 [204,] 160.496827 -1481.573495 [205,] -1021.841361 160.496827 [206,] -839.829927 -1021.841361 [207,] -1061.735897 -839.829927 [208,] -1900.832563 -1061.735897 [209,] -1029.650703 -1900.832563 [210,] -139.270215 -1029.650703 [211,] 127.729785 -139.270215 [212,] 1739.108604 127.729785 [213,] -1184.234238 1739.108604 [214,] -1096.934730 -1184.234238 [215,] -124.391717 -1096.934730 [216,] 247.952847 -124.391717 [217,] -7.271879 247.952847 [218,] 1726.608283 -7.271879 [219,] 1484.825340 1726.608283 [220,] -1015.125087 1484.825340 [221,] -136.270215 -1015.125087 [222,] 332.260570 -136.270215 [223,] -279.024713 332.260570 [224,] 129.498368 -279.024713 [225,] -253.640739 129.498368 [226,] -180.646183 -253.640739 [227,] 52.729785 -180.646183 [228,] -1298.125172 52.729785 [229,] 177.729785 -1298.125172 [230,] -930.761739 177.729785 [231,] 155.492980 -930.761739 [232,] 228.729785 155.492980 [233,] -952.389572 228.729785 [234,] -137.270215 -952.389572 [235,] -139.093865 -137.270215 [236,] -86.270215 -139.093865 [237,] -37.133233 -86.270215 [238,] -139.270215 -37.133233 [239,] -139.270215 -139.270215 [240,] -45.270215 -139.270215 [241,] -403.968451 -45.270215 [242,] -115.270215 -403.968451 [243,] 2404.716804 -115.270215 [244,] -2344.517538 2404.716804 [245,] -139.270215 -2344.517538 [246,] -8.270215 -139.270215 [247,] -695.891470 -8.270215 [248,] -139.270215 -695.891470 [249,] 66.729785 -139.270215 [250,] 625.561964 66.729785 [251,] 27.729785 625.561964 [252,] -694.834735 27.729785 [253,] -1076.587638 -694.834735 [254,] -8060.784634 -1076.587638 [255,] 225.729785 -8060.784634 [256,] 1244.262432 225.729785 [257,] 1908.386400 1244.262432 [258,] -139.270215 1908.386400 [259,] -139.270215 -139.270215 [260,] -139.270215 -139.270215 [261,] 86.729785 -139.270215 [262,] -1109.947060 86.729785 [263,] 1089.257151 -1109.947060 [264,] -139.270215 1089.257151 [265,] -139.270215 -139.270215 [266,] 48.729785 -139.270215 [267,] -649.474595 48.729785 [268,] -1.270215 -649.474595 [269,] -1365.460633 -1.270215 [270,] -139.270215 -1365.460633 [271,] -139.270215 -139.270215 [272,] -14.270215 -139.270215 [273,] -412.769924 -14.270215 [274,] 142.729785 -412.769924 [275,] 39.764643 142.729785 [276,] -658.661062 39.764643 [277,] 1184.729785 -658.661062 [278,] -2612.549884 1184.729785 [279,] -2014.058505 -2612.549884 [280,] -139.270215 -2014.058505 [281,] 109.729785 -139.270215 [282,] 650.616351 109.729785 [283,] 193.729785 650.616351 [284,] -1589.334055 193.729785 [285,] 461.729785 -1589.334055 [286,] 79.057971 461.729785 [287,] -299.019094 79.057971 [288,] 109.729785 -299.019094 [289,] 162.247354 109.729785 [290,] 25.729785 162.247354 [291,] 325.753632 25.729785 [292,] -2168.191760 325.753632 [293,] -86.270215 -2168.191760 [294,] -127.865024 -86.270215 [295,] 1420.496699 -127.865024 [296,] -139.270215 1420.496699 [297,] -139.270215 -139.270215 [298,] -139.270215 -139.270215 [299,] -109.816972 -139.270215 [300,] 24.072522 -109.816972 [301,] -412.072342 24.072522 [302,] -139.330966 -412.072342 [303,] 226.729785 -139.330966 [304,] 1051.433365 226.729785 [305,] -142.631572 1051.433365 [306,] 69.729785 -142.631572 [307,] -423.947404 69.729785 [308,] 923.587931 -423.947404 [309,] -139.270215 923.587931 [310,] 225.729785 -139.270215 [311,] -4465.476016 225.729785 [312,] -90.270215 -4465.476016 [313,] -559.628047 -90.270215 [314,] -2.242662 -559.628047 [315,] -89.365450 -2.242662 [316,] 370.348242 -89.365450 [317,] -1136.768660 370.348242 [318,] -139.746521 -1136.768660 [319,] -139.270215 -139.746521 [320,] -139.270215 -139.270215 [321,] -139.270215 -139.270215 [322,] -139.270215 -139.270215 [323,] -139.270215 -139.270215 [324,] -117.270215 -139.270215 [325,] -281.518718 -117.270215 [326,] -139.270215 -281.518718 [327,] -139.270215 -139.270215 [328,] -139.270215 -139.270215 [329,] -139.270215 -139.270215 [330,] -139.270215 -139.270215 [331,] -139.270215 -139.270215 [332,] -43.270215 -139.270215 [333,] -144.309852 -43.270215 [334,] 173.512840 -144.309852 [335,] 704.043521 173.512840 [336,] -1133.525590 704.043521 [337,] -113.270215 -1133.525590 [338,] -127.184237 -113.270215 [339,] 496.643627 -127.184237 [340,] -1613.230251 496.643627 [341,] -139.270215 -1613.230251 [342,] -139.270215 -139.270215 [343,] 481.729785 -139.270215 [344,] -490.925235 481.729785 [345,] -20.270215 -490.925235 [346,] -361.063762 -20.270215 [347,] -139.270215 -361.063762 [348,] 1455.729785 -139.270215 [349,] 368.585152 1455.729785 [350,] 553.501829 368.585152 [351,] 292.282379 553.501829 [352,] 752.484755 292.282379 [353,] -494.815251 752.484755 [354,] -139.270215 -494.815251 [355,] 374.729785 -139.270215 [356,] -3795.045538 374.729785 [357,] -139.270215 -3795.045538 [358,] -139.270215 -139.270215 [359,] -138.270215 -139.270215 [360,] -139.069613 -138.270215 [361,] -139.270215 -139.069613 [362,] 1623.729785 -139.270215 [363,] -502.779063 1623.729785 [364,] 150.011155 -502.779063 [365,] -139.270215 150.011155 [366,] -139.270215 -139.270215 [367,] -139.270215 -139.270215 [368,] 78.729785 -139.270215 [369,] -131.175020 78.729785 [370,] 308.729785 -131.175020 [371,] 688.326689 308.729785 [372,] -2445.155383 688.326689 [373,] -1439.270720 -2445.155383 [374,] -139.270215 -1439.270720 [375,] -18.270215 -139.270215 [376,] 487.406100 -18.270215 [377,] 2477.986548 487.406100 [378,] 80.964476 2477.986548 [379,] -139.270215 80.964476 [380,] 390.729785 -139.270215 [381,] -2290.634823 390.729785 [382,] 2398.483194 -2290.634823 [383,] -61.270215 2398.483194 [384,] 2006.554552 -61.270215 [385,] -5383.735718 2006.554552 [386,] 1593.866401 -5383.735718 [387,] 85.075320 1593.866401 [388,] 446.233013 85.075320 [389,] 364.247630 446.233013 [390,] 211.575315 364.247630 [391,] -353.489569 211.575315 [392,] -2062.862285 -353.489569 [393,] 4234.482677 -2062.862285 [394,] -6199.111883 4234.482677 [395,] 3712.798981 -6199.111883 [396,] -480.110115 3712.798981 [397,] -6949.088627 -480.110115 [398,] 347.097654 -6949.088627 [399,] 5540.850501 347.097654 [400,] -2058.527317 5540.850501 [401,] 586.628153 -2058.527317 [402,] 815.992878 586.628153 [403,] -3252.172767 815.992878 [404,] -1639.129780 -3252.172767 [405,] -1155.084059 -1639.129780 [406,] -682.175458 -1155.084059 [407,] 161.294988 -682.175458 [408,] 266.453118 161.294988 [409,] -5000.528162 266.453118 [410,] 4958.071007 -5000.528162 [411,] -937.482146 4958.071007 [412,] 1675.626363 -937.482146 [413,] 5109.128002 1675.626363 [414,] -10702.652116 5109.128002 [415,] -5180.946020 -10702.652116 [416,] 11155.125275 -5180.946020 [417,] -443.565059 11155.125275 [418,] -3332.403928 -443.565059 [419,] -3668.241362 -3332.403928 [420,] -2348.678631 -3668.241362 [421,] 5066.650137 -2348.678631 [422,] 3417.018967 5066.650137 [423,] -14401.817777 3417.018967 [424,] 4095.498813 -14401.817777 [425,] -2470.243562 4095.498813 [426,] 24153.403695 -2470.243562 [427,] -12733.619516 24153.403695 [428,] 7888.505936 -12733.619516 [429,] 17046.328143 7888.505936 [430,] 153933.352077 17046.328143 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 3189.149011 18112.030820 2 -5975.240229 3189.149011 3 -5896.762010 -5975.240229 4 -4963.274386 -5896.762010 5 -3586.218801 -4963.274386 6 -3639.496594 -3586.218801 7 -1771.217901 -3639.496594 8 -5120.670902 -1771.217901 9 -7276.418018 -5120.670902 10 3969.352621 -7276.418018 11 -4492.098869 3969.352621 12 -10688.605698 -4492.098869 13 -3692.210291 -10688.605698 14 -3658.988881 -3692.210291 15 -571.565104 -3658.988881 16 8580.303693 -571.565104 17 -2205.676265 8580.303693 18 -7792.585933 -2205.676265 19 749.071535 -7792.585933 20 -3340.109525 749.071535 21 -3344.907924 -3340.109525 22 1462.809853 -3344.907924 23 760.744286 1462.809853 24 2502.076015 760.744286 25 -3164.005247 2502.076015 26 -3829.109271 -3164.005247 27 -2410.225134 -3829.109271 28 -5503.404364 -2410.225134 29 -3004.098847 -5503.404364 30 -1916.099554 -3004.098847 31 -860.802328 -1916.099554 32 -2790.326038 -860.802328 33 -779.212044 -2790.326038 34 8916.934765 -779.212044 35 -407.005211 8916.934765 36 -2080.687204 -407.005211 37 -935.992092 -2080.687204 38 3988.236079 -935.992092 39 346.561017 3988.236079 40 1884.552701 346.561017 41 -220.002535 1884.552701 42 -705.251195 -220.002535 43 1014.684702 -705.251195 44 -1268.559636 1014.684702 45 446.715966 -1268.559636 46 5762.193305 446.715966 47 1607.948969 5762.193305 48 818.275166 1607.948969 49 329.161773 818.275166 50 -121.639573 329.161773 51 1192.597994 -121.639573 52 3825.734858 1192.597994 53 1898.113007 3825.734858 54 3883.404235 1898.113007 55 -1663.508514 3883.404235 56 -3703.968891 -1663.508514 57 -777.589508 -3703.968891 58 656.930989 -777.589508 59 -707.185995 656.930989 60 -745.130968 -707.185995 61 -164.659171 -745.130968 62 -120.864169 -164.659171 63 -715.556170 -120.864169 64 -1514.636577 -715.556170 65 -650.065442 -1514.636577 66 217.786117 -650.065442 67 318.539775 217.786117 68 1989.611071 318.539775 69 19.621278 1989.611071 70 -2214.551765 19.621278 71 -2783.762310 -2214.551765 72 1173.976909 -2783.762310 73 -213.043438 1173.976909 74 1246.830179 -213.043438 75 82.895705 1246.830179 76 1252.997086 82.895705 77 22.872089 1252.997086 78 -676.919991 22.872089 79 1592.900422 -676.919991 80 -1299.502738 1592.900422 81 -530.426761 -1299.502738 82 2538.618480 -530.426761 83 -2052.259268 2538.618480 84 635.543149 -2052.259268 85 -1470.195210 635.543149 86 -656.375725 -1470.195210 87 829.144472 -656.375725 88 1120.037178 829.144472 89 609.702605 1120.037178 90 91.538895 609.702605 91 -504.930538 91.538895 92 -1849.356020 -504.930538 93 1524.364018 -1849.356020 94 -2564.966168 1524.364018 95 -1618.168056 -2564.966168 96 849.475477 -1618.168056 97 -761.256688 849.475477 98 1699.459918 -761.256688 99 -533.611160 1699.459918 100 -328.546992 -533.611160 101 821.087835 -328.546992 102 1179.268373 821.087835 103 1533.402505 1179.268373 104 -1072.452550 1533.402505 105 -304.685718 -1072.452550 106 -199.410741 -304.685718 107 -481.158937 -199.410741 108 561.232612 -481.158937 109 -191.644522 561.232612 110 -1300.599370 -191.644522 111 -985.347261 -1300.599370 112 10.897599 -985.347261 113 -399.340499 10.897599 114 -1661.614585 -399.340499 115 -1532.042011 -1661.614585 116 11.808421 -1532.042011 117 -562.540244 11.808421 118 202.842125 -562.540244 119 -1088.325075 202.842125 120 375.001446 -1088.325075 121 -914.745069 375.001446 122 179.341066 -914.745069 123 185.809164 179.341066 124 -1335.232702 185.809164 125 47.161008 -1335.232702 126 -1061.735897 47.161008 127 -1061.158843 -1061.735897 128 -1245.133403 -1061.158843 129 -1061.735897 -1245.133403 130 -1063.834583 -1061.735897 131 -1725.494028 -1063.834583 132 -1061.735897 -1725.494028 133 -1045.872342 -1061.735897 134 -1082.659339 -1045.872342 135 -1090.258677 -1082.659339 136 -1080.264849 -1090.258677 137 -1059.405046 -1080.264849 138 -1114.409716 -1059.405046 139 -1061.735897 -1114.409716 140 -1231.006356 -1061.735897 141 -1061.735897 -1231.006356 142 -1223.175503 -1061.735897 143 -1058.869179 -1223.175503 144 -1047.115252 -1058.869179 145 -1274.122439 -1047.115252 146 -1061.735897 -1274.122439 147 -1055.892234 -1061.735897 148 -1135.815656 -1055.892234 149 -1050.887300 -1135.815656 150 -1059.353040 -1050.887300 151 -1138.519584 -1059.353040 152 -1062.843185 -1138.519584 153 -1092.243375 -1062.843185 154 -1051.874989 -1092.243375 155 -1061.735897 -1051.874989 156 -1061.735897 -1061.735897 157 -1029.374768 -1061.735897 158 -1203.209376 -1029.374768 159 -1061.735897 -1203.209376 160 -1061.735897 -1061.735897 161 -1061.735897 -1061.735897 162 -1061.735897 -1061.735897 163 -1034.423174 -1061.735897 164 -893.054439 -1034.423174 165 -1061.927612 -893.054439 166 -1061.735897 -1061.927612 167 -2242.074409 -1061.735897 168 -1508.435528 -2242.074409 169 -1061.735897 -1508.435528 170 -1061.735897 -1061.735897 171 -1005.825011 -1061.735897 172 -1061.735897 -1005.825011 173 -982.565463 -1061.735897 174 -1061.735897 -982.565463 175 -1060.752605 -1061.735897 176 -1098.819809 -1060.752605 177 -1069.429333 -1098.819809 178 -1061.735897 -1069.429333 179 -1061.735897 -1061.735897 180 -841.713596 -1061.735897 181 -643.649635 -841.713596 182 -1453.169244 -643.649635 183 -434.857855 -1453.169244 184 -970.539261 -434.857855 185 -1223.462113 -970.539261 186 -1214.789853 -1223.462113 187 -995.191188 -1214.789853 188 -692.406099 -995.191188 189 -1061.735897 -692.406099 190 -1159.573192 -1061.735897 191 -1061.735897 -1159.573192 192 -1026.952296 -1061.735897 193 -1061.735897 -1026.952296 194 -1057.512909 -1061.735897 195 -1061.735897 -1057.512909 196 670.444619 -1061.735897 197 -1061.735897 670.444619 198 -1061.735897 -1061.735897 199 122.450100 -1061.735897 200 -1061.735897 122.450100 201 -1061.735897 -1061.735897 202 -1061.735897 -1061.735897 203 -1481.573495 -1061.735897 204 160.496827 -1481.573495 205 -1021.841361 160.496827 206 -839.829927 -1021.841361 207 -1061.735897 -839.829927 208 -1900.832563 -1061.735897 209 -1029.650703 -1900.832563 210 -139.270215 -1029.650703 211 127.729785 -139.270215 212 1739.108604 127.729785 213 -1184.234238 1739.108604 214 -1096.934730 -1184.234238 215 -124.391717 -1096.934730 216 247.952847 -124.391717 217 -7.271879 247.952847 218 1726.608283 -7.271879 219 1484.825340 1726.608283 220 -1015.125087 1484.825340 221 -136.270215 -1015.125087 222 332.260570 -136.270215 223 -279.024713 332.260570 224 129.498368 -279.024713 225 -253.640739 129.498368 226 -180.646183 -253.640739 227 52.729785 -180.646183 228 -1298.125172 52.729785 229 177.729785 -1298.125172 230 -930.761739 177.729785 231 155.492980 -930.761739 232 228.729785 155.492980 233 -952.389572 228.729785 234 -137.270215 -952.389572 235 -139.093865 -137.270215 236 -86.270215 -139.093865 237 -37.133233 -86.270215 238 -139.270215 -37.133233 239 -139.270215 -139.270215 240 -45.270215 -139.270215 241 -403.968451 -45.270215 242 -115.270215 -403.968451 243 2404.716804 -115.270215 244 -2344.517538 2404.716804 245 -139.270215 -2344.517538 246 -8.270215 -139.270215 247 -695.891470 -8.270215 248 -139.270215 -695.891470 249 66.729785 -139.270215 250 625.561964 66.729785 251 27.729785 625.561964 252 -694.834735 27.729785 253 -1076.587638 -694.834735 254 -8060.784634 -1076.587638 255 225.729785 -8060.784634 256 1244.262432 225.729785 257 1908.386400 1244.262432 258 -139.270215 1908.386400 259 -139.270215 -139.270215 260 -139.270215 -139.270215 261 86.729785 -139.270215 262 -1109.947060 86.729785 263 1089.257151 -1109.947060 264 -139.270215 1089.257151 265 -139.270215 -139.270215 266 48.729785 -139.270215 267 -649.474595 48.729785 268 -1.270215 -649.474595 269 -1365.460633 -1.270215 270 -139.270215 -1365.460633 271 -139.270215 -139.270215 272 -14.270215 -139.270215 273 -412.769924 -14.270215 274 142.729785 -412.769924 275 39.764643 142.729785 276 -658.661062 39.764643 277 1184.729785 -658.661062 278 -2612.549884 1184.729785 279 -2014.058505 -2612.549884 280 -139.270215 -2014.058505 281 109.729785 -139.270215 282 650.616351 109.729785 283 193.729785 650.616351 284 -1589.334055 193.729785 285 461.729785 -1589.334055 286 79.057971 461.729785 287 -299.019094 79.057971 288 109.729785 -299.019094 289 162.247354 109.729785 290 25.729785 162.247354 291 325.753632 25.729785 292 -2168.191760 325.753632 293 -86.270215 -2168.191760 294 -127.865024 -86.270215 295 1420.496699 -127.865024 296 -139.270215 1420.496699 297 -139.270215 -139.270215 298 -139.270215 -139.270215 299 -109.816972 -139.270215 300 24.072522 -109.816972 301 -412.072342 24.072522 302 -139.330966 -412.072342 303 226.729785 -139.330966 304 1051.433365 226.729785 305 -142.631572 1051.433365 306 69.729785 -142.631572 307 -423.947404 69.729785 308 923.587931 -423.947404 309 -139.270215 923.587931 310 225.729785 -139.270215 311 -4465.476016 225.729785 312 -90.270215 -4465.476016 313 -559.628047 -90.270215 314 -2.242662 -559.628047 315 -89.365450 -2.242662 316 370.348242 -89.365450 317 -1136.768660 370.348242 318 -139.746521 -1136.768660 319 -139.270215 -139.746521 320 -139.270215 -139.270215 321 -139.270215 -139.270215 322 -139.270215 -139.270215 323 -139.270215 -139.270215 324 -117.270215 -139.270215 325 -281.518718 -117.270215 326 -139.270215 -281.518718 327 -139.270215 -139.270215 328 -139.270215 -139.270215 329 -139.270215 -139.270215 330 -139.270215 -139.270215 331 -139.270215 -139.270215 332 -43.270215 -139.270215 333 -144.309852 -43.270215 334 173.512840 -144.309852 335 704.043521 173.512840 336 -1133.525590 704.043521 337 -113.270215 -1133.525590 338 -127.184237 -113.270215 339 496.643627 -127.184237 340 -1613.230251 496.643627 341 -139.270215 -1613.230251 342 -139.270215 -139.270215 343 481.729785 -139.270215 344 -490.925235 481.729785 345 -20.270215 -490.925235 346 -361.063762 -20.270215 347 -139.270215 -361.063762 348 1455.729785 -139.270215 349 368.585152 1455.729785 350 553.501829 368.585152 351 292.282379 553.501829 352 752.484755 292.282379 353 -494.815251 752.484755 354 -139.270215 -494.815251 355 374.729785 -139.270215 356 -3795.045538 374.729785 357 -139.270215 -3795.045538 358 -139.270215 -139.270215 359 -138.270215 -139.270215 360 -139.069613 -138.270215 361 -139.270215 -139.069613 362 1623.729785 -139.270215 363 -502.779063 1623.729785 364 150.011155 -502.779063 365 -139.270215 150.011155 366 -139.270215 -139.270215 367 -139.270215 -139.270215 368 78.729785 -139.270215 369 -131.175020 78.729785 370 308.729785 -131.175020 371 688.326689 308.729785 372 -2445.155383 688.326689 373 -1439.270720 -2445.155383 374 -139.270215 -1439.270720 375 -18.270215 -139.270215 376 487.406100 -18.270215 377 2477.986548 487.406100 378 80.964476 2477.986548 379 -139.270215 80.964476 380 390.729785 -139.270215 381 -2290.634823 390.729785 382 2398.483194 -2290.634823 383 -61.270215 2398.483194 384 2006.554552 -61.270215 385 -5383.735718 2006.554552 386 1593.866401 -5383.735718 387 85.075320 1593.866401 388 446.233013 85.075320 389 364.247630 446.233013 390 211.575315 364.247630 391 -353.489569 211.575315 392 -2062.862285 -353.489569 393 4234.482677 -2062.862285 394 -6199.111883 4234.482677 395 3712.798981 -6199.111883 396 -480.110115 3712.798981 397 -6949.088627 -480.110115 398 347.097654 -6949.088627 399 5540.850501 347.097654 400 -2058.527317 5540.850501 401 586.628153 -2058.527317 402 815.992878 586.628153 403 -3252.172767 815.992878 404 -1639.129780 -3252.172767 405 -1155.084059 -1639.129780 406 -682.175458 -1155.084059 407 161.294988 -682.175458 408 266.453118 161.294988 409 -5000.528162 266.453118 410 4958.071007 -5000.528162 411 -937.482146 4958.071007 412 1675.626363 -937.482146 413 5109.128002 1675.626363 414 -10702.652116 5109.128002 415 -5180.946020 -10702.652116 416 11155.125275 -5180.946020 417 -443.565059 11155.125275 418 -3332.403928 -443.565059 419 -3668.241362 -3332.403928 420 -2348.678631 -3668.241362 421 5066.650137 -2348.678631 422 3417.018967 5066.650137 423 -14401.817777 3417.018967 424 4095.498813 -14401.817777 425 -2470.243562 4095.498813 426 24153.403695 -2470.243562 427 -12733.619516 24153.403695 428 7888.505936 -12733.619516 429 17046.328143 7888.505936 430 153933.352077 17046.328143 > plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals') > lines(lowess(z)) > abline(lm(z)) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/70ok61291297725.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/80ok61291297725.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/90ok61291297725.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0)) > plot(mylm, las = 1, sub='Residual Diagnostics') > par(opar) > dev.off() null device 1 > if (n > n25) { + postscript(file="/var/www/html/rcomp/tmp/10bf291291297725.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) + plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint') + grid() + dev.off() + } null device 1 > > #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE) > a<-table.row.end(a) > myeq <- colnames(x)[1] > myeq <- paste(myeq, '[t] = ', sep='') > for (i in 1:k){ + if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '') + myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ') + if (rownames(mysum$coefficients)[i] != '(Intercept)') { + myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='') + if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='') + } + } > myeq <- paste(myeq, ' + e[t]') > a<-table.row.start(a) > a<-table.element(a, myeq) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/11wx0x1291297725.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Variable',header=TRUE) > a<-table.element(a,'Parameter',header=TRUE) > a<-table.element(a,'S.D.',header=TRUE) > a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE) > a<-table.element(a,'2-tail p-value',header=TRUE) > a<-table.element(a,'1-tail p-value',header=TRUE) > a<-table.row.end(a) > for (i in 1:k){ + a<-table.row.start(a) + a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE) + a<-table.element(a,mysum$coefficients[i,1]) + a<-table.element(a, round(mysum$coefficients[i,2],6)) + a<-table.element(a, round(mysum$coefficients[i,3],4)) + a<-table.element(a, round(mysum$coefficients[i,4],6)) + a<-table.element(a, round(mysum$coefficients[i,4]/2,6)) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/12zgzl1291297725.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple R',1,TRUE) > a<-table.element(a, sqrt(mysum$r.squared)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'R-squared',1,TRUE) > a<-table.element(a, mysum$r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Adjusted R-squared',1,TRUE) > a<-table.element(a, mysum$adj.r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (value)',1,TRUE) > a<-table.element(a, mysum$fstatistic[1]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[2]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[3]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'p-value',1,TRUE) > a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3])) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Residual Standard Deviation',1,TRUE) > a<-table.element(a, mysum$sigma) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Sum Squared Residuals',1,TRUE) > a<-table.element(a, sum(myerror*myerror)) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/13v7wc1291297725.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Time or Index', 1, TRUE) > a<-table.element(a, 'Actuals', 1, TRUE) > a<-table.element(a, 'Interpolation
Forecast', 1, TRUE) > a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,i, 1, TRUE) + a<-table.element(a,x[i]) + a<-table.element(a,x[i]-mysum$resid[i]) + a<-table.element(a,mysum$resid[i]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/14hqv01291297725.tab") > if (n > n25) { + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'p-values',header=TRUE) + a<-table.element(a,'Alternative Hypothesis',3,header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'breakpoint index',header=TRUE) + a<-table.element(a,'greater',header=TRUE) + a<-table.element(a,'2-sided',header=TRUE) + a<-table.element(a,'less',header=TRUE) + a<-table.row.end(a) + for (mypoint in kp3:nmkm3) { + a<-table.row.start(a) + a<-table.element(a,mypoint,header=TRUE) + a<-table.element(a,gqarr[mypoint-kp3+1,1]) + a<-table.element(a,gqarr[mypoint-kp3+1,2]) + a<-table.element(a,gqarr[mypoint-kp3+1,3]) + a<-table.row.end(a) + } + a<-table.end(a) + table.save(a,file="/var/www/html/rcomp/tmp/15kqb51291297725.tab") + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'Description',header=TRUE) + a<-table.element(a,'# significant tests',header=TRUE) + a<-table.element(a,'% significant tests',header=TRUE) + a<-table.element(a,'OK/NOK',header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'1% type I error level',header=TRUE) + a<-table.element(a,numsignificant1) + a<-table.element(a,numsignificant1/numgqtests) + if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'5% type I error level',header=TRUE) + a<-table.element(a,numsignificant5) + a<-table.element(a,numsignificant5/numgqtests) + if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'10% type I error level',header=TRUE) + a<-table.element(a,numsignificant10) + a<-table.element(a,numsignificant10/numgqtests) + if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.end(a) + table.save(a,file="/var/www/html/rcomp/tmp/166rat1291297725.tab") + } > > try(system("convert tmp/1menx1291297725.ps tmp/1menx1291297725.png",intern=TRUE)) character(0) > try(system("convert tmp/2menx1291297725.ps tmp/2menx1291297725.png",intern=TRUE)) character(0) > try(system("convert tmp/3f5401291297725.ps tmp/3f5401291297725.png",intern=TRUE)) character(0) > try(system("convert tmp/4f5401291297725.ps tmp/4f5401291297725.png",intern=TRUE)) character(0) > try(system("convert tmp/5f5401291297725.ps tmp/5f5401291297725.png",intern=TRUE)) character(0) > try(system("convert tmp/6pell1291297725.ps tmp/6pell1291297725.png",intern=TRUE)) character(0) > try(system("convert tmp/70ok61291297725.ps tmp/70ok61291297725.png",intern=TRUE)) character(0) > try(system("convert tmp/80ok61291297725.ps tmp/80ok61291297725.png",intern=TRUE)) character(0) > try(system("convert tmp/90ok61291297725.ps tmp/90ok61291297725.png",intern=TRUE)) character(0) > try(system("convert tmp/10bf291291297725.ps tmp/10bf291291297725.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 12.029 2.177 25.622