R version 2.12.0 (2010-10-15) Copyright (C) 2010 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-pc-linux-gnu (32-bit) 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(1 + ,1595 + ,17 + ,60720 + ,319369 + ,0 + ,5565 + ,47 + ,94355 + ,493408 + ,0 + ,601 + ,6 + ,60720 + ,319210 + ,1 + ,188 + ,11 + ,77655 + ,381180 + ,1 + ,7146 + ,105 + ,134028 + ,297978 + ,0 + ,1135 + ,17 + ,62285 + ,290476 + ,1 + ,450 + ,8 + ,59325 + ,292136 + ,1 + ,34 + ,8 + ,60630 + ,314353 + ,1 + ,133 + ,14 + ,65990 + ,339445 + ,1 + ,119 + ,6 + ,59118 + ,303677 + ,0 + ,2053 + ,53 + ,100423 + ,397144 + ,0 + ,4036 + ,63 + ,100269 + ,424898 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,4655 + ,393 + ,60720 + ,341570 + ,1 + ,131 + ,11 + ,61808 + ,308989 + ,1 + ,1766 + ,73 + ,60438 + ,305959 + ,1 + ,312 + ,14 + ,58598 + ,318690 + ,1 + ,448 + ,40 + ,59781 + ,323361 + ,1 + ,115 + ,13 + ,60945 + ,318903 + ,1 + ,60 + ,10 + ,61124 + ,314049 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,364 + ,35 + ,47705 + ,298466 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1442 + ,118 + ,121173 + ,438493 + ,0 + ,1389 + ,9 + ,57530 + ,378049 + ,1 + ,149 + ,5 + ,62041 + ,313332 + ,1 + ,2212 + ,14 + ,60720 + ,319864 + ,0 + ,7489 + ,322 + ,136996 + ,430866 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,402 + ,26 + ,84990 + ,332743 + ,0 + ,7419 + ,29 + ,63255 + ,286963 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,307 + ,15 + ,58856 + ,331955 + ,0 + ,1134 + ,162 + ,146216 + ,527021 + ,0 + ,7561 + ,87 + ,82425 + ,364304 + ,1 + ,5131 + ,71 + ,86111 + ,292154 + ,1 + ,9 + ,6 + ,60835 + ,314987 + ,1 + ,2264 + ,24 + ,55174 + ,272713 + ,1 + ,40949 + ,481 + ,129352 + ,1398893 + ,0 + ,4980 + ,91 + ,113521 + ,409642 + ,1 + ,272 + ,8 + ,112995 + ,429112 + ,1 + ,757 + ,19 + ,45689 + ,382712 + ,1 + ,1424 + ,94 + ,131116 + ,374943 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,378 + ,52 + ,71873 + ,297413 + ,1 + ,1061 + ,5 + ,70415 + ,304252 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,203 + ,9 + ,61938 + ,309596 + ,0 + ,3689 + ,25 + ,101792 + ,377305 + ,1 + ,567 + ,11 + ,103345 + ,386212 + ,0 + ,27664 + ,282 + ,74996 + ,657954 + ,1 + ,1686 + ,40 + ,68696 + ,368186 + ,0 + ,1164 + ,17 + ,57320 + ,269753 + ,1 + ,15824 + ,64 + ,63838 + ,585715 + ,1 + ,6575 + ,159 + ,151352 + ,1071292 + ,0 + ,37597 + ,248 + ,37884 + ,992426 + ,1 + ,314 + ,5 + ,60720 + ,315371 + ,0 + ,1932 + ,36 + ,72954 + ,249898 + ,0 + ,5692 + ,72 + ,64239 + ,357312 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,6533 + ,75 + ,126910 + ,364839 + ,0 + ,2126 + ,138 + ,21001 + ,230621 + ,1 + ,14 + ,1 + ,61016 + ,315877 + ,0 + ,3192 + ,60 + ,59831 + ,125390 + ,1 + ,47 + ,3 + ,60720 + ,314882 + ,0 + ,4891 + ,217 + ,210568 + ,370837 + ,0 + ,531 + ,10 + ,71641 + ,330068 + ,1 + ,2332 + ,4 + ,72680 + ,324385 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,16275 + ,109 + ,47818 + ,1443586 + ,1 + ,9252 + ,48 + ,79510 + ,253537 + ,0 + ,29790 + ,309 + ,81767 + ,4321023 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,2300 + ,29 + ,64873 + ,352108 + ,1 + ,267 + ,32 + ,57640 + ,330059 + ,0 + ,382 + ,21 + ,59365 + ,340968 + ,1 + ,397 + ,28 + ,60646 + ,317736 + ,1 + ,11920 + ,182 + ,31080 + ,209458 + ,0 + ,54660 + ,302 + ,126942 + ,1491348 + ,1 + ,5 + ,2 + ,60720 + ,314887 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,2647 + ,37 + ,62920 + ,353058 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,94 + ,7 + ,60793 + ,314533 + ,1 + ,3 + ,2 + ,60698 + ,315354 + ,1 + ,454 + ,51 + ,63261 + ,302187 + ,1 + ,227 + ,9 + ,76644 + ,336639 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,206 + ,9 + ,53110 + ,296702 + ,1 + ,25830 + ,115 + ,245546 + ,1073089 + ,1 + ,1528 + ,38 + ,60326 + ,146494 + ,1 + ,271 + ,2 + ,60720 + ,325249 + ,1 + ,138 + ,6 + ,69817 + ,331420 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,278 + ,26 + ,61404 + ,314922 + ,1 + ,282 + ,30 + ,62452 + ,320016 + ,1 + ,571 + ,80 + ,41477 + ,280398 + ,1 + ,2253 + ,99 + ,63593 + ,452469 + ,1 + ,290 + ,2 + ,58790 + ,301164 + ,1 + ,78 + ,16 + ,62700 + ,317330 + ,1 + ,20 + ,2 + ,60805 + ,315576 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,18073 + ,180 + ,27284 + ,-7170 + ,1 + ,1866 + ,29 + ,56225 + ,322331 + ,0 + ,42634 + ,156 + ,157214 + ,1629616 + ,1 + ,249 + ,62 + ,54323 + ,292754 + ,1 + ,422 + ,28 + ,57935 + ,318056 + ,1 + ,2675 + ,79 + ,59017 + ,355178 + ,1 + ,965 + ,3 + ,73490 + ,204325 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,621 + ,16 + ,68005 + ,317046 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,365 + ,8 + ,54820 + ,309560 + ,0 + ,3122 + ,34 + ,94670 + ,414462 + ,0 + ,12988 + ,119 + ,82340 + ,857217 + ,0 + ,5336 + ,81 + ,112477 + ,697458 + ,0 + ,3160 + ,108 + ,108094 + ,530670 + ,1 + ,2489 + ,40 + ,67804 + ,238125 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,6984 + ,45 + ,80570 + ,741409 + ,1 + ,2001 + ,52 + ,95551 + ,393343 + ,1 + ,15236 + ,90 + ,77440 + ,372631 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,530 + ,60 + ,73433 + ,317291 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,35624 + ,223 + ,157278 + ,306275 + ,0 + ,1383 + ,11 + ,73221 + ,317892 + ,1 + ,875 + ,14 + ,67000 + ,334280 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,72 + ,3 + ,60398 + ,314210 + ,0 + ,265 + ,31 + ,60720 + ,306948 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,335 + ,33 + ,64175 + ,320398 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,4525 + ,226 + ,93811 + ,501749 + ,0 + ,3045 + ,58 + ,27330 + ,202055 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,638 + ,14 + ,60370 + ,333210 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,607 + ,11 + ,60436 + ,322340 + ,1 + ,1558 + ,63 + ,55637 + ,369448 + ,0 + ,1324 + ,117 + ,67440 + ,291841 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,611 + ,8 + ,59190 + ,296919 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,923 + ,8 + ,58620 + ,309038 + ,1 + ,661 + ,3 + ,65920 + ,246541 + ,0 + ,2397 + ,14 + ,74020 + ,289513 + ,1 + ,366 + ,53 + ,61808 + ,344425 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,135 + ,2 + ,62065 + ,314210 + ,0 + ,1659 + ,29 + ,107577 + ,480382 + ,1 + ,316 + ,9 + ,60505 + ,315009 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,309 + ,11 + ,56535 + ,312878 + ,0 + ,49 + ,8 + ,64107 + ,322031 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,4519 + ,49 + ,102129 + ,597793 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,837 + ,69 + ,61262 + ,315688 + ,1 + ,5119 + ,49 + ,39039 + ,378525 + ,1 + ,1280 + ,117 + ,69465 + ,312378 + ,1 + ,2564 + ,22 + ,130140 + ,403560 + ,0 + ,2045 + ,20 + ,97890 + ,510834 + ,0 + ,2234 + ,16 + ,77200 + ,214215 + ,1 + ,975 + ,23 + ,90534 + ,235133 + ,0 + ,1136 + ,32 + ,48522 + ,343613 + ,1 + ,453 + ,21 + ,81125 + ,365959 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,61 + ,14 + ,60720 + ,314551 + ,1 + ,368 + ,17 + ,61686 + ,303230 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,4901 + ,84 + ,121920 + ,469107 + ,1 + ,540 + ,14 + ,103960 + ,354228 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,2 + ,9 + ,60735 + ,315394 + ,1 + ,36 + ,4 + ,61564 + ,312412 + ,1 + ,776 + ,9 + ,64230 + ,333505 + ,1 + ,84738 + ,588 + ,-26007 + ,223193 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,3 + ,4 + ,60761 + ,315656 + ,1 + ,529 + ,9 + ,63870 + ,296261 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,405 + ,7 + ,60845 + ,336425 + ,1 + ,972 + ,34 + ,71642 + ,359335 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,2099 + ,64 + ,106611 + ,308636 + ,0 + ,3437 + ,47 + ,48022 + ,158492 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,22330 + ,84 + ,79801 + ,711969 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,483 + ,27 + ,60830 + ,306268 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,2239 + ,21 + ,88590 + ,442882 + ,0 + ,2949 + ,41 + ,82903 + ,378509 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,365 + ,83 + ,87192 + ,346611 + ,1 + ,2461 + ,57 + ,55792 + ,314289 + ,0 + ,21950 + ,519 + ,114337 + ,856956 + ,0 + ,3294 + ,23 + ,75832 + ,217193 + ,1 + ,141 + ,6 + ,61630 + ,315366 + ,1 + ,572 + ,16 + ,58580 + ,307930 + ,1 + ,13326 + ,102 + ,165548 + ,702380 + ,1 + ,2284 + ,33 + ,76403 + ,194493 + ,0 + ,10 + ,2 + ,61656 + ,316155 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1414 + ,198 + ,70184 + ,330546 + ,1 + ,1975 + ,35 + ,118881 + ,394510 + ,1 + ,43 + ,3 + ,60887 + ,312846 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,844 + ,78 + ,60925 + ,296139 + ,1 + ,304 + ,15 + ,62969 + ,295580 + ,1 + ,458 + ,11 + ,58625 + ,297765 + ,1 + ,18562 + ,155 + ,102313 + ,377934 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,7123 + ,109 + ,288170 + ,638830 + ,1 + ,622 + ,73 + ,73007 + ,304376 + ,1 + ,174 + ,1 + ,64820 + ,307424 + ,1 + ,2220 + ,22 + ,301670 + ,644190 + ,1 + ,121 + ,24 + ,56178 + ,295370 + ,0 + ,11819 + ,85 + ,106113 + ,574339 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,125 + ,12 + ,60798 + ,310201 + ,1 + ,1182 + ,13 + ,70694 + ,327007 + ,1 + ,1503 + ,39 + ,56364 + ,343466 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,30 + ,4 + ,62045 + ,318098 + ,0 + ,3310 + ,33 + ,75230 + ,448243 + ,1 + ,554 + ,11 + ,79285 + ,325738 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,468 + ,21 + ,60720 + ,312161 + ,1 + ,4917 + ,42 + ,52811 + ,243650 + ,1 + ,3256 + ,66 + ,35250 + ,407159 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,125 + ,13 + ,59734 + ,317698 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,22 + ,2 + ,60722 + ,312502 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,514 + ,52 + ,78780 + ,322378 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,10327 + ,91 + ,88577 + ,640273 + ,1 + ,13718 + ,61 + ,96448 + ,345783 + ,1 + ,3748 + ,9 + ,50350 + ,652925 + ,1 + ,14416 + ,361 + ,49857 + ,439798 + ,0 + ,1526 + ,25 + ,69351 + ,278990 + ,0 + ,666 + ,171 + ,117869 + ,339836 + ,1 + ,2844 + ,67 + ,72683 + ,240897 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,368 + ,17 + ,61167 + ,297141 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,333 + ,13 + ,70811 + ,331323 + ,1 + ,26 + ,1 + ,60896 + ,313880 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,1303 + ,61 + ,69863 + ,309422 + ,1 + ,20 + ,3 + ,60938 + ,315245 + ,1 + ,2384 + ,97 + ,61348 + ,405972 + ,0 + ,203 + ,27 + ,50804 + ,300962 + ,0 + ,71 + ,1 + ,60745 + ,316176 + ,1 + ,53 + ,14 + ,59506 + ,302409 + ,1 + ,562 + ,13 + ,58456 + ,283587 + ,1 + ,622 + ,7 + ,60950 + ,263276 + ,1 + ,645 + ,5 + ,60720 + ,312075 + ,0 + ,1763 + ,24 + ,61600 + ,308336 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,317 + ,9 + ,63915 + ,298700 + ,1 + ,1 + ,6 + ,60719 + ,315372 + ,1 + ,275 + ,4 + ,59500 + ,318745 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,936 + ,30 + ,67939 + ,408881 + ,1 + ,8568 + ,37 + ,32168 + ,786690 + ,0 + ,11528 + ,70 + ,-14545 + ,-83265 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,738 + ,10 + ,60720 + ,321376 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,592 + ,7 + ,64270 + ,276898 + ,1 + ,126 + ,10 + ,60951 + ,315547 + ,1 + ,2 + ,1 + ,60743 + ,315487 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,18014 + ,239 + ,-1710 + ,1405225 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,37704 + ,452 + ,60448 + ,983660 + ,0 + ,63 + ,4 + ,65688 + ,318574 + ,1 + ,1431 + ,34 + ,106885 + ,310768 + ,1 + ,94 + ,3 + ,61360 + ,312887 + ,1 + ,192 + ,4 + ,65276 + ,312339 + ,1 + ,32 + ,1 + ,59988 + ,314964 + ,1 + ,6869 + ,105 + ,117520 + ,379983 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,1 + ,4 + ,60722 + ,315398 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,2328 + ,94 + ,82732 + ,253588 + ,0 + ,209 + ,13 + ,64016 + ,316647 + ,1 + ,28 + ,2 + ,60890 + ,315688 + ,1 + ,176 + ,28 + ,68136 + ,310670 + ,1 + ,1920 + ,15 + ,79420 + ,165404 + ,0 + ,87550 + ,458 + ,153198 + ,4111912 + ,1 + ,520 + ,10 + ,58650 + ,291650 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,1013 + ,2 + ,69770 + ,253468 + ,1 + ,15 + ,8 + ,60831 + ,315688 + ,1 + ,587 + ,12 + ,59595 + ,325699 + ,1 + ,5371 + ,66 + ,87720 + ,446211 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,1012 + ,47 + ,114768 + ,368078 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,876 + ,59 + ,138971 + ,352850 + ,0 + ,162556 + ,1081 + ,213118 + ,6282154 + ,1 + ,43556 + ,247 + ,32648 + ,227132 + ,0 + ,3425 + ,12 + ,83620 + ,283910 + ,1 + ,810 + ,43 + ,74015 + ,236761 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,2365 + ,79 + ,191778 + ,550608 + ,0 + ,1261 + ,33 + ,76114 + ,307528 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,1585 + ,107 + ,93099 + ,355864 + ,1 + ,16189 + ,295 + ,116384 + ,358589 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,10579 + ,92 + ,110309 + ,375195 + ,1 + ,474 + ,20 + ,61977 + ,288985 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,3642 + ,86 + ,90262 + ,458343 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,472 + ,20 + ,80045 + ,269587 + ,1 + ,98 + ,2 + ,61490 + ,315236 + ,1 + ,3999 + ,30 + ,51252 + ,42754 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,621 + ,48 + ,60720 + ,308256 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,30 + ,4 + ,60798 + ,313164 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,746 + ,23 + ,71561 + ,269661 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,4150 + ,150 + ,134759 + ,518365 + ,0 + ,4658 + ,158 + ,156608 + ,233773 + ,1 + ,814 + ,18 + ,39625 + ,301881 + ,1 + ,1002 + ,15 + ,56750 + ,298568 + ,0 + ,496 + ,17 + ,87390 + ,325479 + ,1 + ,389 + ,14 + ,58990 + ,325506 + ,1 + ,12679 + ,110 + ,48020 + ,984885 + ,1 + ,400 + ,152 + ,60720 + ,313267 + ,1 + ,53 + ,6 + ,60349 + ,315793 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,5109 + ,97 + ,67038 + ,215362 + ,1 + ,576 + ,18 + ,113761 + ,314073 + ,1 + ,438 + ,11 + ,58320 + ,298096 + ,1 + ,165 + ,24 + ,62841 + ,325176 + ,1 + ,4069 + ,145 + ,79804 + ,207393 + ,0 + ,23 + ,5 + ,62555 + ,314806 + ,1 + ,1285 + ,107 + ,99489 + ,341340 + ,0 + ,4677 + ,75 + ,90131 + ,426280 + ,0 + ,24811 + ,505 + ,95350 + ,929118 + ,1 + ,167 + ,2 + ,64245 + ,307322 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,4628 + ,62 + ,69159 + ,387475 + ,1 + ,226 + ,5 + ,65745 + ,291787 + ,0 + ,1765 + ,20 + ,77623 + ,247060 + ,1 + ,460 + ,15 + ,63346 + ,329784 + ,1 + ,36 + ,3 + ,60894 + ,315834 + ,1 + ,989 + ,11 + ,58930 + ,304555 + ,0 + ,3055 + ,43 + ,60247 + ,376641 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,13253 + ,947 + ,90829 + ,357760 + ,1 + ,24 + ,8 + ,59818 + ,315637 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,131 + ,6 + ,59661 + ,327071 + ,1 + ,514 + ,15 + ,102725 + ,377516 + ,1 + ,3366 + ,80 + ,108479 + ,299243 + ,0 + ,9327 + ,101 + ,87419 + ,387699 + ,1 + ,384 + ,22 + ,54683 + ,309836 + ,1 + ,17821 + ,373 + ,122844 + ,444477 + ,1 + ,397 + ,10 + ,62710 + ,322327 + ,1 + ,897 + ,5 + ,60720 + ,314913 + ,1 + ,218 + ,21 + ,60720 + ,315553 + ,0 + ,3369 + ,28 + ,87161 + ,688779 + ,1 + ,5702 + ,61 + ,101481 + ,321896 + ,1 + ,8636 + ,855 + ,128294 + ,301607 + ,1 + ,534 + ,7 + ,62620 + ,304485 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,726 + ,3 + ,69980 + ,231861 + ,1 + ,1380 + ,26 + ,60982 + ,347385 + ,1 + ,180 + ,11 + ,59635 + ,316386 + ,0 + ,7285 + ,115 + ,188873 + ,491303 + ,1 + ,880 + ,74 + ,80791 + ,261216 + ,1 + ,1 + ,2 + ,60727 + ,315388 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,96 + ,5 + ,60379 + ,313729 + ,1 + ,1889 + ,45 + ,37527 + ,358649 + ,0 + ,45187 + ,353 + ,234817 + ,1926517 + ,1 + ,288 + ,4 + ,60510 + ,296656 + ,1 + ,1270 + ,26 + ,69206 + ,275311 + ,1 + ,6526 + ,26 + ,55830 + ,-42143 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,226 + ,12 + ,72835 + ,343929 + ,1 + ,694 + ,8 + ,68060 + ,367655 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,249 + ,72 + ,56726 + ,313491) + ,dim=c(5 + ,431) + ,dimnames=list(c('group' + ,'costs' + ,'trades' + ,'dividends' + ,'wealth') + ,1:431)) > y <- array(NA,dim=c(5,431),dimnames=list(c('group','costs','trades','dividends','wealth'),1:431)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par20 = '' > par19 = '' > par18 = '' > par17 = '' > par16 = '' > par15 = '' > par14 = '' > par13 = '' > par12 = '' > par11 = '' > par10 = '' > par9 = '' > par8 = '' > par7 = '' > par6 = '' > par5 = '' > par4 = '' > par3 = 'No Linear Trend' > par2 = 'Do not include Seasonal Dummies' > par1 = '5' > ylab = '' > xlab = '' > main = '' > #'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 > n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test > par1 <- as.numeric(par1) > x <- t(y) > k <- length(x[1,]) > n <- length(x[,1]) > x1 <- cbind(x[,par1], x[,1:k!=par1]) > mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1]) > colnames(x1) <- mycolnames #colnames(x)[par1] > x <- x1 > if (par3 == 'First Differences'){ + x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep=''))) + for (i in 1:n-1) { + for (j in 1:k) { + x2[i,j] <- x[i+1,j] - x[i,j] + } + } + x <- x2 + } > if (par2 == 'Include Monthly Dummies'){ + x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep =''))) + for (i in 1:11){ + x2[seq(i,n,12),i] <- 1 + } + x <- cbind(x, x2) + } > if (par2 == 'Include Quarterly Dummies'){ + x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep =''))) + for (i in 1:3){ + x2[seq(i,n,4),i] <- 1 + } + x <- cbind(x, x2) + } > k <- length(x[1,]) > if (par3 == 'Linear Trend'){ + x <- cbind(x, c(1:n)) + colnames(x)[k+1] <- 't' + } > x wealth group costs trades dividends 1 319369 1 1595 17 60720 2 493408 0 5565 47 94355 3 319210 0 601 6 60720 4 381180 1 188 11 77655 5 297978 1 7146 105 134028 6 290476 0 1135 17 62285 7 292136 1 450 8 59325 8 314353 1 34 8 60630 9 339445 1 133 14 65990 10 303677 1 119 6 59118 11 397144 0 2053 53 100423 12 424898 0 4036 63 100269 13 315380 1 0 0 60720 14 341570 1 4655 393 60720 15 308989 1 131 11 61808 16 305959 1 1766 73 60438 17 318690 1 312 14 58598 18 323361 1 448 40 59781 19 318903 1 115 13 60945 20 314049 1 60 10 61124 21 315380 1 0 0 60720 22 298466 1 364 35 47705 23 315380 1 0 0 60720 24 438493 0 1442 118 121173 25 378049 0 1389 9 57530 26 313332 1 149 5 62041 27 319864 1 2212 14 60720 28 430866 0 7489 322 136996 29 315380 1 0 0 60720 30 332743 1 402 26 84990 31 286963 0 7419 29 63255 32 315380 1 0 0 60720 33 331955 1 307 15 58856 34 527021 0 1134 162 146216 35 364304 0 7561 87 82425 36 292154 1 5131 71 86111 37 314987 1 9 6 60835 38 272713 1 2264 24 55174 39 1398893 1 40949 481 129352 40 409642 0 4980 91 113521 41 429112 1 272 8 112995 42 382712 1 757 19 45689 43 374943 1 1424 94 131116 44 315380 1 0 0 60720 45 297413 0 378 52 71873 46 304252 1 1061 5 70415 47 315380 1 0 0 60720 48 309596 1 203 9 61938 49 377305 0 3689 25 101792 50 386212 1 567 11 103345 51 657954 0 27664 282 74996 52 368186 1 1686 40 68696 53 269753 0 1164 17 57320 54 585715 1 15824 64 63838 55 1071292 1 6575 159 151352 56 992426 0 37597 248 37884 57 315371 1 314 5 60720 58 249898 0 1932 36 72954 59 357312 0 5692 72 64239 60 315380 1 0 0 60720 61 315380 1 0 0 60720 62 364839 0 6533 75 126910 63 230621 0 2126 138 21001 64 315877 1 14 1 61016 65 125390 0 3192 60 59831 66 314882 1 47 3 60720 67 370837 0 4891 217 210568 68 330068 0 531 10 71641 69 324385 1 2332 4 72680 70 315380 1 0 0 60720 71 1443586 0 16275 109 47818 72 253537 1 9252 48 79510 73 4321023 0 29790 309 81767 74 315380 1 0 0 60720 75 352108 0 2300 29 64873 76 330059 1 267 32 57640 77 340968 0 382 21 59365 78 317736 1 397 28 60646 79 209458 1 11920 182 31080 80 1491348 0 54660 302 126942 81 314887 1 5 2 60720 82 315380 1 0 0 60720 83 315380 1 0 0 60720 84 353058 0 2647 37 62920 85 315380 1 0 0 60720 86 315380 0 0 0 60720 87 314533 1 94 7 60793 88 315354 1 3 2 60698 89 302187 1 454 51 63261 90 336639 1 227 9 76644 91 315380 0 0 0 60720 92 296702 1 206 9 53110 93 1073089 1 25830 115 245546 94 146494 1 1528 38 60326 95 325249 1 271 2 60720 96 331420 1 138 6 69817 97 315380 1 0 0 60720 98 314922 1 278 26 61404 99 320016 1 282 30 62452 100 280398 1 571 80 41477 101 452469 1 2253 99 63593 102 301164 1 290 2 58790 103 317330 1 78 16 62700 104 315576 1 20 2 60805 105 315380 1 0 0 60720 106 -7170 1 18073 180 27284 107 322331 1 1866 29 56225 108 1629616 0 42634 156 157214 109 292754 1 249 62 54323 110 318056 1 422 28 57935 111 355178 1 2675 79 59017 112 204325 1 965 3 73490 113 315380 1 0 0 60720 114 317046 1 621 16 68005 115 315380 0 0 0 60720 116 309560 1 365 8 54820 117 414462 0 3122 34 94670 118 857217 0 12988 119 82340 119 697458 0 5336 81 112477 120 530670 0 3160 108 108094 121 238125 1 2489 40 67804 122 315380 1 0 0 60720 123 741409 0 6984 45 80570 124 393343 1 2001 52 95551 125 372631 1 15236 90 77440 126 315380 0 0 0 60720 127 317291 1 530 60 73433 128 315380 1 0 0 60720 129 306275 1 35624 223 157278 130 317892 0 1383 11 73221 131 334280 1 875 14 67000 132 315380 0 0 0 60720 133 315380 0 0 0 60720 134 315380 0 0 0 60720 135 314210 1 72 3 60398 136 306948 0 265 31 60720 137 315380 1 0 0 60720 138 320398 1 335 33 64175 139 315380 0 0 0 60720 140 315380 1 0 0 60720 141 501749 1 4525 226 93811 142 202055 0 3045 58 27330 143 315380 1 0 0 60720 144 315380 1 0 0 60720 145 333210 0 638 14 60370 146 315380 0 0 0 60720 147 322340 1 607 11 60436 148 369448 1 1558 63 55637 149 291841 0 1324 117 67440 150 315380 1 0 0 60720 151 315380 0 0 0 60720 152 296919 1 611 8 59190 153 315380 1 0 0 60720 154 309038 1 923 8 58620 155 246541 1 661 3 65920 156 289513 0 2397 14 74020 157 344425 1 366 53 61808 158 315380 1 0 0 60720 159 314210 1 135 2 62065 160 480382 0 1659 29 107577 161 315009 1 316 9 60505 162 315380 1 0 0 60720 163 312878 1 309 11 56535 164 322031 0 49 8 64107 165 315380 1 0 0 60720 166 597793 0 4519 49 102129 167 315380 1 0 0 60720 168 315688 1 837 69 61262 169 378525 1 5119 49 39039 170 312378 1 1280 117 69465 171 403560 1 2564 22 130140 172 510834 0 2045 20 97890 173 214215 0 2234 16 77200 174 235133 1 975 23 90534 175 343613 0 1136 32 48522 176 365959 1 453 21 81125 177 315380 1 0 0 60720 178 314551 1 61 14 60720 179 303230 1 368 17 61686 180 315380 0 0 0 60720 181 469107 1 4901 84 121920 182 354228 1 540 14 103960 183 315380 1 0 0 60720 184 315380 1 0 0 60720 185 315380 0 0 0 60720 186 315394 1 2 9 60735 187 312412 1 36 4 61564 188 333505 1 776 9 64230 189 223193 1 84738 588 -26007 190 315380 1 0 0 60720 191 315656 1 3 4 60761 192 296261 1 529 9 63870 193 315380 1 0 0 60720 194 336425 1 405 7 60845 195 359335 1 972 34 71642 196 315380 1 0 0 60720 197 308636 1 2099 64 106611 198 158492 0 3437 47 48022 199 315380 0 0 0 60720 200 315380 1 0 0 60720 201 711969 0 22330 84 79801 202 315380 0 0 0 60720 203 315380 0 0 0 60720 204 306268 1 483 27 60830 205 315380 1 0 0 60720 206 442882 1 2239 21 88590 207 378509 0 2949 41 82903 208 315380 1 0 0 60720 209 346611 1 365 83 87192 210 314289 1 2461 57 55792 211 856956 0 21950 519 114337 212 217193 0 3294 23 75832 213 315366 1 141 6 61630 214 307930 1 572 16 58580 215 702380 1 13326 102 165548 216 194493 1 2284 33 76403 217 316155 0 10 2 61656 218 315380 1 0 0 60720 219 330546 0 1414 198 70184 220 394510 1 1975 35 118881 221 312846 1 43 3 60887 222 315380 1 0 0 60720 223 296139 1 844 78 60925 224 295580 1 304 15 62969 225 297765 1 458 11 58625 226 377934 1 18562 155 102313 227 315380 1 0 0 60720 228 638830 1 7123 109 288170 229 304376 1 622 73 73007 230 307424 1 174 1 64820 231 644190 1 2220 22 301670 232 295370 1 121 24 56178 233 574339 0 11819 85 106113 234 315380 1 0 0 60720 235 310201 1 125 12 60798 236 327007 1 1182 13 70694 237 343466 1 1503 39 56364 238 315380 1 0 0 60720 239 315380 1 0 0 60720 240 318098 1 30 4 62045 241 448243 0 3310 33 75230 242 325738 1 554 11 79285 243 315380 0 0 0 60720 244 312161 1 468 21 60720 245 243650 1 4917 42 52811 246 407159 1 3256 66 35250 247 315380 1 0 0 60720 248 317698 1 125 13 59734 249 315380 1 0 0 60720 250 312502 1 22 2 60722 251 315380 0 0 0 60720 252 322378 1 514 52 78780 253 315380 0 0 0 60720 254 315380 1 0 0 60720 255 315380 0 0 0 60720 256 315380 1 0 0 60720 257 315380 1 0 0 60720 258 315380 1 0 0 60720 259 640273 1 10327 91 88577 260 345783 1 13718 61 96448 261 652925 1 3748 9 50350 262 439798 1 14416 361 49857 263 278990 0 1526 25 69351 264 339836 0 666 171 117869 265 240897 1 2844 67 72683 266 315380 1 0 0 60720 267 297141 1 368 17 61167 268 315380 1 0 0 60720 269 331323 1 333 13 70811 270 313880 1 26 1 60896 271 315380 1 0 0 60720 272 315380 1 0 0 60720 273 309422 1 1303 61 69863 274 315245 1 20 3 60938 275 405972 1 2384 97 61348 276 300962 0 203 27 50804 277 316176 0 71 1 60745 278 302409 1 53 14 59506 279 283587 1 562 13 58456 280 263276 1 622 7 60950 281 312075 1 645 5 60720 282 308336 0 1763 24 61600 283 315380 1 0 0 60720 284 298700 1 317 9 63915 285 315372 1 1 6 60719 286 318745 1 275 4 59500 287 315380 1 0 0 60720 288 408881 1 936 30 67939 289 786690 1 8568 37 32168 290 -83265 0 11528 70 -14545 291 315380 0 0 0 60720 292 321376 1 738 10 60720 293 315380 1 0 0 60720 294 276898 1 592 7 64270 295 315547 1 126 10 60951 296 315487 1 2 1 60743 297 315380 1 0 0 60720 298 1405225 0 18014 239 -1710 299 315380 0 0 0 60720 300 315380 0 0 0 60720 301 983660 0 37704 452 60448 302 318574 0 63 4 65688 303 310768 1 1431 34 106885 304 312887 1 94 3 61360 305 312339 1 192 4 65276 306 314964 1 32 1 59988 307 379983 1 6869 105 117520 308 315380 1 0 0 60720 309 315380 1 0 0 60720 310 315398 1 1 4 60722 311 315380 0 0 0 60720 312 253588 1 2328 94 82732 313 316647 0 209 13 64016 314 315688 1 28 2 60890 315 310670 1 176 28 68136 316 165404 1 1920 15 79420 317 4111912 0 87550 458 153198 318 291650 1 520 10 58650 319 315380 0 0 0 60720 320 253468 1 1013 2 69770 321 315688 1 15 8 60831 322 325699 1 587 12 59595 323 446211 1 5371 66 87720 324 315380 1 0 0 60720 325 368078 1 1012 47 114768 326 315380 1 0 0 60720 327 352850 1 876 59 138971 328 6282154 0 162556 1081 213118 329 227132 1 43556 247 32648 330 283910 0 3425 12 83620 331 236761 1 810 43 74015 332 315380 1 0 0 60720 333 315380 0 0 0 60720 334 550608 0 2365 79 191778 335 307528 0 1261 33 76114 336 315380 1 0 0 60720 337 355864 1 1585 107 93099 338 358589 1 16189 295 116384 339 315380 1 0 0 60720 340 315380 0 0 0 60720 341 375195 1 10579 92 110309 342 288985 1 474 20 61977 343 315380 1 0 0 60720 344 458343 1 3642 86 90262 345 315380 1 0 0 60720 346 269587 0 472 20 80045 347 315236 1 98 2 61490 348 42754 1 3999 30 51252 349 315380 1 0 0 60720 350 308256 1 621 48 60720 351 315380 1 0 0 60720 352 315380 1 0 0 60720 353 313164 0 30 4 60798 354 315380 1 0 0 60720 355 269661 1 746 23 71561 356 315380 0 0 0 60720 357 315380 1 0 0 60720 358 315380 0 0 0 60720 359 315380 1 0 0 60720 360 518365 0 4150 150 134759 361 233773 0 4658 158 156608 362 301881 1 814 18 39625 363 298568 1 1002 15 56750 364 325479 0 496 17 87390 365 325506 1 389 14 58990 366 984885 1 12679 110 48020 367 313267 1 400 152 60720 368 315793 1 53 6 60349 369 315380 1 0 0 60720 370 215362 0 5109 97 67038 371 314073 1 576 18 113761 372 298096 1 438 11 58320 373 325176 1 165 24 62841 374 207393 1 4069 145 79804 375 314806 0 23 5 62555 376 341340 1 1285 107 99489 377 426280 0 4677 75 90131 378 929118 0 24811 505 95350 379 307322 1 167 2 64245 380 315380 1 0 0 60720 381 387475 1 4628 62 69159 382 291787 1 226 5 65745 383 247060 0 1765 20 77623 384 329784 1 460 15 63346 385 315834 1 36 3 60894 386 304555 1 989 11 58930 387 376641 0 3055 43 60247 388 315380 0 0 0 60720 389 315380 1 0 0 60720 390 315380 0 0 0 60720 391 315380 1 0 0 60720 392 315380 1 0 0 60720 393 357760 1 13253 947 90829 394 315637 1 24 8 59818 395 315380 1 0 0 60720 396 315380 0 0 0 60720 397 315380 1 0 0 60720 398 315380 1 0 0 60720 399 327071 1 131 6 59661 400 377516 1 514 15 102725 401 299243 1 3366 80 108479 402 387699 0 9327 101 87419 403 309836 1 384 22 54683 404 444477 1 17821 373 122844 405 322327 1 397 10 62710 406 314913 1 897 5 60720 407 315553 1 218 21 60720 408 688779 0 3369 28 87161 409 321896 1 5702 61 101481 410 301607 1 8636 855 128294 411 304485 1 534 7 62620 412 315380 1 0 0 60720 413 315380 0 0 0 60720 414 231861 1 726 3 69980 415 347385 1 1380 26 60982 416 316386 1 180 11 59635 417 491303 0 7285 115 188873 418 261216 1 880 74 80791 419 315388 1 1 2 60727 420 315380 1 0 0 60720 421 313729 1 96 5 60379 422 358649 1 1889 45 37527 423 1926517 0 45187 353 234817 424 296656 1 288 4 60510 425 275311 1 1270 26 69206 426 -42143 1 6526 26 55830 427 315380 0 0 0 60720 428 343929 1 226 12 72835 429 367655 1 694 8 68060 430 315380 1 0 0 60720 431 313491 1 249 72 56726 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) group costs trades dividends 205978.971 -55485.819 29.811 -407.150 1.926 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -2163938 -14344 35201 47931 3195284 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 2.060e+05 3.854e+04 5.345 1.48e-07 *** group -5.549e+04 2.784e+04 -1.993 0.0469 * costs 2.981e+01 1.578e+00 18.889 < 2e-16 *** trades -4.072e+02 1.663e+02 -2.448 0.0148 * dividends 1.926e+00 4.105e-01 4.692 3.65e-06 *** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 253400 on 426 degrees of freedom Multiple R-squared: 0.6569, Adjusted R-squared: 0.6537 F-statistic: 203.9 on 4 and 426 DF, p-value: < 2.2e-16 > if (n > n25) { + kp3 <- k + 3 + nmkm3 <- n - k - 3 + gqarr <- array(NA, dim=c(nmkm3-kp3+1,3)) + numgqtests <- 0 + numsignificant1 <- 0 + numsignificant5 <- 0 + numsignificant10 <- 0 + for (mypoint in kp3:nmkm3) { + j <- 0 + numgqtests <- numgqtests + 1 + for (myalt in c('greater', 'two.sided', 'less')) { + j <- j + 1 + gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value + } + if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1 + } + gqarr + } [,1] [,2] [,3] [1,] 8.484368e-04 1.696874e-03 9.991516e-01 [2,] 5.665972e-04 1.133194e-03 9.994334e-01 [3,] 7.190595e-05 1.438119e-04 9.999281e-01 [4,] 9.710950e-05 1.942190e-04 9.999029e-01 [5,] 4.832291e-05 9.664581e-05 9.999517e-01 [6,] 9.861619e-06 1.972324e-05 9.999901e-01 [7,] 1.611924e-05 3.223847e-05 9.999839e-01 [8,] 3.240916e-06 6.481832e-06 9.999968e-01 [9,] 6.401563e-07 1.280313e-06 9.999994e-01 [10,] 1.197364e-07 2.394728e-07 9.999999e-01 [11,] 2.181095e-08 4.362190e-08 1.000000e+00 [12,] 3.780405e-09 7.560810e-09 1.000000e+00 [13,] 6.273737e-10 1.254747e-09 1.000000e+00 [14,] 1.010434e-10 2.020868e-10 1.000000e+00 [15,] 1.615546e-11 3.231092e-11 1.000000e+00 [16,] 2.459365e-12 4.918730e-12 1.000000e+00 [17,] 4.223568e-13 8.447135e-13 1.000000e+00 [18,] 6.406709e-14 1.281342e-13 1.000000e+00 [19,] 9.077889e-15 1.815578e-14 1.000000e+00 [20,] 1.279272e-15 2.558545e-15 1.000000e+00 [21,] 1.855800e-16 3.711600e-16 1.000000e+00 [22,] 2.469671e-17 4.939342e-17 1.000000e+00 [23,] 3.209765e-18 6.419529e-18 1.000000e+00 [24,] 1.583470e-18 3.166941e-18 1.000000e+00 [25,] 2.086646e-19 4.173293e-19 1.000000e+00 [26,] 3.129255e-20 6.258509e-20 1.000000e+00 [27,] 7.759474e-21 1.551895e-20 1.000000e+00 [28,] 1.146319e-21 2.292638e-21 1.000000e+00 [29,] 1.806038e-22 3.612076e-22 1.000000e+00 [30,] 2.268877e-23 4.537755e-23 1.000000e+00 [31,] 3.258163e-24 6.516326e-24 1.000000e+00 [32,] 1.981634e-13 3.963269e-13 1.000000e+00 [33,] 6.034170e-14 1.206834e-13 1.000000e+00 [34,] 2.678091e-14 5.356181e-14 1.000000e+00 [35,] 1.362866e-14 2.725733e-14 1.000000e+00 [36,] 3.854280e-15 7.708560e-15 1.000000e+00 [37,] 1.057449e-15 2.114899e-15 1.000000e+00 [38,] 2.799644e-16 5.599288e-16 1.000000e+00 [39,] 7.883505e-17 1.576701e-16 1.000000e+00 [40,] 2.089755e-17 4.179509e-17 1.000000e+00 [41,] 5.343678e-18 1.068736e-17 1.000000e+00 [42,] 1.465547e-18 2.931094e-18 1.000000e+00 [43,] 4.083408e-19 8.166816e-19 1.000000e+00 [44,] 5.765247e-18 1.153049e-17 1.000000e+00 [45,] 1.622127e-18 3.244253e-18 1.000000e+00 [46,] 4.383994e-19 8.767988e-19 1.000000e+00 [47,] 1.250605e-19 2.501209e-19 1.000000e+00 [48,] 5.163040e-14 1.032608e-13 1.000000e+00 [49,] 2.802594e-14 5.605188e-14 1.000000e+00 [50,] 9.907984e-15 1.981597e-14 1.000000e+00 [51,] 4.019064e-15 8.038128e-15 1.000000e+00 [52,] 1.414009e-15 2.828018e-15 1.000000e+00 [53,] 4.893971e-16 9.787943e-16 1.000000e+00 [54,] 1.676641e-16 3.353283e-16 1.000000e+00 [55,] 1.897401e-16 3.794802e-16 1.000000e+00 [56,] 7.157465e-17 1.431493e-16 1.000000e+00 [57,] 2.439192e-17 4.878383e-17 1.000000e+00 [58,] 3.525940e-17 7.051880e-17 1.000000e+00 [59,] 1.211744e-17 2.423487e-17 1.000000e+00 [60,] 7.198683e-17 1.439737e-16 1.000000e+00 [61,] 3.299063e-17 6.598125e-17 1.000000e+00 [62,] 1.320711e-17 2.641422e-17 1.000000e+00 [63,] 4.637037e-18 9.274073e-18 1.000000e+00 [64,] 4.366882e-09 8.733765e-09 1.000000e+00 [65,] 1.001800e-08 2.003599e-08 1.000000e+00 [66,] 1.000000e+00 9.921018e-16 4.960509e-16 [67,] 1.000000e+00 1.895110e-15 9.475552e-16 [68,] 1.000000e+00 3.642158e-15 1.821079e-15 [69,] 1.000000e+00 6.757272e-15 3.378636e-15 [70,] 1.000000e+00 1.271375e-14 6.356875e-15 [71,] 1.000000e+00 2.362851e-14 1.181425e-14 [72,] 1.000000e+00 9.671149e-15 4.835575e-15 [73,] 1.000000e+00 8.345993e-16 4.172996e-16 [74,] 1.000000e+00 1.571996e-15 7.859979e-16 [75,] 1.000000e+00 2.937336e-15 1.468668e-15 [76,] 1.000000e+00 5.448409e-15 2.724204e-15 [77,] 1.000000e+00 1.008237e-14 5.041186e-15 [78,] 1.000000e+00 1.842538e-14 9.212692e-15 [79,] 1.000000e+00 3.396701e-14 1.698350e-14 [80,] 1.000000e+00 6.133035e-14 3.066518e-14 [81,] 1.000000e+00 1.096820e-13 5.484099e-14 [82,] 1.000000e+00 1.963150e-13 9.815749e-14 [83,] 1.000000e+00 3.462971e-13 1.731485e-13 [84,] 1.000000e+00 6.153513e-13 3.076757e-13 [85,] 1.000000e+00 1.077475e-12 5.387375e-13 [86,] 1.000000e+00 1.313099e-12 6.565497e-13 [87,] 1.000000e+00 1.836230e-12 9.181151e-13 [88,] 1.000000e+00 3.136701e-12 1.568351e-12 [89,] 1.000000e+00 5.324004e-12 2.662002e-12 [90,] 1.000000e+00 8.989909e-12 4.494954e-12 [91,] 1.000000e+00 1.514019e-11 7.570096e-12 [92,] 1.000000e+00 2.527376e-11 1.263688e-11 [93,] 1.000000e+00 4.088254e-11 2.044127e-11 [94,] 1.000000e+00 6.143946e-11 3.071973e-11 [95,] 1.000000e+00 1.015079e-10 5.075395e-11 [96,] 1.000000e+00 1.651878e-10 8.259389e-11 [97,] 1.000000e+00 2.669176e-10 1.334588e-10 [98,] 1.000000e+00 4.285435e-10 2.142717e-10 [99,] 1.000000e+00 1.176614e-11 5.883070e-12 [100,] 1.000000e+00 1.969451e-11 9.847255e-12 [101,] 1.000000e+00 3.119349e-11 1.559675e-11 [102,] 1.000000e+00 5.087096e-11 2.543548e-11 [103,] 1.000000e+00 8.247451e-11 4.123725e-11 [104,] 1.000000e+00 1.334495e-10 6.672476e-11 [105,] 1.000000e+00 2.006704e-10 1.003352e-10 [106,] 1.000000e+00 3.195677e-10 1.597839e-10 [107,] 1.000000e+00 5.123797e-10 2.561899e-10 [108,] 1.000000e+00 8.196487e-10 4.098244e-10 [109,] 1.000000e+00 1.287450e-09 6.437252e-10 [110,] 1.000000e+00 2.010781e-09 1.005391e-09 [111,] 1.000000e+00 2.762543e-09 1.381271e-09 [112,] 1.000000e+00 3.681272e-09 1.840636e-09 [113,] 1.000000e+00 5.603798e-09 2.801899e-09 [114,] 1.000000e+00 8.074612e-09 4.037306e-09 [115,] 1.000000e+00 1.225828e-08 6.129138e-09 [116,] 1.000000e+00 1.486480e-08 7.432400e-09 [117,] 1.000000e+00 2.263544e-08 1.131772e-08 [118,] 1.000000e+00 1.652429e-08 8.262144e-09 [119,] 1.000000e+00 2.516410e-08 1.258205e-08 [120,] 1.000000e+00 3.794774e-08 1.897387e-08 [121,] 1.000000e+00 5.632909e-08 2.816454e-08 [122,] 1.000000e+00 3.982383e-11 1.991191e-11 [123,] 1.000000e+00 6.242953e-11 3.121476e-11 [124,] 1.000000e+00 9.976004e-11 4.988002e-11 [125,] 1.000000e+00 1.591813e-10 7.959063e-11 [126,] 1.000000e+00 2.527101e-10 1.263551e-10 [127,] 1.000000e+00 3.991536e-10 1.995768e-10 [128,] 1.000000e+00 6.204676e-10 3.102338e-10 [129,] 1.000000e+00 9.594216e-10 4.797108e-10 [130,] 1.000000e+00 1.473770e-09 7.368852e-10 [131,] 1.000000e+00 2.260209e-09 1.130105e-09 [132,] 1.000000e+00 3.479437e-09 1.739719e-09 [133,] 1.000000e+00 5.264166e-09 2.632083e-09 [134,] 1.000000e+00 7.481755e-09 3.740877e-09 [135,] 1.000000e+00 9.381567e-09 4.690784e-09 [136,] 1.000000e+00 1.401844e-08 7.009222e-09 [137,] 1.000000e+00 2.084636e-08 1.042318e-08 [138,] 1.000000e+00 3.120861e-08 1.560431e-08 [139,] 1.000000e+00 4.654986e-08 2.327493e-08 [140,] 1.000000e+00 6.851278e-08 3.425639e-08 [141,] 1.000000e+00 9.736479e-08 4.868239e-08 [142,] 9.999999e-01 1.355086e-07 6.775431e-08 [143,] 9.999999e-01 1.958388e-07 9.791941e-08 [144,] 9.999999e-01 2.852404e-07 1.426202e-07 [145,] 9.999998e-01 4.133511e-07 2.066756e-07 [146,] 9.999997e-01 5.887539e-07 2.943769e-07 [147,] 9.999996e-01 8.439698e-07 4.219849e-07 [148,] 9.999994e-01 1.193783e-06 5.968914e-07 [149,] 9.999992e-01 1.576432e-06 7.882159e-07 [150,] 9.999989e-01 2.163169e-06 1.081585e-06 [151,] 9.999985e-01 3.012143e-06 1.506072e-06 [152,] 9.999979e-01 4.189007e-06 2.094504e-06 [153,] 9.999971e-01 5.770129e-06 2.885064e-06 [154,] 9.999960e-01 7.952795e-06 3.976398e-06 [155,] 9.999946e-01 1.087242e-05 5.436209e-06 [156,] 9.999926e-01 1.480999e-05 7.404995e-06 [157,] 9.999899e-01 2.027990e-05 1.013995e-05 [158,] 9.999863e-01 2.734521e-05 1.367261e-05 [159,] 9.999822e-01 3.559695e-05 1.779848e-05 [160,] 9.999762e-01 4.759168e-05 2.379584e-05 [161,] 9.999683e-01 6.331389e-05 3.165695e-05 [162,] 9.999578e-01 8.433054e-05 4.216527e-05 [163,] 9.999445e-01 1.109703e-04 5.548515e-05 [164,] 9.999285e-01 1.429787e-04 7.148937e-05 [165,] 9.999080e-01 1.840707e-04 9.203536e-05 [166,] 9.998970e-01 2.059166e-04 1.029583e-04 [167,] 9.998708e-01 2.584229e-04 1.292114e-04 [168,] 9.998323e-01 3.353787e-04 1.676894e-04 [169,] 9.997847e-01 4.306850e-04 2.153425e-04 [170,] 9.997235e-01 5.529107e-04 2.764554e-04 [171,] 9.996467e-01 7.066665e-04 3.533332e-04 [172,] 9.995475e-01 9.050856e-04 4.525428e-04 [173,] 9.994216e-01 1.156865e-03 5.784324e-04 [174,] 9.992685e-01 1.462979e-03 7.314896e-04 [175,] 9.990778e-01 1.844470e-03 9.222352e-04 [176,] 9.988463e-01 2.307397e-03 1.153699e-03 [177,] 9.985628e-01 2.874349e-03 1.437175e-03 [178,] 9.982033e-01 3.593447e-03 1.796723e-03 [179,] 9.977823e-01 4.435495e-03 2.217748e-03 [180,] 9.972682e-01 5.463658e-03 2.731829e-03 [181,] 9.966456e-01 6.708709e-03 3.354354e-03 [182,] 1.000000e+00 6.486876e-17 3.243438e-17 [183,] 1.000000e+00 1.145452e-16 5.727261e-17 [184,] 1.000000e+00 2.010375e-16 1.005187e-16 [185,] 1.000000e+00 3.590393e-16 1.795196e-16 [186,] 1.000000e+00 6.254204e-16 3.127102e-16 [187,] 1.000000e+00 1.071454e-15 5.357270e-16 [188,] 1.000000e+00 1.827067e-15 9.135334e-16 [189,] 1.000000e+00 3.140341e-15 1.570170e-15 [190,] 1.000000e+00 5.086919e-15 2.543460e-15 [191,] 1.000000e+00 5.017429e-15 2.508715e-15 [192,] 1.000000e+00 8.608374e-15 4.304187e-15 [193,] 1.000000e+00 1.461258e-14 7.306292e-15 [194,] 1.000000e+00 7.335428e-15 3.667714e-15 [195,] 1.000000e+00 1.258696e-14 6.293480e-15 [196,] 1.000000e+00 2.150346e-14 1.075173e-14 [197,] 1.000000e+00 3.663851e-14 1.831926e-14 [198,] 1.000000e+00 6.147306e-14 3.073653e-14 [199,] 1.000000e+00 9.945364e-14 4.972682e-14 [200,] 1.000000e+00 1.639913e-13 8.199563e-14 [201,] 1.000000e+00 2.716200e-13 1.358100e-13 [202,] 1.000000e+00 4.312368e-13 2.156184e-13 [203,] 1.000000e+00 7.215039e-13 3.607519e-13 [204,] 1.000000e+00 1.186368e-12 5.931840e-13 [205,] 1.000000e+00 1.220945e-12 6.104724e-13 [206,] 1.000000e+00 1.993336e-12 9.966680e-13 [207,] 1.000000e+00 3.268640e-12 1.634320e-12 [208,] 1.000000e+00 4.798555e-12 2.399278e-12 [209,] 1.000000e+00 6.016286e-12 3.008143e-12 [210,] 1.000000e+00 9.757872e-12 4.878936e-12 [211,] 1.000000e+00 1.560175e-11 7.800875e-12 [212,] 1.000000e+00 2.280856e-11 1.140428e-11 [213,] 1.000000e+00 3.648137e-11 1.824068e-11 [214,] 1.000000e+00 5.775923e-11 2.887961e-11 [215,] 1.000000e+00 9.083114e-11 4.541557e-11 [216,] 1.000000e+00 1.429135e-10 7.145675e-11 [217,] 1.000000e+00 2.260035e-10 1.130018e-10 [218,] 1.000000e+00 3.556400e-10 1.778200e-10 [219,] 1.000000e+00 5.220834e-11 2.610417e-11 [220,] 1.000000e+00 8.242160e-11 4.121080e-11 [221,] 1.000000e+00 1.002645e-10 5.013225e-11 [222,] 1.000000e+00 1.580578e-10 7.902889e-11 [223,] 1.000000e+00 2.498562e-10 1.249281e-10 [224,] 1.000000e+00 3.737215e-10 1.868607e-10 [225,] 1.000000e+00 5.816809e-10 2.908404e-10 [226,] 1.000000e+00 7.356029e-10 3.678015e-10 [227,] 1.000000e+00 1.132015e-09 5.660074e-10 [228,] 1.000000e+00 1.739411e-09 8.697056e-10 [229,] 1.000000e+00 2.694508e-09 1.347254e-09 [230,] 1.000000e+00 4.081513e-09 2.040756e-09 [231,] 1.000000e+00 6.169855e-09 3.084927e-09 [232,] 1.000000e+00 9.285235e-09 4.642618e-09 [233,] 1.000000e+00 1.389478e-08 6.947392e-09 [234,] 1.000000e+00 2.091167e-08 1.045583e-08 [235,] 1.000000e+00 3.148867e-08 1.574433e-08 [236,] 1.000000e+00 4.686588e-08 2.343294e-08 [237,] 1.000000e+00 6.936867e-08 3.468434e-08 [238,] 1.000000e+00 8.318332e-08 4.159166e-08 [239,] 9.999999e-01 1.143899e-07 5.719495e-08 [240,] 9.999999e-01 1.667603e-07 8.338016e-08 [241,] 9.999999e-01 2.410357e-07 1.205179e-07 [242,] 9.999998e-01 3.482596e-07 1.741298e-07 [243,] 9.999997e-01 5.018126e-07 2.509063e-07 [244,] 9.999996e-01 7.232469e-07 3.616235e-07 [245,] 9.999995e-01 1.035168e-06 5.175841e-07 [246,] 9.999993e-01 1.478584e-06 7.392921e-07 [247,] 9.999990e-01 2.088956e-06 1.044478e-06 [248,] 9.999985e-01 2.956987e-06 1.478493e-06 [249,] 9.999979e-01 4.140453e-06 2.070227e-06 [250,] 9.999971e-01 5.771556e-06 2.885778e-06 [251,] 9.999960e-01 8.009095e-06 4.004547e-06 [252,] 9.999946e-01 1.081374e-05 5.406868e-06 [253,] 9.999977e-01 4.615971e-06 2.307985e-06 [254,] 9.999982e-01 3.631730e-06 1.815865e-06 [255,] 9.999977e-01 4.605737e-06 2.302868e-06 [256,] 9.999969e-01 6.203360e-06 3.101680e-06 [257,] 9.999959e-01 8.165397e-06 4.082698e-06 [258,] 9.999948e-01 1.046507e-05 5.232537e-06 [259,] 9.999928e-01 1.447155e-05 7.235777e-06 [260,] 9.999899e-01 2.011157e-05 1.005578e-05 [261,] 9.999862e-01 2.755454e-05 1.377727e-05 [262,] 9.999812e-01 3.765978e-05 1.882989e-05 [263,] 9.999744e-01 5.118334e-05 2.559167e-05 [264,] 9.999654e-01 6.916072e-05 3.458036e-05 [265,] 9.999535e-01 9.302223e-05 4.651111e-05 [266,] 9.999371e-01 1.257906e-04 6.289529e-05 [267,] 9.999162e-01 1.675937e-04 8.379683e-05 [268,] 9.998942e-01 2.115774e-04 1.057887e-04 [269,] 9.998600e-01 2.800421e-04 1.400211e-04 [270,] 9.998149e-01 3.702998e-04 1.851499e-04 [271,] 9.997569e-01 4.861448e-04 2.430724e-04 [272,] 9.996796e-01 6.407712e-04 3.203856e-04 [273,] 9.995811e-01 8.377625e-04 4.188813e-04 [274,] 9.994546e-01 1.090767e-03 5.453835e-04 [275,] 9.992996e-01 1.400709e-03 7.003545e-04 [276,] 9.991039e-01 1.792185e-03 8.960926e-04 [277,] 9.988486e-01 2.302798e-03 1.151399e-03 [278,] 9.985426e-01 2.914834e-03 1.457417e-03 [279,] 9.981595e-01 3.681073e-03 1.840536e-03 [280,] 9.976886e-01 4.622843e-03 2.311422e-03 [281,] 9.972482e-01 5.503583e-03 2.751792e-03 [282,] 9.978107e-01 4.378580e-03 2.189290e-03 [283,] 9.996868e-01 6.264513e-04 3.132256e-04 [284,] 9.995868e-01 8.264326e-04 4.132163e-04 [285,] 9.994602e-01 1.079672e-03 5.398359e-04 [286,] 9.993029e-01 1.394196e-03 6.970981e-04 [287,] 9.990934e-01 1.813229e-03 9.066147e-04 [288,] 9.988408e-01 2.318448e-03 1.159224e-03 [289,] 9.985249e-01 2.950238e-03 1.475119e-03 [290,] 9.981314e-01 3.737125e-03 1.868562e-03 [291,] 9.998542e-01 2.915175e-04 1.457587e-04 [292,] 9.998048e-01 3.904402e-04 1.952201e-04 [293,] 9.997399e-01 5.201471e-04 2.600735e-04 [294,] 9.998766e-01 2.468940e-04 1.234470e-04 [295,] 9.998331e-01 3.338266e-04 1.669133e-04 [296,] 9.997783e-01 4.433009e-04 2.216504e-04 [297,] 9.997064e-01 5.871813e-04 2.935906e-04 [298,] 9.996113e-01 7.774578e-04 3.887289e-04 [299,] 9.994920e-01 1.016099e-03 5.080495e-04 [300,] 9.994095e-01 1.181021e-03 5.905106e-04 [301,] 9.992334e-01 1.533102e-03 7.665509e-04 [302,] 9.990100e-01 1.980046e-03 9.900230e-04 [303,] 9.987295e-01 2.541049e-03 1.270525e-03 [304,] 9.983548e-01 3.290451e-03 1.645225e-03 [305,] 9.979449e-01 4.110107e-03 2.055053e-03 [306,] 9.973663e-01 5.267438e-03 2.633719e-03 [307,] 9.966800e-01 6.639974e-03 3.319987e-03 [308,] 9.958227e-01 8.354700e-03 4.177350e-03 [309,] 9.954982e-01 9.003552e-03 4.501776e-03 [310,] 9.999903e-01 1.934470e-05 9.672349e-06 [311,] 9.999860e-01 2.800932e-05 1.400466e-05 [312,] 9.999799e-01 4.010688e-05 2.005344e-05 [313,] 9.999718e-01 5.641322e-05 2.820661e-05 [314,] 9.999607e-01 7.858445e-05 3.929222e-05 [315,] 9.999453e-01 1.093790e-04 5.468950e-05 [316,] 9.999227e-01 1.545355e-04 7.726775e-05 [317,] 9.998938e-01 2.123160e-04 1.061580e-04 [318,] 9.998522e-01 2.956329e-04 1.478164e-04 [319,] 9.997992e-01 4.015649e-04 2.007825e-04 [320,] 9.997250e-01 5.499780e-04 2.749890e-04 [321,] 1.000000e+00 9.631258e-13 4.815629e-13 [322,] 1.000000e+00 1.299596e-20 6.497978e-21 [323,] 1.000000e+00 1.726870e-20 8.634350e-21 [324,] 1.000000e+00 3.829884e-20 1.914942e-20 [325,] 1.000000e+00 9.285787e-20 4.642894e-20 [326,] 1.000000e+00 2.312863e-19 1.156431e-19 [327,] 1.000000e+00 2.987356e-19 1.493678e-19 [328,] 1.000000e+00 7.117202e-19 3.558601e-19 [329,] 1.000000e+00 1.697864e-18 8.489321e-19 [330,] 1.000000e+00 3.749497e-18 1.874748e-18 [331,] 1.000000e+00 2.041115e-19 1.020557e-19 [332,] 1.000000e+00 5.034865e-19 2.517433e-19 [333,] 1.000000e+00 1.269140e-18 6.345702e-19 [334,] 1.000000e+00 4.378779e-19 2.189390e-19 [335,] 1.000000e+00 1.123281e-18 5.616406e-19 [336,] 1.000000e+00 2.783901e-18 1.391951e-18 [337,] 1.000000e+00 5.795603e-18 2.897802e-18 [338,] 1.000000e+00 1.417015e-17 7.085075e-18 [339,] 1.000000e+00 3.297251e-17 1.648626e-17 [340,] 1.000000e+00 7.999336e-17 3.999668e-17 [341,] 1.000000e+00 1.736588e-18 8.682939e-19 [342,] 1.000000e+00 4.473357e-18 2.236679e-18 [343,] 1.000000e+00 1.157456e-17 5.787279e-18 [344,] 1.000000e+00 2.935581e-17 1.467791e-17 [345,] 1.000000e+00 7.386488e-17 3.693244e-17 [346,] 1.000000e+00 1.879568e-16 9.397839e-17 [347,] 1.000000e+00 4.655715e-16 2.327857e-16 [348,] 1.000000e+00 1.118305e-15 5.591527e-16 [349,] 1.000000e+00 2.774310e-15 1.387155e-15 [350,] 1.000000e+00 6.722435e-15 3.361217e-15 [351,] 1.000000e+00 1.639046e-14 8.195231e-15 [352,] 1.000000e+00 3.906479e-14 1.953239e-14 [353,] 1.000000e+00 5.036680e-14 2.518340e-14 [354,] 1.000000e+00 4.458427e-14 2.229214e-14 [355,] 1.000000e+00 1.088889e-13 5.444446e-14 [356,] 1.000000e+00 2.622499e-13 1.311249e-13 [357,] 1.000000e+00 6.275484e-13 3.137742e-13 [358,] 1.000000e+00 1.453913e-12 7.269567e-13 [359,] 1.000000e+00 5.679096e-14 2.839548e-14 [360,] 1.000000e+00 1.118570e-13 5.592850e-14 [361,] 1.000000e+00 2.688117e-13 1.344058e-13 [362,] 1.000000e+00 6.425918e-13 3.212959e-13 [363,] 1.000000e+00 3.013306e-13 1.506653e-13 [364,] 1.000000e+00 7.514672e-13 3.757336e-13 [365,] 1.000000e+00 1.890680e-12 9.453402e-13 [366,] 1.000000e+00 4.330187e-12 2.165093e-12 [367,] 1.000000e+00 4.670571e-12 2.335285e-12 [368,] 1.000000e+00 1.171694e-11 5.858472e-12 [369,] 1.000000e+00 2.832921e-11 1.416460e-11 [370,] 1.000000e+00 6.706326e-11 3.353163e-11 [371,] 1.000000e+00 1.551195e-10 7.755977e-11 [372,] 1.000000e+00 3.685980e-10 1.842990e-10 [373,] 1.000000e+00 8.404425e-10 4.202213e-10 [374,] 1.000000e+00 1.972418e-09 9.862089e-10 [375,] 1.000000e+00 4.590103e-09 2.295052e-09 [376,] 1.000000e+00 5.883611e-09 2.941806e-09 [377,] 1.000000e+00 1.291454e-08 6.457270e-09 [378,] 1.000000e+00 2.823164e-08 1.411582e-08 [379,] 1.000000e+00 6.367452e-08 3.183726e-08 [380,] 9.999999e-01 1.386605e-07 6.933026e-08 [381,] 9.999998e-01 3.009328e-07 1.504664e-07 [382,] 9.999997e-01 6.274968e-07 3.137484e-07 [383,] 9.999993e-01 1.321403e-06 6.607016e-07 [384,] 9.999987e-01 2.684638e-06 1.342319e-06 [385,] 9.999973e-01 5.375235e-06 2.687617e-06 [386,] 9.999948e-01 1.038571e-05 5.192855e-06 [387,] 9.999899e-01 2.011283e-05 1.005642e-05 [388,] 9.999807e-01 3.866451e-05 1.933225e-05 [389,] 9.999625e-01 7.508213e-05 3.754107e-05 [390,] 9.999299e-01 1.401252e-04 7.006262e-05 [391,] 9.998714e-01 2.571669e-04 1.285835e-04 [392,] 9.997735e-01 4.530684e-04 2.265342e-04 [393,] 9.996045e-01 7.909579e-04 3.954789e-04 [394,] 9.993499e-01 1.300221e-03 6.501103e-04 [395,] 9.995763e-01 8.473422e-04 4.236711e-04 [396,] 9.992353e-01 1.529455e-03 7.647274e-04 [397,] 9.996841e-01 6.318893e-04 3.159447e-04 [398,] 9.994032e-01 1.193665e-03 5.968325e-04 [399,] 9.988576e-01 2.284762e-03 1.142381e-03 [400,] 9.979459e-01 4.108281e-03 2.054141e-03 [401,] 9.983773e-01 3.245341e-03 1.622670e-03 [402,] 9.978334e-01 4.333206e-03 2.166603e-03 [403,] 9.969624e-01 6.075206e-03 3.037603e-03 [404,] 9.943040e-01 1.139194e-02 5.695970e-03 [405,] 9.902485e-01 1.950303e-02 9.751513e-03 [406,] 9.824527e-01 3.509458e-02 1.754729e-02 [407,] 9.695200e-01 6.096005e-02 3.048003e-02 [408,] 9.491105e-01 1.017791e-01 5.088953e-02 [409,] 9.195370e-01 1.609259e-01 8.046297e-02 [410,] 9.640072e-01 7.198563e-02 3.599282e-02 [411,] 9.844260e-01 3.114798e-02 1.557399e-02 [412,] 9.682975e-01 6.340502e-02 3.170251e-02 [413,] 9.400735e-01 1.198531e-01 5.992655e-02 [414,] 8.891404e-01 2.217192e-01 1.108596e-01 [415,] 9.810401e-01 3.791974e-02 1.895987e-02 [416,] 9.967602e-01 6.479602e-03 3.239801e-03 > postscript(file="/var/www/rcomp/tmp/1kz5g1293276484.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/rcomp/tmp/2c84j1293276484.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/rcomp/tmp/3c84j1293276484.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/rcomp/tmp/4c84j1293276484.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/rcomp/tmp/5n0m41293276484.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > qqnorm(mysum$resid, main='Residual Normal Q-Q Plot') > qqline(mysum$resid) > grid() > dev.off() null device 1 > (myerror <- as.ts(mysum$resid)) Time Series: Start = 1 End = 431 Frequency = 1 1 2 3 4 5 6 11292.552 -41076.166 -19198.651 79985.494 -280953.361 -62387.577 7 8 9 10 11 12 17215.934 49320.735 63580.147 38208.833 -41888.465 -68881.815 13 14 15 16 17 18 47930.758 95359.943 40017.496 16128.433 51727.086 60651.036 19 20 21 22 23 24 52885.045 48104.426 47930.758 59484.659 47930.758 4172.197 25 26 27 28 29 30 23515.058 40932.200 -7827.369 -131141.638 47930.758 17147.822 31 32 33 34 35 36 -250216.313 47930.758 65051.344 71559.952 -190418.162 -148255.495 37 38 39 40 41 42 49490.851 -41774.594 -25648.977 -126404.718 56121.615 129383.608 43 44 45 46 47 48 -32278.740 47930.758 -37101.259 -11465.171 47930.758 37413.394 49 50 51 52 53 54 -124535.591 24236.170 -402357.714 51398.186 -74411.743 -133413.771 55 56 57 58 59 60 498000.209 -306359.379 40596.813 -139539.338 -112771.362 47930.758 61 62 63 64 65 66 47930.758 -249808.092 -23000.886 47847.411 -266560.840 47253.085 67 68 69 70 71 72 -298183.103 -25660.802 -33992.012 47930.758 704705.338 -306374.089 73 74 75 76 77 78 3195283.878 47930.758 -35584.633 73611.546 17805.179 49994.478 79 80 81 82 83 84 -282147.420 -465658.107 48103.003 47930.758 47930.758 -37960.116 85 86 87 88 89 90 47930.758 -7555.060 46990.954 48672.001 37073.804 35414.901 91 92 93 94 95 96 -7555.060 41434.050 -273562.970 -150276.044 50535.241 44777.497 97 98 99 100 101 102 47930.758 48453.678 53038.424 65563.601 152629.298 29601.307 103 104 105 106 107 108 50256.107 48181.113 47930.758 -675705.948 19719.600 -86633.872 109 110 111 112 113 114 55446.705 54791.004 43429.078 -115267.526 47930.758 23566.412 115 116 117 118 119 120 -7555.060 45851.206 -53093.255 153902.340 148737.922 66254.339 121 122 123 124 125 126 -100883.026 47930.758 190360.572 20323.606 -344582.470 -7555.060 127 128 129 130 131 132 33983.667 47930.758 -1118357.214 -65872.062 34349.869 -7555.060 133 134 135 136 137 138 -7555.060 -7555.060 46456.029 -11265.356 47930.758 49743.121 139 140 141 142 143 144 -7555.060 47930.758 127681.941 -123725.963 47930.758 47930.758 145 146 147 148 149 150 -2370.308 -7555.060 41821.079 90994.114 -35871.190 47930.758 151 152 153 154 155 156 -7555.060 17459.373 47930.758 21375.207 -49407.953 -124797.099 157 158 159 160 161 162 85548.186 47930.758 40959.878 29543.810 42217.914 47930.758 163 164 165 166 167 168 48756.725 -5631.489 47930.758 80331.658 47930.758 50336.225 169 170 171 172 173 174 20183.871 37562.856 -65081.319 63483.018 -200546.755 -109444.008 175 176 177 178 179 180 23336.527 54252.289 47930.758 50983.381 29871.149 -7555.060 181 182 183 184 185 186 -28126.641 -6906.064 47930.758 47930.758 -7555.060 51520.595 187 188 189 190 191 192 43892.484 39825.868 -2163938.266 47930.758 49666.953 10638.634 193 194 195 196 197 198 47930.758 59511.531 55714.984 47930.758 -83722.349 -223309.564 199 200 201 202 203 204 -7555.060 47930.758 -279201.028 -7555.060 -7555.060 35201.158 205 206 207 208 209 210 47930.758 63553.861 -58373.820 47930.758 51085.001 6174.204 211 212 213 214 215 216 -12297.139 -223683.527 44403.489 34065.162 -122718.005 -157816.795 217 218 219 220 221 222 -8066.752 47930.758 27844.604 -29593.037 45014.662 47930.758 223 224 225 226 227 228 34892.012 20843.505 25176.204 -459875.779 47930.758 -234688.381 229 230 231 232 233 234 24439.537 27297.538 -144589.537 42833.808 -153760.016 47930.758 235 236 237 238 239 240 43760.928 10402.487 55479.808 47930.758 47930.758 48830.870 241 242 243 244 245 246 12120.541 10492.987 -7555.060 39310.301 -138046.328 118575.766 247 248 249 250 251 252 47930.758 53714.507 47930.758 45207.361 -7555.060 25991.297 253 254 255 256 257 258 -7555.060 47930.758 -7555.060 47930.758 47930.758 47930.758 259 260 261 262 263 264 48357.958 -374596.864 297382.198 -89503.531 -95882.735 -43408.390 265 266 267 268 269 270 -107098.631 47930.758 24781.823 47930.758 39802.780 45723.816 271 272 273 274 275 276 47930.758 47930.758 10354.181 48001.084 105736.955 2068.080 277 278 279 280 281 282 -8516.655 41418.222 9037.666 -20308.732 27433.327 -59079.502 283 284 285 286 287 288 47930.758 19310.917 50337.774 47076.205 47930.758 111838.134 289 290 291 292 293 294 333879.073 -576390.302 -7555.060 35997.643 47930.758 -12187.230 295 296 297 298 299 300 47968.115 48340.985 47930.758 762830.868 -7555.060 -7555.060 301 302 303 304 305 306 -278718.282 -14179.696 -74418.783 42624.223 32019.062 48377.897 307 308 309 310 311 312 -158893.721 47930.758 47930.758 49543.695 -7555.060 -87387.956 313 314 315 316 317 318 -13574.240 47890.901 35089.843 -189194.453 1187360.002 16757.609 319 320 321 322 323 324 -7555.060 -60797.317 50834.989 47803.348 -6488.100 47930.758 325 326 327 328 329 330 -14508.840 47930.758 -67415.436 1259827.443 -1184133.917 -180351.327 331 332 333 334 335 336 -62936.025 47930.758 -7555.060 -63103.461 -69214.165 47930.758 337 338 339 340 341 342 22362.228 -378580.851 47930.758 -7555.060 -265684.483 13127.108 343 344 345 346 347 348 47930.758 60434.071 47930.758 -96498.844 44196.429 -313458.634 349 350 351 352 353 354 47930.758 41837.252 47930.758 47930.758 -9187.034 47930.758 355 356 357 358 359 360 -31544.333 -7555.060 47930.758 -7555.060 47930.758 -9824.282 361 362 363 364 365 366 -348387.682 58126.423 15002.087 -56691.365 55492.576 408709.063 367 368 369 370 371 372 95780.130 49921.272 47930.758 -232554.023 -65383.903 26690.904 373 374 375 376 377 378 58494.152 -159079.702 -10313.459 4473.442 -62195.598 5446.994 379 380 381 382 383 384 28918.905 47930.758 -8951.678 9957.267 -152906.495 49670.808 385 386 387 388 389 390 48197.857 15549.014 -18948.548 -7555.060 47930.758 -7555.060 391 392 393 394 395 396 47930.758 47930.758 22800.421 52466.884 47930.758 -7555.060 397 398 399 400 401 402 47930.758 47930.758 60199.198 19942.976 -127969.710 -223588.738 403 404 405 406 407 408 51524.780 -322029.866 43281.192 22758.922 50155.084 225880.987 409 410 411 412 413 414 -169212.146 -5335.757 20306.970 47930.758 -7555.060 -73845.864 415 416 417 418 419 420 48877.645 50139.283 -248826.358 -40997.771 48709.764 47930.758 421 422 423 424 425 426 46110.459 97881.580 64892.532 22654.245 -35757.824 -484134.906 427 428 429 430 431 -7555.060 51292.884 68636.060 47930.758 75626.658 > postscript(file="/var/www/rcomp/tmp/6n0m41293276484.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 11292.552 NA 1 -41076.166 11292.552 2 -19198.651 -41076.166 3 79985.494 -19198.651 4 -280953.361 79985.494 5 -62387.577 -280953.361 6 17215.934 -62387.577 7 49320.735 17215.934 8 63580.147 49320.735 9 38208.833 63580.147 10 -41888.465 38208.833 11 -68881.815 -41888.465 12 47930.758 -68881.815 13 95359.943 47930.758 14 40017.496 95359.943 15 16128.433 40017.496 16 51727.086 16128.433 17 60651.036 51727.086 18 52885.045 60651.036 19 48104.426 52885.045 20 47930.758 48104.426 21 59484.659 47930.758 22 47930.758 59484.659 23 4172.197 47930.758 24 23515.058 4172.197 25 40932.200 23515.058 26 -7827.369 40932.200 27 -131141.638 -7827.369 28 47930.758 -131141.638 29 17147.822 47930.758 30 -250216.313 17147.822 31 47930.758 -250216.313 32 65051.344 47930.758 33 71559.952 65051.344 34 -190418.162 71559.952 35 -148255.495 -190418.162 36 49490.851 -148255.495 37 -41774.594 49490.851 38 -25648.977 -41774.594 39 -126404.718 -25648.977 40 56121.615 -126404.718 41 129383.608 56121.615 42 -32278.740 129383.608 43 47930.758 -32278.740 44 -37101.259 47930.758 45 -11465.171 -37101.259 46 47930.758 -11465.171 47 37413.394 47930.758 48 -124535.591 37413.394 49 24236.170 -124535.591 50 -402357.714 24236.170 51 51398.186 -402357.714 52 -74411.743 51398.186 53 -133413.771 -74411.743 54 498000.209 -133413.771 55 -306359.379 498000.209 56 40596.813 -306359.379 57 -139539.338 40596.813 58 -112771.362 -139539.338 59 47930.758 -112771.362 60 47930.758 47930.758 61 -249808.092 47930.758 62 -23000.886 -249808.092 63 47847.411 -23000.886 64 -266560.840 47847.411 65 47253.085 -266560.840 66 -298183.103 47253.085 67 -25660.802 -298183.103 68 -33992.012 -25660.802 69 47930.758 -33992.012 70 704705.338 47930.758 71 -306374.089 704705.338 72 3195283.878 -306374.089 73 47930.758 3195283.878 74 -35584.633 47930.758 75 73611.546 -35584.633 76 17805.179 73611.546 77 49994.478 17805.179 78 -282147.420 49994.478 79 -465658.107 -282147.420 80 48103.003 -465658.107 81 47930.758 48103.003 82 47930.758 47930.758 83 -37960.116 47930.758 84 47930.758 -37960.116 85 -7555.060 47930.758 86 46990.954 -7555.060 87 48672.001 46990.954 88 37073.804 48672.001 89 35414.901 37073.804 90 -7555.060 35414.901 91 41434.050 -7555.060 92 -273562.970 41434.050 93 -150276.044 -273562.970 94 50535.241 -150276.044 95 44777.497 50535.241 96 47930.758 44777.497 97 48453.678 47930.758 98 53038.424 48453.678 99 65563.601 53038.424 100 152629.298 65563.601 101 29601.307 152629.298 102 50256.107 29601.307 103 48181.113 50256.107 104 47930.758 48181.113 105 -675705.948 47930.758 106 19719.600 -675705.948 107 -86633.872 19719.600 108 55446.705 -86633.872 109 54791.004 55446.705 110 43429.078 54791.004 111 -115267.526 43429.078 112 47930.758 -115267.526 113 23566.412 47930.758 114 -7555.060 23566.412 115 45851.206 -7555.060 116 -53093.255 45851.206 117 153902.340 -53093.255 118 148737.922 153902.340 119 66254.339 148737.922 120 -100883.026 66254.339 121 47930.758 -100883.026 122 190360.572 47930.758 123 20323.606 190360.572 124 -344582.470 20323.606 125 -7555.060 -344582.470 126 33983.667 -7555.060 127 47930.758 33983.667 128 -1118357.214 47930.758 129 -65872.062 -1118357.214 130 34349.869 -65872.062 131 -7555.060 34349.869 132 -7555.060 -7555.060 133 -7555.060 -7555.060 134 46456.029 -7555.060 135 -11265.356 46456.029 136 47930.758 -11265.356 137 49743.121 47930.758 138 -7555.060 49743.121 139 47930.758 -7555.060 140 127681.941 47930.758 141 -123725.963 127681.941 142 47930.758 -123725.963 143 47930.758 47930.758 144 -2370.308 47930.758 145 -7555.060 -2370.308 146 41821.079 -7555.060 147 90994.114 41821.079 148 -35871.190 90994.114 149 47930.758 -35871.190 150 -7555.060 47930.758 151 17459.373 -7555.060 152 47930.758 17459.373 153 21375.207 47930.758 154 -49407.953 21375.207 155 -124797.099 -49407.953 156 85548.186 -124797.099 157 47930.758 85548.186 158 40959.878 47930.758 159 29543.810 40959.878 160 42217.914 29543.810 161 47930.758 42217.914 162 48756.725 47930.758 163 -5631.489 48756.725 164 47930.758 -5631.489 165 80331.658 47930.758 166 47930.758 80331.658 167 50336.225 47930.758 168 20183.871 50336.225 169 37562.856 20183.871 170 -65081.319 37562.856 171 63483.018 -65081.319 172 -200546.755 63483.018 173 -109444.008 -200546.755 174 23336.527 -109444.008 175 54252.289 23336.527 176 47930.758 54252.289 177 50983.381 47930.758 178 29871.149 50983.381 179 -7555.060 29871.149 180 -28126.641 -7555.060 181 -6906.064 -28126.641 182 47930.758 -6906.064 183 47930.758 47930.758 184 -7555.060 47930.758 185 51520.595 -7555.060 186 43892.484 51520.595 187 39825.868 43892.484 188 -2163938.266 39825.868 189 47930.758 -2163938.266 190 49666.953 47930.758 191 10638.634 49666.953 192 47930.758 10638.634 193 59511.531 47930.758 194 55714.984 59511.531 195 47930.758 55714.984 196 -83722.349 47930.758 197 -223309.564 -83722.349 198 -7555.060 -223309.564 199 47930.758 -7555.060 200 -279201.028 47930.758 201 -7555.060 -279201.028 202 -7555.060 -7555.060 203 35201.158 -7555.060 204 47930.758 35201.158 205 63553.861 47930.758 206 -58373.820 63553.861 207 47930.758 -58373.820 208 51085.001 47930.758 209 6174.204 51085.001 210 -12297.139 6174.204 211 -223683.527 -12297.139 212 44403.489 -223683.527 213 34065.162 44403.489 214 -122718.005 34065.162 215 -157816.795 -122718.005 216 -8066.752 -157816.795 217 47930.758 -8066.752 218 27844.604 47930.758 219 -29593.037 27844.604 220 45014.662 -29593.037 221 47930.758 45014.662 222 34892.012 47930.758 223 20843.505 34892.012 224 25176.204 20843.505 225 -459875.779 25176.204 226 47930.758 -459875.779 227 -234688.381 47930.758 228 24439.537 -234688.381 229 27297.538 24439.537 230 -144589.537 27297.538 231 42833.808 -144589.537 232 -153760.016 42833.808 233 47930.758 -153760.016 234 43760.928 47930.758 235 10402.487 43760.928 236 55479.808 10402.487 237 47930.758 55479.808 238 47930.758 47930.758 239 48830.870 47930.758 240 12120.541 48830.870 241 10492.987 12120.541 242 -7555.060 10492.987 243 39310.301 -7555.060 244 -138046.328 39310.301 245 118575.766 -138046.328 246 47930.758 118575.766 247 53714.507 47930.758 248 47930.758 53714.507 249 45207.361 47930.758 250 -7555.060 45207.361 251 25991.297 -7555.060 252 -7555.060 25991.297 253 47930.758 -7555.060 254 -7555.060 47930.758 255 47930.758 -7555.060 256 47930.758 47930.758 257 47930.758 47930.758 258 48357.958 47930.758 259 -374596.864 48357.958 260 297382.198 -374596.864 261 -89503.531 297382.198 262 -95882.735 -89503.531 263 -43408.390 -95882.735 264 -107098.631 -43408.390 265 47930.758 -107098.631 266 24781.823 47930.758 267 47930.758 24781.823 268 39802.780 47930.758 269 45723.816 39802.780 270 47930.758 45723.816 271 47930.758 47930.758 272 10354.181 47930.758 273 48001.084 10354.181 274 105736.955 48001.084 275 2068.080 105736.955 276 -8516.655 2068.080 277 41418.222 -8516.655 278 9037.666 41418.222 279 -20308.732 9037.666 280 27433.327 -20308.732 281 -59079.502 27433.327 282 47930.758 -59079.502 283 19310.917 47930.758 284 50337.774 19310.917 285 47076.205 50337.774 286 47930.758 47076.205 287 111838.134 47930.758 288 333879.073 111838.134 289 -576390.302 333879.073 290 -7555.060 -576390.302 291 35997.643 -7555.060 292 47930.758 35997.643 293 -12187.230 47930.758 294 47968.115 -12187.230 295 48340.985 47968.115 296 47930.758 48340.985 297 762830.868 47930.758 298 -7555.060 762830.868 299 -7555.060 -7555.060 300 -278718.282 -7555.060 301 -14179.696 -278718.282 302 -74418.783 -14179.696 303 42624.223 -74418.783 304 32019.062 42624.223 305 48377.897 32019.062 306 -158893.721 48377.897 307 47930.758 -158893.721 308 47930.758 47930.758 309 49543.695 47930.758 310 -7555.060 49543.695 311 -87387.956 -7555.060 312 -13574.240 -87387.956 313 47890.901 -13574.240 314 35089.843 47890.901 315 -189194.453 35089.843 316 1187360.002 -189194.453 317 16757.609 1187360.002 318 -7555.060 16757.609 319 -60797.317 -7555.060 320 50834.989 -60797.317 321 47803.348 50834.989 322 -6488.100 47803.348 323 47930.758 -6488.100 324 -14508.840 47930.758 325 47930.758 -14508.840 326 -67415.436 47930.758 327 1259827.443 -67415.436 328 -1184133.917 1259827.443 329 -180351.327 -1184133.917 330 -62936.025 -180351.327 331 47930.758 -62936.025 332 -7555.060 47930.758 333 -63103.461 -7555.060 334 -69214.165 -63103.461 335 47930.758 -69214.165 336 22362.228 47930.758 337 -378580.851 22362.228 338 47930.758 -378580.851 339 -7555.060 47930.758 340 -265684.483 -7555.060 341 13127.108 -265684.483 342 47930.758 13127.108 343 60434.071 47930.758 344 47930.758 60434.071 345 -96498.844 47930.758 346 44196.429 -96498.844 347 -313458.634 44196.429 348 47930.758 -313458.634 349 41837.252 47930.758 350 47930.758 41837.252 351 47930.758 47930.758 352 -9187.034 47930.758 353 47930.758 -9187.034 354 -31544.333 47930.758 355 -7555.060 -31544.333 356 47930.758 -7555.060 357 -7555.060 47930.758 358 47930.758 -7555.060 359 -9824.282 47930.758 360 -348387.682 -9824.282 361 58126.423 -348387.682 362 15002.087 58126.423 363 -56691.365 15002.087 364 55492.576 -56691.365 365 408709.063 55492.576 366 95780.130 408709.063 367 49921.272 95780.130 368 47930.758 49921.272 369 -232554.023 47930.758 370 -65383.903 -232554.023 371 26690.904 -65383.903 372 58494.152 26690.904 373 -159079.702 58494.152 374 -10313.459 -159079.702 375 4473.442 -10313.459 376 -62195.598 4473.442 377 5446.994 -62195.598 378 28918.905 5446.994 379 47930.758 28918.905 380 -8951.678 47930.758 381 9957.267 -8951.678 382 -152906.495 9957.267 383 49670.808 -152906.495 384 48197.857 49670.808 385 15549.014 48197.857 386 -18948.548 15549.014 387 -7555.060 -18948.548 388 47930.758 -7555.060 389 -7555.060 47930.758 390 47930.758 -7555.060 391 47930.758 47930.758 392 22800.421 47930.758 393 52466.884 22800.421 394 47930.758 52466.884 395 -7555.060 47930.758 396 47930.758 -7555.060 397 47930.758 47930.758 398 60199.198 47930.758 399 19942.976 60199.198 400 -127969.710 19942.976 401 -223588.738 -127969.710 402 51524.780 -223588.738 403 -322029.866 51524.780 404 43281.192 -322029.866 405 22758.922 43281.192 406 50155.084 22758.922 407 225880.987 50155.084 408 -169212.146 225880.987 409 -5335.757 -169212.146 410 20306.970 -5335.757 411 47930.758 20306.970 412 -7555.060 47930.758 413 -73845.864 -7555.060 414 48877.645 -73845.864 415 50139.283 48877.645 416 -248826.358 50139.283 417 -40997.771 -248826.358 418 48709.764 -40997.771 419 47930.758 48709.764 420 46110.459 47930.758 421 97881.580 46110.459 422 64892.532 97881.580 423 22654.245 64892.532 424 -35757.824 22654.245 425 -484134.906 -35757.824 426 -7555.060 -484134.906 427 51292.884 -7555.060 428 68636.060 51292.884 429 47930.758 68636.060 430 75626.658 47930.758 431 NA 75626.658 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] -41076.166 11292.552 [2,] -19198.651 -41076.166 [3,] 79985.494 -19198.651 [4,] -280953.361 79985.494 [5,] -62387.577 -280953.361 [6,] 17215.934 -62387.577 [7,] 49320.735 17215.934 [8,] 63580.147 49320.735 [9,] 38208.833 63580.147 [10,] -41888.465 38208.833 [11,] -68881.815 -41888.465 [12,] 47930.758 -68881.815 [13,] 95359.943 47930.758 [14,] 40017.496 95359.943 [15,] 16128.433 40017.496 [16,] 51727.086 16128.433 [17,] 60651.036 51727.086 [18,] 52885.045 60651.036 [19,] 48104.426 52885.045 [20,] 47930.758 48104.426 [21,] 59484.659 47930.758 [22,] 47930.758 59484.659 [23,] 4172.197 47930.758 [24,] 23515.058 4172.197 [25,] 40932.200 23515.058 [26,] -7827.369 40932.200 [27,] -131141.638 -7827.369 [28,] 47930.758 -131141.638 [29,] 17147.822 47930.758 [30,] -250216.313 17147.822 [31,] 47930.758 -250216.313 [32,] 65051.344 47930.758 [33,] 71559.952 65051.344 [34,] -190418.162 71559.952 [35,] -148255.495 -190418.162 [36,] 49490.851 -148255.495 [37,] -41774.594 49490.851 [38,] -25648.977 -41774.594 [39,] -126404.718 -25648.977 [40,] 56121.615 -126404.718 [41,] 129383.608 56121.615 [42,] -32278.740 129383.608 [43,] 47930.758 -32278.740 [44,] -37101.259 47930.758 [45,] -11465.171 -37101.259 [46,] 47930.758 -11465.171 [47,] 37413.394 47930.758 [48,] -124535.591 37413.394 [49,] 24236.170 -124535.591 [50,] -402357.714 24236.170 [51,] 51398.186 -402357.714 [52,] -74411.743 51398.186 [53,] -133413.771 -74411.743 [54,] 498000.209 -133413.771 [55,] -306359.379 498000.209 [56,] 40596.813 -306359.379 [57,] -139539.338 40596.813 [58,] -112771.362 -139539.338 [59,] 47930.758 -112771.362 [60,] 47930.758 47930.758 [61,] -249808.092 47930.758 [62,] -23000.886 -249808.092 [63,] 47847.411 -23000.886 [64,] -266560.840 47847.411 [65,] 47253.085 -266560.840 [66,] -298183.103 47253.085 [67,] -25660.802 -298183.103 [68,] -33992.012 -25660.802 [69,] 47930.758 -33992.012 [70,] 704705.338 47930.758 [71,] -306374.089 704705.338 [72,] 3195283.878 -306374.089 [73,] 47930.758 3195283.878 [74,] -35584.633 47930.758 [75,] 73611.546 -35584.633 [76,] 17805.179 73611.546 [77,] 49994.478 17805.179 [78,] -282147.420 49994.478 [79,] -465658.107 -282147.420 [80,] 48103.003 -465658.107 [81,] 47930.758 48103.003 [82,] 47930.758 47930.758 [83,] -37960.116 47930.758 [84,] 47930.758 -37960.116 [85,] -7555.060 47930.758 [86,] 46990.954 -7555.060 [87,] 48672.001 46990.954 [88,] 37073.804 48672.001 [89,] 35414.901 37073.804 [90,] -7555.060 35414.901 [91,] 41434.050 -7555.060 [92,] -273562.970 41434.050 [93,] -150276.044 -273562.970 [94,] 50535.241 -150276.044 [95,] 44777.497 50535.241 [96,] 47930.758 44777.497 [97,] 48453.678 47930.758 [98,] 53038.424 48453.678 [99,] 65563.601 53038.424 [100,] 152629.298 65563.601 [101,] 29601.307 152629.298 [102,] 50256.107 29601.307 [103,] 48181.113 50256.107 [104,] 47930.758 48181.113 [105,] -675705.948 47930.758 [106,] 19719.600 -675705.948 [107,] -86633.872 19719.600 [108,] 55446.705 -86633.872 [109,] 54791.004 55446.705 [110,] 43429.078 54791.004 [111,] -115267.526 43429.078 [112,] 47930.758 -115267.526 [113,] 23566.412 47930.758 [114,] -7555.060 23566.412 [115,] 45851.206 -7555.060 [116,] -53093.255 45851.206 [117,] 153902.340 -53093.255 [118,] 148737.922 153902.340 [119,] 66254.339 148737.922 [120,] -100883.026 66254.339 [121,] 47930.758 -100883.026 [122,] 190360.572 47930.758 [123,] 20323.606 190360.572 [124,] -344582.470 20323.606 [125,] -7555.060 -344582.470 [126,] 33983.667 -7555.060 [127,] 47930.758 33983.667 [128,] -1118357.214 47930.758 [129,] -65872.062 -1118357.214 [130,] 34349.869 -65872.062 [131,] -7555.060 34349.869 [132,] -7555.060 -7555.060 [133,] -7555.060 -7555.060 [134,] 46456.029 -7555.060 [135,] -11265.356 46456.029 [136,] 47930.758 -11265.356 [137,] 49743.121 47930.758 [138,] -7555.060 49743.121 [139,] 47930.758 -7555.060 [140,] 127681.941 47930.758 [141,] -123725.963 127681.941 [142,] 47930.758 -123725.963 [143,] 47930.758 47930.758 [144,] -2370.308 47930.758 [145,] -7555.060 -2370.308 [146,] 41821.079 -7555.060 [147,] 90994.114 41821.079 [148,] -35871.190 90994.114 [149,] 47930.758 -35871.190 [150,] -7555.060 47930.758 [151,] 17459.373 -7555.060 [152,] 47930.758 17459.373 [153,] 21375.207 47930.758 [154,] -49407.953 21375.207 [155,] -124797.099 -49407.953 [156,] 85548.186 -124797.099 [157,] 47930.758 85548.186 [158,] 40959.878 47930.758 [159,] 29543.810 40959.878 [160,] 42217.914 29543.810 [161,] 47930.758 42217.914 [162,] 48756.725 47930.758 [163,] -5631.489 48756.725 [164,] 47930.758 -5631.489 [165,] 80331.658 47930.758 [166,] 47930.758 80331.658 [167,] 50336.225 47930.758 [168,] 20183.871 50336.225 [169,] 37562.856 20183.871 [170,] -65081.319 37562.856 [171,] 63483.018 -65081.319 [172,] -200546.755 63483.018 [173,] -109444.008 -200546.755 [174,] 23336.527 -109444.008 [175,] 54252.289 23336.527 [176,] 47930.758 54252.289 [177,] 50983.381 47930.758 [178,] 29871.149 50983.381 [179,] -7555.060 29871.149 [180,] -28126.641 -7555.060 [181,] -6906.064 -28126.641 [182,] 47930.758 -6906.064 [183,] 47930.758 47930.758 [184,] -7555.060 47930.758 [185,] 51520.595 -7555.060 [186,] 43892.484 51520.595 [187,] 39825.868 43892.484 [188,] -2163938.266 39825.868 [189,] 47930.758 -2163938.266 [190,] 49666.953 47930.758 [191,] 10638.634 49666.953 [192,] 47930.758 10638.634 [193,] 59511.531 47930.758 [194,] 55714.984 59511.531 [195,] 47930.758 55714.984 [196,] -83722.349 47930.758 [197,] -223309.564 -83722.349 [198,] -7555.060 -223309.564 [199,] 47930.758 -7555.060 [200,] -279201.028 47930.758 [201,] -7555.060 -279201.028 [202,] -7555.060 -7555.060 [203,] 35201.158 -7555.060 [204,] 47930.758 35201.158 [205,] 63553.861 47930.758 [206,] -58373.820 63553.861 [207,] 47930.758 -58373.820 [208,] 51085.001 47930.758 [209,] 6174.204 51085.001 [210,] -12297.139 6174.204 [211,] -223683.527 -12297.139 [212,] 44403.489 -223683.527 [213,] 34065.162 44403.489 [214,] -122718.005 34065.162 [215,] -157816.795 -122718.005 [216,] -8066.752 -157816.795 [217,] 47930.758 -8066.752 [218,] 27844.604 47930.758 [219,] -29593.037 27844.604 [220,] 45014.662 -29593.037 [221,] 47930.758 45014.662 [222,] 34892.012 47930.758 [223,] 20843.505 34892.012 [224,] 25176.204 20843.505 [225,] -459875.779 25176.204 [226,] 47930.758 -459875.779 [227,] -234688.381 47930.758 [228,] 24439.537 -234688.381 [229,] 27297.538 24439.537 [230,] -144589.537 27297.538 [231,] 42833.808 -144589.537 [232,] -153760.016 42833.808 [233,] 47930.758 -153760.016 [234,] 43760.928 47930.758 [235,] 10402.487 43760.928 [236,] 55479.808 10402.487 [237,] 47930.758 55479.808 [238,] 47930.758 47930.758 [239,] 48830.870 47930.758 [240,] 12120.541 48830.870 [241,] 10492.987 12120.541 [242,] -7555.060 10492.987 [243,] 39310.301 -7555.060 [244,] -138046.328 39310.301 [245,] 118575.766 -138046.328 [246,] 47930.758 118575.766 [247,] 53714.507 47930.758 [248,] 47930.758 53714.507 [249,] 45207.361 47930.758 [250,] -7555.060 45207.361 [251,] 25991.297 -7555.060 [252,] -7555.060 25991.297 [253,] 47930.758 -7555.060 [254,] -7555.060 47930.758 [255,] 47930.758 -7555.060 [256,] 47930.758 47930.758 [257,] 47930.758 47930.758 [258,] 48357.958 47930.758 [259,] -374596.864 48357.958 [260,] 297382.198 -374596.864 [261,] -89503.531 297382.198 [262,] -95882.735 -89503.531 [263,] -43408.390 -95882.735 [264,] -107098.631 -43408.390 [265,] 47930.758 -107098.631 [266,] 24781.823 47930.758 [267,] 47930.758 24781.823 [268,] 39802.780 47930.758 [269,] 45723.816 39802.780 [270,] 47930.758 45723.816 [271,] 47930.758 47930.758 [272,] 10354.181 47930.758 [273,] 48001.084 10354.181 [274,] 105736.955 48001.084 [275,] 2068.080 105736.955 [276,] -8516.655 2068.080 [277,] 41418.222 -8516.655 [278,] 9037.666 41418.222 [279,] -20308.732 9037.666 [280,] 27433.327 -20308.732 [281,] -59079.502 27433.327 [282,] 47930.758 -59079.502 [283,] 19310.917 47930.758 [284,] 50337.774 19310.917 [285,] 47076.205 50337.774 [286,] 47930.758 47076.205 [287,] 111838.134 47930.758 [288,] 333879.073 111838.134 [289,] -576390.302 333879.073 [290,] -7555.060 -576390.302 [291,] 35997.643 -7555.060 [292,] 47930.758 35997.643 [293,] -12187.230 47930.758 [294,] 47968.115 -12187.230 [295,] 48340.985 47968.115 [296,] 47930.758 48340.985 [297,] 762830.868 47930.758 [298,] -7555.060 762830.868 [299,] -7555.060 -7555.060 [300,] -278718.282 -7555.060 [301,] -14179.696 -278718.282 [302,] -74418.783 -14179.696 [303,] 42624.223 -74418.783 [304,] 32019.062 42624.223 [305,] 48377.897 32019.062 [306,] -158893.721 48377.897 [307,] 47930.758 -158893.721 [308,] 47930.758 47930.758 [309,] 49543.695 47930.758 [310,] -7555.060 49543.695 [311,] -87387.956 -7555.060 [312,] -13574.240 -87387.956 [313,] 47890.901 -13574.240 [314,] 35089.843 47890.901 [315,] -189194.453 35089.843 [316,] 1187360.002 -189194.453 [317,] 16757.609 1187360.002 [318,] -7555.060 16757.609 [319,] -60797.317 -7555.060 [320,] 50834.989 -60797.317 [321,] 47803.348 50834.989 [322,] -6488.100 47803.348 [323,] 47930.758 -6488.100 [324,] -14508.840 47930.758 [325,] 47930.758 -14508.840 [326,] -67415.436 47930.758 [327,] 1259827.443 -67415.436 [328,] -1184133.917 1259827.443 [329,] -180351.327 -1184133.917 [330,] -62936.025 -180351.327 [331,] 47930.758 -62936.025 [332,] -7555.060 47930.758 [333,] -63103.461 -7555.060 [334,] -69214.165 -63103.461 [335,] 47930.758 -69214.165 [336,] 22362.228 47930.758 [337,] -378580.851 22362.228 [338,] 47930.758 -378580.851 [339,] -7555.060 47930.758 [340,] -265684.483 -7555.060 [341,] 13127.108 -265684.483 [342,] 47930.758 13127.108 [343,] 60434.071 47930.758 [344,] 47930.758 60434.071 [345,] -96498.844 47930.758 [346,] 44196.429 -96498.844 [347,] -313458.634 44196.429 [348,] 47930.758 -313458.634 [349,] 41837.252 47930.758 [350,] 47930.758 41837.252 [351,] 47930.758 47930.758 [352,] -9187.034 47930.758 [353,] 47930.758 -9187.034 [354,] -31544.333 47930.758 [355,] -7555.060 -31544.333 [356,] 47930.758 -7555.060 [357,] -7555.060 47930.758 [358,] 47930.758 -7555.060 [359,] -9824.282 47930.758 [360,] -348387.682 -9824.282 [361,] 58126.423 -348387.682 [362,] 15002.087 58126.423 [363,] -56691.365 15002.087 [364,] 55492.576 -56691.365 [365,] 408709.063 55492.576 [366,] 95780.130 408709.063 [367,] 49921.272 95780.130 [368,] 47930.758 49921.272 [369,] -232554.023 47930.758 [370,] -65383.903 -232554.023 [371,] 26690.904 -65383.903 [372,] 58494.152 26690.904 [373,] -159079.702 58494.152 [374,] -10313.459 -159079.702 [375,] 4473.442 -10313.459 [376,] -62195.598 4473.442 [377,] 5446.994 -62195.598 [378,] 28918.905 5446.994 [379,] 47930.758 28918.905 [380,] -8951.678 47930.758 [381,] 9957.267 -8951.678 [382,] -152906.495 9957.267 [383,] 49670.808 -152906.495 [384,] 48197.857 49670.808 [385,] 15549.014 48197.857 [386,] -18948.548 15549.014 [387,] -7555.060 -18948.548 [388,] 47930.758 -7555.060 [389,] -7555.060 47930.758 [390,] 47930.758 -7555.060 [391,] 47930.758 47930.758 [392,] 22800.421 47930.758 [393,] 52466.884 22800.421 [394,] 47930.758 52466.884 [395,] -7555.060 47930.758 [396,] 47930.758 -7555.060 [397,] 47930.758 47930.758 [398,] 60199.198 47930.758 [399,] 19942.976 60199.198 [400,] -127969.710 19942.976 [401,] -223588.738 -127969.710 [402,] 51524.780 -223588.738 [403,] -322029.866 51524.780 [404,] 43281.192 -322029.866 [405,] 22758.922 43281.192 [406,] 50155.084 22758.922 [407,] 225880.987 50155.084 [408,] -169212.146 225880.987 [409,] -5335.757 -169212.146 [410,] 20306.970 -5335.757 [411,] 47930.758 20306.970 [412,] -7555.060 47930.758 [413,] -73845.864 -7555.060 [414,] 48877.645 -73845.864 [415,] 50139.283 48877.645 [416,] -248826.358 50139.283 [417,] -40997.771 -248826.358 [418,] 48709.764 -40997.771 [419,] 47930.758 48709.764 [420,] 46110.459 47930.758 [421,] 97881.580 46110.459 [422,] 64892.532 97881.580 [423,] 22654.245 64892.532 [424,] -35757.824 22654.245 [425,] -484134.906 -35757.824 [426,] -7555.060 -484134.906 [427,] 51292.884 -7555.060 [428,] 68636.060 51292.884 [429,] 47930.758 68636.060 [430,] 75626.658 47930.758 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 -41076.166 11292.552 2 -19198.651 -41076.166 3 79985.494 -19198.651 4 -280953.361 79985.494 5 -62387.577 -280953.361 6 17215.934 -62387.577 7 49320.735 17215.934 8 63580.147 49320.735 9 38208.833 63580.147 10 -41888.465 38208.833 11 -68881.815 -41888.465 12 47930.758 -68881.815 13 95359.943 47930.758 14 40017.496 95359.943 15 16128.433 40017.496 16 51727.086 16128.433 17 60651.036 51727.086 18 52885.045 60651.036 19 48104.426 52885.045 20 47930.758 48104.426 21 59484.659 47930.758 22 47930.758 59484.659 23 4172.197 47930.758 24 23515.058 4172.197 25 40932.200 23515.058 26 -7827.369 40932.200 27 -131141.638 -7827.369 28 47930.758 -131141.638 29 17147.822 47930.758 30 -250216.313 17147.822 31 47930.758 -250216.313 32 65051.344 47930.758 33 71559.952 65051.344 34 -190418.162 71559.952 35 -148255.495 -190418.162 36 49490.851 -148255.495 37 -41774.594 49490.851 38 -25648.977 -41774.594 39 -126404.718 -25648.977 40 56121.615 -126404.718 41 129383.608 56121.615 42 -32278.740 129383.608 43 47930.758 -32278.740 44 -37101.259 47930.758 45 -11465.171 -37101.259 46 47930.758 -11465.171 47 37413.394 47930.758 48 -124535.591 37413.394 49 24236.170 -124535.591 50 -402357.714 24236.170 51 51398.186 -402357.714 52 -74411.743 51398.186 53 -133413.771 -74411.743 54 498000.209 -133413.771 55 -306359.379 498000.209 56 40596.813 -306359.379 57 -139539.338 40596.813 58 -112771.362 -139539.338 59 47930.758 -112771.362 60 47930.758 47930.758 61 -249808.092 47930.758 62 -23000.886 -249808.092 63 47847.411 -23000.886 64 -266560.840 47847.411 65 47253.085 -266560.840 66 -298183.103 47253.085 67 -25660.802 -298183.103 68 -33992.012 -25660.802 69 47930.758 -33992.012 70 704705.338 47930.758 71 -306374.089 704705.338 72 3195283.878 -306374.089 73 47930.758 3195283.878 74 -35584.633 47930.758 75 73611.546 -35584.633 76 17805.179 73611.546 77 49994.478 17805.179 78 -282147.420 49994.478 79 -465658.107 -282147.420 80 48103.003 -465658.107 81 47930.758 48103.003 82 47930.758 47930.758 83 -37960.116 47930.758 84 47930.758 -37960.116 85 -7555.060 47930.758 86 46990.954 -7555.060 87 48672.001 46990.954 88 37073.804 48672.001 89 35414.901 37073.804 90 -7555.060 35414.901 91 41434.050 -7555.060 92 -273562.970 41434.050 93 -150276.044 -273562.970 94 50535.241 -150276.044 95 44777.497 50535.241 96 47930.758 44777.497 97 48453.678 47930.758 98 53038.424 48453.678 99 65563.601 53038.424 100 152629.298 65563.601 101 29601.307 152629.298 102 50256.107 29601.307 103 48181.113 50256.107 104 47930.758 48181.113 105 -675705.948 47930.758 106 19719.600 -675705.948 107 -86633.872 19719.600 108 55446.705 -86633.872 109 54791.004 55446.705 110 43429.078 54791.004 111 -115267.526 43429.078 112 47930.758 -115267.526 113 23566.412 47930.758 114 -7555.060 23566.412 115 45851.206 -7555.060 116 -53093.255 45851.206 117 153902.340 -53093.255 118 148737.922 153902.340 119 66254.339 148737.922 120 -100883.026 66254.339 121 47930.758 -100883.026 122 190360.572 47930.758 123 20323.606 190360.572 124 -344582.470 20323.606 125 -7555.060 -344582.470 126 33983.667 -7555.060 127 47930.758 33983.667 128 -1118357.214 47930.758 129 -65872.062 -1118357.214 130 34349.869 -65872.062 131 -7555.060 34349.869 132 -7555.060 -7555.060 133 -7555.060 -7555.060 134 46456.029 -7555.060 135 -11265.356 46456.029 136 47930.758 -11265.356 137 49743.121 47930.758 138 -7555.060 49743.121 139 47930.758 -7555.060 140 127681.941 47930.758 141 -123725.963 127681.941 142 47930.758 -123725.963 143 47930.758 47930.758 144 -2370.308 47930.758 145 -7555.060 -2370.308 146 41821.079 -7555.060 147 90994.114 41821.079 148 -35871.190 90994.114 149 47930.758 -35871.190 150 -7555.060 47930.758 151 17459.373 -7555.060 152 47930.758 17459.373 153 21375.207 47930.758 154 -49407.953 21375.207 155 -124797.099 -49407.953 156 85548.186 -124797.099 157 47930.758 85548.186 158 40959.878 47930.758 159 29543.810 40959.878 160 42217.914 29543.810 161 47930.758 42217.914 162 48756.725 47930.758 163 -5631.489 48756.725 164 47930.758 -5631.489 165 80331.658 47930.758 166 47930.758 80331.658 167 50336.225 47930.758 168 20183.871 50336.225 169 37562.856 20183.871 170 -65081.319 37562.856 171 63483.018 -65081.319 172 -200546.755 63483.018 173 -109444.008 -200546.755 174 23336.527 -109444.008 175 54252.289 23336.527 176 47930.758 54252.289 177 50983.381 47930.758 178 29871.149 50983.381 179 -7555.060 29871.149 180 -28126.641 -7555.060 181 -6906.064 -28126.641 182 47930.758 -6906.064 183 47930.758 47930.758 184 -7555.060 47930.758 185 51520.595 -7555.060 186 43892.484 51520.595 187 39825.868 43892.484 188 -2163938.266 39825.868 189 47930.758 -2163938.266 190 49666.953 47930.758 191 10638.634 49666.953 192 47930.758 10638.634 193 59511.531 47930.758 194 55714.984 59511.531 195 47930.758 55714.984 196 -83722.349 47930.758 197 -223309.564 -83722.349 198 -7555.060 -223309.564 199 47930.758 -7555.060 200 -279201.028 47930.758 201 -7555.060 -279201.028 202 -7555.060 -7555.060 203 35201.158 -7555.060 204 47930.758 35201.158 205 63553.861 47930.758 206 -58373.820 63553.861 207 47930.758 -58373.820 208 51085.001 47930.758 209 6174.204 51085.001 210 -12297.139 6174.204 211 -223683.527 -12297.139 212 44403.489 -223683.527 213 34065.162 44403.489 214 -122718.005 34065.162 215 -157816.795 -122718.005 216 -8066.752 -157816.795 217 47930.758 -8066.752 218 27844.604 47930.758 219 -29593.037 27844.604 220 45014.662 -29593.037 221 47930.758 45014.662 222 34892.012 47930.758 223 20843.505 34892.012 224 25176.204 20843.505 225 -459875.779 25176.204 226 47930.758 -459875.779 227 -234688.381 47930.758 228 24439.537 -234688.381 229 27297.538 24439.537 230 -144589.537 27297.538 231 42833.808 -144589.537 232 -153760.016 42833.808 233 47930.758 -153760.016 234 43760.928 47930.758 235 10402.487 43760.928 236 55479.808 10402.487 237 47930.758 55479.808 238 47930.758 47930.758 239 48830.870 47930.758 240 12120.541 48830.870 241 10492.987 12120.541 242 -7555.060 10492.987 243 39310.301 -7555.060 244 -138046.328 39310.301 245 118575.766 -138046.328 246 47930.758 118575.766 247 53714.507 47930.758 248 47930.758 53714.507 249 45207.361 47930.758 250 -7555.060 45207.361 251 25991.297 -7555.060 252 -7555.060 25991.297 253 47930.758 -7555.060 254 -7555.060 47930.758 255 47930.758 -7555.060 256 47930.758 47930.758 257 47930.758 47930.758 258 48357.958 47930.758 259 -374596.864 48357.958 260 297382.198 -374596.864 261 -89503.531 297382.198 262 -95882.735 -89503.531 263 -43408.390 -95882.735 264 -107098.631 -43408.390 265 47930.758 -107098.631 266 24781.823 47930.758 267 47930.758 24781.823 268 39802.780 47930.758 269 45723.816 39802.780 270 47930.758 45723.816 271 47930.758 47930.758 272 10354.181 47930.758 273 48001.084 10354.181 274 105736.955 48001.084 275 2068.080 105736.955 276 -8516.655 2068.080 277 41418.222 -8516.655 278 9037.666 41418.222 279 -20308.732 9037.666 280 27433.327 -20308.732 281 -59079.502 27433.327 282 47930.758 -59079.502 283 19310.917 47930.758 284 50337.774 19310.917 285 47076.205 50337.774 286 47930.758 47076.205 287 111838.134 47930.758 288 333879.073 111838.134 289 -576390.302 333879.073 290 -7555.060 -576390.302 291 35997.643 -7555.060 292 47930.758 35997.643 293 -12187.230 47930.758 294 47968.115 -12187.230 295 48340.985 47968.115 296 47930.758 48340.985 297 762830.868 47930.758 298 -7555.060 762830.868 299 -7555.060 -7555.060 300 -278718.282 -7555.060 301 -14179.696 -278718.282 302 -74418.783 -14179.696 303 42624.223 -74418.783 304 32019.062 42624.223 305 48377.897 32019.062 306 -158893.721 48377.897 307 47930.758 -158893.721 308 47930.758 47930.758 309 49543.695 47930.758 310 -7555.060 49543.695 311 -87387.956 -7555.060 312 -13574.240 -87387.956 313 47890.901 -13574.240 314 35089.843 47890.901 315 -189194.453 35089.843 316 1187360.002 -189194.453 317 16757.609 1187360.002 318 -7555.060 16757.609 319 -60797.317 -7555.060 320 50834.989 -60797.317 321 47803.348 50834.989 322 -6488.100 47803.348 323 47930.758 -6488.100 324 -14508.840 47930.758 325 47930.758 -14508.840 326 -67415.436 47930.758 327 1259827.443 -67415.436 328 -1184133.917 1259827.443 329 -180351.327 -1184133.917 330 -62936.025 -180351.327 331 47930.758 -62936.025 332 -7555.060 47930.758 333 -63103.461 -7555.060 334 -69214.165 -63103.461 335 47930.758 -69214.165 336 22362.228 47930.758 337 -378580.851 22362.228 338 47930.758 -378580.851 339 -7555.060 47930.758 340 -265684.483 -7555.060 341 13127.108 -265684.483 342 47930.758 13127.108 343 60434.071 47930.758 344 47930.758 60434.071 345 -96498.844 47930.758 346 44196.429 -96498.844 347 -313458.634 44196.429 348 47930.758 -313458.634 349 41837.252 47930.758 350 47930.758 41837.252 351 47930.758 47930.758 352 -9187.034 47930.758 353 47930.758 -9187.034 354 -31544.333 47930.758 355 -7555.060 -31544.333 356 47930.758 -7555.060 357 -7555.060 47930.758 358 47930.758 -7555.060 359 -9824.282 47930.758 360 -348387.682 -9824.282 361 58126.423 -348387.682 362 15002.087 58126.423 363 -56691.365 15002.087 364 55492.576 -56691.365 365 408709.063 55492.576 366 95780.130 408709.063 367 49921.272 95780.130 368 47930.758 49921.272 369 -232554.023 47930.758 370 -65383.903 -232554.023 371 26690.904 -65383.903 372 58494.152 26690.904 373 -159079.702 58494.152 374 -10313.459 -159079.702 375 4473.442 -10313.459 376 -62195.598 4473.442 377 5446.994 -62195.598 378 28918.905 5446.994 379 47930.758 28918.905 380 -8951.678 47930.758 381 9957.267 -8951.678 382 -152906.495 9957.267 383 49670.808 -152906.495 384 48197.857 49670.808 385 15549.014 48197.857 386 -18948.548 15549.014 387 -7555.060 -18948.548 388 47930.758 -7555.060 389 -7555.060 47930.758 390 47930.758 -7555.060 391 47930.758 47930.758 392 22800.421 47930.758 393 52466.884 22800.421 394 47930.758 52466.884 395 -7555.060 47930.758 396 47930.758 -7555.060 397 47930.758 47930.758 398 60199.198 47930.758 399 19942.976 60199.198 400 -127969.710 19942.976 401 -223588.738 -127969.710 402 51524.780 -223588.738 403 -322029.866 51524.780 404 43281.192 -322029.866 405 22758.922 43281.192 406 50155.084 22758.922 407 225880.987 50155.084 408 -169212.146 225880.987 409 -5335.757 -169212.146 410 20306.970 -5335.757 411 47930.758 20306.970 412 -7555.060 47930.758 413 -73845.864 -7555.060 414 48877.645 -73845.864 415 50139.283 48877.645 416 -248826.358 50139.283 417 -40997.771 -248826.358 418 48709.764 -40997.771 419 47930.758 48709.764 420 46110.459 47930.758 421 97881.580 46110.459 422 64892.532 97881.580 423 22654.245 64892.532 424 -35757.824 22654.245 425 -484134.906 -35757.824 426 -7555.060 -484134.906 427 51292.884 -7555.060 428 68636.060 51292.884 429 47930.758 68636.060 430 75626.658 47930.758 > 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/rcomp/tmp/7grlp1293276484.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/rcomp/tmp/8grlp1293276484.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/rcomp/tmp/9r0ks1293276484.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/rcomp/tmp/10r0ks1293276484.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/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/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/rcomp/tmp/11cj1y1293276484.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/rcomp/tmp/12f1hm1293276484.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/rcomp/tmp/13cbfc1293276484.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/rcomp/tmp/14ftdi1293276484.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/rcomp/tmp/150uc61293276484.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/rcomp/tmp/164dtu1293276484.tab") + } > > try(system("convert tmp/1kz5g1293276484.ps tmp/1kz5g1293276484.png",intern=TRUE)) character(0) > try(system("convert tmp/2c84j1293276484.ps tmp/2c84j1293276484.png",intern=TRUE)) character(0) > try(system("convert tmp/3c84j1293276484.ps tmp/3c84j1293276484.png",intern=TRUE)) character(0) > try(system("convert tmp/4c84j1293276484.ps tmp/4c84j1293276484.png",intern=TRUE)) character(0) > try(system("convert tmp/5n0m41293276484.ps tmp/5n0m41293276484.png",intern=TRUE)) character(0) > try(system("convert tmp/6n0m41293276484.ps tmp/6n0m41293276484.png",intern=TRUE)) character(0) > try(system("convert tmp/7grlp1293276484.ps tmp/7grlp1293276484.png",intern=TRUE)) character(0) > try(system("convert tmp/8grlp1293276484.ps tmp/8grlp1293276484.png",intern=TRUE)) character(0) > try(system("convert tmp/9r0ks1293276484.ps tmp/9r0ks1293276484.png",intern=TRUE)) character(0) > try(system("convert tmp/10r0ks1293276484.ps tmp/10r0ks1293276484.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 12.230 1.810 14.224