R version 2.9.0 (2009-04-17) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(1 + ,162556 + ,1081 + ,807 + ,213118 + ,230380558 + ,6282154 + ,1 + ,29790 + ,309 + ,444 + ,81767 + ,25266003 + ,4321023 + ,1 + ,87550 + ,458 + ,412 + ,153198 + ,70164684 + ,4111912 + ,0 + ,84738 + ,588 + ,428 + ,-26007 + ,-15292116 + ,223193 + ,1 + ,54660 + ,302 + ,315 + ,126942 + ,37955658 + ,1491348 + ,1 + ,42634 + ,156 + ,168 + ,157214 + ,24525384 + ,1629616 + ,0 + ,40949 + ,481 + ,263 + ,129352 + ,62218312 + ,1398893 + ,1 + ,45187 + ,353 + ,267 + ,234817 + ,75845891 + ,1926517 + ,1 + ,37704 + ,452 + ,228 + ,60448 + ,27322496 + ,983660 + ,1 + ,16275 + ,109 + ,129 + ,47818 + ,5212162 + ,1443586 + ,0 + ,25830 + ,115 + ,104 + ,245546 + ,28237790 + ,1073089 + ,0 + ,12679 + ,110 + ,122 + ,48020 + ,5282200 + ,984885 + ,1 + ,18014 + ,239 + ,393 + ,-1710 + ,-408690 + ,1405225 + ,0 + ,43556 + ,247 + ,190 + ,32648 + ,8064056 + ,227132 + ,1 + ,24811 + ,505 + ,280 + ,95350 + ,47388950 + ,929118 + ,0 + ,6575 + ,159 + ,63 + ,151352 + ,15589256 + ,1071292 + ,0 + ,7123 + ,109 + ,102 + ,288170 + ,31410530 + ,638830 + ,1 + ,21950 + ,519 + ,265 + ,114337 + ,57397174 + ,856956 + ,1 + ,37597 + ,248 + ,234 + ,37884 + ,9395232 + ,992426 + ,0 + ,17821 + ,373 + ,277 + ,122844 + ,45820812 + ,444477 + ,1 + ,12988 + ,119 + ,73 + ,82340 + ,9798460 + ,857217 + ,1 + ,22330 + ,84 + ,67 + ,79801 + ,6703284 + ,711969 + ,0 + ,13326 + ,102 + ,103 + ,165548 + ,16885896 + ,702380 + ,0 + ,16189 + ,295 + ,290 + ,116384 + ,34333280 + ,358589 + ,0 + ,7146 + ,105 + ,83 + ,134028 + ,14072940 + ,297978 + ,0 + ,15824 + ,64 + ,56 + ,63838 + ,4085632 + ,585715 + ,1 + ,27664 + ,282 + ,236 + ,74996 + ,20023932 + ,657954 + ,0 + ,11920 + ,182 + ,73 + ,31080 + ,4009320 + ,209458 + ,0 + ,8568 + ,37 + ,34 + ,32168 + ,1190216 + ,786690 + ,0 + ,14416 + ,361 + ,139 + ,49857 + ,17998377 + ,439798 + ,1 + ,3369 + ,28 + ,26 + ,87161 + ,2440508 + ,688779 + ,1 + ,11819 + ,85 + ,70 + ,106113 + ,9019605 + ,574339 + ,1 + ,6984 + ,45 + ,40 + ,80570 + ,3545080 + ,741409 + ,1 + ,4519 + ,49 + ,42 + ,102129 + ,5004321 + ,597793 + ,0 + ,2220 + ,22 + ,12 + ,301670 + ,6636740 + ,644190 + ,0 + ,18562 + ,155 + ,211 + ,102313 + ,15858515 + ,377934 + ,0 + ,10327 + ,91 + ,74 + ,88577 + ,8060507 + ,640273 + ,1 + ,5336 + ,81 + ,80 + ,112477 + ,9110637 + ,697458 + ,1 + ,2365 + ,79 + ,83 + ,191778 + ,15150462 + ,550608 + ,0 + ,4069 + ,145 + ,131 + ,79804 + ,11571580 + ,207393 + ,0 + ,8636 + ,855 + ,203 + ,128294 + ,104687904 + ,301607 + ,0 + ,13718 + ,61 + ,56 + ,96448 + ,5883328 + ,345783 + ,0 + ,4525 + ,226 + ,89 + ,93811 + ,21201286 + ,501749 + ,0 + ,6869 + ,105 + ,88 + ,117520 + ,12339600 + ,379983 + ,0 + ,4628 + ,62 + ,39 + ,69159 + ,4287858 + ,387475 + ,1 + ,3689 + ,25 + ,25 + ,101792 + ,2443008 + ,377305 + ,1 + ,4891 + ,217 + ,49 + ,210568 + ,5474768 + ,370837 + ,1 + ,7489 + ,322 + ,149 + ,136996 + ,44112712 + ,430866 + ,0 + ,4901 + ,84 + ,58 + ,121920 + ,10241280 + ,469107 + ,0 + ,2284 + ,33 + ,41 + ,76403 + ,2521299 + ,194493 + ,1 + ,3160 + ,108 + ,90 + ,108094 + ,11674152 + ,530670 + ,1 + ,4150 + ,150 + ,136 + ,134759 + ,20213850 + ,518365 + ,1 + ,7285 + ,115 + ,97 + ,188873 + ,21720395 + ,491303 + ,1 + ,1134 + ,162 + ,63 + ,146216 + ,23686992 + ,527021 + ,1 + ,4658 + ,158 + ,114 + ,156608 + ,24744064 + ,233773 + ,0 + ,2384 + ,97 + ,77 + ,61348 + ,5950756 + ,405972 + ,0 + ,3748 + ,9 + ,6 + ,50350 + ,453150 + ,652925 + ,0 + ,5371 + ,66 + ,47 + ,87720 + ,5789520 + ,446211 + ,0 + ,1285 + ,107 + ,51 + ,99489 + ,10645323 + ,341340 + ,1 + ,9327 + ,101 + ,85 + ,87419 + ,8829319 + ,387699 + ,1 + ,5565 + ,47 + ,43 + ,94355 + ,4434685 + ,493408 + ,0 + ,1528 + ,38 + ,32 + ,60326 + ,2292388 + ,146494 + ,1 + ,3122 + ,34 + ,25 + ,94670 + ,3218780 + ,414462 + ,1 + ,7561 + ,87 + ,77 + ,82425 + ,6923700 + ,364304 + ,0 + ,2675 + ,79 + ,54 + ,59017 + ,4662343 + ,355178 + ,0 + ,13253 + ,947 + ,251 + ,90829 + ,86015063 + ,357760 + ,0 + ,880 + ,74 + ,15 + ,80791 + ,5978534 + ,261216 + ,1 + ,2053 + ,53 + ,44 + ,100423 + ,5322419 + ,397144 + ,0 + ,1424 + ,94 + ,73 + ,131116 + ,12324904 + ,374943 + ,1 + ,4036 + ,63 + ,85 + ,100269 + ,6316947 + ,424898 + ,1 + ,3045 + ,58 + ,49 + ,27330 + ,1585140 + ,202055 + ,0 + ,5119 + ,49 + ,38 + ,39039 + ,1912911 + ,378525 + ,0 + ,1431 + ,34 + ,35 + ,106885 + ,3634090 + ,310768 + ,0 + ,554 + ,11 + ,9 + ,79285 + ,872135 + ,325738 + ,0 + ,1975 + ,35 + ,34 + ,118881 + ,4160835 + ,394510 + ,1 + ,1765 + ,20 + ,20 + ,77623 + ,1319591 + ,247060 + ,0 + ,1012 + ,47 + ,29 + ,114768 + ,5394096 + ,368078 + ,0 + ,810 + ,43 + ,11 + ,74015 + ,3182645 + ,236761 + ,0 + ,1280 + ,117 + ,52 + ,69465 + ,8127405 + ,312378 + ,1 + ,666 + ,171 + ,13 + ,117869 + ,20155599 + ,339836 + ,0 + ,1380 + ,26 + ,29 + ,60982 + ,1585532 + ,347385 + ,1 + ,4677 + ,75 + ,66 + ,90131 + ,6579563 + ,426280 + ,0 + ,876 + ,59 + ,33 + ,138971 + ,8199289 + ,352850 + ,0 + ,814 + ,18 + ,15 + ,39625 + ,713250 + ,301881 + ,0 + ,514 + ,15 + ,15 + ,102725 + ,1540875 + ,377516 + ,1 + ,5692 + ,72 + ,68 + ,64239 + ,4625208 + ,357312 + ,0 + ,3642 + ,86 + ,100 + ,90262 + ,7762532 + ,458343 + ,0 + ,540 + ,14 + ,13 + ,103960 + ,1455440 + ,354228 + ,0 + ,2099 + ,64 + ,45 + ,106611 + ,6823104 + ,308636 + ,0 + ,567 + ,11 + ,14 + ,103345 + ,1136795 + ,386212 + ,0 + ,2001 + ,52 + ,36 + ,95551 + ,4968652 + ,393343 + ,1 + ,2949 + ,41 + ,40 + ,82903 + ,3399023 + ,378509 + ,0 + ,2253 + ,99 + ,68 + ,63593 + ,6295707 + ,452469 + ,1 + ,6533 + ,75 + ,29 + ,126910 + ,9518250 + ,364839 + ,0 + ,1889 + ,45 + ,43 + ,37527 + ,1688715 + ,358649 + ,1 + ,3055 + ,43 + ,30 + ,60247 + ,2590621 + ,376641 + ,0 + ,272 + ,8 + ,9 + ,112995 + ,903960 + ,429112 + ,1 + ,1414 + ,198 + ,22 + ,70184 + ,13896432 + ,330546 + ,0 + ,2564 + ,22 + ,19 + ,130140 + ,2863080 + ,403560 + ,1 + ,1383 + ,11 + ,9 + ,73221 + ,805431 + ,317892 + ,1 + ,1261 + ,33 + ,31 + ,76114 + ,2511762 + ,307528 + ,0 + ,975 + ,23 + ,19 + ,90534 + ,2082282 + ,235133 + ,0 + ,3366 + ,80 + ,55 + ,108479 + ,8678320 + ,299243 + ,0 + ,576 + ,18 + ,8 + ,113761 + ,2047698 + ,314073 + ,0 + ,1686 + ,40 + ,28 + ,68696 + ,1923488 + ,368186 + ,0 + ,746 + ,23 + ,29 + ,71561 + ,1645903 + ,269661 + ,1 + ,3192 + ,60 + ,48 + ,59831 + ,3589860 + ,125390 + ,1 + ,2045 + ,20 + ,16 + ,97890 + ,1957800 + ,510834 + ,0 + ,5702 + ,61 + ,47 + ,101481 + ,5987379 + ,321896 + ,1 + ,1932 + ,36 + ,20 + ,72954 + ,2626344 + ,249898 + ,0 + ,936 + ,30 + ,22 + ,67939 + ,2038170 + ,408881 + ,1 + ,3437 + ,47 + ,33 + ,48022 + ,2257034 + ,158492 + ,0 + ,5131 + ,71 + ,44 + ,86111 + ,6113881 + ,292154 + ,1 + ,2397 + ,14 + ,13 + ,74020 + ,1036280 + ,289513 + ,1 + ,1389 + ,9 + ,6 + ,57530 + ,517770 + ,378049 + ,0 + ,1503 + ,39 + ,35 + ,56364 + ,2198196 + ,343466 + ,0 + ,402 + ,26 + ,8 + ,84990 + ,2209740 + ,332743 + ,0 + ,2239 + ,21 + ,17 + ,88590 + ,1860390 + ,442882 + ,1 + ,2234 + ,16 + ,11 + ,77200 + ,1235200 + ,214215 + ,0 + ,837 + ,69 + ,21 + ,61262 + ,4227078 + ,315688 + ,0 + ,10579 + ,92 + ,92 + ,110309 + ,10148428 + ,375195 + ,0 + ,875 + ,14 + ,12 + ,67000 + ,938000 + ,334280 + ,0 + ,1585 + ,107 + ,112 + ,93099 + ,9589197 + ,355864 + ,1 + ,1659 + ,29 + ,25 + ,107577 + ,3119733 + ,480382 + ,1 + ,2647 + ,37 + ,17 + ,62920 + ,2328040 + ,353058 + ,1 + ,3294 + ,23 + ,23 + ,75832 + ,1744136 + ,217193 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,94 + ,7 + ,10 + ,60793 + ,425551 + ,314533 + ,0 + ,422 + ,28 + ,23 + ,57935 + ,1622180 + ,318056 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,34 + ,8 + ,7 + ,60630 + ,485040 + ,314353 + ,0 + ,1558 + ,63 + ,25 + ,55637 + ,3505131 + ,369448 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,43 + ,3 + ,20 + ,60887 + ,182661 + ,312846 + ,0 + ,645 + ,5 + ,4 + ,60720 + ,0 + ,312075 + ,0 + ,316 + ,9 + ,4 + ,60505 + ,544545 + ,315009 + ,0 + ,115 + ,13 + ,10 + ,60945 + ,792285 + ,318903 + ,0 + ,5 + ,2 + ,1 + ,60720 + ,0 + ,314887 + ,0 + ,897 + ,5 + ,4 + ,60720 + ,0 + ,314913 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,389 + ,14 + ,8 + ,58990 + ,825860 + ,325506 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,1002 + ,15 + ,11 + ,56750 + ,851250 + ,298568 + ,0 + ,36 + ,3 + ,4 + ,60894 + ,182682 + ,315834 + ,0 + ,460 + ,15 + ,15 + ,63346 + ,950190 + ,329784 + ,0 + ,309 + ,11 + ,9 + ,56535 + ,621885 + ,312878 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,9 + ,6 + ,7 + ,60835 + ,365010 + ,314987 + ,0 + ,271 + ,2 + ,2 + ,60720 + ,0 + ,325249 + ,0 + ,14 + ,1 + ,0 + ,61016 + ,61016 + ,315877 + ,0 + ,520 + ,10 + ,7 + ,58650 + ,586500 + ,291650 + ,0 + ,1766 + ,73 + ,46 + ,60438 + ,4411974 + ,305959 + ,0 + ,0 + ,0 + ,5 + ,60720 + ,0 + ,315380 + ,0 + ,458 + ,11 + ,7 + ,58625 + ,644875 + ,297765 + ,0 + ,20 + ,3 + ,2 + ,60938 + ,182814 + ,315245 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,98 + ,2 + ,2 + ,61490 + ,122980 + ,315236 + ,0 + ,405 + ,7 + ,5 + ,60845 + ,425915 + ,336425 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,483 + ,27 + ,7 + ,60830 + ,1642410 + ,306268 + ,0 + ,454 + ,51 + ,24 + ,63261 + ,3226311 + ,302187 + ,0 + ,47 + ,3 + ,1 + ,60720 + ,0 + ,314882 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,757 + ,19 + ,18 + ,45689 + ,868091 + ,382712 + ,0 + ,4655 + ,393 + ,55 + ,60720 + ,0 + ,341570 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,36 + ,4 + ,3 + ,61564 + ,246256 + ,312412 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,203 + ,9 + ,9 + ,61938 + ,557442 + ,309596 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,126 + ,10 + ,8 + ,60951 + ,487608 + ,315547 + ,0 + ,400 + ,152 + ,113 + ,60720 + ,0 + ,313267 + ,1 + ,71 + ,1 + ,0 + ,60745 + ,60745 + ,316176 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,972 + ,34 + ,19 + ,71642 + ,2435828 + ,359335 + ,1 + ,531 + ,10 + ,11 + ,71641 + ,716410 + ,330068 + ,0 + ,2461 + ,57 + ,25 + ,55792 + ,2120096 + ,314289 + ,1 + ,378 + ,52 + ,16 + ,71873 + ,718730 + ,297413 + ,1 + ,23 + ,5 + ,5 + ,62555 + ,312775 + ,314806 + ,1 + ,638 + ,14 + ,11 + ,60370 + ,845180 + ,333210 + ,1 + ,2300 + ,29 + ,23 + ,64873 + ,1037968 + ,352108 + ,0 + ,149 + ,5 + ,6 + ,62041 + ,310205 + ,313332 + ,0 + ,226 + ,5 + ,5 + ,65745 + ,328725 + ,291787 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,275 + ,4 + ,7 + ,59500 + ,238000 + ,318745 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,141 + ,6 + ,7 + ,61630 + ,369780 + ,315366 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,28 + ,2 + ,3 + ,60890 + ,121780 + ,315688 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,4980 + ,91 + ,89 + ,113521 + ,9081680 + ,409642 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,472 + ,20 + ,19 + ,80045 + ,1600900 + ,269587 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,203 + ,27 + ,12 + ,50804 + ,1371708 + ,300962 + ,1 + ,496 + ,17 + ,12 + ,87390 + ,1485630 + ,325479 + ,1 + ,10 + ,2 + ,5 + ,61656 + ,123312 + ,316155 + ,1 + ,63 + ,4 + ,2 + ,65688 + ,262752 + ,318574 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,1136 + ,32 + ,26 + ,48522 + ,1552704 + ,343613 + ,1 + ,265 + ,31 + ,3 + ,60720 + ,0 + ,306948 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,267 + ,32 + ,11 + ,57640 + ,1844480 + ,330059 + ,0 + ,474 + ,20 + ,10 + ,61977 + ,1239540 + ,288985 + ,0 + ,534 + ,7 + ,5 + ,62620 + ,438340 + ,304485 + ,0 + ,0 + ,0 + ,2 + ,60720 + ,0 + ,315380 + ,0 + ,15 + ,8 + ,6 + ,60831 + ,486648 + ,315688 + ,0 + ,397 + ,28 + ,7 + ,60646 + ,1698088 + ,317736 + ,0 + ,0 + ,0 + ,2 + ,60720 + ,0 + ,315380 + ,0 + ,1866 + ,29 + ,28 + ,56225 + ,1124500 + ,322331 + ,0 + ,288 + ,4 + ,3 + ,60510 + ,242040 + ,296656 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,3 + ,2 + ,1 + ,60698 + ,121396 + ,315354 + ,0 + ,468 + ,21 + ,20 + ,60720 + ,0 + ,312161 + ,0 + ,20 + ,2 + ,1 + ,60805 + ,121610 + ,315576 + ,0 + ,278 + ,26 + ,22 + ,61404 + ,1596504 + ,314922 + ,0 + ,61 + ,14 + ,9 + ,60720 + ,0 + ,314551 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,192 + ,4 + ,2 + ,65276 + ,261104 + ,312339 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,317 + ,9 + ,7 + ,63915 + ,575235 + ,298700 + ,0 + ,738 + ,10 + ,9 + ,60720 + ,0 + ,321376 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,368 + ,17 + ,13 + ,61686 + ,1048662 + ,303230 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,2 + ,1 + ,0 + ,60743 + ,60743 + ,315487 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,53 + ,6 + ,6 + ,60349 + ,362094 + ,315793 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,94 + ,3 + ,3 + ,61360 + ,184080 + ,312887 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,24 + ,8 + ,7 + ,59818 + ,478544 + ,315637 + ,0 + ,2332 + ,4 + ,2 + ,72680 + ,290720 + ,324385 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,131 + ,11 + ,15 + ,61808 + ,679888 + ,308989 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,206 + ,9 + ,9 + ,53110 + ,477990 + ,296702 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,167 + ,2 + ,1 + ,64245 + ,128490 + ,307322 + ,0 + ,622 + ,73 + ,38 + ,73007 + ,5329511 + ,304376 + ,0 + ,2328 + ,94 + ,57 + ,82732 + ,7032220 + ,253588 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,365 + ,8 + ,7 + ,54820 + ,438560 + ,309560 + ,0 + ,364 + ,35 + ,26 + ,47705 + ,1669675 + ,298466 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,226 + ,12 + ,13 + ,72835 + ,874020 + ,343929 + ,0 + ,307 + ,15 + ,10 + ,58856 + ,882840 + ,331955 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,188 + ,11 + ,9 + ,77655 + ,854205 + ,381180 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,138 + ,6 + ,26 + ,69817 + ,418902 + ,331420 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,125 + ,12 + ,19 + ,60798 + ,729576 + ,310201 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,282 + ,30 + ,12 + ,62452 + ,1873560 + ,320016 + ,0 + ,335 + ,33 + ,23 + ,64175 + ,2117775 + ,320398 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,1324 + ,117 + ,29 + ,67440 + ,5530080 + ,291841 + ,0 + ,176 + ,28 + ,8 + ,68136 + ,1907808 + ,310670 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,249 + ,72 + ,26 + ,56726 + ,4084272 + ,313491 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,333 + ,13 + ,9 + ,70811 + ,920543 + ,331323 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,601 + ,6 + ,5 + ,60720 + ,0 + ,319210 + ,0 + ,30 + ,4 + ,3 + ,62045 + ,248180 + ,318098 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,249 + ,62 + ,13 + ,54323 + ,3368026 + ,292754 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,165 + ,24 + ,12 + ,62841 + ,1508184 + ,325176 + ,0 + ,453 + ,21 + ,19 + ,81125 + ,1703625 + ,365959 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,53 + ,14 + ,10 + ,59506 + ,833084 + ,302409 + ,1 + ,382 + ,21 + ,9 + ,59365 + ,1246665 + ,340968 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,9 + ,60720 + ,0 + ,315380 + ,1 + ,30 + ,4 + ,4 + ,60798 + ,243192 + ,313164 + ,0 + ,290 + ,2 + ,1 + ,58790 + ,117580 + ,301164 + ,0 + ,0 + ,0 + ,1 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,366 + ,53 + ,14 + ,61808 + ,3275824 + ,344425 + ,0 + ,2 + ,9 + ,12 + ,60735 + ,546615 + ,315394 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,209 + ,13 + ,19 + ,64016 + ,832208 + ,316647 + ,0 + ,384 + ,22 + ,17 + ,54683 + ,1203026 + ,309836 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,365 + ,83 + ,32 + ,87192 + ,7236936 + ,346611 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,49 + ,8 + ,14 + ,64107 + ,512856 + ,322031 + ,0 + ,3 + ,4 + ,8 + ,60761 + ,243044 + ,315656 + ,0 + ,133 + ,14 + ,4 + ,65990 + ,923860 + ,339445 + ,0 + ,32 + ,1 + ,0 + ,59988 + ,59988 + ,314964 + ,0 + ,368 + ,17 + ,20 + ,61167 + ,1039839 + ,297141 + ,0 + ,1 + ,6 + ,5 + ,60719 + ,364314 + ,315372 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,22 + ,2 + ,1 + ,60722 + ,121444 + ,312502 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,96 + ,5 + ,4 + ,60379 + ,301895 + ,313729 + ,0 + ,1 + ,2 + ,1 + ,60727 + ,121454 + ,315388 + ,0 + ,314 + ,5 + ,4 + ,60720 + ,0 + ,315371 + ,0 + ,844 + ,78 + ,20 + ,60925 + ,426475 + ,296139 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,26 + ,1 + ,1 + ,60896 + ,60896 + ,313880 + ,0 + ,125 + ,13 + ,10 + ,59734 + ,776542 + ,317698 + ,0 + ,304 + ,15 + ,12 + ,62969 + ,944535 + ,295580 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,621 + ,48 + ,13 + ,60720 + ,0 + ,308256 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,119 + ,6 + ,3 + ,59118 + ,354708 + ,303677 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,1595 + ,17 + ,10 + ,60720 + ,0 + ,319369 + ,0 + ,312 + ,14 + ,3 + ,58598 + ,820372 + ,318690 + ,0 + ,60 + ,10 + ,7 + ,61124 + ,611240 + ,314049 + ,0 + ,587 + ,12 + ,10 + ,59595 + ,715140 + ,325699 + ,0 + ,135 + ,2 + ,1 + ,62065 + ,124130 + ,314210 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,514 + ,52 + ,15 + ,78780 + ,4096560 + ,322378 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,1 + ,4 + ,4 + ,60722 + ,242888 + ,315398 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,1 + ,1763 + ,24 + ,28 + ,61600 + ,184800 + ,308336 + ,0 + ,180 + ,11 + ,9 + ,59635 + ,655985 + ,316386 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,218 + ,21 + ,7 + ,60720 + ,0 + ,315553 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,448 + ,40 + ,7 + ,59781 + ,2391240 + ,323361 + ,0 + ,227 + ,9 + ,7 + ,76644 + ,689796 + ,336639 + ,0 + ,174 + ,1 + ,3 + ,64820 + ,64820 + ,307424 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,121 + ,24 + ,11 + ,56178 + ,1348272 + ,295370 + ,0 + ,607 + ,11 + ,7 + ,60436 + ,664796 + ,322340 + ,0 + ,2212 + ,14 + ,10 + ,60720 + ,0 + ,319864 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,530 + ,60 + ,18 + ,73433 + ,4405980 + ,317291 + ,0 + ,571 + ,80 + ,14 + ,41477 + ,3318160 + ,280398 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,0 + ,315380 + ,0 + ,78 + ,16 + ,12 + ,62700 + ,1003200 + ,317330 + ,0 + ,2489 + ,40 + ,29 + ,67804 + ,2712160 + ,238125 + ,0 + ,131 + ,6 + ,3 + ,59661 + ,357966 + ,327071 + ,0 + ,923 + ,8 + ,6 + ,58620 + ,468960 + ,309038 + ,0 + ,72 + ,3 + ,3 + ,60398 + ,181194 + ,314210 + ,0 + ,572 + ,16 + ,8 + ,58580 + ,937280 + ,307930 + ,0 + ,397 + ,10 + ,10 + ,62710 + ,627100 + ,322327 + ,0 + ,450 + ,8 + ,6 + ,59325 + ,474600 + ,292136 + ,0 + ,622 + ,7 + ,8 + ,60950 + ,426650 + ,263276 + ,0 + ,694 + ,8 + ,6 + ,68060 + ,544480 + ,367655 + ,1 + ,3425 + ,12 + ,9 + ,83620 + ,1003440 + ,283910 + ,0 + ,562 + ,13 + ,8 + ,58456 + ,759928 + ,283587 + ,0 + ,4917 + ,42 + ,26 + ,52811 + ,2218062 + ,243650 + ,1 + ,1442 + ,118 + ,239 + ,121173 + ,14298414 + ,438493 + ,0 + ,529 + ,9 + ,7 + ,63870 + ,574830 + ,296261 + ,1 + ,2126 + ,138 + ,41 + ,21001 + ,2898138 + ,230621 + ,0 + ,1061 + ,5 + ,3 + ,70415 + ,352075 + ,304252 + ,0 + ,776 + ,9 + ,8 + ,64230 + ,578070 + ,333505 + ,0 + ,611 + ,8 + ,6 + ,59190 + ,473520 + ,296919 + ,1 + ,1526 + ,25 + ,21 + ,69351 + ,1733775 + ,278990 + ,0 + ,592 + ,7 + ,7 + ,64270 + ,449890 + ,276898 + ,0 + ,1182 + ,13 + ,11 + ,70694 + ,919022 + ,327007 + ,0 + ,621 + ,16 + ,11 + ,68005 + ,1088080 + ,317046 + ,0 + ,989 + ,11 + ,12 + ,58930 + ,648230 + ,304555 + ,0 + ,438 + ,11 + ,9 + ,58320 + ,641520 + ,298096 + ,0 + ,726 + ,3 + ,3 + ,69980 + ,209940 + ,231861 + ,0 + ,1303 + ,61 + ,57 + ,69863 + ,4261643 + ,309422 + ,1 + ,7419 + ,29 + ,21 + ,63255 + ,1518120 + ,286963 + ,1 + ,1164 + ,17 + ,15 + ,57320 + ,974440 + ,269753 + ,1 + ,3310 + ,33 + ,32 + ,75230 + ,2482590 + ,448243 + ,0 + ,1920 + ,15 + ,11 + ,79420 + ,555940 + ,165404 + ,0 + ,965 + ,3 + ,2 + ,73490 + ,220470 + ,204325 + ,0 + ,3256 + ,66 + ,23 + ,35250 + ,2326500 + ,407159 + ,1 + ,1135 + ,17 + ,20 + ,62285 + ,1058845 + ,290476 + ,0 + ,1270 + ,26 + ,24 + ,69206 + ,1799356 + ,275311 + ,0 + ,661 + ,3 + ,1 + ,65920 + ,197760 + ,246541 + ,0 + ,1013 + ,2 + ,1 + ,69770 + ,139540 + ,253468 + ,0 + ,2844 + ,67 + ,74 + ,72683 + ,4869761 + ,240897 + ,1 + ,11528 + ,70 + ,68 + ,-14545 + ,-1018150 + ,-83265 + ,0 + ,6526 + ,26 + ,20 + ,55830 + ,1451580 + ,-42143 + ,0 + ,2264 + ,24 + ,20 + ,55174 + ,1324176 + ,272713 + ,1 + ,5109 + ,97 + ,82 + ,67038 + ,6301572 + ,215362 + ,0 + ,3999 + ,30 + ,21 + ,51252 + ,1537560 + ,42754 + ,0 + ,35624 + ,223 + ,244 + ,157278 + ,35072994 + ,306275 + ,0 + ,9252 + ,48 + ,32 + ,79510 + ,3816480 + ,253537 + ,0 + ,15236 + ,90 + ,86 + ,77440 + ,6969600 + ,372631 + ,0 + ,18073 + ,180 + ,69 + ,27284 + ,4911120 + ,-7170) + ,dim=c(7 + ,431) + ,dimnames=list(c('Group' + ,'Costs' + ,'Trades' + ,'Orders' + ,'Dividends' + ,'TrDiv' + ,'Wealth ') + ,1:431)) > y <- array(NA,dim=c(7,431),dimnames=list(c('Group','Costs','Trades','Orders','Dividends','TrDiv','Wealth '),1:431)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par3 = 'No Linear Trend' > par2 = 'Do not include Seasonal Dummies' > par1 = '7' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > library(lattice) > library(lmtest) Loading required package: zoo Attaching package: 'zoo' The following object(s) are masked from package:base : as.Date.numeric > n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test > par1 <- as.numeric(par1) > x <- t(y) > k <- length(x[1,]) > n <- length(x[,1]) > x1 <- cbind(x[,par1], x[,1:k!=par1]) > mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1]) > colnames(x1) <- mycolnames #colnames(x)[par1] > x <- x1 > if (par3 == 'First Differences'){ + x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep=''))) + for (i in 1:n-1) { + for (j in 1:k) { + x2[i,j] <- x[i+1,j] - x[i,j] + } + } + x <- x2 + } > if (par2 == 'Include Monthly Dummies'){ + x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep =''))) + for (i in 1:11){ + x2[seq(i,n,12),i] <- 1 + } + x <- cbind(x, x2) + } > if (par2 == 'Include Quarterly Dummies'){ + x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep =''))) + for (i in 1:3){ + x2[seq(i,n,4),i] <- 1 + } + x <- cbind(x, x2) + } > k <- length(x[1,]) > if (par3 == 'Linear Trend'){ + x <- cbind(x, c(1:n)) + colnames(x)[k+1] <- 't' + } > x Wealth\r Group Costs Trades Orders Dividends TrDiv 1 6282154 1 162556 1081 807 213118 230380558 2 4321023 1 29790 309 444 81767 25266003 3 4111912 1 87550 458 412 153198 70164684 4 223193 0 84738 588 428 -26007 -15292116 5 1491348 1 54660 302 315 126942 37955658 6 1629616 1 42634 156 168 157214 24525384 7 1398893 0 40949 481 263 129352 62218312 8 1926517 1 45187 353 267 234817 75845891 9 983660 1 37704 452 228 60448 27322496 10 1443586 1 16275 109 129 47818 5212162 11 1073089 0 25830 115 104 245546 28237790 12 984885 0 12679 110 122 48020 5282200 13 1405225 1 18014 239 393 -1710 -408690 14 227132 0 43556 247 190 32648 8064056 15 929118 1 24811 505 280 95350 47388950 16 1071292 0 6575 159 63 151352 15589256 17 638830 0 7123 109 102 288170 31410530 18 856956 1 21950 519 265 114337 57397174 19 992426 1 37597 248 234 37884 9395232 20 444477 0 17821 373 277 122844 45820812 21 857217 1 12988 119 73 82340 9798460 22 711969 1 22330 84 67 79801 6703284 23 702380 0 13326 102 103 165548 16885896 24 358589 0 16189 295 290 116384 34333280 25 297978 0 7146 105 83 134028 14072940 26 585715 0 15824 64 56 63838 4085632 27 657954 1 27664 282 236 74996 20023932 28 209458 0 11920 182 73 31080 4009320 29 786690 0 8568 37 34 32168 1190216 30 439798 0 14416 361 139 49857 17998377 31 688779 1 3369 28 26 87161 2440508 32 574339 1 11819 85 70 106113 9019605 33 741409 1 6984 45 40 80570 3545080 34 597793 1 4519 49 42 102129 5004321 35 644190 0 2220 22 12 301670 6636740 36 377934 0 18562 155 211 102313 15858515 37 640273 0 10327 91 74 88577 8060507 38 697458 1 5336 81 80 112477 9110637 39 550608 1 2365 79 83 191778 15150462 40 207393 0 4069 145 131 79804 11571580 41 301607 0 8636 855 203 128294 104687904 42 345783 0 13718 61 56 96448 5883328 43 501749 0 4525 226 89 93811 21201286 44 379983 0 6869 105 88 117520 12339600 45 387475 0 4628 62 39 69159 4287858 46 377305 1 3689 25 25 101792 2443008 47 370837 1 4891 217 49 210568 5474768 48 430866 1 7489 322 149 136996 44112712 49 469107 0 4901 84 58 121920 10241280 50 194493 0 2284 33 41 76403 2521299 51 530670 1 3160 108 90 108094 11674152 52 518365 1 4150 150 136 134759 20213850 53 491303 1 7285 115 97 188873 21720395 54 527021 1 1134 162 63 146216 23686992 55 233773 1 4658 158 114 156608 24744064 56 405972 0 2384 97 77 61348 5950756 57 652925 0 3748 9 6 50350 453150 58 446211 0 5371 66 47 87720 5789520 59 341340 0 1285 107 51 99489 10645323 60 387699 1 9327 101 85 87419 8829319 61 493408 1 5565 47 43 94355 4434685 62 146494 0 1528 38 32 60326 2292388 63 414462 1 3122 34 25 94670 3218780 64 364304 1 7561 87 77 82425 6923700 65 355178 0 2675 79 54 59017 4662343 66 357760 0 13253 947 251 90829 86015063 67 261216 0 880 74 15 80791 5978534 68 397144 1 2053 53 44 100423 5322419 69 374943 0 1424 94 73 131116 12324904 70 424898 1 4036 63 85 100269 6316947 71 202055 1 3045 58 49 27330 1585140 72 378525 0 5119 49 38 39039 1912911 73 310768 0 1431 34 35 106885 3634090 74 325738 0 554 11 9 79285 872135 75 394510 0 1975 35 34 118881 4160835 76 247060 1 1765 20 20 77623 1319591 77 368078 0 1012 47 29 114768 5394096 78 236761 0 810 43 11 74015 3182645 79 312378 0 1280 117 52 69465 8127405 80 339836 1 666 171 13 117869 20155599 81 347385 0 1380 26 29 60982 1585532 82 426280 1 4677 75 66 90131 6579563 83 352850 0 876 59 33 138971 8199289 84 301881 0 814 18 15 39625 713250 85 377516 0 514 15 15 102725 1540875 86 357312 1 5692 72 68 64239 4625208 87 458343 0 3642 86 100 90262 7762532 88 354228 0 540 14 13 103960 1455440 89 308636 0 2099 64 45 106611 6823104 90 386212 0 567 11 14 103345 1136795 91 393343 0 2001 52 36 95551 4968652 92 378509 1 2949 41 40 82903 3399023 93 452469 0 2253 99 68 63593 6295707 94 364839 1 6533 75 29 126910 9518250 95 358649 0 1889 45 43 37527 1688715 96 376641 1 3055 43 30 60247 2590621 97 429112 0 272 8 9 112995 903960 98 330546 1 1414 198 22 70184 13896432 99 403560 0 2564 22 19 130140 2863080 100 317892 1 1383 11 9 73221 805431 101 307528 1 1261 33 31 76114 2511762 102 235133 0 975 23 19 90534 2082282 103 299243 0 3366 80 55 108479 8678320 104 314073 0 576 18 8 113761 2047698 105 368186 0 1686 40 28 68696 1923488 106 269661 0 746 23 29 71561 1645903 107 125390 1 3192 60 48 59831 3589860 108 510834 1 2045 20 16 97890 1957800 109 321896 0 5702 61 47 101481 5987379 110 249898 1 1932 36 20 72954 2626344 111 408881 0 936 30 22 67939 2038170 112 158492 1 3437 47 33 48022 2257034 113 292154 0 5131 71 44 86111 6113881 114 289513 1 2397 14 13 74020 1036280 115 378049 1 1389 9 6 57530 517770 116 343466 0 1503 39 35 56364 2198196 117 332743 0 402 26 8 84990 2209740 118 442882 0 2239 21 17 88590 1860390 119 214215 1 2234 16 11 77200 1235200 120 315688 0 837 69 21 61262 4227078 121 375195 0 10579 92 92 110309 10148428 122 334280 0 875 14 12 67000 938000 123 355864 0 1585 107 112 93099 9589197 124 480382 1 1659 29 25 107577 3119733 125 353058 1 2647 37 17 62920 2328040 126 217193 1 3294 23 23 75832 1744136 127 315380 0 0 0 0 60720 0 128 314533 0 94 7 10 60793 425551 129 318056 0 422 28 23 57935 1622180 130 315380 0 0 0 0 60720 0 131 314353 0 34 8 7 60630 485040 132 369448 0 1558 63 25 55637 3505131 133 315380 0 0 0 0 60720 0 134 312846 0 43 3 20 60887 182661 135 312075 0 645 5 4 60720 0 136 315009 0 316 9 4 60505 544545 137 318903 0 115 13 10 60945 792285 138 314887 0 5 2 1 60720 0 139 314913 0 897 5 4 60720 0 140 315380 0 0 0 0 60720 0 141 325506 0 389 14 8 58990 825860 142 315380 0 0 0 0 60720 0 143 298568 0 1002 15 11 56750 851250 144 315834 0 36 3 4 60894 182682 145 329784 0 460 15 15 63346 950190 146 312878 0 309 11 9 56535 621885 147 315380 0 0 0 0 60720 0 148 314987 0 9 6 7 60835 365010 149 325249 0 271 2 2 60720 0 150 315877 0 14 1 0 61016 61016 151 291650 0 520 10 7 58650 586500 152 305959 0 1766 73 46 60438 4411974 153 315380 0 0 0 5 60720 0 154 297765 0 458 11 7 58625 644875 155 315245 0 20 3 2 60938 182814 156 315380 0 0 0 0 60720 0 157 315380 0 0 0 0 60720 0 158 315236 0 98 2 2 61490 122980 159 336425 0 405 7 5 60845 425915 160 315380 0 0 0 0 60720 0 161 315380 0 0 0 0 60720 0 162 315380 0 0 0 0 60720 0 163 315380 0 0 0 0 60720 0 164 306268 0 483 27 7 60830 1642410 165 302187 0 454 51 24 63261 3226311 166 314882 0 47 3 1 60720 0 167 315380 0 0 0 0 60720 0 168 382712 0 757 19 18 45689 868091 169 341570 0 4655 393 55 60720 0 170 315380 0 0 0 0 60720 0 171 315380 0 0 0 0 60720 0 172 312412 0 36 4 3 61564 246256 173 315380 0 0 0 0 60720 0 174 309596 0 203 9 9 61938 557442 175 315380 0 0 0 0 60720 0 176 315547 0 126 10 8 60951 487608 177 313267 0 400 152 113 60720 0 178 316176 1 71 1 0 60745 60745 179 315380 0 0 0 0 60720 0 180 315380 0 0 0 0 60720 0 181 359335 0 972 34 19 71642 2435828 182 330068 1 531 10 11 71641 716410 183 314289 0 2461 57 25 55792 2120096 184 297413 1 378 52 16 71873 718730 185 314806 1 23 5 5 62555 312775 186 333210 1 638 14 11 60370 845180 187 352108 1 2300 29 23 64873 1037968 188 313332 0 149 5 6 62041 310205 189 291787 0 226 5 5 65745 328725 190 315380 0 0 0 0 60720 0 191 318745 0 275 4 7 59500 238000 192 315380 0 0 0 0 60720 0 193 315366 0 141 6 7 61630 369780 194 315380 0 0 0 0 60720 0 195 315688 0 28 2 3 60890 121780 196 315380 1 0 0 0 60720 0 197 409642 1 4980 91 89 113521 9081680 198 315380 1 0 0 0 60720 0 199 315380 1 0 0 0 60720 0 200 269587 1 472 20 19 80045 1600900 201 315380 1 0 0 0 60720 0 202 315380 1 0 0 0 60720 0 203 315380 1 0 0 0 60720 0 204 300962 1 203 27 12 50804 1371708 205 325479 1 496 17 12 87390 1485630 206 316155 1 10 2 5 61656 123312 207 318574 1 63 4 2 65688 262752 208 315380 1 0 0 0 60720 0 209 343613 1 1136 32 26 48522 1552704 210 306948 1 265 31 3 60720 0 211 315380 0 0 0 0 60720 0 212 315380 0 0 0 0 60720 0 213 330059 0 267 32 11 57640 1844480 214 288985 0 474 20 10 61977 1239540 215 304485 0 534 7 5 62620 438340 216 315380 0 0 0 2 60720 0 217 315688 0 15 8 6 60831 486648 218 317736 0 397 28 7 60646 1698088 219 315380 0 0 0 2 60720 0 220 322331 0 1866 29 28 56225 1124500 221 296656 0 288 4 3 60510 242040 222 315380 0 0 0 0 60720 0 223 315354 0 3 2 1 60698 121396 224 312161 0 468 21 20 60720 0 225 315576 0 20 2 1 60805 121610 226 314922 0 278 26 22 61404 1596504 227 314551 0 61 14 9 60720 0 228 315380 0 0 0 0 60720 0 229 312339 0 192 4 2 65276 261104 230 315380 0 0 0 0 60720 0 231 298700 0 317 9 7 63915 575235 232 321376 0 738 10 9 60720 0 233 315380 0 0 0 0 60720 0 234 303230 0 368 17 13 61686 1048662 235 315380 0 0 0 0 60720 0 236 315487 0 2 1 0 60743 60743 237 315380 0 0 0 0 60720 0 238 315793 0 53 6 6 60349 362094 239 315380 0 0 0 0 60720 0 240 315380 0 0 0 0 60720 0 241 315380 0 0 0 0 60720 0 242 312887 0 94 3 3 61360 184080 243 315380 0 0 0 0 60720 0 244 315637 0 24 8 7 59818 478544 245 324385 0 2332 4 2 72680 290720 246 315380 1 0 0 0 60720 0 247 315380 1 0 0 0 60720 0 248 308989 0 131 11 15 61808 679888 249 315380 1 0 0 0 60720 0 250 315380 1 0 0 0 60720 0 251 296702 0 206 9 9 53110 477990 252 315380 1 0 0 0 60720 0 253 307322 0 167 2 1 64245 128490 254 304376 0 622 73 38 73007 5329511 255 253588 0 2328 94 57 82732 7032220 256 315380 0 0 0 0 60720 0 257 309560 0 365 8 7 54820 438560 258 298466 0 364 35 26 47705 1669675 259 315380 1 0 0 0 60720 0 260 315380 0 0 0 0 60720 0 261 315380 1 0 0 0 60720 0 262 315380 1 0 0 0 60720 0 263 343929 0 226 12 13 72835 874020 264 331955 0 307 15 10 58856 882840 265 315380 1 0 0 0 60720 0 266 315380 0 0 0 0 60720 0 267 315380 1 0 0 0 60720 0 268 381180 0 188 11 9 77655 854205 269 315380 1 0 0 0 60720 0 270 331420 0 138 6 26 69817 418902 271 315380 1 0 0 0 60720 0 272 315380 1 0 0 0 60720 0 273 315380 1 0 0 0 60720 0 274 310201 0 125 12 19 60798 729576 275 315380 0 0 0 0 60720 0 276 320016 0 282 30 12 62452 1873560 277 320398 0 335 33 23 64175 2117775 278 315380 0 0 0 0 60720 0 279 291841 1 1324 117 29 67440 5530080 280 310670 0 176 28 8 68136 1907808 281 315380 1 0 0 0 60720 0 282 315380 0 0 0 0 60720 0 283 313491 0 249 72 26 56726 4084272 284 315380 0 0 0 0 60720 0 285 331323 0 333 13 9 70811 920543 286 315380 0 0 0 0 60720 0 287 319210 1 601 6 5 60720 0 288 318098 0 30 4 3 62045 248180 289 315380 1 0 0 0 60720 0 290 292754 0 249 62 13 54323 3368026 291 315380 1 0 0 0 60720 0 292 325176 0 165 24 12 62841 1508184 293 365959 0 453 21 19 81125 1703625 294 315380 0 0 0 0 60720 0 295 302409 0 53 14 10 59506 833084 296 340968 1 382 21 9 59365 1246665 297 315380 0 0 0 0 60720 0 298 315380 1 0 0 0 60720 0 299 315380 0 0 0 0 60720 0 300 315380 0 0 0 9 60720 0 301 313164 1 30 4 4 60798 243192 302 301164 0 290 2 1 58790 117580 303 315380 0 0 0 1 60720 0 304 315380 1 0 0 0 60720 0 305 344425 0 366 53 14 61808 3275824 306 315394 0 2 9 12 60735 546615 307 315380 1 0 0 0 60720 0 308 316647 1 209 13 19 64016 832208 309 309836 0 384 22 17 54683 1203026 310 315380 1 0 0 0 60720 0 311 315380 1 0 0 0 60720 0 312 346611 0 365 83 32 87192 7236936 313 315380 1 0 0 0 60720 0 314 322031 1 49 8 14 64107 512856 315 315656 0 3 4 8 60761 243044 316 339445 0 133 14 4 65990 923860 317 314964 0 32 1 0 59988 59988 318 297141 0 368 17 20 61167 1039839 319 315372 0 1 6 5 60719 364314 320 315380 0 0 0 0 60720 0 321 315380 0 0 0 0 60720 0 322 315380 0 0 0 0 60720 0 323 315380 0 0 0 0 60720 0 324 315380 1 0 0 0 60720 0 325 315380 0 0 0 0 60720 0 326 312502 0 22 2 1 60722 121444 327 315380 1 0 0 0 60720 0 328 315380 1 0 0 0 60720 0 329 315380 0 0 0 0 60720 0 330 315380 0 0 0 0 60720 0 331 315380 0 0 0 0 60720 0 332 315380 0 0 0 0 60720 0 333 315380 0 0 0 0 60720 0 334 313729 0 96 5 4 60379 301895 335 315388 0 1 2 1 60727 121454 336 315371 0 314 5 4 60720 0 337 296139 0 844 78 20 60925 426475 338 315380 0 0 0 0 60720 0 339 313880 0 26 1 1 60896 60896 340 317698 0 125 13 10 59734 776542 341 295580 0 304 15 12 62969 944535 342 315380 0 0 0 0 60720 0 343 315380 0 0 0 0 60720 0 344 315380 0 0 0 0 60720 0 345 308256 0 621 48 13 60720 0 346 315380 0 0 0 0 60720 0 347 303677 0 119 6 3 59118 354708 348 315380 0 0 0 0 60720 0 349 315380 0 0 0 0 60720 0 350 319369 0 1595 17 10 60720 0 351 318690 0 312 14 3 58598 820372 352 314049 0 60 10 7 61124 611240 353 325699 0 587 12 10 59595 715140 354 314210 0 135 2 1 62065 124130 355 315380 0 0 0 0 60720 0 356 315380 0 0 0 0 60720 0 357 322378 0 514 52 15 78780 4096560 358 315380 0 0 0 0 60720 0 359 315380 0 0 0 0 60720 0 360 315380 0 0 0 0 60720 0 361 315398 0 1 4 4 60722 242888 362 315380 0 0 0 0 60720 0 363 315380 0 0 0 0 60720 0 364 308336 1 1763 24 28 61600 184800 365 316386 0 180 11 9 59635 655985 366 315380 0 0 0 0 60720 0 367 315380 0 0 0 0 60720 0 368 315380 0 0 0 0 60720 0 369 315380 0 0 0 0 60720 0 370 315553 0 218 21 7 60720 0 371 315380 0 0 0 0 60720 0 372 323361 0 448 40 7 59781 2391240 373 336639 0 227 9 7 76644 689796 374 307424 0 174 1 3 64820 64820 375 315380 0 0 0 0 60720 0 376 315380 0 0 0 0 60720 0 377 295370 0 121 24 11 56178 1348272 378 322340 0 607 11 7 60436 664796 379 319864 0 2212 14 10 60720 0 380 315380 0 0 0 0 60720 0 381 315380 0 0 0 0 60720 0 382 317291 0 530 60 18 73433 4405980 383 280398 0 571 80 14 41477 3318160 384 315380 0 0 0 0 60720 0 385 317330 0 78 16 12 62700 1003200 386 238125 0 2489 40 29 67804 2712160 387 327071 0 131 6 3 59661 357966 388 309038 0 923 8 6 58620 468960 389 314210 0 72 3 3 60398 181194 390 307930 0 572 16 8 58580 937280 391 322327 0 397 10 10 62710 627100 392 292136 0 450 8 6 59325 474600 393 263276 0 622 7 8 60950 426650 394 367655 0 694 8 6 68060 544480 395 283910 1 3425 12 9 83620 1003440 396 283587 0 562 13 8 58456 759928 397 243650 0 4917 42 26 52811 2218062 398 438493 1 1442 118 239 121173 14298414 399 296261 0 529 9 7 63870 574830 400 230621 1 2126 138 41 21001 2898138 401 304252 0 1061 5 3 70415 352075 402 333505 0 776 9 8 64230 578070 403 296919 0 611 8 6 59190 473520 404 278990 1 1526 25 21 69351 1733775 405 276898 0 592 7 7 64270 449890 406 327007 0 1182 13 11 70694 919022 407 317046 0 621 16 11 68005 1088080 408 304555 0 989 11 12 58930 648230 409 298096 0 438 11 9 58320 641520 410 231861 0 726 3 3 69980 209940 411 309422 0 1303 61 57 69863 4261643 412 286963 1 7419 29 21 63255 1518120 413 269753 1 1164 17 15 57320 974440 414 448243 1 3310 33 32 75230 2482590 415 165404 0 1920 15 11 79420 555940 416 204325 0 965 3 2 73490 220470 417 407159 0 3256 66 23 35250 2326500 418 290476 1 1135 17 20 62285 1058845 419 275311 0 1270 26 24 69206 1799356 420 246541 0 661 3 1 65920 197760 421 253468 0 1013 2 1 69770 139540 422 240897 0 2844 67 74 72683 4869761 423 -83265 1 11528 70 68 -14545 -1018150 424 -42143 0 6526 26 20 55830 1451580 425 272713 0 2264 24 20 55174 1324176 426 215362 1 5109 97 82 67038 6301572 427 42754 0 3999 30 21 51252 1537560 428 306275 0 35624 223 244 157278 35072994 429 253537 0 9252 48 32 79510 3816480 430 372631 0 15236 90 86 77440 6969600 431 -7170 0 18073 180 69 27284 4911120 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Group Costs Trades Orders Dividends 3.217e+05 4.216e+04 1.408e+01 -2.725e+03 3.065e+03 -5.856e-01 TrDiv 1.778e-02 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -1188448 -15430 19435 29268 2617776 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 3.217e+05 2.993e+04 10.748 < 2e-16 *** Group 4.216e+04 2.370e+04 1.779 0.0759 . Costs 1.408e+01 1.940e+00 7.259 1.88e-12 *** Trades -2.725e+03 2.277e+02 -11.968 < 2e-16 *** Orders 3.065e+03 3.867e+02 7.926 2.02e-14 *** Dividends -5.856e-01 4.011e-01 -1.460 0.1450 TrDiv 1.778e-02 1.557e-03 11.416 < 2e-16 *** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 214100 on 424 degrees of freedom Multiple R-squared: 0.7562, Adjusted R-squared: 0.7528 F-statistic: 219.2 on 6 and 424 DF, p-value: < 2.2e-16 > if (n > n25) { + kp3 <- k + 3 + nmkm3 <- n - k - 3 + gqarr <- array(NA, dim=c(nmkm3-kp3+1,3)) + numgqtests <- 0 + numsignificant1 <- 0 + numsignificant5 <- 0 + numsignificant10 <- 0 + for (mypoint in kp3:nmkm3) { + j <- 0 + numgqtests <- numgqtests + 1 + for (myalt in c('greater', 'two.sided', 'less')) { + j <- j + 1 + gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value + } + if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1 + } + gqarr + } [,1] [,2] [,3] [1,] 1.0000000 3.908873e-42 1.954437e-42 [2,] 1.0000000 1.044946e-43 5.224728e-44 [3,] 1.0000000 2.145338e-45 1.072669e-45 [4,] 1.0000000 1.390369e-77 6.951846e-78 [5,] 1.0000000 9.215240e-81 4.607620e-81 [6,] 1.0000000 6.074299e-83 3.037149e-83 [7,] 1.0000000 2.638137e-94 1.319068e-94 [8,] 1.0000000 3.102184e-97 1.551092e-97 [9,] 1.0000000 1.055356e-99 5.276782e-100 [10,] 1.0000000 1.122485e-104 5.612427e-105 [11,] 1.0000000 1.189060e-110 5.945298e-111 [12,] 1.0000000 4.807750e-115 2.403875e-115 [13,] 1.0000000 8.390606e-117 4.195303e-117 [14,] 1.0000000 2.660184e-117 1.330092e-117 [15,] 1.0000000 1.186422e-124 5.932112e-125 [16,] 1.0000000 3.104100e-124 1.552050e-124 [17,] 1.0000000 1.138221e-125 5.691105e-126 [18,] 1.0000000 9.545324e-131 4.772662e-131 [19,] 1.0000000 6.637554e-131 3.318777e-131 [20,] 1.0000000 2.294628e-141 1.147314e-141 [21,] 1.0000000 3.684278e-145 1.842139e-145 [22,] 1.0000000 6.822887e-149 3.411443e-149 [23,] 1.0000000 8.966636e-151 4.483318e-151 [24,] 1.0000000 3.247504e-158 1.623752e-158 [25,] 1.0000000 3.243316e-160 1.621658e-160 [26,] 1.0000000 1.026129e-160 5.130646e-161 [27,] 1.0000000 3.089878e-164 1.544939e-164 [28,] 1.0000000 1.301769e-169 6.508847e-170 [29,] 1.0000000 4.695807e-175 2.347903e-175 [30,] 1.0000000 3.984533e-175 1.992267e-175 [31,] 1.0000000 3.919484e-176 1.959742e-176 [32,] 1.0000000 9.203015e-178 4.601508e-178 [33,] 1.0000000 3.881099e-177 1.940550e-177 [34,] 1.0000000 4.157070e-178 2.078535e-178 [35,] 1.0000000 3.561892e-177 1.780946e-177 [36,] 1.0000000 9.179370e-177 4.589685e-177 [37,] 1.0000000 2.859478e-176 1.429739e-176 [38,] 1.0000000 9.835408e-177 4.917704e-177 [39,] 1.0000000 2.389990e-177 1.194995e-177 [40,] 1.0000000 6.622025e-177 3.311012e-177 [41,] 1.0000000 2.786805e-177 1.393402e-177 [42,] 1.0000000 4.739549e-178 2.369775e-178 [43,] 1.0000000 7.221390e-179 3.610695e-179 [44,] 1.0000000 3.636777e-179 1.818389e-179 [45,] 1.0000000 1.009147e-178 5.045737e-179 [46,] 1.0000000 1.072704e-183 5.363522e-184 [47,] 1.0000000 2.443819e-183 1.221909e-183 [48,] 1.0000000 8.313775e-196 4.156888e-196 [49,] 1.0000000 2.982411e-196 1.491206e-196 [50,] 1.0000000 2.097453e-195 1.048727e-195 [51,] 1.0000000 1.061922e-195 5.309611e-196 [52,] 1.0000000 2.564566e-197 1.282283e-197 [53,] 1.0000000 6.161342e-199 3.080671e-199 [54,] 1.0000000 1.766825e-198 8.834127e-199 [55,] 1.0000000 2.406412e-198 1.203206e-198 [56,] 1.0000000 1.271599e-197 6.357996e-198 [57,] 1.0000000 1.195076e-198 5.975380e-199 [58,] 1.0000000 2.632198e-198 1.316099e-198 [59,] 1.0000000 1.915897e-197 9.579484e-198 [60,] 1.0000000 8.405656e-197 4.202828e-197 [61,] 1.0000000 1.539951e-196 7.699756e-197 [62,] 1.0000000 3.866107e-196 1.933054e-196 [63,] 1.0000000 4.167390e-197 2.083695e-197 [64,] 1.0000000 1.307778e-196 6.538890e-197 [65,] 1.0000000 1.117219e-195 5.586095e-196 [66,] 1.0000000 1.033454e-194 5.167271e-195 [67,] 1.0000000 1.564942e-194 7.824709e-195 [68,] 1.0000000 1.079778e-193 5.398891e-194 [69,] 1.0000000 1.303994e-193 6.519971e-194 [70,] 1.0000000 1.017860e-192 5.089298e-193 [71,] 1.0000000 8.121289e-193 4.060645e-193 [72,] 1.0000000 5.123164e-192 2.561582e-192 [73,] 1.0000000 5.731890e-192 2.865945e-192 [74,] 1.0000000 3.135882e-192 1.567941e-192 [75,] 1.0000000 2.381538e-191 1.190769e-191 [76,] 1.0000000 1.892175e-190 9.460876e-191 [77,] 1.0000000 2.074694e-190 1.037347e-190 [78,] 1.0000000 1.230158e-190 6.150791e-191 [79,] 1.0000000 8.278442e-190 4.139221e-190 [80,] 1.0000000 1.177325e-189 5.886625e-190 [81,] 1.0000000 9.189858e-189 4.594929e-189 [82,] 1.0000000 6.997138e-188 3.498569e-188 [83,] 1.0000000 3.132277e-187 1.566138e-187 [84,] 1.0000000 1.744358e-188 8.721792e-189 [85,] 1.0000000 1.202912e-187 6.014562e-188 [86,] 1.0000000 4.958419e-188 2.479210e-188 [87,] 1.0000000 4.372690e-188 2.186345e-188 [88,] 1.0000000 2.439241e-187 1.219621e-187 [89,] 1.0000000 6.159480e-187 3.079740e-187 [90,] 1.0000000 4.680896e-186 2.340448e-186 [91,] 1.0000000 4.144283e-185 2.072141e-185 [92,] 1.0000000 2.724700e-184 1.362350e-184 [93,] 1.0000000 2.513352e-185 1.256676e-185 [94,] 1.0000000 9.711445e-186 4.855722e-186 [95,] 1.0000000 1.863421e-186 9.317103e-187 [96,] 1.0000000 6.935786e-186 3.467893e-186 [97,] 1.0000000 2.343780e-185 1.171890e-185 [98,] 1.0000000 3.870440e-188 1.935220e-188 [99,] 1.0000000 8.919708e-190 4.459854e-190 [100,] 1.0000000 5.935222e-189 2.967611e-189 [101,] 1.0000000 7.442567e-189 3.721284e-189 [102,] 1.0000000 5.301103e-189 2.650551e-189 [103,] 1.0000000 1.180028e-189 5.900140e-190 [104,] 1.0000000 6.496618e-189 3.248309e-189 [105,] 1.0000000 4.546480e-188 2.273240e-188 [106,] 1.0000000 3.638798e-188 1.819399e-188 [107,] 1.0000000 1.374292e-187 6.871458e-188 [108,] 1.0000000 7.320454e-187 3.660227e-187 [109,] 1.0000000 2.781333e-187 1.390667e-187 [110,] 1.0000000 2.654477e-188 1.327239e-188 [111,] 1.0000000 2.067562e-187 1.033781e-187 [112,] 1.0000000 5.752587e-187 2.876293e-187 [113,] 1.0000000 4.714060e-186 2.357030e-186 [114,] 1.0000000 2.202281e-185 1.101140e-185 [115,] 1.0000000 1.740570e-185 8.702851e-186 [116,] 1.0000000 5.710215e-185 2.855107e-185 [117,] 1.0000000 1.815285e-185 9.076425e-186 [118,] 1.0000000 1.622450e-184 8.112252e-185 [119,] 1.0000000 1.509894e-183 7.549469e-184 [120,] 1.0000000 1.354318e-182 6.771590e-183 [121,] 1.0000000 1.200635e-181 6.003176e-182 [122,] 1.0000000 1.090444e-180 5.452220e-181 [123,] 1.0000000 1.599948e-180 7.999739e-181 [124,] 1.0000000 1.418869e-179 7.094347e-180 [125,] 1.0000000 1.317322e-178 6.586611e-179 [126,] 1.0000000 1.145655e-177 5.728273e-178 [127,] 1.0000000 1.002756e-176 5.013780e-177 [128,] 1.0000000 8.831197e-176 4.415599e-176 [129,] 1.0000000 7.678882e-175 3.839441e-175 [130,] 1.0000000 6.237804e-174 3.118902e-174 [131,] 1.0000000 5.392747e-173 2.696373e-173 [132,] 1.0000000 4.076088e-172 2.038044e-172 [133,] 1.0000000 3.499756e-171 1.749878e-171 [134,] 1.0000000 3.110761e-170 1.555380e-170 [135,] 1.0000000 2.673970e-169 1.336985e-169 [136,] 1.0000000 2.098609e-168 1.049304e-168 [137,] 1.0000000 1.724587e-167 8.622933e-168 [138,] 1.0000000 1.449684e-166 7.248418e-167 [139,] 1.0000000 1.238711e-165 6.193557e-166 [140,] 1.0000000 9.131562e-165 4.565781e-165 [141,] 1.0000000 7.573076e-164 3.786538e-164 [142,] 1.0000000 6.302734e-163 3.151367e-163 [143,] 1.0000000 5.412046e-162 2.706023e-162 [144,] 1.0000000 4.483365e-161 2.241683e-161 [145,] 1.0000000 3.755662e-160 1.877831e-160 [146,] 1.0000000 3.067906e-159 1.533953e-159 [147,] 1.0000000 2.480127e-158 1.240064e-158 [148,] 1.0000000 1.996510e-157 9.982548e-158 [149,] 1.0000000 1.612870e-156 8.064349e-157 [150,] 1.0000000 9.553520e-156 4.776760e-156 [151,] 1.0000000 7.597415e-155 3.798707e-155 [152,] 1.0000000 6.015653e-154 3.007827e-154 [153,] 1.0000000 4.742465e-153 2.371232e-153 [154,] 1.0000000 3.722376e-152 1.861188e-152 [155,] 1.0000000 2.857934e-151 1.428967e-151 [156,] 1.0000000 1.893190e-150 9.465952e-151 [157,] 1.0000000 1.462255e-149 7.311275e-150 [158,] 1.0000000 1.129286e-148 5.646429e-149 [159,] 1.0000000 1.521012e-149 7.605062e-150 [160,] 1.0000000 4.066609e-154 2.033305e-154 [161,] 1.0000000 3.332437e-153 1.666218e-153 [162,] 1.0000000 2.718575e-152 1.359288e-152 [163,] 1.0000000 2.238306e-151 1.119153e-151 [164,] 1.0000000 1.809575e-150 9.047876e-151 [165,] 1.0000000 1.487498e-149 7.437492e-150 [166,] 1.0000000 1.191570e-148 5.957848e-149 [167,] 1.0000000 9.534674e-148 4.767337e-148 [168,] 1.0000000 5.968665e-148 2.984332e-148 [169,] 1.0000000 4.844685e-147 2.422342e-147 [170,] 1.0000000 3.860345e-146 1.930173e-146 [171,] 1.0000000 3.062135e-145 1.531067e-145 [172,] 1.0000000 1.552456e-144 7.762282e-145 [173,] 1.0000000 1.224564e-143 6.122821e-144 [174,] 1.0000000 4.691638e-143 2.345819e-143 [175,] 1.0000000 3.631547e-142 1.815774e-142 [176,] 1.0000000 2.824773e-141 1.412386e-141 [177,] 1.0000000 1.859185e-140 9.295926e-141 [178,] 1.0000000 3.590788e-140 1.795394e-140 [179,] 1.0000000 2.809971e-139 1.404986e-139 [180,] 1.0000000 1.845239e-138 9.226197e-139 [181,] 1.0000000 1.414216e-137 7.071081e-138 [182,] 1.0000000 1.030555e-136 5.152775e-137 [183,] 1.0000000 7.826531e-136 3.913266e-136 [184,] 1.0000000 5.949083e-135 2.974542e-135 [185,] 1.0000000 4.476732e-134 2.238366e-134 [186,] 1.0000000 3.358439e-133 1.679220e-133 [187,] 1.0000000 2.533453e-132 1.266726e-132 [188,] 1.0000000 9.281480e-132 4.640740e-132 [189,] 1.0000000 6.968238e-131 3.484119e-131 [190,] 1.0000000 5.207917e-130 2.603959e-130 [191,] 1.0000000 7.281771e-130 3.640886e-130 [192,] 1.0000000 5.432854e-129 2.716427e-129 [193,] 1.0000000 4.034670e-128 2.017335e-128 [194,] 1.0000000 2.982435e-127 1.491218e-127 [195,] 1.0000000 2.184677e-126 1.092338e-126 [196,] 1.0000000 1.231170e-125 6.155851e-126 [197,] 1.0000000 8.941063e-125 4.470531e-125 [198,] 1.0000000 6.432923e-124 3.216462e-124 [199,] 1.0000000 4.647641e-123 2.323821e-123 [200,] 1.0000000 1.060691e-122 5.303455e-123 [201,] 1.0000000 7.377447e-122 3.688723e-122 [202,] 1.0000000 5.217086e-121 2.608543e-121 [203,] 1.0000000 3.672792e-120 1.836396e-120 [204,] 1.0000000 2.162860e-119 1.081430e-119 [205,] 1.0000000 1.346674e-118 6.733372e-119 [206,] 1.0000000 9.478394e-118 4.739197e-118 [207,] 1.0000000 6.572324e-117 3.286162e-117 [208,] 1.0000000 4.544835e-116 2.272418e-116 [209,] 1.0000000 3.058214e-115 1.529107e-115 [210,] 1.0000000 2.091731e-114 1.045866e-114 [211,] 1.0000000 7.033129e-114 3.516565e-114 [212,] 1.0000000 4.702885e-113 2.351442e-113 [213,] 1.0000000 3.178368e-112 1.589184e-112 [214,] 1.0000000 2.139505e-111 1.069752e-111 [215,] 1.0000000 1.220844e-110 6.104218e-111 [216,] 1.0000000 8.140874e-110 4.070437e-110 [217,] 1.0000000 5.419023e-109 2.709511e-109 [218,] 1.0000000 3.366314e-108 1.683157e-108 [219,] 1.0000000 2.214464e-107 1.107232e-107 [220,] 1.0000000 1.461147e-106 7.305736e-107 [221,] 1.0000000 9.524529e-106 4.762264e-106 [222,] 1.0000000 5.988260e-105 2.994130e-105 [223,] 1.0000000 3.130107e-104 1.565053e-104 [224,] 1.0000000 2.014392e-103 1.007196e-103 [225,] 1.0000000 1.303920e-102 6.519602e-103 [226,] 1.0000000 8.313912e-102 4.156956e-102 [227,] 1.0000000 5.275470e-101 2.637735e-101 [228,] 1.0000000 3.332327e-100 1.666163e-100 [229,] 1.0000000 2.087629e-99 1.043814e-99 [230,] 1.0000000 1.306239e-98 6.531196e-99 [231,] 1.0000000 8.134663e-98 4.067331e-98 [232,] 1.0000000 5.041957e-97 2.520979e-97 [233,] 1.0000000 3.136886e-96 1.568443e-96 [234,] 1.0000000 1.925868e-95 9.629338e-96 [235,] 1.0000000 1.171682e-94 5.858409e-95 [236,] 1.0000000 6.229775e-94 3.114888e-94 [237,] 1.0000000 3.818448e-93 1.909224e-93 [238,] 1.0000000 2.329318e-92 1.164659e-92 [239,] 1.0000000 1.417569e-91 7.087847e-92 [240,] 1.0000000 8.564235e-91 4.282117e-91 [241,] 1.0000000 5.149227e-90 2.574613e-90 [242,] 1.0000000 3.086889e-89 1.543444e-89 [243,] 1.0000000 1.838023e-88 9.190116e-89 [244,] 1.0000000 1.081284e-87 5.406419e-88 [245,] 1.0000000 4.510494e-87 2.255247e-87 [246,] 1.0000000 2.888843e-87 1.444422e-87 [247,] 1.0000000 1.691541e-86 8.457707e-87 [248,] 1.0000000 9.506920e-86 4.753460e-86 [249,] 1.0000000 5.245745e-85 2.622872e-85 [250,] 1.0000000 3.070208e-84 1.535104e-84 [251,] 1.0000000 1.762620e-83 8.813099e-84 [252,] 1.0000000 1.021562e-82 5.107811e-83 [253,] 1.0000000 5.891581e-82 2.945790e-82 [254,] 1.0000000 3.125666e-81 1.562833e-81 [255,] 1.0000000 1.463538e-80 7.317690e-81 [256,] 1.0000000 8.329938e-80 4.164969e-80 [257,] 1.0000000 4.646510e-79 2.323255e-79 [258,] 1.0000000 2.618481e-78 1.309240e-78 [259,] 1.0000000 8.576041e-78 4.288020e-78 [260,] 1.0000000 4.801570e-77 2.400785e-77 [261,] 1.0000000 2.438888e-76 1.219444e-76 [262,] 1.0000000 1.352480e-75 6.762402e-76 [263,] 1.0000000 7.462130e-75 3.731065e-75 [264,] 1.0000000 4.096175e-74 2.048087e-74 [265,] 1.0000000 2.214539e-73 1.107270e-73 [266,] 1.0000000 1.184138e-72 5.920690e-73 [267,] 1.0000000 6.333426e-72 3.166713e-72 [268,] 1.0000000 3.390794e-71 1.695397e-71 [269,] 1.0000000 1.785530e-70 8.927650e-71 [270,] 1.0000000 6.201449e-70 3.100724e-70 [271,] 1.0000000 2.995341e-69 1.497670e-69 [272,] 1.0000000 1.581789e-68 7.908944e-69 [273,] 1.0000000 8.174506e-68 4.087253e-68 [274,] 1.0000000 4.158217e-67 2.079109e-67 [275,] 1.0000000 2.127296e-66 1.063648e-66 [276,] 1.0000000 1.080600e-65 5.403002e-66 [277,] 1.0000000 5.470083e-65 2.735042e-65 [278,] 1.0000000 2.698814e-64 1.349407e-64 [279,] 1.0000000 1.350834e-63 6.754172e-64 [280,] 1.0000000 6.848623e-63 3.424312e-63 [281,] 1.0000000 2.968749e-62 1.484375e-62 [282,] 1.0000000 1.490526e-61 7.452632e-62 [283,] 1.0000000 7.306993e-61 3.653496e-61 [284,] 1.0000000 3.128628e-60 1.564314e-60 [285,] 1.0000000 1.521587e-59 7.607937e-60 [286,] 1.0000000 7.474342e-59 3.737171e-59 [287,] 1.0000000 3.293209e-58 1.646604e-58 [288,] 1.0000000 1.577055e-57 7.885275e-58 [289,] 1.0000000 7.642456e-57 3.821228e-57 [290,] 1.0000000 3.619388e-56 1.809694e-56 [291,] 1.0000000 1.672293e-55 8.361466e-56 [292,] 1.0000000 7.941475e-55 3.970737e-55 [293,] 1.0000000 3.768594e-54 1.884297e-54 [294,] 1.0000000 1.742584e-53 8.712922e-54 [295,] 1.0000000 8.166366e-53 4.083183e-53 [296,] 1.0000000 3.502691e-52 1.751345e-52 [297,] 1.0000000 1.592927e-51 7.964636e-52 [298,] 1.0000000 7.342702e-51 3.671351e-51 [299,] 1.0000000 3.323503e-50 1.661752e-50 [300,] 1.0000000 1.448429e-49 7.242144e-50 [301,] 1.0000000 6.563139e-49 3.281570e-49 [302,] 1.0000000 2.956061e-48 1.478030e-48 [303,] 1.0000000 6.519018e-48 3.259509e-48 [304,] 1.0000000 2.919433e-47 1.459717e-47 [305,] 1.0000000 1.292070e-46 6.460349e-47 [306,] 1.0000000 5.576111e-46 2.788056e-46 [307,] 1.0000000 2.288733e-45 1.144366e-45 [308,] 1.0000000 9.816873e-45 4.908437e-45 [309,] 1.0000000 4.252258e-44 2.126129e-44 [310,] 1.0000000 1.813617e-43 9.068084e-44 [311,] 1.0000000 7.648296e-43 3.824148e-43 [312,] 1.0000000 3.204841e-42 1.602421e-42 [313,] 1.0000000 1.334292e-41 6.671458e-42 [314,] 1.0000000 5.519206e-41 2.759603e-41 [315,] 1.0000000 2.313664e-40 1.156832e-40 [316,] 1.0000000 9.448908e-40 4.724454e-40 [317,] 1.0000000 3.874082e-39 1.937041e-39 [318,] 1.0000000 1.592963e-38 7.964815e-39 [319,] 1.0000000 6.505751e-38 3.252875e-38 [320,] 1.0000000 2.588829e-37 1.294415e-37 [321,] 1.0000000 1.023169e-36 5.115845e-37 [322,] 1.0000000 4.016102e-36 2.008051e-36 [323,] 1.0000000 1.565488e-35 7.827442e-36 [324,] 1.0000000 6.059778e-35 3.029889e-35 [325,] 1.0000000 2.339080e-34 1.169540e-34 [326,] 1.0000000 8.952955e-34 4.476477e-34 [327,] 1.0000000 3.284060e-33 1.642030e-33 [328,] 1.0000000 1.128166e-32 5.640828e-33 [329,] 1.0000000 4.225197e-32 2.112598e-32 [330,] 1.0000000 1.579129e-31 7.895643e-32 [331,] 1.0000000 5.779044e-31 2.889522e-31 [332,] 1.0000000 2.138975e-30 1.069488e-30 [333,] 1.0000000 7.783977e-30 3.891988e-30 [334,] 1.0000000 2.811085e-29 1.405542e-29 [335,] 1.0000000 1.007367e-28 5.036837e-29 [336,] 1.0000000 3.170151e-28 1.585076e-28 [337,] 1.0000000 1.121482e-27 5.607408e-28 [338,] 1.0000000 4.026490e-27 2.013245e-27 [339,] 1.0000000 1.402266e-26 7.011330e-27 [340,] 1.0000000 4.843943e-26 2.421971e-26 [341,] 1.0000000 1.277061e-25 6.385307e-26 [342,] 1.0000000 4.305731e-25 2.152866e-25 [343,] 1.0000000 1.474220e-24 7.371099e-25 [344,] 1.0000000 4.487394e-24 2.243697e-24 [345,] 1.0000000 1.503227e-23 7.516136e-24 [346,] 1.0000000 4.951910e-23 2.475955e-23 [347,] 1.0000000 1.616816e-22 8.084078e-23 [348,] 1.0000000 4.367584e-22 2.183792e-22 [349,] 1.0000000 1.407272e-21 7.036359e-22 [350,] 1.0000000 4.492721e-21 2.246361e-21 [351,] 1.0000000 1.420939e-20 7.104693e-21 [352,] 1.0000000 4.461536e-20 2.230768e-20 [353,] 1.0000000 1.384384e-19 6.921920e-20 [354,] 1.0000000 4.253691e-19 2.126845e-19 [355,] 1.0000000 1.183508e-18 5.917541e-19 [356,] 1.0000000 3.546460e-18 1.773230e-18 [357,] 1.0000000 1.060546e-17 5.302731e-18 [358,] 1.0000000 3.138536e-17 1.569268e-17 [359,] 1.0000000 9.189677e-17 4.594838e-17 [360,] 1.0000000 2.661682e-16 1.330841e-16 [361,] 1.0000000 7.332853e-16 3.666426e-16 [362,] 1.0000000 2.080129e-15 1.040065e-15 [363,] 1.0000000 6.049051e-15 3.024526e-15 [364,] 1.0000000 1.724903e-14 8.624513e-15 [365,] 1.0000000 4.908203e-14 2.454102e-14 [366,] 1.0000000 1.333673e-13 6.668367e-14 [367,] 1.0000000 3.579477e-13 1.789738e-13 [368,] 1.0000000 1.000170e-12 5.000849e-13 [369,] 1.0000000 2.498746e-12 1.249373e-12 [370,] 1.0000000 4.887412e-12 2.443706e-12 [371,] 1.0000000 1.252272e-11 6.261358e-12 [372,] 1.0000000 3.163682e-11 1.581841e-11 [373,] 1.0000000 7.443512e-11 3.721756e-11 [374,] 1.0000000 1.913755e-10 9.568776e-11 [375,] 1.0000000 4.700726e-10 2.350363e-10 [376,] 1.0000000 1.175939e-09 5.879697e-10 [377,] 1.0000000 2.553986e-09 1.276993e-09 [378,] 1.0000000 5.768896e-09 2.884448e-09 [379,] 1.0000000 1.310138e-08 6.550692e-09 [380,] 1.0000000 3.016517e-08 1.508258e-08 [381,] 1.0000000 7.001368e-08 3.500684e-08 [382,] 0.9999999 1.521280e-07 7.606399e-08 [383,] 0.9999998 3.528164e-07 1.764082e-07 [384,] 0.9999996 8.262271e-07 4.131135e-07 [385,] 0.9999993 1.347806e-06 6.739028e-07 [386,] 0.9999985 2.914259e-06 1.457130e-06 [387,] 0.9999968 6.495370e-06 3.247685e-06 [388,] 0.9999932 1.356252e-05 6.781259e-06 [389,] 0.9999881 2.378738e-05 1.189369e-05 [390,] 0.9999747 5.050905e-05 2.525452e-05 [391,] 0.9999476 1.047618e-04 5.238090e-05 [392,] 0.9998931 2.138492e-04 1.069246e-04 [393,] 0.9998151 3.698316e-04 1.849158e-04 [394,] 0.9996518 6.963006e-04 3.481503e-04 [395,] 0.9993570 1.285937e-03 6.429687e-04 [396,] 0.9987532 2.493687e-03 1.246843e-03 [397,] 0.9978238 4.352468e-03 2.176234e-03 [398,] 0.9961317 7.736677e-03 3.868338e-03 [399,] 0.9939039 1.219229e-02 6.096145e-03 [400,] 0.9902459 1.950828e-02 9.754138e-03 [401,] 0.9827845 3.443101e-02 1.721550e-02 [402,] 0.9709392 5.812166e-02 2.906083e-02 [403,] 0.9537908 9.241842e-02 4.620921e-02 [404,] 0.9250910 1.498180e-01 7.490902e-02 [405,] 0.9238618 1.522764e-01 7.613821e-02 [406,] 0.8936620 2.126759e-01 1.063380e-01 [407,] 0.8401800 3.196400e-01 1.598200e-01 [408,] 0.9696029 6.079416e-02 3.039708e-02 [409,] 0.9434311 1.131379e-01 5.656895e-02 [410,] 0.8923107 2.153785e-01 1.076893e-01 [411,] 0.8148940 3.702120e-01 1.851060e-01 [412,] 0.6874003 6.251993e-01 3.125997e-01 > postscript(file="/var/www/html/rcomp/tmp/1lsvn1291394141.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index') > points(x[,1]-mysum$resid) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/2lsvn1291394141.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/3w1u81291394141.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/4w1u81291394141.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals') > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5w1u81291394141.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 131691.2321 2617776.3344 1343261.9381 -744280.2782 -384882.4697 6 7 8 9 10 231829.6052 -24799.1285 -140602.1496 171701.7832 687631.7910 11 12 13 14 15 24217.3337 344779.5263 240841.5619 -741217.9665 -52513.3362 16 17 18 19 20 708793.0810 -188311.0806 -167054.1605 -86916.7805 -703145.6635 21 22 23 24 25 285123.3410 -15121.8885 -47851.8834 -818053.2764 -264218.2790 26 27 28 29 30 8784.4069 -362208.6164 -60845.5727 338692.8381 182219.7224 31 32 33 34 35 281788.1807 -36995.6012 263445.7065 145992.6598 373117.5424 36 37 38 39 40 -651360.7695 102982.2074 137960.6316 -42632.1754 -336872.4361 41 42 43 44 45 -219519.1761 -222531.9421 137566.0210 -172492.4321 14357.7386 46 47 48 49 50 -30779.2236 405337.5305 -321471.9194 18934.8649 -195140.8815 51 52 53 54 55 -3391.6077 -192350.0646 -234497.9226 60187.4370 -462598.5773 56 57 58 59 60 9233.5964 306052.0936 33188.9153 5900.7235 -198482.0380 61 62 63 64 65 23938.1649 -196629.2342 20928.2372 -179696.4179 -2681.0677 66 67 68 69 70 185234.2138 23891.0866 -21817.5531 -76648.7399 -138163.7547 71 72 73 74 75 -208938.8859 -9289.8180 -47672.1817 29587.7757 31861.7253 76 77 78 79 80 -126419.0457 42686.2541 -26071.1896 28375.2746 103544.1770 81 82 83 84 85 -4214.0980 -65473.3759 14124.1325 -17642.2383 76278.7647 86 87 88 89 90 -143459.0879 -71856.1787 58271.2104 -64948.5081 83938.2662 91 92 93 94 95 42507.8170 -49580.7887 85793.1811 -70344.4744 -6814.0289 96 97 98 99 100 -15737.3701 147930.4943 213063.0816 72826.2490 -34459.7025 101 102 103 104 105 -79216.1575 -79815.2978 -111110.0897 39046.6756 52007.6985 106 107 108 109 110 -76065.7715 -295762.5098 146195.4603 -104872.2276 -108290.3137 111 112 113 114 115 91917.1948 -238787.2513 -101375.4457 -84835.3560 25282.2672 116 117 118 119 120 -6420.5109 62241.0108 113622.0525 -147931.6276 66650.5828 121 122 123 124 125 -242477.7777 24226.5214 -155730.5335 103139.2098 -3850.1514 126 127 128 129 130 -187428.7167 29268.1262 8003.2527 1350.2192 29268.1262 131 132 133 134 135 19435.3611 91186.7223 29268.1262 -30144.0315 18248.3820 136 137 138 139 140 26910.1165 21999.4275 31090.4241 17538.2510 29268.1262 141 142 143 144 145 31858.1740 29268.1262 -11943.4238 21985.8429 16747.6812 146 147 148 149 150 11303.9014 29268.1262 17224.4191 34642.2003 31382.0373 151 152 153 154 155 -7622.8240 -25650.9295 13943.2545 1038.1341 27775.4888 156 157 158 159 160 29268.1262 29268.1262 25329.7892 40865.2598 29268.1262 161 162 163 164 165 29268.1262 29268.1262 29268.1262 36353.2418 19251.7140 166 167 168 169 170 33219.4047 29268.1262 58320.1911 892399.6433 29268.1262 171 172 173 174 175 29268.1262 23616.3821 29268.1262 8373.1483 29268.1262 176 177 178 179 180 21862.0309 89432.1141 -11439.4325 29268.1262 29268.1262 181 182 183 184 185 57060.1061 -18485.6793 31673.1233 50247.3889 -19977.3842 186 187 188 189 190 -14838.0182 -16030.6275 15618.3139 -2106.0886 29268.1262 191 192 193 194 195 13262.5146 29268.1262 16125.6181 29268.1262 23372.3899 196 197 198 199 200 -12895.9151 -144048.2459 -12895.9151 -12895.9151 -86204.1344 201 202 203 204 205 -12895.9151 -12895.9151 -12895.9151 -23558.2587 -11021.0744 206 207 208 209 210 -23779.8406 -7579.1442 -12895.9151 -27880.4863 50231.4022 211 212 213 214 215 29268.1262 29268.1262 59092.2205 -1242.1168 7927.4901 216 217 218 219 220 23138.1775 24191.9708 50659.9525 23138.1775 -19460.1087 221 222 223 224 225 3769.9938 29268.1262 29414.7303 15392.3051 29456.2253 226 227 228 229 230 345.8293 38150.1854 29268.1262 26321.6280 29268.1262 231 232 233 234 235 2843.3897 24541.7601 29268.1262 347.1197 29268.1262 236 237 238 239 240 31005.9853 29268.1262 20243.1219 29268.1262 29268.1262 241 242 243 244 245 29268.1262 21535.2133 29268.1262 20500.1432 12045.8516 246 247 248 249 250 -12895.9151 -12895.9151 -6412.0358 -12895.9151 -12895.9151 251 252 253 254 255 -8320.2368 -12895.9151 21024.5792 4443.0129 -95940.8109 256 257 258 259 260 29268.1262 7405.9368 -14375.4965 -12895.9151 29268.1262 261 262 263 264 265 -12895.9151 -12895.9151 39051.9120 34965.7507 -12895.9151 266 267 268 269 270 29268.1262 -12895.9151 89547.2459 -12895.9151 -22091.7116 271 272 273 274 275 -12895.9151 -12895.9151 -12895.9151 -16124.8231 29268.1262 276 277 278 279 280 42623.2566 13388.0309 29268.1262 80534.1352 44298.5890 281 282 283 284 285 -12895.9151 29268.1262 65466.0576 29268.1262 37912.3754 286 287 288 289 290 29268.1262 -16500.7829 29634.3242 -12895.9151 68645.4134 291 292 293 294 295 -12895.9151 39801.4098 54131.1182 29268.1262 7535.8105 296 297 298 299 300 14006.3083 29268.1262 -12895.9151 29268.1262 1683.3572 301 302 303 304 305 -21170.2426 10134.3494 26203.1519 -12895.9151 97098.1523 306 307 308 309 310 7294.2924 -12895.9151 -50240.2563 1249.7813 -12895.9151 311 312 313 314 315 -12895.9151 70339.2670 -12895.9151 -35175.0997 11587.0227 316 317 318 319 320 64018.5109 29631.8963 -27343.7769 23796.4518 29268.1262 321 322 323 324 325 29268.1262 29268.1262 29268.1262 -12895.9151 29268.1262 326 327 328 329 330 26308.4131 -12895.9151 -12895.9151 29268.1262 29268.1262 331 332 333 334 335 29268.1262 29268.1262 29268.1262 22065.9821 29492.8409 336 337 338 339 340 26204.8240 141959.3090 29268.1262 26082.9678 20224.3425 341 342 343 344 345 -6185.1728 29268.1262 29268.1262 29268.1262 104371.9698 346 347 348 349 350 29268.1262 15803.2182 29268.1262 29268.1262 26480.6749 351 352 353 354 355 41319.2054 22261.8646 20005.2039 27164.0740 29268.1262 356 357 358 359 360 29268.1262 62525.8537 29268.1262 29268.1262 29268.1262 361 362 363 364 365 23597.0153 29268.1262 29268.1262 -67943.7223 17837.3355 366 367 368 369 370 29268.1262 29268.1262 29268.1262 29268.1262 62148.9428 371 372 373 374 375 29268.1262 75442.7081 47466.9676 13641.2611 29268.1262 376 377 378 379 380 29268.1262 12620.8166 24221.5966 10112.3784 29268.1262 381 382 383 384 385 29268.1262 61189.8980 91110.8221 29268.1262 20271.8496 386 387 388 389 390 -106966.4255 39288.3148 3777.1431 22655.9444 14935.5040 391 392 393 394 395 17246.8195 -6152.4952 -44485.5778 69803.8539 -91897.9029 396 397 398 399 400 -14362.6556 -120978.1743 -539798.0202 -2599.6976 48066.8103 401 402 403 404 405 7051.7420 28254.8104 -3696.2117 -92768.8288 -25845.1997 406 407 408 409 410 15471.1292 17005.1009 -14854.1938 -4598.1795 -63801.2624 411 412 413 414 415 -73896.2632 -156603.4022 -93868.7126 29585.1004 -139508.1245 416 417 418 419 420 -89769.1764 128309.4869 -86655.2669 -58399.6203 -44237.0641 421 422 423 424 425 -41702.0976 -209029.6603 -617473.3856 -439248.1597 -67953.7316 426 427 428 429 430 -280137.2562 -315144.1159 -1188447.8307 -186944.8306 -260412.3078 431 -375551.4045 > postscript(file="/var/www/html/rcomp/tmp/66tbb1291394141.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 131691.2321 NA 1 2617776.3344 131691.2321 2 1343261.9381 2617776.3344 3 -744280.2782 1343261.9381 4 -384882.4697 -744280.2782 5 231829.6052 -384882.4697 6 -24799.1285 231829.6052 7 -140602.1496 -24799.1285 8 171701.7832 -140602.1496 9 687631.7910 171701.7832 10 24217.3337 687631.7910 11 344779.5263 24217.3337 12 240841.5619 344779.5263 13 -741217.9665 240841.5619 14 -52513.3362 -741217.9665 15 708793.0810 -52513.3362 16 -188311.0806 708793.0810 17 -167054.1605 -188311.0806 18 -86916.7805 -167054.1605 19 -703145.6635 -86916.7805 20 285123.3410 -703145.6635 21 -15121.8885 285123.3410 22 -47851.8834 -15121.8885 23 -818053.2764 -47851.8834 24 -264218.2790 -818053.2764 25 8784.4069 -264218.2790 26 -362208.6164 8784.4069 27 -60845.5727 -362208.6164 28 338692.8381 -60845.5727 29 182219.7224 338692.8381 30 281788.1807 182219.7224 31 -36995.6012 281788.1807 32 263445.7065 -36995.6012 33 145992.6598 263445.7065 34 373117.5424 145992.6598 35 -651360.7695 373117.5424 36 102982.2074 -651360.7695 37 137960.6316 102982.2074 38 -42632.1754 137960.6316 39 -336872.4361 -42632.1754 40 -219519.1761 -336872.4361 41 -222531.9421 -219519.1761 42 137566.0210 -222531.9421 43 -172492.4321 137566.0210 44 14357.7386 -172492.4321 45 -30779.2236 14357.7386 46 405337.5305 -30779.2236 47 -321471.9194 405337.5305 48 18934.8649 -321471.9194 49 -195140.8815 18934.8649 50 -3391.6077 -195140.8815 51 -192350.0646 -3391.6077 52 -234497.9226 -192350.0646 53 60187.4370 -234497.9226 54 -462598.5773 60187.4370 55 9233.5964 -462598.5773 56 306052.0936 9233.5964 57 33188.9153 306052.0936 58 5900.7235 33188.9153 59 -198482.0380 5900.7235 60 23938.1649 -198482.0380 61 -196629.2342 23938.1649 62 20928.2372 -196629.2342 63 -179696.4179 20928.2372 64 -2681.0677 -179696.4179 65 185234.2138 -2681.0677 66 23891.0866 185234.2138 67 -21817.5531 23891.0866 68 -76648.7399 -21817.5531 69 -138163.7547 -76648.7399 70 -208938.8859 -138163.7547 71 -9289.8180 -208938.8859 72 -47672.1817 -9289.8180 73 29587.7757 -47672.1817 74 31861.7253 29587.7757 75 -126419.0457 31861.7253 76 42686.2541 -126419.0457 77 -26071.1896 42686.2541 78 28375.2746 -26071.1896 79 103544.1770 28375.2746 80 -4214.0980 103544.1770 81 -65473.3759 -4214.0980 82 14124.1325 -65473.3759 83 -17642.2383 14124.1325 84 76278.7647 -17642.2383 85 -143459.0879 76278.7647 86 -71856.1787 -143459.0879 87 58271.2104 -71856.1787 88 -64948.5081 58271.2104 89 83938.2662 -64948.5081 90 42507.8170 83938.2662 91 -49580.7887 42507.8170 92 85793.1811 -49580.7887 93 -70344.4744 85793.1811 94 -6814.0289 -70344.4744 95 -15737.3701 -6814.0289 96 147930.4943 -15737.3701 97 213063.0816 147930.4943 98 72826.2490 213063.0816 99 -34459.7025 72826.2490 100 -79216.1575 -34459.7025 101 -79815.2978 -79216.1575 102 -111110.0897 -79815.2978 103 39046.6756 -111110.0897 104 52007.6985 39046.6756 105 -76065.7715 52007.6985 106 -295762.5098 -76065.7715 107 146195.4603 -295762.5098 108 -104872.2276 146195.4603 109 -108290.3137 -104872.2276 110 91917.1948 -108290.3137 111 -238787.2513 91917.1948 112 -101375.4457 -238787.2513 113 -84835.3560 -101375.4457 114 25282.2672 -84835.3560 115 -6420.5109 25282.2672 116 62241.0108 -6420.5109 117 113622.0525 62241.0108 118 -147931.6276 113622.0525 119 66650.5828 -147931.6276 120 -242477.7777 66650.5828 121 24226.5214 -242477.7777 122 -155730.5335 24226.5214 123 103139.2098 -155730.5335 124 -3850.1514 103139.2098 125 -187428.7167 -3850.1514 126 29268.1262 -187428.7167 127 8003.2527 29268.1262 128 1350.2192 8003.2527 129 29268.1262 1350.2192 130 19435.3611 29268.1262 131 91186.7223 19435.3611 132 29268.1262 91186.7223 133 -30144.0315 29268.1262 134 18248.3820 -30144.0315 135 26910.1165 18248.3820 136 21999.4275 26910.1165 137 31090.4241 21999.4275 138 17538.2510 31090.4241 139 29268.1262 17538.2510 140 31858.1740 29268.1262 141 29268.1262 31858.1740 142 -11943.4238 29268.1262 143 21985.8429 -11943.4238 144 16747.6812 21985.8429 145 11303.9014 16747.6812 146 29268.1262 11303.9014 147 17224.4191 29268.1262 148 34642.2003 17224.4191 149 31382.0373 34642.2003 150 -7622.8240 31382.0373 151 -25650.9295 -7622.8240 152 13943.2545 -25650.9295 153 1038.1341 13943.2545 154 27775.4888 1038.1341 155 29268.1262 27775.4888 156 29268.1262 29268.1262 157 25329.7892 29268.1262 158 40865.2598 25329.7892 159 29268.1262 40865.2598 160 29268.1262 29268.1262 161 29268.1262 29268.1262 162 29268.1262 29268.1262 163 36353.2418 29268.1262 164 19251.7140 36353.2418 165 33219.4047 19251.7140 166 29268.1262 33219.4047 167 58320.1911 29268.1262 168 892399.6433 58320.1911 169 29268.1262 892399.6433 170 29268.1262 29268.1262 171 23616.3821 29268.1262 172 29268.1262 23616.3821 173 8373.1483 29268.1262 174 29268.1262 8373.1483 175 21862.0309 29268.1262 176 89432.1141 21862.0309 177 -11439.4325 89432.1141 178 29268.1262 -11439.4325 179 29268.1262 29268.1262 180 57060.1061 29268.1262 181 -18485.6793 57060.1061 182 31673.1233 -18485.6793 183 50247.3889 31673.1233 184 -19977.3842 50247.3889 185 -14838.0182 -19977.3842 186 -16030.6275 -14838.0182 187 15618.3139 -16030.6275 188 -2106.0886 15618.3139 189 29268.1262 -2106.0886 190 13262.5146 29268.1262 191 29268.1262 13262.5146 192 16125.6181 29268.1262 193 29268.1262 16125.6181 194 23372.3899 29268.1262 195 -12895.9151 23372.3899 196 -144048.2459 -12895.9151 197 -12895.9151 -144048.2459 198 -12895.9151 -12895.9151 199 -86204.1344 -12895.9151 200 -12895.9151 -86204.1344 201 -12895.9151 -12895.9151 202 -12895.9151 -12895.9151 203 -23558.2587 -12895.9151 204 -11021.0744 -23558.2587 205 -23779.8406 -11021.0744 206 -7579.1442 -23779.8406 207 -12895.9151 -7579.1442 208 -27880.4863 -12895.9151 209 50231.4022 -27880.4863 210 29268.1262 50231.4022 211 29268.1262 29268.1262 212 59092.2205 29268.1262 213 -1242.1168 59092.2205 214 7927.4901 -1242.1168 215 23138.1775 7927.4901 216 24191.9708 23138.1775 217 50659.9525 24191.9708 218 23138.1775 50659.9525 219 -19460.1087 23138.1775 220 3769.9938 -19460.1087 221 29268.1262 3769.9938 222 29414.7303 29268.1262 223 15392.3051 29414.7303 224 29456.2253 15392.3051 225 345.8293 29456.2253 226 38150.1854 345.8293 227 29268.1262 38150.1854 228 26321.6280 29268.1262 229 29268.1262 26321.6280 230 2843.3897 29268.1262 231 24541.7601 2843.3897 232 29268.1262 24541.7601 233 347.1197 29268.1262 234 29268.1262 347.1197 235 31005.9853 29268.1262 236 29268.1262 31005.9853 237 20243.1219 29268.1262 238 29268.1262 20243.1219 239 29268.1262 29268.1262 240 29268.1262 29268.1262 241 21535.2133 29268.1262 242 29268.1262 21535.2133 243 20500.1432 29268.1262 244 12045.8516 20500.1432 245 -12895.9151 12045.8516 246 -12895.9151 -12895.9151 247 -6412.0358 -12895.9151 248 -12895.9151 -6412.0358 249 -12895.9151 -12895.9151 250 -8320.2368 -12895.9151 251 -12895.9151 -8320.2368 252 21024.5792 -12895.9151 253 4443.0129 21024.5792 254 -95940.8109 4443.0129 255 29268.1262 -95940.8109 256 7405.9368 29268.1262 257 -14375.4965 7405.9368 258 -12895.9151 -14375.4965 259 29268.1262 -12895.9151 260 -12895.9151 29268.1262 261 -12895.9151 -12895.9151 262 39051.9120 -12895.9151 263 34965.7507 39051.9120 264 -12895.9151 34965.7507 265 29268.1262 -12895.9151 266 -12895.9151 29268.1262 267 89547.2459 -12895.9151 268 -12895.9151 89547.2459 269 -22091.7116 -12895.9151 270 -12895.9151 -22091.7116 271 -12895.9151 -12895.9151 272 -12895.9151 -12895.9151 273 -16124.8231 -12895.9151 274 29268.1262 -16124.8231 275 42623.2566 29268.1262 276 13388.0309 42623.2566 277 29268.1262 13388.0309 278 80534.1352 29268.1262 279 44298.5890 80534.1352 280 -12895.9151 44298.5890 281 29268.1262 -12895.9151 282 65466.0576 29268.1262 283 29268.1262 65466.0576 284 37912.3754 29268.1262 285 29268.1262 37912.3754 286 -16500.7829 29268.1262 287 29634.3242 -16500.7829 288 -12895.9151 29634.3242 289 68645.4134 -12895.9151 290 -12895.9151 68645.4134 291 39801.4098 -12895.9151 292 54131.1182 39801.4098 293 29268.1262 54131.1182 294 7535.8105 29268.1262 295 14006.3083 7535.8105 296 29268.1262 14006.3083 297 -12895.9151 29268.1262 298 29268.1262 -12895.9151 299 1683.3572 29268.1262 300 -21170.2426 1683.3572 301 10134.3494 -21170.2426 302 26203.1519 10134.3494 303 -12895.9151 26203.1519 304 97098.1523 -12895.9151 305 7294.2924 97098.1523 306 -12895.9151 7294.2924 307 -50240.2563 -12895.9151 308 1249.7813 -50240.2563 309 -12895.9151 1249.7813 310 -12895.9151 -12895.9151 311 70339.2670 -12895.9151 312 -12895.9151 70339.2670 313 -35175.0997 -12895.9151 314 11587.0227 -35175.0997 315 64018.5109 11587.0227 316 29631.8963 64018.5109 317 -27343.7769 29631.8963 318 23796.4518 -27343.7769 319 29268.1262 23796.4518 320 29268.1262 29268.1262 321 29268.1262 29268.1262 322 29268.1262 29268.1262 323 -12895.9151 29268.1262 324 29268.1262 -12895.9151 325 26308.4131 29268.1262 326 -12895.9151 26308.4131 327 -12895.9151 -12895.9151 328 29268.1262 -12895.9151 329 29268.1262 29268.1262 330 29268.1262 29268.1262 331 29268.1262 29268.1262 332 29268.1262 29268.1262 333 22065.9821 29268.1262 334 29492.8409 22065.9821 335 26204.8240 29492.8409 336 141959.3090 26204.8240 337 29268.1262 141959.3090 338 26082.9678 29268.1262 339 20224.3425 26082.9678 340 -6185.1728 20224.3425 341 29268.1262 -6185.1728 342 29268.1262 29268.1262 343 29268.1262 29268.1262 344 104371.9698 29268.1262 345 29268.1262 104371.9698 346 15803.2182 29268.1262 347 29268.1262 15803.2182 348 29268.1262 29268.1262 349 26480.6749 29268.1262 350 41319.2054 26480.6749 351 22261.8646 41319.2054 352 20005.2039 22261.8646 353 27164.0740 20005.2039 354 29268.1262 27164.0740 355 29268.1262 29268.1262 356 62525.8537 29268.1262 357 29268.1262 62525.8537 358 29268.1262 29268.1262 359 29268.1262 29268.1262 360 23597.0153 29268.1262 361 29268.1262 23597.0153 362 29268.1262 29268.1262 363 -67943.7223 29268.1262 364 17837.3355 -67943.7223 365 29268.1262 17837.3355 366 29268.1262 29268.1262 367 29268.1262 29268.1262 368 29268.1262 29268.1262 369 62148.9428 29268.1262 370 29268.1262 62148.9428 371 75442.7081 29268.1262 372 47466.9676 75442.7081 373 13641.2611 47466.9676 374 29268.1262 13641.2611 375 29268.1262 29268.1262 376 12620.8166 29268.1262 377 24221.5966 12620.8166 378 10112.3784 24221.5966 379 29268.1262 10112.3784 380 29268.1262 29268.1262 381 61189.8980 29268.1262 382 91110.8221 61189.8980 383 29268.1262 91110.8221 384 20271.8496 29268.1262 385 -106966.4255 20271.8496 386 39288.3148 -106966.4255 387 3777.1431 39288.3148 388 22655.9444 3777.1431 389 14935.5040 22655.9444 390 17246.8195 14935.5040 391 -6152.4952 17246.8195 392 -44485.5778 -6152.4952 393 69803.8539 -44485.5778 394 -91897.9029 69803.8539 395 -14362.6556 -91897.9029 396 -120978.1743 -14362.6556 397 -539798.0202 -120978.1743 398 -2599.6976 -539798.0202 399 48066.8103 -2599.6976 400 7051.7420 48066.8103 401 28254.8104 7051.7420 402 -3696.2117 28254.8104 403 -92768.8288 -3696.2117 404 -25845.1997 -92768.8288 405 15471.1292 -25845.1997 406 17005.1009 15471.1292 407 -14854.1938 17005.1009 408 -4598.1795 -14854.1938 409 -63801.2624 -4598.1795 410 -73896.2632 -63801.2624 411 -156603.4022 -73896.2632 412 -93868.7126 -156603.4022 413 29585.1004 -93868.7126 414 -139508.1245 29585.1004 415 -89769.1764 -139508.1245 416 128309.4869 -89769.1764 417 -86655.2669 128309.4869 418 -58399.6203 -86655.2669 419 -44237.0641 -58399.6203 420 -41702.0976 -44237.0641 421 -209029.6603 -41702.0976 422 -617473.3856 -209029.6603 423 -439248.1597 -617473.3856 424 -67953.7316 -439248.1597 425 -280137.2562 -67953.7316 426 -315144.1159 -280137.2562 427 -1188447.8307 -315144.1159 428 -186944.8306 -1188447.8307 429 -260412.3078 -186944.8306 430 -375551.4045 -260412.3078 431 NA -375551.4045 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 2617776.3344 131691.2321 [2,] 1343261.9381 2617776.3344 [3,] -744280.2782 1343261.9381 [4,] -384882.4697 -744280.2782 [5,] 231829.6052 -384882.4697 [6,] -24799.1285 231829.6052 [7,] -140602.1496 -24799.1285 [8,] 171701.7832 -140602.1496 [9,] 687631.7910 171701.7832 [10,] 24217.3337 687631.7910 [11,] 344779.5263 24217.3337 [12,] 240841.5619 344779.5263 [13,] -741217.9665 240841.5619 [14,] -52513.3362 -741217.9665 [15,] 708793.0810 -52513.3362 [16,] -188311.0806 708793.0810 [17,] -167054.1605 -188311.0806 [18,] -86916.7805 -167054.1605 [19,] -703145.6635 -86916.7805 [20,] 285123.3410 -703145.6635 [21,] -15121.8885 285123.3410 [22,] -47851.8834 -15121.8885 [23,] -818053.2764 -47851.8834 [24,] -264218.2790 -818053.2764 [25,] 8784.4069 -264218.2790 [26,] -362208.6164 8784.4069 [27,] -60845.5727 -362208.6164 [28,] 338692.8381 -60845.5727 [29,] 182219.7224 338692.8381 [30,] 281788.1807 182219.7224 [31,] -36995.6012 281788.1807 [32,] 263445.7065 -36995.6012 [33,] 145992.6598 263445.7065 [34,] 373117.5424 145992.6598 [35,] -651360.7695 373117.5424 [36,] 102982.2074 -651360.7695 [37,] 137960.6316 102982.2074 [38,] -42632.1754 137960.6316 [39,] -336872.4361 -42632.1754 [40,] -219519.1761 -336872.4361 [41,] -222531.9421 -219519.1761 [42,] 137566.0210 -222531.9421 [43,] -172492.4321 137566.0210 [44,] 14357.7386 -172492.4321 [45,] -30779.2236 14357.7386 [46,] 405337.5305 -30779.2236 [47,] -321471.9194 405337.5305 [48,] 18934.8649 -321471.9194 [49,] -195140.8815 18934.8649 [50,] -3391.6077 -195140.8815 [51,] -192350.0646 -3391.6077 [52,] -234497.9226 -192350.0646 [53,] 60187.4370 -234497.9226 [54,] -462598.5773 60187.4370 [55,] 9233.5964 -462598.5773 [56,] 306052.0936 9233.5964 [57,] 33188.9153 306052.0936 [58,] 5900.7235 33188.9153 [59,] -198482.0380 5900.7235 [60,] 23938.1649 -198482.0380 [61,] -196629.2342 23938.1649 [62,] 20928.2372 -196629.2342 [63,] -179696.4179 20928.2372 [64,] -2681.0677 -179696.4179 [65,] 185234.2138 -2681.0677 [66,] 23891.0866 185234.2138 [67,] -21817.5531 23891.0866 [68,] -76648.7399 -21817.5531 [69,] -138163.7547 -76648.7399 [70,] -208938.8859 -138163.7547 [71,] -9289.8180 -208938.8859 [72,] -47672.1817 -9289.8180 [73,] 29587.7757 -47672.1817 [74,] 31861.7253 29587.7757 [75,] -126419.0457 31861.7253 [76,] 42686.2541 -126419.0457 [77,] -26071.1896 42686.2541 [78,] 28375.2746 -26071.1896 [79,] 103544.1770 28375.2746 [80,] -4214.0980 103544.1770 [81,] -65473.3759 -4214.0980 [82,] 14124.1325 -65473.3759 [83,] -17642.2383 14124.1325 [84,] 76278.7647 -17642.2383 [85,] -143459.0879 76278.7647 [86,] -71856.1787 -143459.0879 [87,] 58271.2104 -71856.1787 [88,] -64948.5081 58271.2104 [89,] 83938.2662 -64948.5081 [90,] 42507.8170 83938.2662 [91,] -49580.7887 42507.8170 [92,] 85793.1811 -49580.7887 [93,] -70344.4744 85793.1811 [94,] -6814.0289 -70344.4744 [95,] -15737.3701 -6814.0289 [96,] 147930.4943 -15737.3701 [97,] 213063.0816 147930.4943 [98,] 72826.2490 213063.0816 [99,] -34459.7025 72826.2490 [100,] -79216.1575 -34459.7025 [101,] -79815.2978 -79216.1575 [102,] -111110.0897 -79815.2978 [103,] 39046.6756 -111110.0897 [104,] 52007.6985 39046.6756 [105,] -76065.7715 52007.6985 [106,] -295762.5098 -76065.7715 [107,] 146195.4603 -295762.5098 [108,] -104872.2276 146195.4603 [109,] -108290.3137 -104872.2276 [110,] 91917.1948 -108290.3137 [111,] -238787.2513 91917.1948 [112,] -101375.4457 -238787.2513 [113,] -84835.3560 -101375.4457 [114,] 25282.2672 -84835.3560 [115,] -6420.5109 25282.2672 [116,] 62241.0108 -6420.5109 [117,] 113622.0525 62241.0108 [118,] -147931.6276 113622.0525 [119,] 66650.5828 -147931.6276 [120,] -242477.7777 66650.5828 [121,] 24226.5214 -242477.7777 [122,] -155730.5335 24226.5214 [123,] 103139.2098 -155730.5335 [124,] -3850.1514 103139.2098 [125,] -187428.7167 -3850.1514 [126,] 29268.1262 -187428.7167 [127,] 8003.2527 29268.1262 [128,] 1350.2192 8003.2527 [129,] 29268.1262 1350.2192 [130,] 19435.3611 29268.1262 [131,] 91186.7223 19435.3611 [132,] 29268.1262 91186.7223 [133,] -30144.0315 29268.1262 [134,] 18248.3820 -30144.0315 [135,] 26910.1165 18248.3820 [136,] 21999.4275 26910.1165 [137,] 31090.4241 21999.4275 [138,] 17538.2510 31090.4241 [139,] 29268.1262 17538.2510 [140,] 31858.1740 29268.1262 [141,] 29268.1262 31858.1740 [142,] -11943.4238 29268.1262 [143,] 21985.8429 -11943.4238 [144,] 16747.6812 21985.8429 [145,] 11303.9014 16747.6812 [146,] 29268.1262 11303.9014 [147,] 17224.4191 29268.1262 [148,] 34642.2003 17224.4191 [149,] 31382.0373 34642.2003 [150,] -7622.8240 31382.0373 [151,] -25650.9295 -7622.8240 [152,] 13943.2545 -25650.9295 [153,] 1038.1341 13943.2545 [154,] 27775.4888 1038.1341 [155,] 29268.1262 27775.4888 [156,] 29268.1262 29268.1262 [157,] 25329.7892 29268.1262 [158,] 40865.2598 25329.7892 [159,] 29268.1262 40865.2598 [160,] 29268.1262 29268.1262 [161,] 29268.1262 29268.1262 [162,] 29268.1262 29268.1262 [163,] 36353.2418 29268.1262 [164,] 19251.7140 36353.2418 [165,] 33219.4047 19251.7140 [166,] 29268.1262 33219.4047 [167,] 58320.1911 29268.1262 [168,] 892399.6433 58320.1911 [169,] 29268.1262 892399.6433 [170,] 29268.1262 29268.1262 [171,] 23616.3821 29268.1262 [172,] 29268.1262 23616.3821 [173,] 8373.1483 29268.1262 [174,] 29268.1262 8373.1483 [175,] 21862.0309 29268.1262 [176,] 89432.1141 21862.0309 [177,] -11439.4325 89432.1141 [178,] 29268.1262 -11439.4325 [179,] 29268.1262 29268.1262 [180,] 57060.1061 29268.1262 [181,] -18485.6793 57060.1061 [182,] 31673.1233 -18485.6793 [183,] 50247.3889 31673.1233 [184,] -19977.3842 50247.3889 [185,] -14838.0182 -19977.3842 [186,] -16030.6275 -14838.0182 [187,] 15618.3139 -16030.6275 [188,] -2106.0886 15618.3139 [189,] 29268.1262 -2106.0886 [190,] 13262.5146 29268.1262 [191,] 29268.1262 13262.5146 [192,] 16125.6181 29268.1262 [193,] 29268.1262 16125.6181 [194,] 23372.3899 29268.1262 [195,] -12895.9151 23372.3899 [196,] -144048.2459 -12895.9151 [197,] -12895.9151 -144048.2459 [198,] -12895.9151 -12895.9151 [199,] -86204.1344 -12895.9151 [200,] -12895.9151 -86204.1344 [201,] -12895.9151 -12895.9151 [202,] -12895.9151 -12895.9151 [203,] -23558.2587 -12895.9151 [204,] -11021.0744 -23558.2587 [205,] -23779.8406 -11021.0744 [206,] -7579.1442 -23779.8406 [207,] -12895.9151 -7579.1442 [208,] -27880.4863 -12895.9151 [209,] 50231.4022 -27880.4863 [210,] 29268.1262 50231.4022 [211,] 29268.1262 29268.1262 [212,] 59092.2205 29268.1262 [213,] -1242.1168 59092.2205 [214,] 7927.4901 -1242.1168 [215,] 23138.1775 7927.4901 [216,] 24191.9708 23138.1775 [217,] 50659.9525 24191.9708 [218,] 23138.1775 50659.9525 [219,] -19460.1087 23138.1775 [220,] 3769.9938 -19460.1087 [221,] 29268.1262 3769.9938 [222,] 29414.7303 29268.1262 [223,] 15392.3051 29414.7303 [224,] 29456.2253 15392.3051 [225,] 345.8293 29456.2253 [226,] 38150.1854 345.8293 [227,] 29268.1262 38150.1854 [228,] 26321.6280 29268.1262 [229,] 29268.1262 26321.6280 [230,] 2843.3897 29268.1262 [231,] 24541.7601 2843.3897 [232,] 29268.1262 24541.7601 [233,] 347.1197 29268.1262 [234,] 29268.1262 347.1197 [235,] 31005.9853 29268.1262 [236,] 29268.1262 31005.9853 [237,] 20243.1219 29268.1262 [238,] 29268.1262 20243.1219 [239,] 29268.1262 29268.1262 [240,] 29268.1262 29268.1262 [241,] 21535.2133 29268.1262 [242,] 29268.1262 21535.2133 [243,] 20500.1432 29268.1262 [244,] 12045.8516 20500.1432 [245,] -12895.9151 12045.8516 [246,] -12895.9151 -12895.9151 [247,] -6412.0358 -12895.9151 [248,] -12895.9151 -6412.0358 [249,] -12895.9151 -12895.9151 [250,] -8320.2368 -12895.9151 [251,] -12895.9151 -8320.2368 [252,] 21024.5792 -12895.9151 [253,] 4443.0129 21024.5792 [254,] -95940.8109 4443.0129 [255,] 29268.1262 -95940.8109 [256,] 7405.9368 29268.1262 [257,] -14375.4965 7405.9368 [258,] -12895.9151 -14375.4965 [259,] 29268.1262 -12895.9151 [260,] -12895.9151 29268.1262 [261,] -12895.9151 -12895.9151 [262,] 39051.9120 -12895.9151 [263,] 34965.7507 39051.9120 [264,] -12895.9151 34965.7507 [265,] 29268.1262 -12895.9151 [266,] -12895.9151 29268.1262 [267,] 89547.2459 -12895.9151 [268,] -12895.9151 89547.2459 [269,] -22091.7116 -12895.9151 [270,] -12895.9151 -22091.7116 [271,] -12895.9151 -12895.9151 [272,] -12895.9151 -12895.9151 [273,] -16124.8231 -12895.9151 [274,] 29268.1262 -16124.8231 [275,] 42623.2566 29268.1262 [276,] 13388.0309 42623.2566 [277,] 29268.1262 13388.0309 [278,] 80534.1352 29268.1262 [279,] 44298.5890 80534.1352 [280,] -12895.9151 44298.5890 [281,] 29268.1262 -12895.9151 [282,] 65466.0576 29268.1262 [283,] 29268.1262 65466.0576 [284,] 37912.3754 29268.1262 [285,] 29268.1262 37912.3754 [286,] -16500.7829 29268.1262 [287,] 29634.3242 -16500.7829 [288,] -12895.9151 29634.3242 [289,] 68645.4134 -12895.9151 [290,] -12895.9151 68645.4134 [291,] 39801.4098 -12895.9151 [292,] 54131.1182 39801.4098 [293,] 29268.1262 54131.1182 [294,] 7535.8105 29268.1262 [295,] 14006.3083 7535.8105 [296,] 29268.1262 14006.3083 [297,] -12895.9151 29268.1262 [298,] 29268.1262 -12895.9151 [299,] 1683.3572 29268.1262 [300,] -21170.2426 1683.3572 [301,] 10134.3494 -21170.2426 [302,] 26203.1519 10134.3494 [303,] -12895.9151 26203.1519 [304,] 97098.1523 -12895.9151 [305,] 7294.2924 97098.1523 [306,] -12895.9151 7294.2924 [307,] -50240.2563 -12895.9151 [308,] 1249.7813 -50240.2563 [309,] -12895.9151 1249.7813 [310,] -12895.9151 -12895.9151 [311,] 70339.2670 -12895.9151 [312,] -12895.9151 70339.2670 [313,] -35175.0997 -12895.9151 [314,] 11587.0227 -35175.0997 [315,] 64018.5109 11587.0227 [316,] 29631.8963 64018.5109 [317,] -27343.7769 29631.8963 [318,] 23796.4518 -27343.7769 [319,] 29268.1262 23796.4518 [320,] 29268.1262 29268.1262 [321,] 29268.1262 29268.1262 [322,] 29268.1262 29268.1262 [323,] -12895.9151 29268.1262 [324,] 29268.1262 -12895.9151 [325,] 26308.4131 29268.1262 [326,] -12895.9151 26308.4131 [327,] -12895.9151 -12895.9151 [328,] 29268.1262 -12895.9151 [329,] 29268.1262 29268.1262 [330,] 29268.1262 29268.1262 [331,] 29268.1262 29268.1262 [332,] 29268.1262 29268.1262 [333,] 22065.9821 29268.1262 [334,] 29492.8409 22065.9821 [335,] 26204.8240 29492.8409 [336,] 141959.3090 26204.8240 [337,] 29268.1262 141959.3090 [338,] 26082.9678 29268.1262 [339,] 20224.3425 26082.9678 [340,] -6185.1728 20224.3425 [341,] 29268.1262 -6185.1728 [342,] 29268.1262 29268.1262 [343,] 29268.1262 29268.1262 [344,] 104371.9698 29268.1262 [345,] 29268.1262 104371.9698 [346,] 15803.2182 29268.1262 [347,] 29268.1262 15803.2182 [348,] 29268.1262 29268.1262 [349,] 26480.6749 29268.1262 [350,] 41319.2054 26480.6749 [351,] 22261.8646 41319.2054 [352,] 20005.2039 22261.8646 [353,] 27164.0740 20005.2039 [354,] 29268.1262 27164.0740 [355,] 29268.1262 29268.1262 [356,] 62525.8537 29268.1262 [357,] 29268.1262 62525.8537 [358,] 29268.1262 29268.1262 [359,] 29268.1262 29268.1262 [360,] 23597.0153 29268.1262 [361,] 29268.1262 23597.0153 [362,] 29268.1262 29268.1262 [363,] -67943.7223 29268.1262 [364,] 17837.3355 -67943.7223 [365,] 29268.1262 17837.3355 [366,] 29268.1262 29268.1262 [367,] 29268.1262 29268.1262 [368,] 29268.1262 29268.1262 [369,] 62148.9428 29268.1262 [370,] 29268.1262 62148.9428 [371,] 75442.7081 29268.1262 [372,] 47466.9676 75442.7081 [373,] 13641.2611 47466.9676 [374,] 29268.1262 13641.2611 [375,] 29268.1262 29268.1262 [376,] 12620.8166 29268.1262 [377,] 24221.5966 12620.8166 [378,] 10112.3784 24221.5966 [379,] 29268.1262 10112.3784 [380,] 29268.1262 29268.1262 [381,] 61189.8980 29268.1262 [382,] 91110.8221 61189.8980 [383,] 29268.1262 91110.8221 [384,] 20271.8496 29268.1262 [385,] -106966.4255 20271.8496 [386,] 39288.3148 -106966.4255 [387,] 3777.1431 39288.3148 [388,] 22655.9444 3777.1431 [389,] 14935.5040 22655.9444 [390,] 17246.8195 14935.5040 [391,] -6152.4952 17246.8195 [392,] -44485.5778 -6152.4952 [393,] 69803.8539 -44485.5778 [394,] -91897.9029 69803.8539 [395,] -14362.6556 -91897.9029 [396,] -120978.1743 -14362.6556 [397,] -539798.0202 -120978.1743 [398,] -2599.6976 -539798.0202 [399,] 48066.8103 -2599.6976 [400,] 7051.7420 48066.8103 [401,] 28254.8104 7051.7420 [402,] -3696.2117 28254.8104 [403,] -92768.8288 -3696.2117 [404,] -25845.1997 -92768.8288 [405,] 15471.1292 -25845.1997 [406,] 17005.1009 15471.1292 [407,] -14854.1938 17005.1009 [408,] -4598.1795 -14854.1938 [409,] -63801.2624 -4598.1795 [410,] -73896.2632 -63801.2624 [411,] -156603.4022 -73896.2632 [412,] -93868.7126 -156603.4022 [413,] 29585.1004 -93868.7126 [414,] -139508.1245 29585.1004 [415,] -89769.1764 -139508.1245 [416,] 128309.4869 -89769.1764 [417,] -86655.2669 128309.4869 [418,] -58399.6203 -86655.2669 [419,] -44237.0641 -58399.6203 [420,] -41702.0976 -44237.0641 [421,] -209029.6603 -41702.0976 [422,] -617473.3856 -209029.6603 [423,] -439248.1597 -617473.3856 [424,] -67953.7316 -439248.1597 [425,] -280137.2562 -67953.7316 [426,] -315144.1159 -280137.2562 [427,] -1188447.8307 -315144.1159 [428,] -186944.8306 -1188447.8307 [429,] -260412.3078 -186944.8306 [430,] -375551.4045 -260412.3078 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 2617776.3344 131691.2321 2 1343261.9381 2617776.3344 3 -744280.2782 1343261.9381 4 -384882.4697 -744280.2782 5 231829.6052 -384882.4697 6 -24799.1285 231829.6052 7 -140602.1496 -24799.1285 8 171701.7832 -140602.1496 9 687631.7910 171701.7832 10 24217.3337 687631.7910 11 344779.5263 24217.3337 12 240841.5619 344779.5263 13 -741217.9665 240841.5619 14 -52513.3362 -741217.9665 15 708793.0810 -52513.3362 16 -188311.0806 708793.0810 17 -167054.1605 -188311.0806 18 -86916.7805 -167054.1605 19 -703145.6635 -86916.7805 20 285123.3410 -703145.6635 21 -15121.8885 285123.3410 22 -47851.8834 -15121.8885 23 -818053.2764 -47851.8834 24 -264218.2790 -818053.2764 25 8784.4069 -264218.2790 26 -362208.6164 8784.4069 27 -60845.5727 -362208.6164 28 338692.8381 -60845.5727 29 182219.7224 338692.8381 30 281788.1807 182219.7224 31 -36995.6012 281788.1807 32 263445.7065 -36995.6012 33 145992.6598 263445.7065 34 373117.5424 145992.6598 35 -651360.7695 373117.5424 36 102982.2074 -651360.7695 37 137960.6316 102982.2074 38 -42632.1754 137960.6316 39 -336872.4361 -42632.1754 40 -219519.1761 -336872.4361 41 -222531.9421 -219519.1761 42 137566.0210 -222531.9421 43 -172492.4321 137566.0210 44 14357.7386 -172492.4321 45 -30779.2236 14357.7386 46 405337.5305 -30779.2236 47 -321471.9194 405337.5305 48 18934.8649 -321471.9194 49 -195140.8815 18934.8649 50 -3391.6077 -195140.8815 51 -192350.0646 -3391.6077 52 -234497.9226 -192350.0646 53 60187.4370 -234497.9226 54 -462598.5773 60187.4370 55 9233.5964 -462598.5773 56 306052.0936 9233.5964 57 33188.9153 306052.0936 58 5900.7235 33188.9153 59 -198482.0380 5900.7235 60 23938.1649 -198482.0380 61 -196629.2342 23938.1649 62 20928.2372 -196629.2342 63 -179696.4179 20928.2372 64 -2681.0677 -179696.4179 65 185234.2138 -2681.0677 66 23891.0866 185234.2138 67 -21817.5531 23891.0866 68 -76648.7399 -21817.5531 69 -138163.7547 -76648.7399 70 -208938.8859 -138163.7547 71 -9289.8180 -208938.8859 72 -47672.1817 -9289.8180 73 29587.7757 -47672.1817 74 31861.7253 29587.7757 75 -126419.0457 31861.7253 76 42686.2541 -126419.0457 77 -26071.1896 42686.2541 78 28375.2746 -26071.1896 79 103544.1770 28375.2746 80 -4214.0980 103544.1770 81 -65473.3759 -4214.0980 82 14124.1325 -65473.3759 83 -17642.2383 14124.1325 84 76278.7647 -17642.2383 85 -143459.0879 76278.7647 86 -71856.1787 -143459.0879 87 58271.2104 -71856.1787 88 -64948.5081 58271.2104 89 83938.2662 -64948.5081 90 42507.8170 83938.2662 91 -49580.7887 42507.8170 92 85793.1811 -49580.7887 93 -70344.4744 85793.1811 94 -6814.0289 -70344.4744 95 -15737.3701 -6814.0289 96 147930.4943 -15737.3701 97 213063.0816 147930.4943 98 72826.2490 213063.0816 99 -34459.7025 72826.2490 100 -79216.1575 -34459.7025 101 -79815.2978 -79216.1575 102 -111110.0897 -79815.2978 103 39046.6756 -111110.0897 104 52007.6985 39046.6756 105 -76065.7715 52007.6985 106 -295762.5098 -76065.7715 107 146195.4603 -295762.5098 108 -104872.2276 146195.4603 109 -108290.3137 -104872.2276 110 91917.1948 -108290.3137 111 -238787.2513 91917.1948 112 -101375.4457 -238787.2513 113 -84835.3560 -101375.4457 114 25282.2672 -84835.3560 115 -6420.5109 25282.2672 116 62241.0108 -6420.5109 117 113622.0525 62241.0108 118 -147931.6276 113622.0525 119 66650.5828 -147931.6276 120 -242477.7777 66650.5828 121 24226.5214 -242477.7777 122 -155730.5335 24226.5214 123 103139.2098 -155730.5335 124 -3850.1514 103139.2098 125 -187428.7167 -3850.1514 126 29268.1262 -187428.7167 127 8003.2527 29268.1262 128 1350.2192 8003.2527 129 29268.1262 1350.2192 130 19435.3611 29268.1262 131 91186.7223 19435.3611 132 29268.1262 91186.7223 133 -30144.0315 29268.1262 134 18248.3820 -30144.0315 135 26910.1165 18248.3820 136 21999.4275 26910.1165 137 31090.4241 21999.4275 138 17538.2510 31090.4241 139 29268.1262 17538.2510 140 31858.1740 29268.1262 141 29268.1262 31858.1740 142 -11943.4238 29268.1262 143 21985.8429 -11943.4238 144 16747.6812 21985.8429 145 11303.9014 16747.6812 146 29268.1262 11303.9014 147 17224.4191 29268.1262 148 34642.2003 17224.4191 149 31382.0373 34642.2003 150 -7622.8240 31382.0373 151 -25650.9295 -7622.8240 152 13943.2545 -25650.9295 153 1038.1341 13943.2545 154 27775.4888 1038.1341 155 29268.1262 27775.4888 156 29268.1262 29268.1262 157 25329.7892 29268.1262 158 40865.2598 25329.7892 159 29268.1262 40865.2598 160 29268.1262 29268.1262 161 29268.1262 29268.1262 162 29268.1262 29268.1262 163 36353.2418 29268.1262 164 19251.7140 36353.2418 165 33219.4047 19251.7140 166 29268.1262 33219.4047 167 58320.1911 29268.1262 168 892399.6433 58320.1911 169 29268.1262 892399.6433 170 29268.1262 29268.1262 171 23616.3821 29268.1262 172 29268.1262 23616.3821 173 8373.1483 29268.1262 174 29268.1262 8373.1483 175 21862.0309 29268.1262 176 89432.1141 21862.0309 177 -11439.4325 89432.1141 178 29268.1262 -11439.4325 179 29268.1262 29268.1262 180 57060.1061 29268.1262 181 -18485.6793 57060.1061 182 31673.1233 -18485.6793 183 50247.3889 31673.1233 184 -19977.3842 50247.3889 185 -14838.0182 -19977.3842 186 -16030.6275 -14838.0182 187 15618.3139 -16030.6275 188 -2106.0886 15618.3139 189 29268.1262 -2106.0886 190 13262.5146 29268.1262 191 29268.1262 13262.5146 192 16125.6181 29268.1262 193 29268.1262 16125.6181 194 23372.3899 29268.1262 195 -12895.9151 23372.3899 196 -144048.2459 -12895.9151 197 -12895.9151 -144048.2459 198 -12895.9151 -12895.9151 199 -86204.1344 -12895.9151 200 -12895.9151 -86204.1344 201 -12895.9151 -12895.9151 202 -12895.9151 -12895.9151 203 -23558.2587 -12895.9151 204 -11021.0744 -23558.2587 205 -23779.8406 -11021.0744 206 -7579.1442 -23779.8406 207 -12895.9151 -7579.1442 208 -27880.4863 -12895.9151 209 50231.4022 -27880.4863 210 29268.1262 50231.4022 211 29268.1262 29268.1262 212 59092.2205 29268.1262 213 -1242.1168 59092.2205 214 7927.4901 -1242.1168 215 23138.1775 7927.4901 216 24191.9708 23138.1775 217 50659.9525 24191.9708 218 23138.1775 50659.9525 219 -19460.1087 23138.1775 220 3769.9938 -19460.1087 221 29268.1262 3769.9938 222 29414.7303 29268.1262 223 15392.3051 29414.7303 224 29456.2253 15392.3051 225 345.8293 29456.2253 226 38150.1854 345.8293 227 29268.1262 38150.1854 228 26321.6280 29268.1262 229 29268.1262 26321.6280 230 2843.3897 29268.1262 231 24541.7601 2843.3897 232 29268.1262 24541.7601 233 347.1197 29268.1262 234 29268.1262 347.1197 235 31005.9853 29268.1262 236 29268.1262 31005.9853 237 20243.1219 29268.1262 238 29268.1262 20243.1219 239 29268.1262 29268.1262 240 29268.1262 29268.1262 241 21535.2133 29268.1262 242 29268.1262 21535.2133 243 20500.1432 29268.1262 244 12045.8516 20500.1432 245 -12895.9151 12045.8516 246 -12895.9151 -12895.9151 247 -6412.0358 -12895.9151 248 -12895.9151 -6412.0358 249 -12895.9151 -12895.9151 250 -8320.2368 -12895.9151 251 -12895.9151 -8320.2368 252 21024.5792 -12895.9151 253 4443.0129 21024.5792 254 -95940.8109 4443.0129 255 29268.1262 -95940.8109 256 7405.9368 29268.1262 257 -14375.4965 7405.9368 258 -12895.9151 -14375.4965 259 29268.1262 -12895.9151 260 -12895.9151 29268.1262 261 -12895.9151 -12895.9151 262 39051.9120 -12895.9151 263 34965.7507 39051.9120 264 -12895.9151 34965.7507 265 29268.1262 -12895.9151 266 -12895.9151 29268.1262 267 89547.2459 -12895.9151 268 -12895.9151 89547.2459 269 -22091.7116 -12895.9151 270 -12895.9151 -22091.7116 271 -12895.9151 -12895.9151 272 -12895.9151 -12895.9151 273 -16124.8231 -12895.9151 274 29268.1262 -16124.8231 275 42623.2566 29268.1262 276 13388.0309 42623.2566 277 29268.1262 13388.0309 278 80534.1352 29268.1262 279 44298.5890 80534.1352 280 -12895.9151 44298.5890 281 29268.1262 -12895.9151 282 65466.0576 29268.1262 283 29268.1262 65466.0576 284 37912.3754 29268.1262 285 29268.1262 37912.3754 286 -16500.7829 29268.1262 287 29634.3242 -16500.7829 288 -12895.9151 29634.3242 289 68645.4134 -12895.9151 290 -12895.9151 68645.4134 291 39801.4098 -12895.9151 292 54131.1182 39801.4098 293 29268.1262 54131.1182 294 7535.8105 29268.1262 295 14006.3083 7535.8105 296 29268.1262 14006.3083 297 -12895.9151 29268.1262 298 29268.1262 -12895.9151 299 1683.3572 29268.1262 300 -21170.2426 1683.3572 301 10134.3494 -21170.2426 302 26203.1519 10134.3494 303 -12895.9151 26203.1519 304 97098.1523 -12895.9151 305 7294.2924 97098.1523 306 -12895.9151 7294.2924 307 -50240.2563 -12895.9151 308 1249.7813 -50240.2563 309 -12895.9151 1249.7813 310 -12895.9151 -12895.9151 311 70339.2670 -12895.9151 312 -12895.9151 70339.2670 313 -35175.0997 -12895.9151 314 11587.0227 -35175.0997 315 64018.5109 11587.0227 316 29631.8963 64018.5109 317 -27343.7769 29631.8963 318 23796.4518 -27343.7769 319 29268.1262 23796.4518 320 29268.1262 29268.1262 321 29268.1262 29268.1262 322 29268.1262 29268.1262 323 -12895.9151 29268.1262 324 29268.1262 -12895.9151 325 26308.4131 29268.1262 326 -12895.9151 26308.4131 327 -12895.9151 -12895.9151 328 29268.1262 -12895.9151 329 29268.1262 29268.1262 330 29268.1262 29268.1262 331 29268.1262 29268.1262 332 29268.1262 29268.1262 333 22065.9821 29268.1262 334 29492.8409 22065.9821 335 26204.8240 29492.8409 336 141959.3090 26204.8240 337 29268.1262 141959.3090 338 26082.9678 29268.1262 339 20224.3425 26082.9678 340 -6185.1728 20224.3425 341 29268.1262 -6185.1728 342 29268.1262 29268.1262 343 29268.1262 29268.1262 344 104371.9698 29268.1262 345 29268.1262 104371.9698 346 15803.2182 29268.1262 347 29268.1262 15803.2182 348 29268.1262 29268.1262 349 26480.6749 29268.1262 350 41319.2054 26480.6749 351 22261.8646 41319.2054 352 20005.2039 22261.8646 353 27164.0740 20005.2039 354 29268.1262 27164.0740 355 29268.1262 29268.1262 356 62525.8537 29268.1262 357 29268.1262 62525.8537 358 29268.1262 29268.1262 359 29268.1262 29268.1262 360 23597.0153 29268.1262 361 29268.1262 23597.0153 362 29268.1262 29268.1262 363 -67943.7223 29268.1262 364 17837.3355 -67943.7223 365 29268.1262 17837.3355 366 29268.1262 29268.1262 367 29268.1262 29268.1262 368 29268.1262 29268.1262 369 62148.9428 29268.1262 370 29268.1262 62148.9428 371 75442.7081 29268.1262 372 47466.9676 75442.7081 373 13641.2611 47466.9676 374 29268.1262 13641.2611 375 29268.1262 29268.1262 376 12620.8166 29268.1262 377 24221.5966 12620.8166 378 10112.3784 24221.5966 379 29268.1262 10112.3784 380 29268.1262 29268.1262 381 61189.8980 29268.1262 382 91110.8221 61189.8980 383 29268.1262 91110.8221 384 20271.8496 29268.1262 385 -106966.4255 20271.8496 386 39288.3148 -106966.4255 387 3777.1431 39288.3148 388 22655.9444 3777.1431 389 14935.5040 22655.9444 390 17246.8195 14935.5040 391 -6152.4952 17246.8195 392 -44485.5778 -6152.4952 393 69803.8539 -44485.5778 394 -91897.9029 69803.8539 395 -14362.6556 -91897.9029 396 -120978.1743 -14362.6556 397 -539798.0202 -120978.1743 398 -2599.6976 -539798.0202 399 48066.8103 -2599.6976 400 7051.7420 48066.8103 401 28254.8104 7051.7420 402 -3696.2117 28254.8104 403 -92768.8288 -3696.2117 404 -25845.1997 -92768.8288 405 15471.1292 -25845.1997 406 17005.1009 15471.1292 407 -14854.1938 17005.1009 408 -4598.1795 -14854.1938 409 -63801.2624 -4598.1795 410 -73896.2632 -63801.2624 411 -156603.4022 -73896.2632 412 -93868.7126 -156603.4022 413 29585.1004 -93868.7126 414 -139508.1245 29585.1004 415 -89769.1764 -139508.1245 416 128309.4869 -89769.1764 417 -86655.2669 128309.4869 418 -58399.6203 -86655.2669 419 -44237.0641 -58399.6203 420 -41702.0976 -44237.0641 421 -209029.6603 -41702.0976 422 -617473.3856 -209029.6603 423 -439248.1597 -617473.3856 424 -67953.7316 -439248.1597 425 -280137.2562 -67953.7316 426 -315144.1159 -280137.2562 427 -1188447.8307 -315144.1159 428 -186944.8306 -1188447.8307 429 -260412.3078 -186944.8306 430 -375551.4045 -260412.3078 > plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals') > lines(lowess(z)) > abline(lm(z)) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7h2ae1291394141.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/8h2ae1291394141.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/9h2ae1291394141.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0)) > plot(mylm, las = 1, sub='Residual Diagnostics') > par(opar) > dev.off() null device 1 > if (n > n25) { + postscript(file="/var/www/html/rcomp/tmp/10absh1291394141.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) + plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint') + grid() + dev.off() + } null device 1 > > #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE) > a<-table.row.end(a) > myeq <- colnames(x)[1] > myeq <- paste(myeq, '[t] = ', sep='') > for (i in 1:k){ + if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '') + myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ') + if (rownames(mysum$coefficients)[i] != '(Intercept)') { + myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='') + if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='') + } + } > myeq <- paste(myeq, ' + e[t]') > a<-table.row.start(a) > a<-table.element(a, myeq) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/11duqn1291394141.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Variable',header=TRUE) > a<-table.element(a,'Parameter',header=TRUE) > a<-table.element(a,'S.D.',header=TRUE) > a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE) > a<-table.element(a,'2-tail p-value',header=TRUE) > a<-table.element(a,'1-tail p-value',header=TRUE) > a<-table.row.end(a) > for (i in 1:k){ + a<-table.row.start(a) + a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE) + a<-table.element(a,mysum$coefficients[i,1]) + a<-table.element(a, round(mysum$coefficients[i,2],6)) + a<-table.element(a, round(mysum$coefficients[i,3],4)) + a<-table.element(a, round(mysum$coefficients[i,4],6)) + a<-table.element(a, round(mysum$coefficients[i,4]/2,6)) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/12zu7t1291394141.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple R',1,TRUE) > a<-table.element(a, sqrt(mysum$r.squared)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'R-squared',1,TRUE) > a<-table.element(a, mysum$r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Adjusted R-squared',1,TRUE) > a<-table.element(a, mysum$adj.r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (value)',1,TRUE) > a<-table.element(a, mysum$fstatistic[1]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[2]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[3]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'p-value',1,TRUE) > a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3])) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Residual Standard Deviation',1,TRUE) > a<-table.element(a, mysum$sigma) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Sum Squared Residuals',1,TRUE) > a<-table.element(a, sum(myerror*myerror)) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/13vm521291394141.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Time or Index', 1, TRUE) > a<-table.element(a, 'Actuals', 1, TRUE) > a<-table.element(a, 'Interpolation
Forecast', 1, TRUE) > a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,i, 1, TRUE) + a<-table.element(a,x[i]) + a<-table.element(a,x[i]-mysum$resid[i]) + a<-table.element(a,mysum$resid[i]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/14y4381291394141.tab") > if (n > n25) { + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'p-values',header=TRUE) + a<-table.element(a,'Alternative Hypothesis',3,header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'breakpoint index',header=TRUE) + a<-table.element(a,'greater',header=TRUE) + a<-table.element(a,'2-sided',header=TRUE) + a<-table.element(a,'less',header=TRUE) + a<-table.row.end(a) + for (mypoint in kp3:nmkm3) { + a<-table.row.start(a) + a<-table.element(a,mypoint,header=TRUE) + a<-table.element(a,gqarr[mypoint-kp3+1,1]) + a<-table.element(a,gqarr[mypoint-kp3+1,2]) + a<-table.element(a,gqarr[mypoint-kp3+1,3]) + a<-table.row.end(a) + } + a<-table.end(a) + table.save(a,file="/var/www/html/rcomp/tmp/15252d1291394141.tab") + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'Description',header=TRUE) + a<-table.element(a,'# significant tests',header=TRUE) + a<-table.element(a,'% significant tests',header=TRUE) + a<-table.element(a,'OK/NOK',header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'1% type I error level',header=TRUE) + a<-table.element(a,numsignificant1) + a<-table.element(a,numsignificant1/numgqtests) + if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'5% type I error level',header=TRUE) + a<-table.element(a,numsignificant5) + a<-table.element(a,numsignificant5/numgqtests) + if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'10% type I error level',header=TRUE) + a<-table.element(a,numsignificant10) + a<-table.element(a,numsignificant10/numgqtests) + if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.end(a) + table.save(a,file="/var/www/html/rcomp/tmp/16ffz41291394141.tab") + } > > try(system("convert tmp/1lsvn1291394141.ps tmp/1lsvn1291394141.png",intern=TRUE)) character(0) > try(system("convert tmp/2lsvn1291394141.ps tmp/2lsvn1291394141.png",intern=TRUE)) character(0) > try(system("convert tmp/3w1u81291394141.ps tmp/3w1u81291394141.png",intern=TRUE)) character(0) > try(system("convert tmp/4w1u81291394141.ps tmp/4w1u81291394141.png",intern=TRUE)) character(0) > try(system("convert tmp/5w1u81291394141.ps tmp/5w1u81291394141.png",intern=TRUE)) character(0) > try(system("convert tmp/66tbb1291394141.ps tmp/66tbb1291394141.png",intern=TRUE)) character(0) > try(system("convert tmp/7h2ae1291394141.ps tmp/7h2ae1291394141.png",intern=TRUE)) character(0) > try(system("convert tmp/8h2ae1291394141.ps tmp/8h2ae1291394141.png",intern=TRUE)) character(0) > try(system("convert tmp/9h2ae1291394141.ps tmp/9h2ae1291394141.png",intern=TRUE)) character(0) > try(system("convert tmp/10absh1291394141.ps tmp/10absh1291394141.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 13.199 2.131 27.349