R version 2.8.0 (2008-10-20) Copyright (C) 2008 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. Natural language support but running in an English locale R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(162556 + ,1081 + ,807 + ,213118 + ,6282154 + ,1 + ,29790 + ,309 + ,444 + ,81767 + ,4321023 + ,1 + ,87550 + ,458 + ,412 + ,153198 + ,4111912 + ,1 + ,84738 + ,588 + ,428 + ,-26007 + ,223193 + ,1 + ,54660 + ,302 + ,315 + ,126942 + ,1491348 + ,1 + ,42634 + ,156 + ,168 + ,157214 + ,1629616 + ,1 + ,40949 + ,481 + ,263 + ,129352 + ,1398893 + ,1 + ,45187 + ,353 + ,267 + ,234817 + ,1926517 + ,1 + ,37704 + ,452 + ,228 + ,60448 + ,983660 + ,1 + ,16275 + ,109 + ,129 + ,47818 + ,1443586 + ,1 + ,25830 + ,115 + ,104 + ,245546 + ,1073089 + ,1 + ,12679 + ,110 + ,122 + ,48020 + ,984885 + ,1 + ,18014 + ,239 + ,393 + ,-1710 + ,1405225 + ,1 + ,43556 + ,247 + ,190 + ,32648 + ,227132 + ,1 + ,24811 + ,505 + ,280 + ,95350 + ,929118 + ,1 + ,6575 + ,159 + ,63 + ,151352 + ,1071292 + ,1 + ,7123 + ,109 + ,102 + ,288170 + ,638830 + ,1 + ,21950 + ,519 + ,265 + ,114337 + ,856956 + ,1 + ,37597 + ,248 + ,234 + ,37884 + ,992426 + ,1 + ,17821 + ,373 + ,277 + ,122844 + ,444477 + ,1 + ,12988 + ,119 + ,73 + ,82340 + ,857217 + ,1 + ,22330 + ,84 + ,67 + ,79801 + ,711969 + ,1 + ,13326 + ,102 + ,103 + ,165548 + ,702380 + ,1 + ,16189 + ,295 + ,290 + ,116384 + ,358589 + ,1 + ,7146 + ,105 + ,83 + ,134028 + ,297978 + ,1 + ,15824 + ,64 + ,56 + ,63838 + ,585715 + ,1 + ,27664 + ,282 + ,236 + ,74996 + ,657954 + ,1 + ,11920 + ,182 + ,73 + ,31080 + ,209458 + ,1 + ,8568 + ,37 + ,34 + ,32168 + ,786690 + ,1 + ,14416 + ,361 + ,139 + ,49857 + ,439798 + ,1 + ,3369 + ,28 + ,26 + ,87161 + ,688779 + ,1 + ,11819 + ,85 + ,70 + ,106113 + ,574339 + ,1 + ,6984 + ,45 + ,40 + ,80570 + ,741409 + ,1 + ,4519 + ,49 + ,42 + ,102129 + ,597793 + ,1 + ,2220 + ,22 + ,12 + ,301670 + ,644190 + ,1 + ,18562 + ,155 + ,211 + ,102313 + ,377934 + ,1 + ,10327 + ,91 + ,74 + ,88577 + ,640273 + ,1 + ,5336 + ,81 + ,80 + ,112477 + ,697458 + ,1 + ,2365 + ,79 + ,83 + ,191778 + ,550608 + ,1 + ,4069 + ,145 + ,131 + ,79804 + ,207393 + ,1 + ,8636 + ,855 + ,203 + ,128294 + ,301607 + ,1 + ,13718 + ,61 + ,56 + ,96448 + ,345783 + ,1 + ,4525 + ,226 + ,89 + ,93811 + ,501749 + ,1 + ,6869 + ,105 + ,88 + ,117520 + ,379983 + ,1 + ,4628 + ,62 + ,39 + ,69159 + ,387475 + ,1 + ,3689 + ,25 + ,25 + ,101792 + ,377305 + ,1 + ,4891 + ,217 + ,49 + ,210568 + ,370837 + ,1 + ,7489 + ,322 + ,149 + ,136996 + ,430866 + ,1 + ,4901 + ,84 + ,58 + ,121920 + ,469107 + ,1 + ,2284 + ,33 + ,41 + ,76403 + ,194493 + ,1 + ,3160 + ,108 + ,90 + ,108094 + ,530670 + ,1 + ,4150 + ,150 + ,136 + ,134759 + ,518365 + ,1 + ,7285 + ,115 + ,97 + ,188873 + ,491303 + ,1 + ,1134 + ,162 + ,63 + ,146216 + ,527021 + ,1 + ,4658 + ,158 + ,114 + ,156608 + ,233773 + ,1 + ,2384 + ,97 + ,77 + ,61348 + ,405972 + ,1 + ,3748 + ,9 + ,6 + ,50350 + ,652925 + ,1 + ,5371 + ,66 + ,47 + ,87720 + ,446211 + ,1 + ,1285 + ,107 + ,51 + ,99489 + ,341340 + ,1 + ,9327 + ,101 + ,85 + ,87419 + ,387699 + ,1 + ,5565 + ,47 + ,43 + ,94355 + ,493408 + ,1 + ,1528 + ,38 + ,32 + ,60326 + ,146494 + ,1 + ,3122 + ,34 + ,25 + ,94670 + ,414462 + ,1 + ,7561 + ,87 + ,77 + ,82425 + ,364304 + ,1 + ,2675 + ,79 + ,54 + ,59017 + ,355178 + ,1 + ,13253 + ,947 + ,251 + ,90829 + ,357760 + ,1 + ,880 + ,74 + ,15 + ,80791 + ,261216 + ,1 + ,2053 + ,53 + ,44 + ,100423 + ,397144 + ,1 + ,1424 + ,94 + ,73 + ,131116 + ,374943 + ,1 + ,4036 + ,63 + ,85 + ,100269 + ,424898 + ,1 + ,3045 + ,58 + ,49 + ,27330 + ,202055 + ,1 + ,5119 + ,49 + ,38 + ,39039 + ,378525 + ,1 + ,1431 + ,34 + ,35 + ,106885 + ,310768 + ,1 + ,554 + ,11 + ,9 + ,79285 + ,325738 + ,1 + ,1975 + ,35 + ,34 + ,118881 + ,394510 + ,1 + ,1765 + ,20 + ,20 + ,77623 + ,247060 + ,1 + ,1012 + ,47 + ,29 + ,114768 + ,368078 + ,1 + ,810 + ,43 + ,11 + ,74015 + ,236761 + ,1 + ,1280 + ,117 + ,52 + ,69465 + ,312378 + ,1 + ,666 + ,171 + ,13 + ,117869 + ,339836 + ,1 + ,1380 + ,26 + ,29 + ,60982 + ,347385 + ,1 + ,4677 + ,75 + ,66 + ,90131 + ,426280 + ,1 + ,876 + ,59 + ,33 + ,138971 + ,352850 + ,1 + ,814 + ,18 + ,15 + ,39625 + ,301881 + ,1 + ,514 + ,15 + ,15 + ,102725 + ,377516 + ,1 + ,5692 + ,72 + ,68 + ,64239 + ,357312 + ,1 + ,3642 + ,86 + ,100 + ,90262 + ,458343 + ,1 + ,540 + ,14 + ,13 + ,103960 + ,354228 + ,1 + ,2099 + ,64 + ,45 + ,106611 + ,308636 + ,1 + ,567 + ,11 + ,14 + ,103345 + ,386212 + ,1 + ,2001 + ,52 + ,36 + ,95551 + ,393343 + ,1 + ,2949 + ,41 + ,40 + ,82903 + ,378509 + ,1 + ,2253 + ,99 + ,68 + ,63593 + ,452469 + ,1 + ,6533 + ,75 + ,29 + ,126910 + ,364839 + ,1 + ,1889 + ,45 + ,43 + ,37527 + ,358649 + ,1 + ,3055 + ,43 + ,30 + ,60247 + ,376641 + ,1 + ,272 + ,8 + ,9 + ,112995 + ,429112 + ,1 + ,1414 + ,198 + ,22 + ,70184 + ,330546 + ,1 + ,2564 + ,22 + ,19 + ,130140 + ,403560 + ,1 + ,1383 + ,11 + ,9 + ,73221 + ,317892 + ,1 + ,1261 + ,33 + ,31 + ,76114 + ,307528 + ,1 + ,975 + ,23 + ,19 + ,90534 + ,235133 + ,1 + ,3366 + ,80 + ,55 + ,108479 + ,299243 + ,1 + ,576 + ,18 + ,8 + ,113761 + ,314073 + ,1 + ,1686 + ,40 + ,28 + ,68696 + ,368186 + ,1 + ,746 + ,23 + ,29 + ,71561 + ,269661 + ,1 + ,3192 + ,60 + ,48 + ,59831 + ,125390 + ,1 + ,2045 + ,20 + ,16 + ,97890 + ,510834 + ,1 + ,5702 + ,61 + ,47 + ,101481 + ,321896 + ,1 + ,1932 + ,36 + ,20 + ,72954 + ,249898 + ,1 + ,936 + ,30 + ,22 + ,67939 + ,408881 + ,1 + ,3437 + ,47 + ,33 + ,48022 + ,158492 + ,1 + ,5131 + ,71 + ,44 + ,86111 + ,292154 + ,1 + ,2397 + ,14 + ,13 + ,74020 + ,289513 + ,1 + ,1389 + ,9 + ,6 + ,57530 + ,378049 + ,1 + ,1503 + ,39 + ,35 + ,56364 + ,343466 + ,1 + ,402 + ,26 + ,8 + ,84990 + ,332743 + ,1 + ,2239 + ,21 + ,17 + ,88590 + ,442882 + ,1 + ,2234 + ,16 + ,11 + ,77200 + ,214215 + ,1 + ,837 + ,69 + ,21 + ,61262 + ,315688 + ,1 + ,10579 + ,92 + ,92 + ,110309 + ,375195 + ,1 + ,875 + ,14 + ,12 + ,67000 + ,334280 + ,1 + ,1585 + ,107 + ,112 + ,93099 + ,355864 + ,1 + ,1659 + ,29 + ,25 + ,107577 + ,480382 + ,1 + ,2647 + ,37 + ,17 + ,62920 + ,353058 + ,1 + ,3294 + ,23 + ,23 + ,75832 + ,217193 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,94 + ,7 + ,10 + ,60793 + ,314533 + ,1 + ,422 + ,28 + ,23 + ,57935 + ,318056 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,34 + ,8 + ,7 + ,60630 + ,314353 + ,1 + ,1558 + ,63 + ,25 + ,55637 + ,369448 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,43 + ,3 + ,20 + ,60887 + ,312846 + ,1 + ,645 + ,5 + ,4 + ,60720 + ,312075 + ,1 + ,316 + ,9 + ,4 + ,60505 + ,315009 + ,1 + ,115 + ,13 + ,10 + ,60945 + ,318903 + ,1 + ,5 + ,2 + ,1 + ,60720 + ,314887 + ,1 + ,897 + ,5 + ,4 + ,60720 + ,314913 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,389 + ,14 + ,8 + ,58990 + ,325506 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,1002 + ,15 + ,11 + ,56750 + ,298568 + ,1 + ,36 + ,3 + ,4 + ,60894 + ,315834 + ,1 + ,460 + ,15 + ,15 + ,63346 + ,329784 + ,1 + ,309 + ,11 + ,9 + ,56535 + ,312878 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,9 + ,6 + ,7 + ,60835 + ,314987 + ,1 + ,271 + ,2 + ,2 + ,60720 + ,325249 + ,1 + ,14 + ,1 + ,0 + ,61016 + ,315877 + ,1 + ,520 + ,10 + ,7 + ,58650 + ,291650 + ,1 + ,1766 + ,73 + ,46 + ,60438 + ,305959 + ,1 + ,0 + ,0 + ,5 + ,60720 + ,315380 + ,1 + ,458 + ,11 + ,7 + ,58625 + ,297765 + ,1 + ,20 + ,3 + ,2 + ,60938 + ,315245 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,98 + ,2 + ,2 + ,61490 + ,315236 + ,1 + ,405 + ,7 + ,5 + ,60845 + ,336425 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,483 + ,27 + ,7 + ,60830 + ,306268 + ,1 + ,454 + ,51 + ,24 + ,63261 + ,302187 + ,1 + ,47 + ,3 + ,1 + ,60720 + ,314882 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,757 + ,19 + ,18 + ,45689 + ,382712 + ,1 + ,4655 + ,393 + ,55 + ,60720 + ,341570 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,36 + ,4 + ,3 + ,61564 + ,312412 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,203 + ,9 + ,9 + ,61938 + ,309596 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,126 + ,10 + ,8 + ,60951 + ,315547 + ,1 + ,400 + ,152 + ,113 + ,60720 + ,313267 + ,1 + ,71 + ,1 + ,0 + ,60745 + ,316176 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,972 + ,34 + ,19 + ,71642 + ,359335 + ,1 + ,531 + ,10 + ,11 + ,71641 + ,330068 + ,1 + ,2461 + ,57 + ,25 + ,55792 + ,314289 + ,1 + ,378 + ,52 + ,16 + ,71873 + ,297413 + ,1 + ,23 + ,5 + ,5 + ,62555 + ,314806 + ,1 + ,638 + ,14 + ,11 + ,60370 + ,333210 + ,1 + ,2300 + ,29 + ,23 + ,64873 + ,352108 + ,1 + ,149 + ,5 + ,6 + ,62041 + ,313332 + ,1 + ,226 + ,5 + ,5 + ,65745 + ,291787 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,275 + ,4 + ,7 + ,59500 + ,318745 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,141 + ,6 + ,7 + ,61630 + ,315366 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,28 + ,2 + ,3 + ,60890 + ,315688 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,4980 + ,91 + ,89 + ,113521 + ,409642 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,472 + ,20 + ,19 + ,80045 + ,269587 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,203 + ,27 + ,12 + ,50804 + ,300962 + ,1 + ,496 + ,17 + ,12 + ,87390 + ,325479 + ,1 + ,10 + ,2 + ,5 + ,61656 + ,316155 + ,1 + ,63 + ,4 + ,2 + ,65688 + ,318574 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,1136 + ,32 + ,26 + ,48522 + ,343613 + ,1 + ,265 + ,31 + ,3 + ,60720 + ,306948 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,267 + ,32 + ,11 + ,57640 + ,330059 + ,1 + ,474 + ,20 + ,10 + ,61977 + ,288985 + ,1 + ,534 + ,7 + ,5 + ,62620 + ,304485 + ,1 + ,0 + ,0 + ,2 + ,60720 + ,315380 + ,1 + ,15 + ,8 + ,6 + ,60831 + ,315688 + ,1 + ,397 + ,28 + ,7 + ,60646 + ,317736 + ,1 + ,0 + ,0 + ,2 + ,60720 + ,315380 + ,1 + ,1866 + ,29 + ,28 + ,56225 + ,322331 + ,1 + ,288 + ,4 + ,3 + ,60510 + ,296656 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,3 + ,2 + ,1 + ,60698 + ,315354 + ,1 + ,468 + ,21 + ,20 + ,60720 + ,312161 + ,1 + ,20 + ,2 + ,1 + ,60805 + ,315576 + ,1 + ,278 + ,26 + ,22 + ,61404 + ,314922 + ,1 + ,61 + ,14 + ,9 + ,60720 + ,314551 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,192 + ,4 + ,2 + ,65276 + ,312339 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,317 + ,9 + ,7 + ,63915 + ,298700 + ,1 + ,738 + ,10 + ,9 + ,60720 + ,321376 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,368 + ,17 + ,13 + ,61686 + ,303230 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,2 + ,1 + ,0 + ,60743 + ,315487 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,53 + ,6 + ,6 + ,60349 + ,315793 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,94 + ,3 + ,3 + ,61360 + ,312887 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,24 + ,8 + ,7 + ,59818 + ,315637 + ,1 + ,2332 + ,4 + ,2 + ,72680 + ,324385 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,131 + ,11 + ,15 + ,61808 + ,308989 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,206 + ,9 + ,9 + ,53110 + ,296702 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,167 + ,2 + ,1 + ,64245 + ,307322 + ,1 + ,622 + ,73 + ,38 + ,73007 + ,304376 + ,1 + ,2328 + ,94 + ,57 + ,82732 + ,253588 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,365 + ,8 + ,7 + ,54820 + ,309560 + ,1 + ,364 + ,35 + ,26 + ,47705 + ,298466 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,226 + ,12 + ,13 + ,72835 + ,343929 + ,1 + ,307 + ,15 + ,10 + ,58856 + ,331955 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,188 + ,11 + ,9 + ,77655 + ,381180 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,138 + ,6 + ,26 + ,69817 + ,331420 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,125 + ,12 + ,19 + ,60798 + ,310201 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,282 + ,30 + ,12 + ,62452 + ,320016 + ,1 + ,335 + ,33 + ,23 + ,64175 + ,320398 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,1324 + ,117 + ,29 + ,67440 + ,291841 + ,1 + ,176 + ,28 + ,8 + ,68136 + ,310670 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,249 + ,72 + ,26 + ,56726 + ,313491 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,333 + ,13 + ,9 + ,70811 + ,331323 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,601 + ,6 + ,5 + ,60720 + ,319210 + ,1 + ,30 + ,4 + ,3 + ,62045 + ,318098 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,249 + ,62 + ,13 + ,54323 + ,292754 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,165 + ,24 + ,12 + ,62841 + ,325176 + ,1 + ,453 + ,21 + ,19 + ,81125 + ,365959 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,53 + ,14 + ,10 + ,59506 + ,302409 + ,1 + ,382 + ,21 + ,9 + ,59365 + ,340968 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,9 + ,60720 + ,315380 + ,1 + ,30 + ,4 + ,4 + ,60798 + ,313164 + ,1 + ,290 + ,2 + ,1 + ,58790 + ,301164 + ,1 + ,0 + ,0 + ,1 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,366 + ,53 + ,14 + ,61808 + ,344425 + ,1 + ,2 + ,9 + ,12 + ,60735 + ,315394 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,209 + ,13 + ,19 + ,64016 + ,316647 + ,1 + ,384 + ,22 + ,17 + ,54683 + ,309836 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,365 + ,83 + ,32 + ,87192 + ,346611 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,49 + ,8 + ,14 + ,64107 + ,322031 + ,1 + ,3 + ,4 + ,8 + ,60761 + ,315656 + ,1 + ,133 + ,14 + ,4 + ,65990 + ,339445 + ,1 + ,32 + ,1 + ,0 + ,59988 + ,314964 + ,1 + ,368 + ,17 + ,20 + ,61167 + ,297141 + ,1 + ,1 + ,6 + ,5 + ,60719 + ,315372 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,22 + ,2 + ,1 + ,60722 + ,312502 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,96 + ,5 + ,4 + ,60379 + ,313729 + ,1 + ,1 + ,2 + ,1 + ,60727 + ,315388 + ,1 + ,314 + ,5 + ,4 + ,60720 + ,315371 + ,1 + ,844 + ,78 + ,20 + ,60925 + ,296139 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,26 + ,1 + ,1 + ,60896 + ,313880 + ,1 + ,125 + ,13 + ,10 + ,59734 + ,317698 + ,1 + ,304 + ,15 + ,12 + ,62969 + ,295580 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,621 + ,48 + ,13 + ,60720 + ,308256 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,119 + ,6 + ,3 + ,59118 + ,303677 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,1595 + ,17 + ,10 + ,60720 + ,319369 + ,1 + ,312 + ,14 + ,3 + ,58598 + ,318690 + ,1 + ,60 + ,10 + ,7 + ,61124 + ,314049 + ,1 + ,587 + ,12 + ,10 + ,59595 + ,325699 + ,1 + ,135 + ,2 + ,1 + ,62065 + ,314210 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,514 + ,52 + ,15 + ,78780 + ,322378 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,1 + ,4 + ,4 + ,60722 + ,315398 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,1763 + ,24 + ,28 + ,61600 + ,308336 + ,1 + ,180 + ,11 + ,9 + ,59635 + ,316386 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,218 + ,21 + ,7 + ,60720 + ,315553 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,448 + ,40 + ,7 + ,59781 + ,323361 + ,1 + ,227 + ,9 + ,7 + ,76644 + ,336639 + ,1 + ,174 + ,1 + ,3 + ,64820 + ,307424 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,121 + ,24 + ,11 + ,56178 + ,295370 + ,1 + ,607 + ,11 + ,7 + ,60436 + ,322340 + ,1 + ,2212 + ,14 + ,10 + ,60720 + ,319864 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,530 + ,60 + ,18 + ,73433 + ,317291 + ,1 + ,571 + ,80 + ,14 + ,41477 + ,280398 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,78 + ,16 + ,12 + ,62700 + ,317330 + ,1 + ,2489 + ,40 + ,29 + ,67804 + ,238125 + ,1 + ,131 + ,6 + ,3 + ,59661 + ,327071 + ,1 + ,923 + ,8 + ,6 + ,58620 + ,309038 + ,1 + ,72 + ,3 + ,3 + ,60398 + ,314210 + ,1 + ,572 + ,16 + ,8 + ,58580 + ,307930 + ,1 + ,397 + ,10 + ,10 + ,62710 + ,322327 + ,1 + ,450 + ,8 + ,6 + ,59325 + ,292136 + ,1 + ,622 + ,7 + ,8 + ,60950 + ,263276 + ,1 + ,694 + ,8 + ,6 + ,68060 + ,367655 + ,1 + ,3425 + ,12 + ,9 + ,83620 + ,283910 + ,1 + ,562 + ,13 + ,8 + ,58456 + ,283587 + ,1 + ,4917 + ,42 + ,26 + ,52811 + ,243650 + ,1 + ,1442 + ,118 + ,239 + ,121173 + ,438493 + ,1 + ,529 + ,9 + ,7 + ,63870 + ,296261 + ,1 + ,2126 + ,138 + ,41 + ,21001 + ,230621 + ,1 + ,1061 + ,5 + ,3 + ,70415 + ,304252 + ,1 + ,776 + ,9 + ,8 + ,64230 + ,333505 + ,1 + ,611 + ,8 + ,6 + ,59190 + ,296919 + ,1 + ,1526 + ,25 + ,21 + ,69351 + ,278990 + ,1 + ,592 + ,7 + ,7 + ,64270 + ,276898 + ,1 + ,1182 + ,13 + ,11 + ,70694 + ,327007 + ,1 + ,621 + ,16 + ,11 + ,68005 + ,317046 + ,1 + ,989 + ,11 + ,12 + ,58930 + ,304555 + ,1 + ,438 + ,11 + ,9 + ,58320 + ,298096 + ,1 + ,726 + ,3 + ,3 + ,69980 + ,231861 + ,1 + ,1303 + ,61 + ,57 + ,69863 + ,309422 + ,1 + ,7419 + ,29 + ,21 + ,63255 + ,286963 + ,1 + ,1164 + ,17 + ,15 + ,57320 + ,269753 + ,1 + ,3310 + ,33 + ,32 + ,75230 + ,448243 + ,1 + ,1920 + ,15 + ,11 + ,79420 + ,165404 + ,1 + ,965 + ,3 + ,2 + ,73490 + ,204325 + ,1 + ,3256 + ,66 + ,23 + ,35250 + ,407159 + ,1 + ,1135 + ,17 + ,20 + ,62285 + ,290476 + ,1 + ,1270 + ,26 + ,24 + ,69206 + ,275311 + ,1 + ,661 + ,3 + ,1 + ,65920 + ,246541 + ,1 + ,1013 + ,2 + ,1 + ,69770 + ,253468 + ,1 + ,2844 + ,67 + ,74 + ,72683 + ,240897 + ,1 + ,11528 + ,70 + ,68 + ,-14545 + ,-83265 + ,1 + ,6526 + ,26 + ,20 + ,55830 + ,-42143 + ,1 + ,2264 + ,24 + ,20 + ,55174 + ,272713 + ,1 + ,5109 + ,97 + ,82 + ,67038 + ,215362 + ,1 + ,3999 + ,30 + ,21 + ,51252 + ,42754 + ,1 + ,35624 + ,223 + ,244 + ,157278 + ,306275 + ,1 + ,9252 + ,48 + ,32 + ,79510 + ,253537 + ,1 + ,15236 + ,90 + ,86 + ,77440 + ,372631 + ,1 + ,18073 + ,180 + ,69 + ,27284 + ,-7170 + ,1) + ,dim=c(6 + ,431) + ,dimnames=list(c('Costs' + ,'Trades' + ,'Orders' + ,'Dividends' + ,'Wealth' + ,'Group') + ,1:431)) > y <- array(NA,dim=c(6,431),dimnames=list(c('Costs','Trades','Orders','Dividends','Wealth','Group'),1:431)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par3 = 'No Linear Trend' > par2 = 'Do not include Seasonal Dummies' > par1 = '5' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > library(lattice) > library(lmtest) Loading required package: zoo Attaching package: 'zoo' The following object(s) are masked from package:base : as.Date.numeric > n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test > par1 <- as.numeric(par1) > x <- t(y) > k <- length(x[1,]) > n <- length(x[,1]) > x1 <- cbind(x[,par1], x[,1:k!=par1]) > mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1]) > colnames(x1) <- mycolnames #colnames(x)[par1] > x <- x1 > if (par3 == 'First Differences'){ + x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep=''))) + for (i in 1:n-1) { + for (j in 1:k) { + x2[i,j] <- x[i+1,j] - x[i,j] + } + } + x <- x2 + } > if (par2 == 'Include Monthly Dummies'){ + x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep =''))) + for (i in 1:11){ + x2[seq(i,n,12),i] <- 1 + } + x <- cbind(x, x2) + } > if (par2 == 'Include Quarterly Dummies'){ + x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep =''))) + for (i in 1:3){ + x2[seq(i,n,4),i] <- 1 + } + x <- cbind(x, x2) + } > k <- length(x[1,]) > if (par3 == 'Linear Trend'){ + x <- cbind(x, c(1:n)) + colnames(x)[k+1] <- 't' + } > x Wealth Costs Trades Orders Dividends Group 1 6282154 162556 1081 807 213118 1 2 4321023 29790 309 444 81767 1 3 4111912 87550 458 412 153198 1 4 223193 84738 588 428 -26007 1 5 1491348 54660 302 315 126942 1 6 1629616 42634 156 168 157214 1 7 1398893 40949 481 263 129352 1 8 1926517 45187 353 267 234817 1 9 983660 37704 452 228 60448 1 10 1443586 16275 109 129 47818 1 11 1073089 25830 115 104 245546 1 12 984885 12679 110 122 48020 1 13 1405225 18014 239 393 -1710 1 14 227132 43556 247 190 32648 1 15 929118 24811 505 280 95350 1 16 1071292 6575 159 63 151352 1 17 638830 7123 109 102 288170 1 18 856956 21950 519 265 114337 1 19 992426 37597 248 234 37884 1 20 444477 17821 373 277 122844 1 21 857217 12988 119 73 82340 1 22 711969 22330 84 67 79801 1 23 702380 13326 102 103 165548 1 24 358589 16189 295 290 116384 1 25 297978 7146 105 83 134028 1 26 585715 15824 64 56 63838 1 27 657954 27664 282 236 74996 1 28 209458 11920 182 73 31080 1 29 786690 8568 37 34 32168 1 30 439798 14416 361 139 49857 1 31 688779 3369 28 26 87161 1 32 574339 11819 85 70 106113 1 33 741409 6984 45 40 80570 1 34 597793 4519 49 42 102129 1 35 644190 2220 22 12 301670 1 36 377934 18562 155 211 102313 1 37 640273 10327 91 74 88577 1 38 697458 5336 81 80 112477 1 39 550608 2365 79 83 191778 1 40 207393 4069 145 131 79804 1 41 301607 8636 855 203 128294 1 42 345783 13718 61 56 96448 1 43 501749 4525 226 89 93811 1 44 379983 6869 105 88 117520 1 45 387475 4628 62 39 69159 1 46 377305 3689 25 25 101792 1 47 370837 4891 217 49 210568 1 48 430866 7489 322 149 136996 1 49 469107 4901 84 58 121920 1 50 194493 2284 33 41 76403 1 51 530670 3160 108 90 108094 1 52 518365 4150 150 136 134759 1 53 491303 7285 115 97 188873 1 54 527021 1134 162 63 146216 1 55 233773 4658 158 114 156608 1 56 405972 2384 97 77 61348 1 57 652925 3748 9 6 50350 1 58 446211 5371 66 47 87720 1 59 341340 1285 107 51 99489 1 60 387699 9327 101 85 87419 1 61 493408 5565 47 43 94355 1 62 146494 1528 38 32 60326 1 63 414462 3122 34 25 94670 1 64 364304 7561 87 77 82425 1 65 355178 2675 79 54 59017 1 66 357760 13253 947 251 90829 1 67 261216 880 74 15 80791 1 68 397144 2053 53 44 100423 1 69 374943 1424 94 73 131116 1 70 424898 4036 63 85 100269 1 71 202055 3045 58 49 27330 1 72 378525 5119 49 38 39039 1 73 310768 1431 34 35 106885 1 74 325738 554 11 9 79285 1 75 394510 1975 35 34 118881 1 76 247060 1765 20 20 77623 1 77 368078 1012 47 29 114768 1 78 236761 810 43 11 74015 1 79 312378 1280 117 52 69465 1 80 339836 666 171 13 117869 1 81 347385 1380 26 29 60982 1 82 426280 4677 75 66 90131 1 83 352850 876 59 33 138971 1 84 301881 814 18 15 39625 1 85 377516 514 15 15 102725 1 86 357312 5692 72 68 64239 1 87 458343 3642 86 100 90262 1 88 354228 540 14 13 103960 1 89 308636 2099 64 45 106611 1 90 386212 567 11 14 103345 1 91 393343 2001 52 36 95551 1 92 378509 2949 41 40 82903 1 93 452469 2253 99 68 63593 1 94 364839 6533 75 29 126910 1 95 358649 1889 45 43 37527 1 96 376641 3055 43 30 60247 1 97 429112 272 8 9 112995 1 98 330546 1414 198 22 70184 1 99 403560 2564 22 19 130140 1 100 317892 1383 11 9 73221 1 101 307528 1261 33 31 76114 1 102 235133 975 23 19 90534 1 103 299243 3366 80 55 108479 1 104 314073 576 18 8 113761 1 105 368186 1686 40 28 68696 1 106 269661 746 23 29 71561 1 107 125390 3192 60 48 59831 1 108 510834 2045 20 16 97890 1 109 321896 5702 61 47 101481 1 110 249898 1932 36 20 72954 1 111 408881 936 30 22 67939 1 112 158492 3437 47 33 48022 1 113 292154 5131 71 44 86111 1 114 289513 2397 14 13 74020 1 115 378049 1389 9 6 57530 1 116 343466 1503 39 35 56364 1 117 332743 402 26 8 84990 1 118 442882 2239 21 17 88590 1 119 214215 2234 16 11 77200 1 120 315688 837 69 21 61262 1 121 375195 10579 92 92 110309 1 122 334280 875 14 12 67000 1 123 355864 1585 107 112 93099 1 124 480382 1659 29 25 107577 1 125 353058 2647 37 17 62920 1 126 217193 3294 23 23 75832 1 127 315380 0 0 0 60720 1 128 314533 94 7 10 60793 1 129 318056 422 28 23 57935 1 130 315380 0 0 0 60720 1 131 314353 34 8 7 60630 1 132 369448 1558 63 25 55637 1 133 315380 0 0 0 60720 1 134 312846 43 3 20 60887 1 135 312075 645 5 4 60720 1 136 315009 316 9 4 60505 1 137 318903 115 13 10 60945 1 138 314887 5 2 1 60720 1 139 314913 897 5 4 60720 1 140 315380 0 0 0 60720 1 141 325506 389 14 8 58990 1 142 315380 0 0 0 60720 1 143 298568 1002 15 11 56750 1 144 315834 36 3 4 60894 1 145 329784 460 15 15 63346 1 146 312878 309 11 9 56535 1 147 315380 0 0 0 60720 1 148 314987 9 6 7 60835 1 149 325249 271 2 2 60720 1 150 315877 14 1 0 61016 1 151 291650 520 10 7 58650 1 152 305959 1766 73 46 60438 1 153 315380 0 0 5 60720 1 154 297765 458 11 7 58625 1 155 315245 20 3 2 60938 1 156 315380 0 0 0 60720 1 157 315380 0 0 0 60720 1 158 315236 98 2 2 61490 1 159 336425 405 7 5 60845 1 160 315380 0 0 0 60720 1 161 315380 0 0 0 60720 1 162 315380 0 0 0 60720 1 163 315380 0 0 0 60720 1 164 306268 483 27 7 60830 1 165 302187 454 51 24 63261 1 166 314882 47 3 1 60720 1 167 315380 0 0 0 60720 1 168 382712 757 19 18 45689 1 169 341570 4655 393 55 60720 1 170 315380 0 0 0 60720 1 171 315380 0 0 0 60720 1 172 312412 36 4 3 61564 1 173 315380 0 0 0 60720 1 174 309596 203 9 9 61938 1 175 315380 0 0 0 60720 1 176 315547 126 10 8 60951 1 177 313267 400 152 113 60720 1 178 316176 71 1 0 60745 1 179 315380 0 0 0 60720 1 180 315380 0 0 0 60720 1 181 359335 972 34 19 71642 1 182 330068 531 10 11 71641 1 183 314289 2461 57 25 55792 1 184 297413 378 52 16 71873 1 185 314806 23 5 5 62555 1 186 333210 638 14 11 60370 1 187 352108 2300 29 23 64873 1 188 313332 149 5 6 62041 1 189 291787 226 5 5 65745 1 190 315380 0 0 0 60720 1 191 318745 275 4 7 59500 1 192 315380 0 0 0 60720 1 193 315366 141 6 7 61630 1 194 315380 0 0 0 60720 1 195 315688 28 2 3 60890 1 196 315380 0 0 0 60720 1 197 409642 4980 91 89 113521 1 198 315380 0 0 0 60720 1 199 315380 0 0 0 60720 1 200 269587 472 20 19 80045 1 201 315380 0 0 0 60720 1 202 315380 0 0 0 60720 1 203 315380 0 0 0 60720 1 204 300962 203 27 12 50804 1 205 325479 496 17 12 87390 1 206 316155 10 2 5 61656 1 207 318574 63 4 2 65688 1 208 315380 0 0 0 60720 1 209 343613 1136 32 26 48522 1 210 306948 265 31 3 60720 1 211 315380 0 0 0 60720 1 212 315380 0 0 0 60720 1 213 330059 267 32 11 57640 1 214 288985 474 20 10 61977 1 215 304485 534 7 5 62620 1 216 315380 0 0 2 60720 1 217 315688 15 8 6 60831 1 218 317736 397 28 7 60646 1 219 315380 0 0 2 60720 1 220 322331 1866 29 28 56225 1 221 296656 288 4 3 60510 1 222 315380 0 0 0 60720 1 223 315354 3 2 1 60698 1 224 312161 468 21 20 60720 1 225 315576 20 2 1 60805 1 226 314922 278 26 22 61404 1 227 314551 61 14 9 60720 1 228 315380 0 0 0 60720 1 229 312339 192 4 2 65276 1 230 315380 0 0 0 60720 1 231 298700 317 9 7 63915 1 232 321376 738 10 9 60720 1 233 315380 0 0 0 60720 1 234 303230 368 17 13 61686 1 235 315380 0 0 0 60720 1 236 315487 2 1 0 60743 1 237 315380 0 0 0 60720 1 238 315793 53 6 6 60349 1 239 315380 0 0 0 60720 1 240 315380 0 0 0 60720 1 241 315380 0 0 0 60720 1 242 312887 94 3 3 61360 1 243 315380 0 0 0 60720 1 244 315637 24 8 7 59818 1 245 324385 2332 4 2 72680 1 246 315380 0 0 0 60720 1 247 315380 0 0 0 60720 1 248 308989 131 11 15 61808 1 249 315380 0 0 0 60720 1 250 315380 0 0 0 60720 1 251 296702 206 9 9 53110 1 252 315380 0 0 0 60720 1 253 307322 167 2 1 64245 1 254 304376 622 73 38 73007 1 255 253588 2328 94 57 82732 1 256 315380 0 0 0 60720 1 257 309560 365 8 7 54820 1 258 298466 364 35 26 47705 1 259 315380 0 0 0 60720 1 260 315380 0 0 0 60720 1 261 315380 0 0 0 60720 1 262 315380 0 0 0 60720 1 263 343929 226 12 13 72835 1 264 331955 307 15 10 58856 1 265 315380 0 0 0 60720 1 266 315380 0 0 0 60720 1 267 315380 0 0 0 60720 1 268 381180 188 11 9 77655 1 269 315380 0 0 0 60720 1 270 331420 138 6 26 69817 1 271 315380 0 0 0 60720 1 272 315380 0 0 0 60720 1 273 315380 0 0 0 60720 1 274 310201 125 12 19 60798 1 275 315380 0 0 0 60720 1 276 320016 282 30 12 62452 1 277 320398 335 33 23 64175 1 278 315380 0 0 0 60720 1 279 291841 1324 117 29 67440 1 280 310670 176 28 8 68136 1 281 315380 0 0 0 60720 1 282 315380 0 0 0 60720 1 283 313491 249 72 26 56726 1 284 315380 0 0 0 60720 1 285 331323 333 13 9 70811 1 286 315380 0 0 0 60720 1 287 319210 601 6 5 60720 1 288 318098 30 4 3 62045 1 289 315380 0 0 0 60720 1 290 292754 249 62 13 54323 1 291 315380 0 0 0 60720 1 292 325176 165 24 12 62841 1 293 365959 453 21 19 81125 1 294 315380 0 0 0 60720 1 295 302409 53 14 10 59506 1 296 340968 382 21 9 59365 1 297 315380 0 0 0 60720 1 298 315380 0 0 0 60720 1 299 315380 0 0 0 60720 1 300 315380 0 0 9 60720 1 301 313164 30 4 4 60798 1 302 301164 290 2 1 58790 1 303 315380 0 0 1 60720 1 304 315380 0 0 0 60720 1 305 344425 366 53 14 61808 1 306 315394 2 9 12 60735 1 307 315380 0 0 0 60720 1 308 316647 209 13 19 64016 1 309 309836 384 22 17 54683 1 310 315380 0 0 0 60720 1 311 315380 0 0 0 60720 1 312 346611 365 83 32 87192 1 313 315380 0 0 0 60720 1 314 322031 49 8 14 64107 1 315 315656 3 4 8 60761 1 316 339445 133 14 4 65990 1 317 314964 32 1 0 59988 1 318 297141 368 17 20 61167 1 319 315372 1 6 5 60719 1 320 315380 0 0 0 60720 1 321 315380 0 0 0 60720 1 322 315380 0 0 0 60720 1 323 315380 0 0 0 60720 1 324 315380 0 0 0 60720 1 325 315380 0 0 0 60720 1 326 312502 22 2 1 60722 1 327 315380 0 0 0 60720 1 328 315380 0 0 0 60720 1 329 315380 0 0 0 60720 1 330 315380 0 0 0 60720 1 331 315380 0 0 0 60720 1 332 315380 0 0 0 60720 1 333 315380 0 0 0 60720 1 334 313729 96 5 4 60379 1 335 315388 1 2 1 60727 1 336 315371 314 5 4 60720 1 337 296139 844 78 20 60925 1 338 315380 0 0 0 60720 1 339 313880 26 1 1 60896 1 340 317698 125 13 10 59734 1 341 295580 304 15 12 62969 1 342 315380 0 0 0 60720 1 343 315380 0 0 0 60720 1 344 315380 0 0 0 60720 1 345 308256 621 48 13 60720 1 346 315380 0 0 0 60720 1 347 303677 119 6 3 59118 1 348 315380 0 0 0 60720 1 349 315380 0 0 0 60720 1 350 319369 1595 17 10 60720 1 351 318690 312 14 3 58598 1 352 314049 60 10 7 61124 1 353 325699 587 12 10 59595 1 354 314210 135 2 1 62065 1 355 315380 0 0 0 60720 1 356 315380 0 0 0 60720 1 357 322378 514 52 15 78780 1 358 315380 0 0 0 60720 1 359 315380 0 0 0 60720 1 360 315380 0 0 0 60720 1 361 315398 1 4 4 60722 1 362 315380 0 0 0 60720 1 363 315380 0 0 0 60720 1 364 308336 1763 24 28 61600 1 365 316386 180 11 9 59635 1 366 315380 0 0 0 60720 1 367 315380 0 0 0 60720 1 368 315380 0 0 0 60720 1 369 315380 0 0 0 60720 1 370 315553 218 21 7 60720 1 371 315380 0 0 0 60720 1 372 323361 448 40 7 59781 1 373 336639 227 9 7 76644 1 374 307424 174 1 3 64820 1 375 315380 0 0 0 60720 1 376 315380 0 0 0 60720 1 377 295370 121 24 11 56178 1 378 322340 607 11 7 60436 1 379 319864 2212 14 10 60720 1 380 315380 0 0 0 60720 1 381 315380 0 0 0 60720 1 382 317291 530 60 18 73433 1 383 280398 571 80 14 41477 1 384 315380 0 0 0 60720 1 385 317330 78 16 12 62700 1 386 238125 2489 40 29 67804 1 387 327071 131 6 3 59661 1 388 309038 923 8 6 58620 1 389 314210 72 3 3 60398 1 390 307930 572 16 8 58580 1 391 322327 397 10 10 62710 1 392 292136 450 8 6 59325 1 393 263276 622 7 8 60950 1 394 367655 694 8 6 68060 1 395 283910 3425 12 9 83620 1 396 283587 562 13 8 58456 1 397 243650 4917 42 26 52811 1 398 438493 1442 118 239 121173 1 399 296261 529 9 7 63870 1 400 230621 2126 138 41 21001 1 401 304252 1061 5 3 70415 1 402 333505 776 9 8 64230 1 403 296919 611 8 6 59190 1 404 278990 1526 25 21 69351 1 405 276898 592 7 7 64270 1 406 327007 1182 13 11 70694 1 407 317046 621 16 11 68005 1 408 304555 989 11 12 58930 1 409 298096 438 11 9 58320 1 410 231861 726 3 3 69980 1 411 309422 1303 61 57 69863 1 412 286963 7419 29 21 63255 1 413 269753 1164 17 15 57320 1 414 448243 3310 33 32 75230 1 415 165404 1920 15 11 79420 1 416 204325 965 3 2 73490 1 417 407159 3256 66 23 35250 1 418 290476 1135 17 20 62285 1 419 275311 1270 26 24 69206 1 420 246541 661 3 1 65920 1 421 253468 1013 2 1 69770 1 422 240897 2844 67 74 72683 1 423 -83265 11528 70 68 -14545 1 424 -42143 6526 26 20 55830 1 425 272713 2264 24 20 55174 1 426 215362 5109 97 82 67038 1 427 42754 3999 30 21 51252 1 428 306275 35624 223 244 157278 1 429 253537 9252 48 32 79510 1 430 372631 15236 90 86 77440 1 431 -7170 18073 180 69 27284 1 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Costs Trades Orders Dividends Group 161391.54 21.71 -1177.26 2571.27 1.71 NA > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -2141698 -11603 39962 50154 2595180 Coefficients: (1 not defined because of singularities) Estimate Std. Error t value Pr(>|t|) (Intercept) 1.614e+05 2.957e+04 5.457 8.23e-08 *** Costs 2.171e+01 2.080e+00 10.435 < 2e-16 *** Trades -1.177e+03 2.072e+02 -5.683 2.46e-08 *** Orders 2.571e+03 4.353e+02 5.907 7.13e-09 *** Dividends 1.710e+00 3.979e-01 4.297 2.14e-05 *** Group NA NA NA NA --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 244700 on 426 degrees of freedom Multiple R-squared: 0.6799, Adjusted R-squared: 0.6769 F-statistic: 226.2 on 4 and 426 DF, p-value: < 2.2e-16 > if (n > n25) { + kp3 <- k + 3 + nmkm3 <- n - k - 3 + gqarr <- array(NA, dim=c(nmkm3-kp3+1,3)) + numgqtests <- 0 + numsignificant1 <- 0 + numsignificant5 <- 0 + numsignificant10 <- 0 + for (mypoint in kp3:nmkm3) { + j <- 0 + numgqtests <- numgqtests + 1 + for (myalt in c('greater', 'two.sided', 'less')) { + j <- j + 1 + gqarr[mypoint-kp3+1,j] <- gqtest(mylm, point=mypoint, alternative=myalt)$p.value + } + if (gqarr[mypoint-kp3+1,2] < 0.01) numsignificant1 <- numsignificant1 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.05) numsignificant5 <- numsignificant5 + 1 + if (gqarr[mypoint-kp3+1,2] < 0.10) numsignificant10 <- numsignificant10 + 1 + } + gqarr + } [,1] [,2] [,3] [1,] 1.0000000 2.515822e-49 1.257911e-49 [2,] 1.0000000 1.627750e-67 8.138748e-68 [3,] 1.0000000 2.679774e-71 1.339887e-71 [4,] 1.0000000 6.279151e-76 3.139576e-76 [5,] 1.0000000 3.494986e-100 1.747493e-100 [6,] 1.0000000 1.990757e-102 9.953783e-103 [7,] 1.0000000 1.939438e-104 9.697188e-105 [8,] 1.0000000 7.733144e-114 3.866572e-114 [9,] 1.0000000 1.871358e-125 9.356792e-126 [10,] 1.0000000 1.976999e-127 9.884993e-128 [11,] 1.0000000 7.486700e-132 3.743350e-132 [12,] 1.0000000 7.166766e-138 3.583383e-138 [13,] 1.0000000 1.349846e-146 6.749230e-147 [14,] 1.0000000 3.517484e-149 1.758742e-149 [15,] 1.0000000 4.911581e-150 2.455791e-150 [16,] 1.0000000 3.176874e-158 1.588437e-158 [17,] 1.0000000 6.409976e-158 3.204988e-158 [18,] 1.0000000 1.010392e-159 5.051959e-160 [19,] 1.0000000 1.485710e-162 7.428551e-163 [20,] 1.0000000 6.738537e-163 3.369269e-163 [21,] 1.0000000 3.716340e-176 1.858170e-176 [22,] 1.0000000 5.937188e-178 2.968594e-178 [23,] 1.0000000 2.334157e-183 1.167079e-183 [24,] 1.0000000 1.332271e-184 6.661354e-185 [25,] 1.0000000 1.389609e-193 6.948047e-194 [26,] 1.0000000 7.856499e-196 3.928249e-196 [27,] 1.0000000 2.122262e-195 1.061131e-195 [28,] 1.0000000 1.021994e-198 5.109968e-199 [29,] 1.0000000 3.933539e-204 1.966769e-204 [30,] 1.0000000 1.057799e-209 5.288994e-210 [31,] 1.0000000 6.061923e-209 3.030962e-209 [32,] 1.0000000 1.006082e-209 5.030410e-210 [33,] 1.0000000 1.654319e-211 8.271594e-212 [34,] 1.0000000 8.719383e-211 4.359692e-211 [35,] 1.0000000 1.507764e-211 7.538820e-212 [36,] 1.0000000 1.473099e-210 7.365495e-211 [37,] 1.0000000 3.246371e-210 1.623186e-210 [38,] 1.0000000 3.126874e-209 1.563437e-209 [39,] 1.0000000 2.062624e-209 1.031312e-209 [40,] 1.0000000 2.170153e-208 1.085076e-208 [41,] 1.0000000 8.262369e-208 4.131184e-208 [42,] 1.0000000 4.387353e-208 2.193677e-208 [43,] 1.0000000 1.743049e-208 8.715246e-209 [44,] 1.0000000 3.699452e-208 1.849726e-208 [45,] 1.0000000 1.306157e-207 6.530783e-208 [46,] 1.0000000 2.762831e-207 1.381416e-207 [47,] 1.0000000 1.265950e-210 6.329751e-211 [48,] 1.0000000 2.850700e-210 1.425350e-210 [49,] 1.0000000 1.221471e-223 6.107355e-224 [50,] 1.0000000 5.133548e-224 2.566774e-224 [51,] 1.0000000 4.181609e-223 2.090804e-223 [52,] 1.0000000 9.841470e-223 4.920735e-223 [53,] 1.0000000 2.292071e-224 1.146035e-224 [54,] 1.0000000 6.980374e-226 3.490187e-226 [55,] 1.0000000 2.031426e-225 1.015713e-225 [56,] 1.0000000 1.033551e-224 5.167753e-225 [57,] 1.0000000 5.395132e-224 2.697566e-224 [58,] 1.0000000 1.073250e-223 5.366252e-224 [59,] 1.0000000 2.527631e-223 1.263816e-223 [60,] 1.0000000 2.414734e-222 1.207367e-222 [61,] 1.0000000 1.278074e-221 6.390371e-222 [62,] 1.0000000 6.543775e-221 3.271888e-221 [63,] 1.0000000 4.512225e-220 2.256112e-220 [64,] 1.0000000 3.661140e-221 1.830570e-221 [65,] 1.0000000 1.228519e-220 6.142597e-221 [66,] 1.0000000 1.034342e-219 5.171712e-220 [67,] 1.0000000 1.231993e-218 6.159964e-219 [68,] 1.0000000 2.837171e-218 1.418585e-218 [69,] 1.0000000 2.497805e-217 1.248902e-217 [70,] 1.0000000 3.128873e-217 1.564437e-217 [71,] 1.0000000 2.790149e-216 1.395075e-216 [72,] 1.0000000 5.051394e-216 2.525697e-216 [73,] 1.0000000 3.102935e-215 1.551467e-215 [74,] 1.0000000 5.099177e-215 2.549589e-215 [75,] 1.0000000 4.489274e-215 2.244637e-215 [76,] 1.0000000 2.385875e-214 1.192937e-214 [77,] 1.0000000 2.337370e-213 1.168685e-213 [78,] 1.0000000 5.901226e-213 2.950613e-213 [79,] 1.0000000 4.201245e-213 2.100622e-213 [80,] 1.0000000 3.101607e-212 1.550804e-212 [81,] 1.0000000 5.986371e-212 2.993186e-212 [82,] 1.0000000 5.878307e-211 2.939153e-211 [83,] 1.0000000 5.248667e-210 2.624334e-210 [84,] 1.0000000 3.381231e-209 1.690615e-209 [85,] 1.0000000 1.238189e-210 6.190946e-211 [86,] 1.0000000 1.184694e-209 5.923469e-210 [87,] 1.0000000 4.035286e-210 2.017643e-210 [88,] 1.0000000 3.436661e-210 1.718331e-210 [89,] 1.0000000 3.109763e-209 1.554882e-209 [90,] 1.0000000 1.108313e-208 5.541567e-209 [91,] 1.0000000 9.101125e-208 4.550562e-208 [92,] 1.0000000 8.690072e-207 4.345036e-207 [93,] 1.0000000 7.668926e-206 3.834463e-206 [94,] 1.0000000 5.594414e-207 2.797207e-207 [95,] 1.0000000 4.036421e-207 2.018210e-207 [96,] 1.0000000 5.126546e-208 2.563273e-208 [97,] 1.0000000 2.232424e-207 1.116212e-207 [98,] 1.0000000 7.987905e-207 3.993952e-207 [99,] 1.0000000 4.665068e-209 2.332534e-209 [100,] 1.0000000 2.267070e-210 1.133535e-210 [101,] 1.0000000 1.502716e-209 7.513578e-210 [102,] 1.0000000 2.679747e-209 1.339873e-209 [103,] 1.0000000 1.895993e-209 9.479967e-210 [104,] 1.0000000 1.186114e-209 5.930570e-210 [105,] 1.0000000 8.135705e-209 4.067853e-209 [106,] 1.0000000 6.374511e-208 3.187256e-208 [107,] 1.0000000 3.801434e-208 1.900717e-208 [108,] 1.0000000 1.447623e-207 7.238117e-208 [109,] 1.0000000 8.916738e-207 4.458369e-207 [110,] 1.0000000 7.665735e-207 3.832867e-207 [111,] 1.0000000 7.502641e-208 3.751320e-208 [112,] 1.0000000 6.575116e-207 3.287558e-207 [113,] 1.0000000 2.238028e-206 1.119014e-206 [114,] 1.0000000 2.000593e-205 1.000296e-205 [115,] 1.0000000 1.370455e-204 6.852274e-205 [116,] 1.0000000 2.230039e-204 1.115020e-204 [117,] 1.0000000 6.196529e-204 3.098264e-204 [118,] 1.0000000 1.936332e-204 9.681661e-205 [119,] 1.0000000 1.712943e-203 8.564714e-204 [120,] 1.0000000 1.640898e-202 8.204489e-203 [121,] 1.0000000 1.503011e-201 7.515054e-202 [122,] 1.0000000 1.337827e-200 6.689137e-201 [123,] 1.0000000 1.246052e-199 6.230261e-200 [124,] 1.0000000 1.076344e-199 5.381718e-200 [125,] 1.0000000 9.724060e-199 4.862030e-199 [126,] 1.0000000 1.010576e-197 5.052882e-198 [127,] 1.0000000 9.484489e-197 4.742245e-197 [128,] 1.0000000 8.632108e-196 4.316054e-196 [129,] 1.0000000 7.996162e-195 3.998081e-195 [130,] 1.0000000 7.265632e-194 3.632816e-194 [131,] 1.0000000 6.576576e-193 3.288288e-193 [132,] 1.0000000 5.943831e-192 2.971916e-192 [133,] 1.0000000 4.622609e-191 2.311304e-191 [134,] 1.0000000 4.178456e-190 2.089228e-190 [135,] 1.0000000 4.056805e-189 2.028402e-189 [136,] 1.0000000 3.732752e-188 1.866376e-188 [137,] 1.0000000 3.276373e-187 1.638186e-187 [138,] 1.0000000 2.839709e-186 1.419855e-186 [139,] 1.0000000 2.549433e-185 1.274716e-185 [140,] 1.0000000 2.365056e-184 1.182528e-184 [141,] 1.0000000 1.921846e-183 9.609232e-184 [142,] 1.0000000 1.716516e-182 8.582579e-183 [143,] 1.0000000 1.592595e-181 7.962976e-182 [144,] 1.0000000 1.571027e-180 7.855134e-181 [145,] 1.0000000 1.433169e-179 7.165846e-180 [146,] 1.0000000 1.333599e-178 6.667994e-179 [147,] 1.0000000 1.187450e-177 5.937252e-178 [148,] 1.0000000 1.045833e-176 5.229166e-177 [149,] 1.0000000 9.192722e-176 4.596361e-176 [150,] 1.0000000 8.195012e-175 4.097506e-175 [151,] 1.0000000 5.395219e-174 2.697610e-174 [152,] 1.0000000 4.713848e-173 2.356924e-173 [153,] 1.0000000 4.108832e-172 2.054416e-172 [154,] 1.0000000 3.572746e-171 1.786373e-171 [155,] 1.0000000 3.098797e-170 1.549398e-170 [156,] 1.0000000 2.724890e-169 1.362445e-169 [157,] 1.0000000 2.327582e-168 1.163791e-168 [158,] 1.0000000 2.005881e-167 1.002941e-167 [159,] 1.0000000 1.720128e-166 8.600638e-167 [160,] 1.0000000 1.807692e-167 9.038461e-168 [161,] 1.0000000 2.333231e-167 1.166616e-167 [162,] 1.0000000 2.043677e-166 1.021838e-166 [163,] 1.0000000 1.784568e-165 8.922841e-166 [164,] 1.0000000 1.586625e-164 7.933123e-165 [165,] 1.0000000 1.376561e-163 6.882804e-164 [166,] 1.0000000 1.238313e-162 6.191564e-163 [167,] 1.0000000 1.066957e-161 5.334785e-162 [168,] 1.0000000 9.307649e-161 4.653824e-161 [169,] 1.0000000 8.118268e-160 4.059134e-160 [170,] 1.0000000 6.870795e-159 3.435398e-159 [171,] 1.0000000 5.832917e-158 2.916458e-158 [172,] 1.0000000 4.935101e-157 2.467550e-157 [173,] 1.0000000 3.059000e-156 1.529500e-156 [174,] 1.0000000 2.690050e-155 1.345025e-155 [175,] 1.0000000 1.417271e-154 7.086354e-155 [176,] 1.0000000 8.619810e-154 4.309905e-154 [177,] 1.0000000 7.360859e-153 3.680430e-153 [178,] 1.0000000 4.718115e-152 2.359057e-152 [179,] 1.0000000 1.695661e-151 8.478306e-152 [180,] 1.0000000 1.443964e-150 7.219821e-151 [181,] 1.0000000 9.742892e-150 4.871446e-150 [182,] 1.0000000 8.008541e-149 4.004271e-149 [183,] 1.0000000 6.323211e-148 3.161605e-148 [184,] 1.0000000 5.159072e-147 2.579536e-147 [185,] 1.0000000 4.282630e-146 2.141315e-146 [186,] 1.0000000 3.467637e-145 1.733818e-145 [187,] 1.0000000 2.817132e-144 1.408566e-144 [188,] 1.0000000 2.263843e-143 1.131922e-143 [189,] 1.0000000 1.452887e-142 7.264433e-143 [190,] 1.0000000 1.159114e-141 5.795569e-142 [191,] 1.0000000 9.212557e-141 4.606279e-141 [192,] 1.0000000 9.468306e-141 4.734153e-141 [193,] 1.0000000 7.538689e-140 3.769344e-140 [194,] 1.0000000 5.979466e-139 2.989733e-139 [195,] 1.0000000 4.724574e-138 2.362287e-138 [196,] 1.0000000 3.402539e-137 1.701269e-137 [197,] 1.0000000 1.466479e-136 7.332393e-137 [198,] 1.0000000 1.167973e-135 5.839864e-136 [199,] 1.0000000 9.240698e-135 4.620349e-135 [200,] 1.0000000 7.183267e-134 3.591633e-134 [201,] 1.0000000 1.015072e-133 5.075362e-134 [202,] 1.0000000 7.836063e-133 3.918032e-133 [203,] 1.0000000 6.064210e-132 3.032105e-132 [204,] 1.0000000 4.674326e-131 2.337163e-131 [205,] 1.0000000 2.621341e-130 1.310671e-130 [206,] 1.0000000 1.876647e-129 9.383237e-130 [207,] 1.0000000 1.465672e-128 7.328358e-129 [208,] 1.0000000 1.119005e-127 5.595026e-128 [209,] 1.0000000 8.508765e-127 4.254383e-127 [210,] 1.0000000 6.115204e-126 3.057602e-126 [211,] 1.0000000 4.611770e-125 2.305885e-125 [212,] 1.0000000 2.130462e-124 1.065231e-124 [213,] 1.0000000 1.594754e-123 7.973768e-124 [214,] 1.0000000 1.184377e-122 5.921884e-123 [215,] 1.0000000 8.759047e-122 4.379523e-122 [216,] 1.0000000 6.555731e-121 3.277865e-121 [217,] 1.0000000 4.805955e-120 2.402977e-120 [218,] 1.0000000 3.567349e-119 1.783675e-119 [219,] 1.0000000 2.610558e-118 1.305279e-118 [220,] 1.0000000 1.889603e-117 9.448014e-118 [221,] 1.0000000 1.374224e-116 6.871118e-117 [222,] 1.0000000 9.864988e-116 4.932494e-116 [223,] 1.0000000 6.808898e-115 3.404449e-115 [224,] 1.0000000 4.483612e-114 2.241806e-114 [225,] 1.0000000 3.179390e-113 1.589695e-113 [226,] 1.0000000 2.301258e-112 1.150629e-112 [227,] 1.0000000 1.617789e-111 8.088943e-112 [228,] 1.0000000 1.131079e-110 5.655397e-111 [229,] 1.0000000 7.883336e-110 3.941668e-110 [230,] 1.0000000 5.458119e-109 2.729060e-109 [231,] 1.0000000 3.771026e-108 1.885513e-108 [232,] 1.0000000 2.594068e-107 1.297034e-107 [233,] 1.0000000 1.776650e-106 8.883252e-107 [234,] 1.0000000 1.229232e-105 6.146162e-106 [235,] 1.0000000 8.344773e-105 4.172387e-105 [236,] 1.0000000 5.594065e-104 2.797032e-104 [237,] 1.0000000 3.850963e-103 1.925481e-103 [238,] 1.0000000 2.578790e-102 1.289395e-102 [239,] 1.0000000 1.719195e-101 8.595975e-102 [240,] 1.0000000 1.173783e-100 5.868914e-101 [241,] 1.0000000 7.754472e-100 3.877236e-100 [242,] 1.0000000 5.099892e-99 2.549946e-99 [243,] 1.0000000 3.322920e-98 1.661460e-98 [244,] 1.0000000 2.165451e-97 1.082726e-97 [245,] 1.0000000 1.410003e-96 7.050015e-97 [246,] 1.0000000 7.596848e-96 3.798424e-96 [247,] 1.0000000 5.326289e-96 2.663145e-96 [248,] 1.0000000 3.452565e-95 1.726283e-95 [249,] 1.0000000 2.035622e-94 1.017811e-94 [250,] 1.0000000 1.100907e-93 5.504537e-94 [251,] 1.0000000 7.044470e-93 3.522235e-93 [252,] 1.0000000 4.486706e-92 2.243353e-92 [253,] 1.0000000 2.844336e-91 1.422168e-91 [254,] 1.0000000 1.794735e-90 8.973677e-91 [255,] 1.0000000 1.149242e-89 5.746208e-90 [256,] 1.0000000 5.724165e-89 2.862083e-89 [257,] 1.0000000 3.564094e-88 1.782047e-88 [258,] 1.0000000 2.208623e-87 1.104312e-87 [259,] 1.0000000 1.362139e-86 6.810697e-87 [260,] 1.0000000 7.000318e-86 3.500159e-86 [261,] 1.0000000 4.279122e-85 2.139561e-85 [262,] 1.0000000 2.681062e-84 1.340531e-84 [263,] 1.0000000 1.622884e-83 8.114419e-84 [264,] 1.0000000 9.775862e-83 4.887931e-83 [265,] 1.0000000 5.860043e-82 2.930022e-82 [266,] 1.0000000 3.585564e-81 1.792782e-81 [267,] 1.0000000 2.128082e-80 1.064041e-80 [268,] 1.0000000 1.257618e-79 6.288088e-80 [269,] 1.0000000 7.551075e-79 3.775537e-79 [270,] 1.0000000 4.416042e-78 2.208021e-78 [271,] 1.0000000 2.038218e-77 1.019109e-77 [272,] 1.0000000 1.124086e-76 5.620431e-77 [273,] 1.0000000 6.487233e-76 3.243617e-76 [274,] 1.0000000 3.724942e-75 1.862471e-75 [275,] 1.0000000 2.099107e-74 1.049554e-74 [276,] 1.0000000 1.193721e-73 5.968606e-74 [277,] 1.0000000 6.929445e-73 3.464723e-73 [278,] 1.0000000 3.900938e-72 1.950469e-72 [279,] 1.0000000 2.079466e-71 1.039733e-71 [280,] 1.0000000 1.165304e-70 5.826519e-71 [281,] 1.0000000 6.458805e-70 3.229403e-70 [282,] 1.0000000 3.592780e-69 1.796390e-69 [283,] 1.0000000 1.971206e-68 9.856032e-69 [284,] 1.0000000 1.065811e-67 5.329054e-68 [285,] 1.0000000 5.907912e-67 2.953956e-67 [286,] 1.0000000 3.190741e-66 1.595371e-66 [287,] 1.0000000 1.761667e-65 8.808334e-66 [288,] 1.0000000 6.948213e-65 3.474107e-65 [289,] 1.0000000 3.699138e-64 1.849569e-64 [290,] 1.0000000 1.958651e-63 9.793254e-64 [291,] 1.0000000 1.031406e-62 5.157028e-63 [292,] 1.0000000 5.411044e-62 2.705522e-62 [293,] 1.0000000 2.850996e-61 1.425498e-61 [294,] 1.0000000 1.510693e-60 7.553465e-61 [295,] 1.0000000 7.782702e-60 3.891351e-60 [296,] 1.0000000 3.985314e-59 1.992657e-59 [297,] 1.0000000 1.678718e-58 8.393590e-59 [298,] 1.0000000 8.572076e-58 4.286038e-58 [299,] 1.0000000 4.323331e-57 2.161666e-57 [300,] 1.0000000 2.226068e-56 1.113034e-56 [301,] 1.0000000 1.024739e-55 5.123697e-56 [302,] 1.0000000 5.081213e-55 2.540607e-55 [303,] 1.0000000 2.504779e-54 1.252390e-54 [304,] 1.0000000 9.347875e-54 4.673937e-54 [305,] 1.0000000 4.563244e-53 2.281622e-53 [306,] 1.0000000 2.248587e-52 1.124294e-52 [307,] 1.0000000 1.086372e-51 5.431859e-52 [308,] 1.0000000 4.993620e-51 2.496810e-51 [309,] 1.0000000 2.359901e-50 1.179951e-50 [310,] 1.0000000 1.154833e-49 5.774165e-50 [311,] 1.0000000 5.465493e-49 2.732747e-49 [312,] 1.0000000 2.558362e-48 1.279181e-48 [313,] 1.0000000 1.190063e-47 5.950317e-48 [314,] 1.0000000 5.500897e-47 2.750449e-47 [315,] 1.0000000 2.526569e-46 1.263285e-46 [316,] 1.0000000 1.153033e-45 5.765166e-46 [317,] 1.0000000 5.228094e-45 2.614047e-45 [318,] 1.0000000 2.385251e-44 1.192625e-44 [319,] 1.0000000 1.067523e-43 5.337617e-44 [320,] 1.0000000 4.746170e-43 2.373085e-43 [321,] 1.0000000 2.096079e-42 1.048040e-42 [322,] 1.0000000 9.194868e-42 4.597434e-42 [323,] 1.0000000 4.006188e-41 2.003094e-41 [324,] 1.0000000 1.733560e-40 8.667799e-41 [325,] 1.0000000 7.449740e-40 3.724870e-40 [326,] 1.0000000 3.191734e-39 1.595867e-39 [327,] 1.0000000 1.355464e-38 6.777319e-39 [328,] 1.0000000 5.654078e-38 2.827039e-38 [329,] 1.0000000 2.354331e-37 1.177166e-37 [330,] 1.0000000 9.773622e-37 4.886811e-37 [331,] 1.0000000 4.063498e-36 2.031749e-36 [332,] 1.0000000 1.634854e-35 8.174271e-36 [333,] 1.0000000 6.872209e-35 3.436105e-35 [334,] 1.0000000 2.773839e-34 1.386920e-34 [335,] 1.0000000 1.111154e-33 5.555768e-34 [336,] 1.0000000 4.417116e-33 2.208558e-33 [337,] 1.0000000 1.814076e-32 9.070378e-33 [338,] 1.0000000 7.106111e-32 3.553055e-32 [339,] 1.0000000 2.840853e-31 1.420426e-31 [340,] 1.0000000 1.095593e-30 5.477964e-31 [341,] 1.0000000 4.191259e-30 2.095630e-30 [342,] 1.0000000 1.448692e-29 7.243459e-30 [343,] 1.0000000 5.223625e-29 2.611812e-29 [344,] 1.0000000 1.990353e-28 9.951764e-29 [345,] 1.0000000 6.629771e-28 3.314885e-28 [346,] 1.0000000 2.484012e-27 1.242006e-27 [347,] 1.0000000 9.051021e-27 4.525510e-27 [348,] 1.0000000 3.269008e-26 1.634504e-26 [349,] 1.0000000 1.042356e-25 5.211780e-26 [350,] 1.0000000 3.713500e-25 1.856750e-25 [351,] 1.0000000 1.310946e-24 6.554730e-25 [352,] 1.0000000 4.585217e-24 2.292608e-24 [353,] 1.0000000 1.597244e-23 7.986220e-24 [354,] 1.0000000 5.482298e-23 2.741149e-23 [355,] 1.0000000 1.863518e-22 9.317588e-23 [356,] 1.0000000 6.147958e-22 3.073979e-22 [357,] 1.0000000 2.024518e-21 1.012259e-21 [358,] 1.0000000 6.684095e-21 3.342048e-21 [359,] 1.0000000 2.184017e-20 1.092009e-20 [360,] 1.0000000 7.061144e-20 3.530572e-20 [361,] 1.0000000 2.258433e-19 1.129217e-19 [362,] 1.0000000 7.403530e-19 3.701765e-19 [363,] 1.0000000 2.318897e-18 1.159449e-18 [364,] 1.0000000 7.372697e-18 3.686349e-18 [365,] 1.0000000 2.444477e-17 1.222239e-17 [366,] 1.0000000 7.929586e-17 3.964793e-17 [367,] 1.0000000 2.384904e-16 1.192452e-16 [368,] 1.0000000 7.085874e-16 3.542937e-16 [369,] 1.0000000 2.211903e-15 1.105951e-15 [370,] 1.0000000 6.098458e-15 3.049229e-15 [371,] 1.0000000 1.517084e-14 7.585421e-15 [372,] 1.0000000 4.292680e-14 2.146340e-14 [373,] 1.0000000 1.197973e-13 5.989865e-14 [374,] 1.0000000 3.403124e-13 1.701562e-13 [375,] 1.0000000 1.001212e-12 5.006062e-13 [376,] 1.0000000 2.708015e-12 1.354007e-12 [377,] 1.0000000 7.585821e-12 3.792910e-12 [378,] 1.0000000 1.828540e-11 9.142700e-12 [379,] 1.0000000 4.461269e-11 2.230635e-11 [380,] 1.0000000 1.116559e-10 5.582793e-11 [381,] 1.0000000 2.851250e-10 1.425625e-10 [382,] 1.0000000 7.300734e-10 3.650367e-10 [383,] 1.0000000 1.798406e-09 8.992030e-10 [384,] 1.0000000 4.706593e-09 2.353296e-09 [385,] 1.0000000 1.268747e-08 6.343734e-09 [386,] 1.0000000 2.417259e-08 1.208630e-08 [387,] 1.0000000 6.157254e-08 3.078627e-08 [388,] 0.9999999 1.559475e-07 7.797375e-08 [389,] 0.9999998 3.770931e-07 1.885465e-07 [390,] 0.9999996 7.592570e-07 3.796285e-07 [391,] 0.9999991 1.851682e-06 9.258412e-07 [392,] 0.9999978 4.472660e-06 2.236330e-06 [393,] 0.9999947 1.056566e-05 5.282829e-06 [394,] 0.9999894 2.118376e-05 1.059188e-05 [395,] 0.9999772 4.558703e-05 2.279352e-05 [396,] 0.9999480 1.040854e-04 5.204271e-05 [397,] 0.9998842 2.315912e-04 1.157956e-04 [398,] 0.9997635 4.730296e-04 2.365148e-04 [399,] 0.9995157 9.686675e-04 4.843337e-04 [400,] 0.9990980 1.804095e-03 9.020473e-04 [401,] 0.9983062 3.387649e-03 1.693824e-03 [402,] 0.9966843 6.631458e-03 3.315729e-03 [403,] 0.9935348 1.293033e-02 6.465164e-03 [404,] 0.9904776 1.904479e-02 9.522396e-03 [405,] 0.9827781 3.444375e-02 1.722188e-02 [406,] 0.9839789 3.204219e-02 1.602110e-02 [407,] 0.9757122 4.857554e-02 2.428777e-02 [408,] 0.9575051 8.498971e-02 4.249486e-02 [409,] 0.9825818 3.483645e-02 1.741822e-02 [410,] 0.9669843 6.603144e-02 3.301572e-02 [411,] 0.9337141 1.325719e-01 6.628593e-02 [412,] 0.8725443 2.549114e-01 1.274557e-01 [413,] 0.7716606 4.566788e-01 2.283394e-01 [414,] 0.6149266 7.701468e-01 3.850734e-01 > postscript(file="/var/www/html/freestat/rcomp/tmp/1ozo91291387065.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index') > points(x[,1]-mysum$resid) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/2h95u1291387065.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/3h95u1291387065.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/4h95u1291387065.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals') > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/5si4x1291387065.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 1424761.6868 2595179.6118 1267617.2969 -2141697.5330 -528231.2923 6 7 8 9 10 25456.4857 17299.3218 111588.9644 -153799.4052 643713.0005 11 12 13 14 15 -201008.0995 281913.3278 126520.7317 -1133466.6604 -59425.3187 16 17 18 19 20 533527.5639 -303941.3313 -46892.0417 -359713.0979 -587011.4116 21 22 23 24 25 225435.3369 -144065.7527 -176182.0570 -751673.8328 -337555.4713 26 27 28 29 30 -97035.4039 -507115.2213 -237311.9801 340409.8527 -52245.6483 31 32 33 34 35 271304.7452 -105029.4459 240738.2350 113337.6690 -86227.9244 36 37 38 39 40 -721470.9502 20062.5278 117533.4460 -110492.7566 -344942.1108 41 42 43 44 45 217921.9000 -350542.1364 118912.6452 -234163.3975 -19947.9048 46 47 48 49 50 -73096.9677 -127351.1621 -131428.4812 -57421.8150 -213711.4726 51 52 53 54 55 11555.7090 -136675.0603 -265261.6904 119697.4669 -403672.9294 56 57 58 59 60 4120.1427 319229.0769 -24944.3045 -23249.7993 -225332.2544 61 62 63 64 65 -5388.1923 -188777.0203 -855.7970 -197758.0913 -11056.7879 66 67 68 69 70 222792.0615 -8890.0679 -31289.6474 -118621.0042 -139973.6658 71 72 73 74 75 -129892.0407 -783.8036 -114438.6411 6545.1960 -59272.4390 76 77 78 79 80 -113270.6585 -30780.3731 -46447.9092 8440.9018 130307.0392 81 82 83 84 85 7792.0332 -72190.1240 -80602.3230 37677.6293 8389.2028 86 87 88 89 90 -127591.6345 -92354.3584 -13609.9017 -120999.5997 12736.5529 91 92 93 94 95 -6237.3608 -43258.9267 75118.1069 -141682.7554 34484.9749 96 97 98 99 100 19382.1451 54863.6374 194966.3872 -58998.6777 -8928.9193 101 102 103 104 105 -52259.7886 -124021.9760 -167970.6595 -53741.3637 27811.1527 106 107 108 109 110 -77790.0753 -260401.3271 120051.7478 -185863.8722 -87237.8808 111 112 113 114 115 89738.5192 -189159.5324 -157439.2056 -67441.9175 83289.6700 116 117 118 119 120 8976.4324 27324.4455 62397.4712 -137142.0463 58597.2943 121 122 123 124 125 -332755.1871 24944.1972 -161159.3664 68868.4749 26449.0078 126 127 128 129 130 -177452.1811 50153.6322 29669.0915 22254.3072 50153.6322 131 132 133 134 135 39961.5367 88974.5340 50153.6322 -1493.2029 28446.6257 136 137 138 139 140 43600.0447 40386.7984 49335.3241 25813.6061 50153.6322 141 142 143 144 145 50704.1245 50153.6322 7751.6224 42775.1877 29170.0055 146 147 148 149 150 37908.1013 50153.6322 38433.2163 51351.0848 51017.7681 151 152 153 154 155 12447.7211 -29464.3939 37297.2609 21128.7762 47600.8600 156 157 158 159 160 50153.6322 50153.6322 43777.2371 57576.6068 50153.6322 161 162 163 164 165 50153.6322 50153.6322 50153.6322 44154.4795 21088.4805 166 167 168 169 170 49595.7466 50153.6322 102839.7560 296524.4666 50153.6322 171 172 173 174 175 50153.6322 41955.9810 50153.6322 25333.4324 50153.6322 176 177 178 179 180 38392.4951 -72254.1418 50542.7018 50153.6322 50153.6322 181 182 183 184 185 45501.4529 18126.4240 6882.4338 24984.8894 38972.2575 186 187 188 189 190 42928.5309 4847.0502 33070.4442 6090.9569 50153.6322 191 192 193 194 195 36344.6628 50153.6322 34586.9467 50153.6322 44203.7263 196 197 198 199 200 50153.6322 -175707.8336 50153.6322 50153.6322 -64242.6081 201 202 203 204 205 50153.6322 50153.6322 50153.6322 49216.0774 -6964.9035 206 207 208 209 210 38609.0592 43050.7884 50153.6322 45402.0907 64749.5854 211 212 213 214 215 50153.6322 50153.6322 73691.2135 9150.7984 19800.6099 216 217 218 219 220 45011.0837 43936.5865 58981.4835 45011.0837 -13575.4133 221 222 223 224 225 22531.3644 50153.6322 49883.3662 10071.1222 49553.3131 226 227 228 229 230 16531.1631 41340.4565 50153.6322 34719.6920 50153.6322 231 232 233 234 235 13724.2079 28758.4829 50153.6322 14949.1276 50153.6322 236 237 238 239 240 51355.1391 50153.6322 41686.3222 50153.6322 50153.6322 241 242 243 244 245 50153.6322 40343.3713 50153.6322 42851.2092 -12355.8354 246 247 248 249 250 50153.6322 50153.6322 13438.7611 50153.6322 50153.6322 251 252 253 254 255 27470.7084 50153.6322 32225.2796 -7134.2484 -135722.2839 256 257 258 259 260 50153.6322 37917.8429 21944.4115 50153.6322 50153.6322 261 262 263 264 265 50153.6322 50153.6322 33779.2527 55197.2354 50153.6322 266 267 268 269 270 50153.6322 50153.6322 72720.5979 50153.6322 -12148.3921 271 272 273 274 275 50153.6322 50153.6322 50153.6322 7400.3458 50153.6322 276 277 278 279 280 50167.9578 21700.6342 50153.6322 49550.8234 41333.8592 281 282 283 284 285 50153.6322 50153.6322 67598.2433 50153.6322 33773.7574 286 287 288 289 290 50153.6322 35142.8678 46949.7047 50153.6322 72624.4917 291 292 293 294 295 50153.6322 50139.3065 31872.2839 50153.6322 28876.8779 296 297 298 299 300 71346.3640 50153.6322 50153.6322 50153.6322 27012.1639 301 302 303 304 305 41576.8749 32725.2768 47582.3580 50153.6322 95788.9728 306 307 308 309 310 29838.6008 50153.6322 7696.9595 28784.5077 50153.6322 311 312 313 314 315 50153.6322 43623.3622 50153.6322 23369.0833 34433.2307 316 317 318 319 320 68515.6647 51471.9224 -8251.2712 44332.8149 50153.6322 321 322 323 324 325 50153.6322 50153.6322 50153.6322 50153.6322 50153.6322 326 327 328 329 330 46577.8273 50153.6322 50153.6322 50153.6322 50153.6322 331 332 333 334 335 50153.6322 50153.6322 50153.6322 42602.7629 49911.1953 336 337 338 339 340 38928.7664 52639.2174 50153.6322 46394.1763 41035.5768 341 342 343 344 345 6711.3425 50153.6322 50153.6322 50153.6322 52629.3460 346 347 348 349 350 50153.6322 37956.3424 50153.6322 50153.6322 13815.2135 351 352 353 354 355 59086.5397 40602.8151 38066.8135 43535.9424 50153.6322 356 357 358 359 360 50153.6322 37757.1676 50153.6322 50153.6322 50153.6322 361 362 363 364 365 44570.4409 50153.6322 50153.6322 -40412.1089 38915.5571 366 367 368 369 370 50153.6322 50153.6322 50153.6322 50153.6322 52317.2864 371 372 373 374 375 50153.6322 79105.5636 31849.7933 24872.2149 50153.6322 376 377 378 379 380 50153.6322 35253.9660 39372.0092 -2616.8774 50153.6322 381 382 383 384 385 50153.6322 43170.7401 93864.5602 50153.6322 35005.1569 386 387 388 389 390 -120729.1955 60161.2553 21354.4899 43789.0774 32210.7648 391 392 393 394 395 31138.4345 13515.9146 -28176.9279 68800.2067 -103849.1978 396 397 398 399 400 4765.1390 -132209.7684 -437035.9785 6759.5568 44199.6986 401 402 403 404 405 -2415.6533 35454.1935 15034.3991 -58691.2392 -17009.7401 406 407 408 409 410 6083.1605 16431.8199 3012.6169 17273.0037 -69144.3131 411 412 413 414 415 -74477.9193 -163524.0366 -33485.7746 42911.0306 -184109.5739 416 417 418 419 420 -105300.1328 133358.5964 -33479.9909 -63100.9831 -40967.7465 421 422 423 424 425 -49443.7944 -217929.1146 -562499.6734 -461505.9634 -55352.9789 426 427 428 429 430 -268237.2948 -311780.3800 -1262344.5769 -270458.3062 -367143.4164 431 -573102.0748 > postscript(file="/var/www/html/freestat/rcomp/tmp/6si4x1291387065.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 1424761.6868 NA 1 2595179.6118 1424761.6868 2 1267617.2969 2595179.6118 3 -2141697.5330 1267617.2969 4 -528231.2923 -2141697.5330 5 25456.4857 -528231.2923 6 17299.3218 25456.4857 7 111588.9644 17299.3218 8 -153799.4052 111588.9644 9 643713.0005 -153799.4052 10 -201008.0995 643713.0005 11 281913.3278 -201008.0995 12 126520.7317 281913.3278 13 -1133466.6604 126520.7317 14 -59425.3187 -1133466.6604 15 533527.5639 -59425.3187 16 -303941.3313 533527.5639 17 -46892.0417 -303941.3313 18 -359713.0979 -46892.0417 19 -587011.4116 -359713.0979 20 225435.3369 -587011.4116 21 -144065.7527 225435.3369 22 -176182.0570 -144065.7527 23 -751673.8328 -176182.0570 24 -337555.4713 -751673.8328 25 -97035.4039 -337555.4713 26 -507115.2213 -97035.4039 27 -237311.9801 -507115.2213 28 340409.8527 -237311.9801 29 -52245.6483 340409.8527 30 271304.7452 -52245.6483 31 -105029.4459 271304.7452 32 240738.2350 -105029.4459 33 113337.6690 240738.2350 34 -86227.9244 113337.6690 35 -721470.9502 -86227.9244 36 20062.5278 -721470.9502 37 117533.4460 20062.5278 38 -110492.7566 117533.4460 39 -344942.1108 -110492.7566 40 217921.9000 -344942.1108 41 -350542.1364 217921.9000 42 118912.6452 -350542.1364 43 -234163.3975 118912.6452 44 -19947.9048 -234163.3975 45 -73096.9677 -19947.9048 46 -127351.1621 -73096.9677 47 -131428.4812 -127351.1621 48 -57421.8150 -131428.4812 49 -213711.4726 -57421.8150 50 11555.7090 -213711.4726 51 -136675.0603 11555.7090 52 -265261.6904 -136675.0603 53 119697.4669 -265261.6904 54 -403672.9294 119697.4669 55 4120.1427 -403672.9294 56 319229.0769 4120.1427 57 -24944.3045 319229.0769 58 -23249.7993 -24944.3045 59 -225332.2544 -23249.7993 60 -5388.1923 -225332.2544 61 -188777.0203 -5388.1923 62 -855.7970 -188777.0203 63 -197758.0913 -855.7970 64 -11056.7879 -197758.0913 65 222792.0615 -11056.7879 66 -8890.0679 222792.0615 67 -31289.6474 -8890.0679 68 -118621.0042 -31289.6474 69 -139973.6658 -118621.0042 70 -129892.0407 -139973.6658 71 -783.8036 -129892.0407 72 -114438.6411 -783.8036 73 6545.1960 -114438.6411 74 -59272.4390 6545.1960 75 -113270.6585 -59272.4390 76 -30780.3731 -113270.6585 77 -46447.9092 -30780.3731 78 8440.9018 -46447.9092 79 130307.0392 8440.9018 80 7792.0332 130307.0392 81 -72190.1240 7792.0332 82 -80602.3230 -72190.1240 83 37677.6293 -80602.3230 84 8389.2028 37677.6293 85 -127591.6345 8389.2028 86 -92354.3584 -127591.6345 87 -13609.9017 -92354.3584 88 -120999.5997 -13609.9017 89 12736.5529 -120999.5997 90 -6237.3608 12736.5529 91 -43258.9267 -6237.3608 92 75118.1069 -43258.9267 93 -141682.7554 75118.1069 94 34484.9749 -141682.7554 95 19382.1451 34484.9749 96 54863.6374 19382.1451 97 194966.3872 54863.6374 98 -58998.6777 194966.3872 99 -8928.9193 -58998.6777 100 -52259.7886 -8928.9193 101 -124021.9760 -52259.7886 102 -167970.6595 -124021.9760 103 -53741.3637 -167970.6595 104 27811.1527 -53741.3637 105 -77790.0753 27811.1527 106 -260401.3271 -77790.0753 107 120051.7478 -260401.3271 108 -185863.8722 120051.7478 109 -87237.8808 -185863.8722 110 89738.5192 -87237.8808 111 -189159.5324 89738.5192 112 -157439.2056 -189159.5324 113 -67441.9175 -157439.2056 114 83289.6700 -67441.9175 115 8976.4324 83289.6700 116 27324.4455 8976.4324 117 62397.4712 27324.4455 118 -137142.0463 62397.4712 119 58597.2943 -137142.0463 120 -332755.1871 58597.2943 121 24944.1972 -332755.1871 122 -161159.3664 24944.1972 123 68868.4749 -161159.3664 124 26449.0078 68868.4749 125 -177452.1811 26449.0078 126 50153.6322 -177452.1811 127 29669.0915 50153.6322 128 22254.3072 29669.0915 129 50153.6322 22254.3072 130 39961.5367 50153.6322 131 88974.5340 39961.5367 132 50153.6322 88974.5340 133 -1493.2029 50153.6322 134 28446.6257 -1493.2029 135 43600.0447 28446.6257 136 40386.7984 43600.0447 137 49335.3241 40386.7984 138 25813.6061 49335.3241 139 50153.6322 25813.6061 140 50704.1245 50153.6322 141 50153.6322 50704.1245 142 7751.6224 50153.6322 143 42775.1877 7751.6224 144 29170.0055 42775.1877 145 37908.1013 29170.0055 146 50153.6322 37908.1013 147 38433.2163 50153.6322 148 51351.0848 38433.2163 149 51017.7681 51351.0848 150 12447.7211 51017.7681 151 -29464.3939 12447.7211 152 37297.2609 -29464.3939 153 21128.7762 37297.2609 154 47600.8600 21128.7762 155 50153.6322 47600.8600 156 50153.6322 50153.6322 157 43777.2371 50153.6322 158 57576.6068 43777.2371 159 50153.6322 57576.6068 160 50153.6322 50153.6322 161 50153.6322 50153.6322 162 50153.6322 50153.6322 163 44154.4795 50153.6322 164 21088.4805 44154.4795 165 49595.7466 21088.4805 166 50153.6322 49595.7466 167 102839.7560 50153.6322 168 296524.4666 102839.7560 169 50153.6322 296524.4666 170 50153.6322 50153.6322 171 41955.9810 50153.6322 172 50153.6322 41955.9810 173 25333.4324 50153.6322 174 50153.6322 25333.4324 175 38392.4951 50153.6322 176 -72254.1418 38392.4951 177 50542.7018 -72254.1418 178 50153.6322 50542.7018 179 50153.6322 50153.6322 180 45501.4529 50153.6322 181 18126.4240 45501.4529 182 6882.4338 18126.4240 183 24984.8894 6882.4338 184 38972.2575 24984.8894 185 42928.5309 38972.2575 186 4847.0502 42928.5309 187 33070.4442 4847.0502 188 6090.9569 33070.4442 189 50153.6322 6090.9569 190 36344.6628 50153.6322 191 50153.6322 36344.6628 192 34586.9467 50153.6322 193 50153.6322 34586.9467 194 44203.7263 50153.6322 195 50153.6322 44203.7263 196 -175707.8336 50153.6322 197 50153.6322 -175707.8336 198 50153.6322 50153.6322 199 -64242.6081 50153.6322 200 50153.6322 -64242.6081 201 50153.6322 50153.6322 202 50153.6322 50153.6322 203 49216.0774 50153.6322 204 -6964.9035 49216.0774 205 38609.0592 -6964.9035 206 43050.7884 38609.0592 207 50153.6322 43050.7884 208 45402.0907 50153.6322 209 64749.5854 45402.0907 210 50153.6322 64749.5854 211 50153.6322 50153.6322 212 73691.2135 50153.6322 213 9150.7984 73691.2135 214 19800.6099 9150.7984 215 45011.0837 19800.6099 216 43936.5865 45011.0837 217 58981.4835 43936.5865 218 45011.0837 58981.4835 219 -13575.4133 45011.0837 220 22531.3644 -13575.4133 221 50153.6322 22531.3644 222 49883.3662 50153.6322 223 10071.1222 49883.3662 224 49553.3131 10071.1222 225 16531.1631 49553.3131 226 41340.4565 16531.1631 227 50153.6322 41340.4565 228 34719.6920 50153.6322 229 50153.6322 34719.6920 230 13724.2079 50153.6322 231 28758.4829 13724.2079 232 50153.6322 28758.4829 233 14949.1276 50153.6322 234 50153.6322 14949.1276 235 51355.1391 50153.6322 236 50153.6322 51355.1391 237 41686.3222 50153.6322 238 50153.6322 41686.3222 239 50153.6322 50153.6322 240 50153.6322 50153.6322 241 40343.3713 50153.6322 242 50153.6322 40343.3713 243 42851.2092 50153.6322 244 -12355.8354 42851.2092 245 50153.6322 -12355.8354 246 50153.6322 50153.6322 247 13438.7611 50153.6322 248 50153.6322 13438.7611 249 50153.6322 50153.6322 250 27470.7084 50153.6322 251 50153.6322 27470.7084 252 32225.2796 50153.6322 253 -7134.2484 32225.2796 254 -135722.2839 -7134.2484 255 50153.6322 -135722.2839 256 37917.8429 50153.6322 257 21944.4115 37917.8429 258 50153.6322 21944.4115 259 50153.6322 50153.6322 260 50153.6322 50153.6322 261 50153.6322 50153.6322 262 33779.2527 50153.6322 263 55197.2354 33779.2527 264 50153.6322 55197.2354 265 50153.6322 50153.6322 266 50153.6322 50153.6322 267 72720.5979 50153.6322 268 50153.6322 72720.5979 269 -12148.3921 50153.6322 270 50153.6322 -12148.3921 271 50153.6322 50153.6322 272 50153.6322 50153.6322 273 7400.3458 50153.6322 274 50153.6322 7400.3458 275 50167.9578 50153.6322 276 21700.6342 50167.9578 277 50153.6322 21700.6342 278 49550.8234 50153.6322 279 41333.8592 49550.8234 280 50153.6322 41333.8592 281 50153.6322 50153.6322 282 67598.2433 50153.6322 283 50153.6322 67598.2433 284 33773.7574 50153.6322 285 50153.6322 33773.7574 286 35142.8678 50153.6322 287 46949.7047 35142.8678 288 50153.6322 46949.7047 289 72624.4917 50153.6322 290 50153.6322 72624.4917 291 50139.3065 50153.6322 292 31872.2839 50139.3065 293 50153.6322 31872.2839 294 28876.8779 50153.6322 295 71346.3640 28876.8779 296 50153.6322 71346.3640 297 50153.6322 50153.6322 298 50153.6322 50153.6322 299 27012.1639 50153.6322 300 41576.8749 27012.1639 301 32725.2768 41576.8749 302 47582.3580 32725.2768 303 50153.6322 47582.3580 304 95788.9728 50153.6322 305 29838.6008 95788.9728 306 50153.6322 29838.6008 307 7696.9595 50153.6322 308 28784.5077 7696.9595 309 50153.6322 28784.5077 310 50153.6322 50153.6322 311 43623.3622 50153.6322 312 50153.6322 43623.3622 313 23369.0833 50153.6322 314 34433.2307 23369.0833 315 68515.6647 34433.2307 316 51471.9224 68515.6647 317 -8251.2712 51471.9224 318 44332.8149 -8251.2712 319 50153.6322 44332.8149 320 50153.6322 50153.6322 321 50153.6322 50153.6322 322 50153.6322 50153.6322 323 50153.6322 50153.6322 324 50153.6322 50153.6322 325 46577.8273 50153.6322 326 50153.6322 46577.8273 327 50153.6322 50153.6322 328 50153.6322 50153.6322 329 50153.6322 50153.6322 330 50153.6322 50153.6322 331 50153.6322 50153.6322 332 50153.6322 50153.6322 333 42602.7629 50153.6322 334 49911.1953 42602.7629 335 38928.7664 49911.1953 336 52639.2174 38928.7664 337 50153.6322 52639.2174 338 46394.1763 50153.6322 339 41035.5768 46394.1763 340 6711.3425 41035.5768 341 50153.6322 6711.3425 342 50153.6322 50153.6322 343 50153.6322 50153.6322 344 52629.3460 50153.6322 345 50153.6322 52629.3460 346 37956.3424 50153.6322 347 50153.6322 37956.3424 348 50153.6322 50153.6322 349 13815.2135 50153.6322 350 59086.5397 13815.2135 351 40602.8151 59086.5397 352 38066.8135 40602.8151 353 43535.9424 38066.8135 354 50153.6322 43535.9424 355 50153.6322 50153.6322 356 37757.1676 50153.6322 357 50153.6322 37757.1676 358 50153.6322 50153.6322 359 50153.6322 50153.6322 360 44570.4409 50153.6322 361 50153.6322 44570.4409 362 50153.6322 50153.6322 363 -40412.1089 50153.6322 364 38915.5571 -40412.1089 365 50153.6322 38915.5571 366 50153.6322 50153.6322 367 50153.6322 50153.6322 368 50153.6322 50153.6322 369 52317.2864 50153.6322 370 50153.6322 52317.2864 371 79105.5636 50153.6322 372 31849.7933 79105.5636 373 24872.2149 31849.7933 374 50153.6322 24872.2149 375 50153.6322 50153.6322 376 35253.9660 50153.6322 377 39372.0092 35253.9660 378 -2616.8774 39372.0092 379 50153.6322 -2616.8774 380 50153.6322 50153.6322 381 43170.7401 50153.6322 382 93864.5602 43170.7401 383 50153.6322 93864.5602 384 35005.1569 50153.6322 385 -120729.1955 35005.1569 386 60161.2553 -120729.1955 387 21354.4899 60161.2553 388 43789.0774 21354.4899 389 32210.7648 43789.0774 390 31138.4345 32210.7648 391 13515.9146 31138.4345 392 -28176.9279 13515.9146 393 68800.2067 -28176.9279 394 -103849.1978 68800.2067 395 4765.1390 -103849.1978 396 -132209.7684 4765.1390 397 -437035.9785 -132209.7684 398 6759.5568 -437035.9785 399 44199.6986 6759.5568 400 -2415.6533 44199.6986 401 35454.1935 -2415.6533 402 15034.3991 35454.1935 403 -58691.2392 15034.3991 404 -17009.7401 -58691.2392 405 6083.1605 -17009.7401 406 16431.8199 6083.1605 407 3012.6169 16431.8199 408 17273.0037 3012.6169 409 -69144.3131 17273.0037 410 -74477.9193 -69144.3131 411 -163524.0366 -74477.9193 412 -33485.7746 -163524.0366 413 42911.0306 -33485.7746 414 -184109.5739 42911.0306 415 -105300.1328 -184109.5739 416 133358.5964 -105300.1328 417 -33479.9909 133358.5964 418 -63100.9831 -33479.9909 419 -40967.7465 -63100.9831 420 -49443.7944 -40967.7465 421 -217929.1146 -49443.7944 422 -562499.6734 -217929.1146 423 -461505.9634 -562499.6734 424 -55352.9789 -461505.9634 425 -268237.2948 -55352.9789 426 -311780.3800 -268237.2948 427 -1262344.5769 -311780.3800 428 -270458.3062 -1262344.5769 429 -367143.4164 -270458.3062 430 -573102.0748 -367143.4164 431 NA -573102.0748 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 2595179.6118 1424761.6868 [2,] 1267617.2969 2595179.6118 [3,] -2141697.5330 1267617.2969 [4,] -528231.2923 -2141697.5330 [5,] 25456.4857 -528231.2923 [6,] 17299.3218 25456.4857 [7,] 111588.9644 17299.3218 [8,] -153799.4052 111588.9644 [9,] 643713.0005 -153799.4052 [10,] -201008.0995 643713.0005 [11,] 281913.3278 -201008.0995 [12,] 126520.7317 281913.3278 [13,] -1133466.6604 126520.7317 [14,] -59425.3187 -1133466.6604 [15,] 533527.5639 -59425.3187 [16,] -303941.3313 533527.5639 [17,] -46892.0417 -303941.3313 [18,] -359713.0979 -46892.0417 [19,] -587011.4116 -359713.0979 [20,] 225435.3369 -587011.4116 [21,] -144065.7527 225435.3369 [22,] -176182.0570 -144065.7527 [23,] -751673.8328 -176182.0570 [24,] -337555.4713 -751673.8328 [25,] -97035.4039 -337555.4713 [26,] -507115.2213 -97035.4039 [27,] -237311.9801 -507115.2213 [28,] 340409.8527 -237311.9801 [29,] -52245.6483 340409.8527 [30,] 271304.7452 -52245.6483 [31,] -105029.4459 271304.7452 [32,] 240738.2350 -105029.4459 [33,] 113337.6690 240738.2350 [34,] -86227.9244 113337.6690 [35,] -721470.9502 -86227.9244 [36,] 20062.5278 -721470.9502 [37,] 117533.4460 20062.5278 [38,] -110492.7566 117533.4460 [39,] -344942.1108 -110492.7566 [40,] 217921.9000 -344942.1108 [41,] -350542.1364 217921.9000 [42,] 118912.6452 -350542.1364 [43,] -234163.3975 118912.6452 [44,] -19947.9048 -234163.3975 [45,] -73096.9677 -19947.9048 [46,] -127351.1621 -73096.9677 [47,] -131428.4812 -127351.1621 [48,] -57421.8150 -131428.4812 [49,] -213711.4726 -57421.8150 [50,] 11555.7090 -213711.4726 [51,] -136675.0603 11555.7090 [52,] -265261.6904 -136675.0603 [53,] 119697.4669 -265261.6904 [54,] -403672.9294 119697.4669 [55,] 4120.1427 -403672.9294 [56,] 319229.0769 4120.1427 [57,] -24944.3045 319229.0769 [58,] -23249.7993 -24944.3045 [59,] -225332.2544 -23249.7993 [60,] -5388.1923 -225332.2544 [61,] -188777.0203 -5388.1923 [62,] -855.7970 -188777.0203 [63,] -197758.0913 -855.7970 [64,] -11056.7879 -197758.0913 [65,] 222792.0615 -11056.7879 [66,] -8890.0679 222792.0615 [67,] -31289.6474 -8890.0679 [68,] -118621.0042 -31289.6474 [69,] -139973.6658 -118621.0042 [70,] -129892.0407 -139973.6658 [71,] -783.8036 -129892.0407 [72,] -114438.6411 -783.8036 [73,] 6545.1960 -114438.6411 [74,] -59272.4390 6545.1960 [75,] -113270.6585 -59272.4390 [76,] -30780.3731 -113270.6585 [77,] -46447.9092 -30780.3731 [78,] 8440.9018 -46447.9092 [79,] 130307.0392 8440.9018 [80,] 7792.0332 130307.0392 [81,] -72190.1240 7792.0332 [82,] -80602.3230 -72190.1240 [83,] 37677.6293 -80602.3230 [84,] 8389.2028 37677.6293 [85,] -127591.6345 8389.2028 [86,] -92354.3584 -127591.6345 [87,] -13609.9017 -92354.3584 [88,] -120999.5997 -13609.9017 [89,] 12736.5529 -120999.5997 [90,] -6237.3608 12736.5529 [91,] -43258.9267 -6237.3608 [92,] 75118.1069 -43258.9267 [93,] -141682.7554 75118.1069 [94,] 34484.9749 -141682.7554 [95,] 19382.1451 34484.9749 [96,] 54863.6374 19382.1451 [97,] 194966.3872 54863.6374 [98,] -58998.6777 194966.3872 [99,] -8928.9193 -58998.6777 [100,] -52259.7886 -8928.9193 [101,] -124021.9760 -52259.7886 [102,] -167970.6595 -124021.9760 [103,] -53741.3637 -167970.6595 [104,] 27811.1527 -53741.3637 [105,] -77790.0753 27811.1527 [106,] -260401.3271 -77790.0753 [107,] 120051.7478 -260401.3271 [108,] -185863.8722 120051.7478 [109,] -87237.8808 -185863.8722 [110,] 89738.5192 -87237.8808 [111,] -189159.5324 89738.5192 [112,] -157439.2056 -189159.5324 [113,] -67441.9175 -157439.2056 [114,] 83289.6700 -67441.9175 [115,] 8976.4324 83289.6700 [116,] 27324.4455 8976.4324 [117,] 62397.4712 27324.4455 [118,] -137142.0463 62397.4712 [119,] 58597.2943 -137142.0463 [120,] -332755.1871 58597.2943 [121,] 24944.1972 -332755.1871 [122,] -161159.3664 24944.1972 [123,] 68868.4749 -161159.3664 [124,] 26449.0078 68868.4749 [125,] -177452.1811 26449.0078 [126,] 50153.6322 -177452.1811 [127,] 29669.0915 50153.6322 [128,] 22254.3072 29669.0915 [129,] 50153.6322 22254.3072 [130,] 39961.5367 50153.6322 [131,] 88974.5340 39961.5367 [132,] 50153.6322 88974.5340 [133,] -1493.2029 50153.6322 [134,] 28446.6257 -1493.2029 [135,] 43600.0447 28446.6257 [136,] 40386.7984 43600.0447 [137,] 49335.3241 40386.7984 [138,] 25813.6061 49335.3241 [139,] 50153.6322 25813.6061 [140,] 50704.1245 50153.6322 [141,] 50153.6322 50704.1245 [142,] 7751.6224 50153.6322 [143,] 42775.1877 7751.6224 [144,] 29170.0055 42775.1877 [145,] 37908.1013 29170.0055 [146,] 50153.6322 37908.1013 [147,] 38433.2163 50153.6322 [148,] 51351.0848 38433.2163 [149,] 51017.7681 51351.0848 [150,] 12447.7211 51017.7681 [151,] -29464.3939 12447.7211 [152,] 37297.2609 -29464.3939 [153,] 21128.7762 37297.2609 [154,] 47600.8600 21128.7762 [155,] 50153.6322 47600.8600 [156,] 50153.6322 50153.6322 [157,] 43777.2371 50153.6322 [158,] 57576.6068 43777.2371 [159,] 50153.6322 57576.6068 [160,] 50153.6322 50153.6322 [161,] 50153.6322 50153.6322 [162,] 50153.6322 50153.6322 [163,] 44154.4795 50153.6322 [164,] 21088.4805 44154.4795 [165,] 49595.7466 21088.4805 [166,] 50153.6322 49595.7466 [167,] 102839.7560 50153.6322 [168,] 296524.4666 102839.7560 [169,] 50153.6322 296524.4666 [170,] 50153.6322 50153.6322 [171,] 41955.9810 50153.6322 [172,] 50153.6322 41955.9810 [173,] 25333.4324 50153.6322 [174,] 50153.6322 25333.4324 [175,] 38392.4951 50153.6322 [176,] -72254.1418 38392.4951 [177,] 50542.7018 -72254.1418 [178,] 50153.6322 50542.7018 [179,] 50153.6322 50153.6322 [180,] 45501.4529 50153.6322 [181,] 18126.4240 45501.4529 [182,] 6882.4338 18126.4240 [183,] 24984.8894 6882.4338 [184,] 38972.2575 24984.8894 [185,] 42928.5309 38972.2575 [186,] 4847.0502 42928.5309 [187,] 33070.4442 4847.0502 [188,] 6090.9569 33070.4442 [189,] 50153.6322 6090.9569 [190,] 36344.6628 50153.6322 [191,] 50153.6322 36344.6628 [192,] 34586.9467 50153.6322 [193,] 50153.6322 34586.9467 [194,] 44203.7263 50153.6322 [195,] 50153.6322 44203.7263 [196,] -175707.8336 50153.6322 [197,] 50153.6322 -175707.8336 [198,] 50153.6322 50153.6322 [199,] -64242.6081 50153.6322 [200,] 50153.6322 -64242.6081 [201,] 50153.6322 50153.6322 [202,] 50153.6322 50153.6322 [203,] 49216.0774 50153.6322 [204,] -6964.9035 49216.0774 [205,] 38609.0592 -6964.9035 [206,] 43050.7884 38609.0592 [207,] 50153.6322 43050.7884 [208,] 45402.0907 50153.6322 [209,] 64749.5854 45402.0907 [210,] 50153.6322 64749.5854 [211,] 50153.6322 50153.6322 [212,] 73691.2135 50153.6322 [213,] 9150.7984 73691.2135 [214,] 19800.6099 9150.7984 [215,] 45011.0837 19800.6099 [216,] 43936.5865 45011.0837 [217,] 58981.4835 43936.5865 [218,] 45011.0837 58981.4835 [219,] -13575.4133 45011.0837 [220,] 22531.3644 -13575.4133 [221,] 50153.6322 22531.3644 [222,] 49883.3662 50153.6322 [223,] 10071.1222 49883.3662 [224,] 49553.3131 10071.1222 [225,] 16531.1631 49553.3131 [226,] 41340.4565 16531.1631 [227,] 50153.6322 41340.4565 [228,] 34719.6920 50153.6322 [229,] 50153.6322 34719.6920 [230,] 13724.2079 50153.6322 [231,] 28758.4829 13724.2079 [232,] 50153.6322 28758.4829 [233,] 14949.1276 50153.6322 [234,] 50153.6322 14949.1276 [235,] 51355.1391 50153.6322 [236,] 50153.6322 51355.1391 [237,] 41686.3222 50153.6322 [238,] 50153.6322 41686.3222 [239,] 50153.6322 50153.6322 [240,] 50153.6322 50153.6322 [241,] 40343.3713 50153.6322 [242,] 50153.6322 40343.3713 [243,] 42851.2092 50153.6322 [244,] -12355.8354 42851.2092 [245,] 50153.6322 -12355.8354 [246,] 50153.6322 50153.6322 [247,] 13438.7611 50153.6322 [248,] 50153.6322 13438.7611 [249,] 50153.6322 50153.6322 [250,] 27470.7084 50153.6322 [251,] 50153.6322 27470.7084 [252,] 32225.2796 50153.6322 [253,] -7134.2484 32225.2796 [254,] -135722.2839 -7134.2484 [255,] 50153.6322 -135722.2839 [256,] 37917.8429 50153.6322 [257,] 21944.4115 37917.8429 [258,] 50153.6322 21944.4115 [259,] 50153.6322 50153.6322 [260,] 50153.6322 50153.6322 [261,] 50153.6322 50153.6322 [262,] 33779.2527 50153.6322 [263,] 55197.2354 33779.2527 [264,] 50153.6322 55197.2354 [265,] 50153.6322 50153.6322 [266,] 50153.6322 50153.6322 [267,] 72720.5979 50153.6322 [268,] 50153.6322 72720.5979 [269,] -12148.3921 50153.6322 [270,] 50153.6322 -12148.3921 [271,] 50153.6322 50153.6322 [272,] 50153.6322 50153.6322 [273,] 7400.3458 50153.6322 [274,] 50153.6322 7400.3458 [275,] 50167.9578 50153.6322 [276,] 21700.6342 50167.9578 [277,] 50153.6322 21700.6342 [278,] 49550.8234 50153.6322 [279,] 41333.8592 49550.8234 [280,] 50153.6322 41333.8592 [281,] 50153.6322 50153.6322 [282,] 67598.2433 50153.6322 [283,] 50153.6322 67598.2433 [284,] 33773.7574 50153.6322 [285,] 50153.6322 33773.7574 [286,] 35142.8678 50153.6322 [287,] 46949.7047 35142.8678 [288,] 50153.6322 46949.7047 [289,] 72624.4917 50153.6322 [290,] 50153.6322 72624.4917 [291,] 50139.3065 50153.6322 [292,] 31872.2839 50139.3065 [293,] 50153.6322 31872.2839 [294,] 28876.8779 50153.6322 [295,] 71346.3640 28876.8779 [296,] 50153.6322 71346.3640 [297,] 50153.6322 50153.6322 [298,] 50153.6322 50153.6322 [299,] 27012.1639 50153.6322 [300,] 41576.8749 27012.1639 [301,] 32725.2768 41576.8749 [302,] 47582.3580 32725.2768 [303,] 50153.6322 47582.3580 [304,] 95788.9728 50153.6322 [305,] 29838.6008 95788.9728 [306,] 50153.6322 29838.6008 [307,] 7696.9595 50153.6322 [308,] 28784.5077 7696.9595 [309,] 50153.6322 28784.5077 [310,] 50153.6322 50153.6322 [311,] 43623.3622 50153.6322 [312,] 50153.6322 43623.3622 [313,] 23369.0833 50153.6322 [314,] 34433.2307 23369.0833 [315,] 68515.6647 34433.2307 [316,] 51471.9224 68515.6647 [317,] -8251.2712 51471.9224 [318,] 44332.8149 -8251.2712 [319,] 50153.6322 44332.8149 [320,] 50153.6322 50153.6322 [321,] 50153.6322 50153.6322 [322,] 50153.6322 50153.6322 [323,] 50153.6322 50153.6322 [324,] 50153.6322 50153.6322 [325,] 46577.8273 50153.6322 [326,] 50153.6322 46577.8273 [327,] 50153.6322 50153.6322 [328,] 50153.6322 50153.6322 [329,] 50153.6322 50153.6322 [330,] 50153.6322 50153.6322 [331,] 50153.6322 50153.6322 [332,] 50153.6322 50153.6322 [333,] 42602.7629 50153.6322 [334,] 49911.1953 42602.7629 [335,] 38928.7664 49911.1953 [336,] 52639.2174 38928.7664 [337,] 50153.6322 52639.2174 [338,] 46394.1763 50153.6322 [339,] 41035.5768 46394.1763 [340,] 6711.3425 41035.5768 [341,] 50153.6322 6711.3425 [342,] 50153.6322 50153.6322 [343,] 50153.6322 50153.6322 [344,] 52629.3460 50153.6322 [345,] 50153.6322 52629.3460 [346,] 37956.3424 50153.6322 [347,] 50153.6322 37956.3424 [348,] 50153.6322 50153.6322 [349,] 13815.2135 50153.6322 [350,] 59086.5397 13815.2135 [351,] 40602.8151 59086.5397 [352,] 38066.8135 40602.8151 [353,] 43535.9424 38066.8135 [354,] 50153.6322 43535.9424 [355,] 50153.6322 50153.6322 [356,] 37757.1676 50153.6322 [357,] 50153.6322 37757.1676 [358,] 50153.6322 50153.6322 [359,] 50153.6322 50153.6322 [360,] 44570.4409 50153.6322 [361,] 50153.6322 44570.4409 [362,] 50153.6322 50153.6322 [363,] -40412.1089 50153.6322 [364,] 38915.5571 -40412.1089 [365,] 50153.6322 38915.5571 [366,] 50153.6322 50153.6322 [367,] 50153.6322 50153.6322 [368,] 50153.6322 50153.6322 [369,] 52317.2864 50153.6322 [370,] 50153.6322 52317.2864 [371,] 79105.5636 50153.6322 [372,] 31849.7933 79105.5636 [373,] 24872.2149 31849.7933 [374,] 50153.6322 24872.2149 [375,] 50153.6322 50153.6322 [376,] 35253.9660 50153.6322 [377,] 39372.0092 35253.9660 [378,] -2616.8774 39372.0092 [379,] 50153.6322 -2616.8774 [380,] 50153.6322 50153.6322 [381,] 43170.7401 50153.6322 [382,] 93864.5602 43170.7401 [383,] 50153.6322 93864.5602 [384,] 35005.1569 50153.6322 [385,] -120729.1955 35005.1569 [386,] 60161.2553 -120729.1955 [387,] 21354.4899 60161.2553 [388,] 43789.0774 21354.4899 [389,] 32210.7648 43789.0774 [390,] 31138.4345 32210.7648 [391,] 13515.9146 31138.4345 [392,] -28176.9279 13515.9146 [393,] 68800.2067 -28176.9279 [394,] -103849.1978 68800.2067 [395,] 4765.1390 -103849.1978 [396,] -132209.7684 4765.1390 [397,] -437035.9785 -132209.7684 [398,] 6759.5568 -437035.9785 [399,] 44199.6986 6759.5568 [400,] -2415.6533 44199.6986 [401,] 35454.1935 -2415.6533 [402,] 15034.3991 35454.1935 [403,] -58691.2392 15034.3991 [404,] -17009.7401 -58691.2392 [405,] 6083.1605 -17009.7401 [406,] 16431.8199 6083.1605 [407,] 3012.6169 16431.8199 [408,] 17273.0037 3012.6169 [409,] -69144.3131 17273.0037 [410,] -74477.9193 -69144.3131 [411,] -163524.0366 -74477.9193 [412,] -33485.7746 -163524.0366 [413,] 42911.0306 -33485.7746 [414,] -184109.5739 42911.0306 [415,] -105300.1328 -184109.5739 [416,] 133358.5964 -105300.1328 [417,] -33479.9909 133358.5964 [418,] -63100.9831 -33479.9909 [419,] -40967.7465 -63100.9831 [420,] -49443.7944 -40967.7465 [421,] -217929.1146 -49443.7944 [422,] -562499.6734 -217929.1146 [423,] -461505.9634 -562499.6734 [424,] -55352.9789 -461505.9634 [425,] -268237.2948 -55352.9789 [426,] -311780.3800 -268237.2948 [427,] -1262344.5769 -311780.3800 [428,] -270458.3062 -1262344.5769 [429,] -367143.4164 -270458.3062 [430,] -573102.0748 -367143.4164 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 2595179.6118 1424761.6868 2 1267617.2969 2595179.6118 3 -2141697.5330 1267617.2969 4 -528231.2923 -2141697.5330 5 25456.4857 -528231.2923 6 17299.3218 25456.4857 7 111588.9644 17299.3218 8 -153799.4052 111588.9644 9 643713.0005 -153799.4052 10 -201008.0995 643713.0005 11 281913.3278 -201008.0995 12 126520.7317 281913.3278 13 -1133466.6604 126520.7317 14 -59425.3187 -1133466.6604 15 533527.5639 -59425.3187 16 -303941.3313 533527.5639 17 -46892.0417 -303941.3313 18 -359713.0979 -46892.0417 19 -587011.4116 -359713.0979 20 225435.3369 -587011.4116 21 -144065.7527 225435.3369 22 -176182.0570 -144065.7527 23 -751673.8328 -176182.0570 24 -337555.4713 -751673.8328 25 -97035.4039 -337555.4713 26 -507115.2213 -97035.4039 27 -237311.9801 -507115.2213 28 340409.8527 -237311.9801 29 -52245.6483 340409.8527 30 271304.7452 -52245.6483 31 -105029.4459 271304.7452 32 240738.2350 -105029.4459 33 113337.6690 240738.2350 34 -86227.9244 113337.6690 35 -721470.9502 -86227.9244 36 20062.5278 -721470.9502 37 117533.4460 20062.5278 38 -110492.7566 117533.4460 39 -344942.1108 -110492.7566 40 217921.9000 -344942.1108 41 -350542.1364 217921.9000 42 118912.6452 -350542.1364 43 -234163.3975 118912.6452 44 -19947.9048 -234163.3975 45 -73096.9677 -19947.9048 46 -127351.1621 -73096.9677 47 -131428.4812 -127351.1621 48 -57421.8150 -131428.4812 49 -213711.4726 -57421.8150 50 11555.7090 -213711.4726 51 -136675.0603 11555.7090 52 -265261.6904 -136675.0603 53 119697.4669 -265261.6904 54 -403672.9294 119697.4669 55 4120.1427 -403672.9294 56 319229.0769 4120.1427 57 -24944.3045 319229.0769 58 -23249.7993 -24944.3045 59 -225332.2544 -23249.7993 60 -5388.1923 -225332.2544 61 -188777.0203 -5388.1923 62 -855.7970 -188777.0203 63 -197758.0913 -855.7970 64 -11056.7879 -197758.0913 65 222792.0615 -11056.7879 66 -8890.0679 222792.0615 67 -31289.6474 -8890.0679 68 -118621.0042 -31289.6474 69 -139973.6658 -118621.0042 70 -129892.0407 -139973.6658 71 -783.8036 -129892.0407 72 -114438.6411 -783.8036 73 6545.1960 -114438.6411 74 -59272.4390 6545.1960 75 -113270.6585 -59272.4390 76 -30780.3731 -113270.6585 77 -46447.9092 -30780.3731 78 8440.9018 -46447.9092 79 130307.0392 8440.9018 80 7792.0332 130307.0392 81 -72190.1240 7792.0332 82 -80602.3230 -72190.1240 83 37677.6293 -80602.3230 84 8389.2028 37677.6293 85 -127591.6345 8389.2028 86 -92354.3584 -127591.6345 87 -13609.9017 -92354.3584 88 -120999.5997 -13609.9017 89 12736.5529 -120999.5997 90 -6237.3608 12736.5529 91 -43258.9267 -6237.3608 92 75118.1069 -43258.9267 93 -141682.7554 75118.1069 94 34484.9749 -141682.7554 95 19382.1451 34484.9749 96 54863.6374 19382.1451 97 194966.3872 54863.6374 98 -58998.6777 194966.3872 99 -8928.9193 -58998.6777 100 -52259.7886 -8928.9193 101 -124021.9760 -52259.7886 102 -167970.6595 -124021.9760 103 -53741.3637 -167970.6595 104 27811.1527 -53741.3637 105 -77790.0753 27811.1527 106 -260401.3271 -77790.0753 107 120051.7478 -260401.3271 108 -185863.8722 120051.7478 109 -87237.8808 -185863.8722 110 89738.5192 -87237.8808 111 -189159.5324 89738.5192 112 -157439.2056 -189159.5324 113 -67441.9175 -157439.2056 114 83289.6700 -67441.9175 115 8976.4324 83289.6700 116 27324.4455 8976.4324 117 62397.4712 27324.4455 118 -137142.0463 62397.4712 119 58597.2943 -137142.0463 120 -332755.1871 58597.2943 121 24944.1972 -332755.1871 122 -161159.3664 24944.1972 123 68868.4749 -161159.3664 124 26449.0078 68868.4749 125 -177452.1811 26449.0078 126 50153.6322 -177452.1811 127 29669.0915 50153.6322 128 22254.3072 29669.0915 129 50153.6322 22254.3072 130 39961.5367 50153.6322 131 88974.5340 39961.5367 132 50153.6322 88974.5340 133 -1493.2029 50153.6322 134 28446.6257 -1493.2029 135 43600.0447 28446.6257 136 40386.7984 43600.0447 137 49335.3241 40386.7984 138 25813.6061 49335.3241 139 50153.6322 25813.6061 140 50704.1245 50153.6322 141 50153.6322 50704.1245 142 7751.6224 50153.6322 143 42775.1877 7751.6224 144 29170.0055 42775.1877 145 37908.1013 29170.0055 146 50153.6322 37908.1013 147 38433.2163 50153.6322 148 51351.0848 38433.2163 149 51017.7681 51351.0848 150 12447.7211 51017.7681 151 -29464.3939 12447.7211 152 37297.2609 -29464.3939 153 21128.7762 37297.2609 154 47600.8600 21128.7762 155 50153.6322 47600.8600 156 50153.6322 50153.6322 157 43777.2371 50153.6322 158 57576.6068 43777.2371 159 50153.6322 57576.6068 160 50153.6322 50153.6322 161 50153.6322 50153.6322 162 50153.6322 50153.6322 163 44154.4795 50153.6322 164 21088.4805 44154.4795 165 49595.7466 21088.4805 166 50153.6322 49595.7466 167 102839.7560 50153.6322 168 296524.4666 102839.7560 169 50153.6322 296524.4666 170 50153.6322 50153.6322 171 41955.9810 50153.6322 172 50153.6322 41955.9810 173 25333.4324 50153.6322 174 50153.6322 25333.4324 175 38392.4951 50153.6322 176 -72254.1418 38392.4951 177 50542.7018 -72254.1418 178 50153.6322 50542.7018 179 50153.6322 50153.6322 180 45501.4529 50153.6322 181 18126.4240 45501.4529 182 6882.4338 18126.4240 183 24984.8894 6882.4338 184 38972.2575 24984.8894 185 42928.5309 38972.2575 186 4847.0502 42928.5309 187 33070.4442 4847.0502 188 6090.9569 33070.4442 189 50153.6322 6090.9569 190 36344.6628 50153.6322 191 50153.6322 36344.6628 192 34586.9467 50153.6322 193 50153.6322 34586.9467 194 44203.7263 50153.6322 195 50153.6322 44203.7263 196 -175707.8336 50153.6322 197 50153.6322 -175707.8336 198 50153.6322 50153.6322 199 -64242.6081 50153.6322 200 50153.6322 -64242.6081 201 50153.6322 50153.6322 202 50153.6322 50153.6322 203 49216.0774 50153.6322 204 -6964.9035 49216.0774 205 38609.0592 -6964.9035 206 43050.7884 38609.0592 207 50153.6322 43050.7884 208 45402.0907 50153.6322 209 64749.5854 45402.0907 210 50153.6322 64749.5854 211 50153.6322 50153.6322 212 73691.2135 50153.6322 213 9150.7984 73691.2135 214 19800.6099 9150.7984 215 45011.0837 19800.6099 216 43936.5865 45011.0837 217 58981.4835 43936.5865 218 45011.0837 58981.4835 219 -13575.4133 45011.0837 220 22531.3644 -13575.4133 221 50153.6322 22531.3644 222 49883.3662 50153.6322 223 10071.1222 49883.3662 224 49553.3131 10071.1222 225 16531.1631 49553.3131 226 41340.4565 16531.1631 227 50153.6322 41340.4565 228 34719.6920 50153.6322 229 50153.6322 34719.6920 230 13724.2079 50153.6322 231 28758.4829 13724.2079 232 50153.6322 28758.4829 233 14949.1276 50153.6322 234 50153.6322 14949.1276 235 51355.1391 50153.6322 236 50153.6322 51355.1391 237 41686.3222 50153.6322 238 50153.6322 41686.3222 239 50153.6322 50153.6322 240 50153.6322 50153.6322 241 40343.3713 50153.6322 242 50153.6322 40343.3713 243 42851.2092 50153.6322 244 -12355.8354 42851.2092 245 50153.6322 -12355.8354 246 50153.6322 50153.6322 247 13438.7611 50153.6322 248 50153.6322 13438.7611 249 50153.6322 50153.6322 250 27470.7084 50153.6322 251 50153.6322 27470.7084 252 32225.2796 50153.6322 253 -7134.2484 32225.2796 254 -135722.2839 -7134.2484 255 50153.6322 -135722.2839 256 37917.8429 50153.6322 257 21944.4115 37917.8429 258 50153.6322 21944.4115 259 50153.6322 50153.6322 260 50153.6322 50153.6322 261 50153.6322 50153.6322 262 33779.2527 50153.6322 263 55197.2354 33779.2527 264 50153.6322 55197.2354 265 50153.6322 50153.6322 266 50153.6322 50153.6322 267 72720.5979 50153.6322 268 50153.6322 72720.5979 269 -12148.3921 50153.6322 270 50153.6322 -12148.3921 271 50153.6322 50153.6322 272 50153.6322 50153.6322 273 7400.3458 50153.6322 274 50153.6322 7400.3458 275 50167.9578 50153.6322 276 21700.6342 50167.9578 277 50153.6322 21700.6342 278 49550.8234 50153.6322 279 41333.8592 49550.8234 280 50153.6322 41333.8592 281 50153.6322 50153.6322 282 67598.2433 50153.6322 283 50153.6322 67598.2433 284 33773.7574 50153.6322 285 50153.6322 33773.7574 286 35142.8678 50153.6322 287 46949.7047 35142.8678 288 50153.6322 46949.7047 289 72624.4917 50153.6322 290 50153.6322 72624.4917 291 50139.3065 50153.6322 292 31872.2839 50139.3065 293 50153.6322 31872.2839 294 28876.8779 50153.6322 295 71346.3640 28876.8779 296 50153.6322 71346.3640 297 50153.6322 50153.6322 298 50153.6322 50153.6322 299 27012.1639 50153.6322 300 41576.8749 27012.1639 301 32725.2768 41576.8749 302 47582.3580 32725.2768 303 50153.6322 47582.3580 304 95788.9728 50153.6322 305 29838.6008 95788.9728 306 50153.6322 29838.6008 307 7696.9595 50153.6322 308 28784.5077 7696.9595 309 50153.6322 28784.5077 310 50153.6322 50153.6322 311 43623.3622 50153.6322 312 50153.6322 43623.3622 313 23369.0833 50153.6322 314 34433.2307 23369.0833 315 68515.6647 34433.2307 316 51471.9224 68515.6647 317 -8251.2712 51471.9224 318 44332.8149 -8251.2712 319 50153.6322 44332.8149 320 50153.6322 50153.6322 321 50153.6322 50153.6322 322 50153.6322 50153.6322 323 50153.6322 50153.6322 324 50153.6322 50153.6322 325 46577.8273 50153.6322 326 50153.6322 46577.8273 327 50153.6322 50153.6322 328 50153.6322 50153.6322 329 50153.6322 50153.6322 330 50153.6322 50153.6322 331 50153.6322 50153.6322 332 50153.6322 50153.6322 333 42602.7629 50153.6322 334 49911.1953 42602.7629 335 38928.7664 49911.1953 336 52639.2174 38928.7664 337 50153.6322 52639.2174 338 46394.1763 50153.6322 339 41035.5768 46394.1763 340 6711.3425 41035.5768 341 50153.6322 6711.3425 342 50153.6322 50153.6322 343 50153.6322 50153.6322 344 52629.3460 50153.6322 345 50153.6322 52629.3460 346 37956.3424 50153.6322 347 50153.6322 37956.3424 348 50153.6322 50153.6322 349 13815.2135 50153.6322 350 59086.5397 13815.2135 351 40602.8151 59086.5397 352 38066.8135 40602.8151 353 43535.9424 38066.8135 354 50153.6322 43535.9424 355 50153.6322 50153.6322 356 37757.1676 50153.6322 357 50153.6322 37757.1676 358 50153.6322 50153.6322 359 50153.6322 50153.6322 360 44570.4409 50153.6322 361 50153.6322 44570.4409 362 50153.6322 50153.6322 363 -40412.1089 50153.6322 364 38915.5571 -40412.1089 365 50153.6322 38915.5571 366 50153.6322 50153.6322 367 50153.6322 50153.6322 368 50153.6322 50153.6322 369 52317.2864 50153.6322 370 50153.6322 52317.2864 371 79105.5636 50153.6322 372 31849.7933 79105.5636 373 24872.2149 31849.7933 374 50153.6322 24872.2149 375 50153.6322 50153.6322 376 35253.9660 50153.6322 377 39372.0092 35253.9660 378 -2616.8774 39372.0092 379 50153.6322 -2616.8774 380 50153.6322 50153.6322 381 43170.7401 50153.6322 382 93864.5602 43170.7401 383 50153.6322 93864.5602 384 35005.1569 50153.6322 385 -120729.1955 35005.1569 386 60161.2553 -120729.1955 387 21354.4899 60161.2553 388 43789.0774 21354.4899 389 32210.7648 43789.0774 390 31138.4345 32210.7648 391 13515.9146 31138.4345 392 -28176.9279 13515.9146 393 68800.2067 -28176.9279 394 -103849.1978 68800.2067 395 4765.1390 -103849.1978 396 -132209.7684 4765.1390 397 -437035.9785 -132209.7684 398 6759.5568 -437035.9785 399 44199.6986 6759.5568 400 -2415.6533 44199.6986 401 35454.1935 -2415.6533 402 15034.3991 35454.1935 403 -58691.2392 15034.3991 404 -17009.7401 -58691.2392 405 6083.1605 -17009.7401 406 16431.8199 6083.1605 407 3012.6169 16431.8199 408 17273.0037 3012.6169 409 -69144.3131 17273.0037 410 -74477.9193 -69144.3131 411 -163524.0366 -74477.9193 412 -33485.7746 -163524.0366 413 42911.0306 -33485.7746 414 -184109.5739 42911.0306 415 -105300.1328 -184109.5739 416 133358.5964 -105300.1328 417 -33479.9909 133358.5964 418 -63100.9831 -33479.9909 419 -40967.7465 -63100.9831 420 -49443.7944 -40967.7465 421 -217929.1146 -49443.7944 422 -562499.6734 -217929.1146 423 -461505.9634 -562499.6734 424 -55352.9789 -461505.9634 425 -268237.2948 -55352.9789 426 -311780.3800 -268237.2948 427 -1262344.5769 -311780.3800 428 -270458.3062 -1262344.5769 429 -367143.4164 -270458.3062 430 -573102.0748 -367143.4164 > plot(z,main=paste('Residual Lag plot, lowess, and regression line'), ylab='values of Residuals', xlab='lagged values of Residuals') > lines(lowess(z)) > abline(lm(z)) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/729l01291387065.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/829l01291387065.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/9di331291387065.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0)) > plot(mylm, las = 1, sub='Residual Diagnostics') > par(opar) > dev.off() null device 1 > if (n > n25) { + postscript(file="/var/www/html/freestat/rcomp/tmp/10di331291387065.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) + plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint') + grid() + dev.off() + } null device 1 > > #Note: the /var/www/html/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/freestat/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE) > a<-table.row.end(a) > myeq <- colnames(x)[1] > myeq <- paste(myeq, '[t] = ', sep='') > for (i in 1:k){ + if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '') + myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ') + if (rownames(mysum$coefficients)[i] != '(Intercept)') { + myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='') + if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='') + } + } Error: subscript out of bounds Execution halted