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(1 + ,162556 + ,1081 + ,213118 + ,6282154 + ,807 + ,1 + ,29790 + ,309 + ,81767 + ,4321023 + ,444 + ,1 + ,87550 + ,458 + ,153198 + ,4111912 + ,412 + ,0 + ,84738 + ,588 + ,-26007 + ,223193 + ,428 + ,1 + ,54660 + ,302 + ,126942 + ,1491348 + ,315 + ,1 + ,42634 + ,156 + ,157214 + ,1629616 + ,168 + ,0 + ,40949 + ,481 + ,129352 + ,1398893 + ,263 + ,1 + ,45187 + ,353 + ,234817 + ,1926517 + ,267 + ,1 + ,37704 + ,452 + ,60448 + ,983660 + ,228 + ,1 + ,16275 + ,109 + ,47818 + ,1443586 + ,129 + ,0 + ,25830 + ,115 + ,245546 + ,1073089 + ,104 + ,0 + ,12679 + ,110 + ,48020 + ,984885 + ,122 + ,1 + ,18014 + ,239 + ,-1710 + ,1405225 + ,393 + ,0 + ,43556 + ,247 + ,32648 + ,227132 + ,190 + ,1 + ,24811 + ,505 + ,95350 + ,929118 + ,280 + ,0 + ,6575 + ,159 + ,151352 + ,1071292 + ,63 + ,0 + ,7123 + ,109 + ,288170 + ,638830 + ,102 + ,1 + ,21950 + ,519 + ,114337 + ,856956 + ,265 + ,1 + ,37597 + ,248 + ,37884 + ,992426 + ,234 + ,0 + ,17821 + ,373 + ,122844 + ,444477 + ,277 + ,1 + ,12988 + ,119 + ,82340 + ,857217 + ,73 + ,1 + ,22330 + ,84 + ,79801 + ,711969 + ,67 + ,0 + ,13326 + ,102 + ,165548 + ,702380 + ,103 + ,0 + ,16189 + ,295 + ,116384 + ,358589 + ,290 + ,0 + ,7146 + ,105 + ,134028 + ,297978 + ,83 + ,0 + ,15824 + ,64 + ,63838 + ,585715 + ,56 + ,1 + ,27664 + ,282 + ,74996 + ,657954 + ,236 + ,0 + ,11920 + ,182 + ,31080 + ,209458 + ,73 + ,0 + ,8568 + ,37 + ,32168 + ,786690 + ,34 + ,0 + ,14416 + ,361 + ,49857 + ,439798 + ,139 + ,1 + ,3369 + ,28 + ,87161 + ,688779 + ,26 + ,1 + ,11819 + ,85 + ,106113 + ,574339 + ,70 + ,1 + ,6984 + ,45 + ,80570 + ,741409 + ,40 + ,1 + ,4519 + ,49 + ,102129 + ,597793 + ,42 + ,0 + ,2220 + ,22 + ,301670 + ,644190 + ,12 + ,0 + ,18562 + ,155 + ,102313 + ,377934 + ,211 + ,0 + ,10327 + ,91 + ,88577 + ,640273 + ,74 + ,1 + ,5336 + ,81 + ,112477 + ,697458 + ,80 + ,1 + ,2365 + ,79 + ,191778 + ,550608 + ,83 + ,0 + ,4069 + ,145 + ,79804 + ,207393 + ,131 + ,0 + ,8636 + ,855 + ,128294 + ,301607 + ,203 + ,0 + ,13718 + ,61 + ,96448 + ,345783 + ,56 + ,0 + ,4525 + ,226 + ,93811 + ,501749 + ,89 + ,0 + ,6869 + ,105 + ,117520 + ,379983 + ,88 + ,0 + ,4628 + ,62 + ,69159 + ,387475 + ,39 + ,1 + ,3689 + ,25 + ,101792 + ,377305 + ,25 + ,1 + ,4891 + ,217 + ,210568 + ,370837 + ,49 + ,1 + ,7489 + ,322 + ,136996 + ,430866 + ,149 + ,0 + ,4901 + ,84 + ,121920 + ,469107 + ,58 + ,0 + ,2284 + ,33 + ,76403 + ,194493 + ,41 + ,1 + ,3160 + ,108 + ,108094 + ,530670 + ,90 + ,1 + ,4150 + ,150 + ,134759 + ,518365 + ,136 + ,1 + ,7285 + ,115 + ,188873 + ,491303 + ,97 + ,1 + ,1134 + ,162 + ,146216 + ,527021 + ,63 + ,1 + ,4658 + ,158 + ,156608 + ,233773 + ,114 + ,0 + ,2384 + ,97 + ,61348 + ,405972 + ,77 + ,0 + ,3748 + ,9 + ,50350 + ,652925 + ,6 + ,0 + ,5371 + ,66 + ,87720 + ,446211 + ,47 + ,0 + ,1285 + ,107 + ,99489 + ,341340 + ,51 + ,1 + ,9327 + ,101 + ,87419 + ,387699 + ,85 + ,1 + ,5565 + ,47 + ,94355 + ,493408 + ,43 + ,0 + ,1528 + ,38 + ,60326 + ,146494 + ,32 + ,1 + ,3122 + ,34 + ,94670 + ,414462 + ,25 + ,1 + ,7561 + ,87 + ,82425 + ,364304 + ,77 + ,0 + ,2675 + ,79 + ,59017 + ,355178 + ,54 + ,0 + ,13253 + ,947 + ,90829 + ,357760 + ,251 + ,0 + ,880 + ,74 + ,80791 + ,261216 + ,15 + ,1 + ,2053 + ,53 + ,100423 + ,397144 + ,44 + ,0 + ,1424 + ,94 + ,131116 + ,374943 + ,73 + ,1 + ,4036 + ,63 + ,100269 + ,424898 + ,85 + ,1 + ,3045 + ,58 + ,27330 + ,202055 + ,49 + ,0 + ,5119 + ,49 + ,39039 + ,378525 + ,38 + ,0 + ,1431 + ,34 + ,106885 + ,310768 + ,35 + ,0 + ,554 + ,11 + ,79285 + ,325738 + ,9 + ,0 + ,1975 + ,35 + ,118881 + ,394510 + ,34 + ,1 + ,1765 + ,20 + ,77623 + ,247060 + ,20 + ,0 + ,1012 + ,47 + ,114768 + ,368078 + ,29 + ,0 + ,810 + ,43 + ,74015 + ,236761 + ,11 + ,0 + ,1280 + ,117 + ,69465 + ,312378 + ,52 + ,1 + ,666 + ,171 + ,117869 + ,339836 + ,13 + ,0 + ,1380 + ,26 + ,60982 + ,347385 + ,29 + ,1 + ,4677 + ,75 + ,90131 + ,426280 + ,66 + ,0 + ,876 + ,59 + ,138971 + ,352850 + ,33 + ,0 + ,814 + ,18 + ,39625 + ,301881 + ,15 + ,0 + ,514 + ,15 + ,102725 + ,377516 + ,15 + ,1 + ,5692 + ,72 + ,64239 + ,357312 + ,68 + ,0 + ,3642 + ,86 + ,90262 + ,458343 + ,100 + ,0 + ,540 + ,14 + ,103960 + ,354228 + ,13 + ,0 + ,2099 + ,64 + ,106611 + ,308636 + ,45 + ,0 + ,567 + ,11 + ,103345 + ,386212 + ,14 + ,0 + ,2001 + ,52 + ,95551 + ,393343 + ,36 + ,1 + ,2949 + ,41 + ,82903 + ,378509 + ,40 + ,0 + ,2253 + ,99 + ,63593 + ,452469 + ,68 + ,1 + ,6533 + ,75 + ,126910 + ,364839 + ,29 + ,0 + ,1889 + ,45 + ,37527 + ,358649 + ,43 + ,1 + ,3055 + ,43 + ,60247 + ,376641 + ,30 + ,0 + ,272 + ,8 + ,112995 + ,429112 + ,9 + ,1 + ,1414 + ,198 + ,70184 + ,330546 + ,22 + ,0 + ,2564 + ,22 + ,130140 + ,403560 + ,19 + ,1 + ,1383 + ,11 + ,73221 + ,317892 + ,9 + ,1 + ,1261 + ,33 + ,76114 + ,307528 + ,31 + ,0 + ,975 + ,23 + ,90534 + ,235133 + ,19 + ,0 + ,3366 + ,80 + ,108479 + ,299243 + ,55 + ,0 + ,576 + ,18 + ,113761 + ,314073 + ,8 + ,0 + ,1686 + ,40 + ,68696 + ,368186 + ,28 + ,0 + ,746 + ,23 + ,71561 + ,269661 + ,29 + ,1 + ,3192 + ,60 + ,59831 + ,125390 + ,48 + ,1 + ,2045 + ,20 + ,97890 + ,510834 + ,16 + ,0 + ,5702 + ,61 + ,101481 + ,321896 + ,47 + ,1 + ,1932 + ,36 + ,72954 + ,249898 + ,20 + ,0 + ,936 + ,30 + ,67939 + ,408881 + ,22 + ,1 + ,3437 + ,47 + ,48022 + ,158492 + ,33 + ,0 + ,5131 + ,71 + ,86111 + ,292154 + ,44 + ,1 + ,2397 + ,14 + ,74020 + ,289513 + ,13 + ,1 + ,1389 + ,9 + ,57530 + ,378049 + ,6 + ,0 + ,1503 + ,39 + ,56364 + ,343466 + ,35 + ,0 + ,402 + ,26 + ,84990 + ,332743 + ,8 + ,0 + ,2239 + ,21 + ,88590 + ,442882 + ,17 + ,1 + ,2234 + ,16 + ,77200 + ,214215 + ,11 + ,0 + ,837 + ,69 + ,61262 + ,315688 + ,21 + ,0 + ,10579 + ,92 + ,110309 + ,375195 + ,92 + ,0 + ,875 + ,14 + ,67000 + ,334280 + ,12 + ,0 + ,1585 + ,107 + ,93099 + ,355864 + ,112 + ,1 + ,1659 + ,29 + ,107577 + ,480382 + ,25 + ,1 + ,2647 + ,37 + ,62920 + ,353058 + ,17 + ,1 + ,3294 + ,23 + ,75832 + ,217193 + ,23 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,94 + ,7 + ,60793 + ,314533 + ,10 + ,0 + ,422 + ,28 + ,57935 + ,318056 + ,23 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,34 + ,8 + ,60630 + ,314353 + ,7 + ,0 + ,1558 + ,63 + ,55637 + ,369448 + ,25 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,43 + ,3 + ,60887 + ,312846 + ,20 + ,0 + ,645 + ,5 + ,60720 + ,312075 + ,4 + ,0 + ,316 + ,9 + ,60505 + ,315009 + ,4 + ,0 + ,115 + ,13 + ,60945 + ,318903 + ,10 + ,0 + ,5 + ,2 + ,60720 + ,314887 + ,1 + ,0 + ,897 + ,5 + ,60720 + ,314913 + ,4 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,389 + ,14 + ,58990 + ,325506 + ,8 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,1002 + ,15 + ,56750 + ,298568 + ,11 + ,0 + ,36 + ,3 + ,60894 + ,315834 + ,4 + ,0 + ,460 + ,15 + ,63346 + ,329784 + ,15 + ,0 + ,309 + ,11 + ,56535 + ,312878 + ,9 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,9 + ,6 + ,60835 + ,314987 + ,7 + ,0 + ,271 + ,2 + ,60720 + ,325249 + ,2 + ,0 + ,14 + ,1 + ,61016 + ,315877 + ,0 + ,0 + ,520 + ,10 + ,58650 + ,291650 + ,7 + ,0 + ,1766 + ,73 + ,60438 + ,305959 + ,46 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,5 + ,0 + ,458 + ,11 + ,58625 + ,297765 + ,7 + ,0 + ,20 + ,3 + ,60938 + ,315245 + ,2 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,98 + ,2 + ,61490 + ,315236 + ,2 + ,0 + ,405 + ,7 + ,60845 + ,336425 + ,5 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,483 + ,27 + ,60830 + ,306268 + ,7 + ,0 + ,454 + ,51 + ,63261 + ,302187 + ,24 + ,0 + ,47 + ,3 + ,60720 + ,314882 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,757 + ,19 + ,45689 + ,382712 + ,18 + ,0 + ,4655 + ,393 + ,60720 + ,341570 + ,55 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,36 + ,4 + ,61564 + ,312412 + ,3 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,203 + ,9 + ,61938 + ,309596 + ,9 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,126 + ,10 + ,60951 + ,315547 + ,8 + ,0 + ,400 + ,152 + ,60720 + ,313267 + ,113 + ,1 + ,71 + ,1 + ,60745 + ,316176 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,972 + ,34 + ,71642 + ,359335 + ,19 + ,1 + ,531 + ,10 + ,71641 + ,330068 + ,11 + ,0 + ,2461 + ,57 + ,55792 + ,314289 + ,25 + ,1 + ,378 + ,52 + ,71873 + ,297413 + ,16 + ,1 + ,23 + ,5 + ,62555 + ,314806 + ,5 + ,1 + ,638 + ,14 + ,60370 + ,333210 + ,11 + ,1 + ,2300 + ,29 + ,64873 + ,352108 + ,23 + ,0 + ,149 + ,5 + ,62041 + ,313332 + ,6 + ,0 + ,226 + ,5 + ,65745 + ,291787 + ,5 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,275 + ,4 + ,59500 + ,318745 + ,7 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,141 + ,6 + ,61630 + ,315366 + ,7 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,28 + ,2 + ,60890 + ,315688 + ,3 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,4980 + ,91 + ,113521 + ,409642 + ,89 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,472 + ,20 + ,80045 + ,269587 + ,19 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,203 + ,27 + ,50804 + ,300962 + ,12 + ,1 + ,496 + ,17 + ,87390 + ,325479 + ,12 + ,1 + ,10 + ,2 + ,61656 + ,316155 + ,5 + ,1 + ,63 + ,4 + ,65688 + ,318574 + ,2 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,1136 + ,32 + ,48522 + ,343613 + ,26 + ,1 + ,265 + ,31 + ,60720 + ,306948 + ,3 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,267 + ,32 + ,57640 + ,330059 + ,11 + ,0 + ,474 + ,20 + ,61977 + ,288985 + ,10 + ,0 + ,534 + ,7 + ,62620 + ,304485 + ,5 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,2 + ,0 + ,15 + ,8 + ,60831 + ,315688 + ,6 + ,0 + ,397 + ,28 + ,60646 + ,317736 + ,7 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,2 + ,0 + ,1866 + ,29 + ,56225 + ,322331 + ,28 + ,0 + ,288 + ,4 + ,60510 + ,296656 + ,3 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,3 + ,2 + ,60698 + ,315354 + ,1 + ,0 + ,468 + ,21 + ,60720 + ,312161 + ,20 + ,0 + ,20 + ,2 + ,60805 + ,315576 + ,1 + ,0 + ,278 + ,26 + ,61404 + ,314922 + ,22 + ,0 + ,61 + ,14 + ,60720 + ,314551 + ,9 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,192 + ,4 + ,65276 + ,312339 + ,2 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,317 + ,9 + ,63915 + ,298700 + ,7 + ,0 + ,738 + ,10 + ,60720 + ,321376 + ,9 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,368 + ,17 + ,61686 + ,303230 + ,13 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,2 + ,1 + ,60743 + ,315487 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,53 + ,6 + ,60349 + ,315793 + ,6 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,94 + ,3 + ,61360 + ,312887 + ,3 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,24 + ,8 + ,59818 + ,315637 + ,7 + ,0 + ,2332 + ,4 + ,72680 + ,324385 + ,2 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,131 + ,11 + ,61808 + ,308989 + ,15 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,206 + ,9 + ,53110 + ,296702 + ,9 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,167 + ,2 + ,64245 + ,307322 + ,1 + ,0 + ,622 + ,73 + ,73007 + ,304376 + ,38 + ,0 + ,2328 + ,94 + ,82732 + ,253588 + ,57 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,365 + ,8 + ,54820 + ,309560 + ,7 + ,0 + ,364 + ,35 + ,47705 + ,298466 + ,26 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,226 + ,12 + ,72835 + ,343929 + ,13 + ,0 + ,307 + ,15 + ,58856 + ,331955 + ,10 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,188 + ,11 + ,77655 + ,381180 + ,9 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,138 + ,6 + ,69817 + ,331420 + ,26 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,125 + ,12 + ,60798 + ,310201 + ,19 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,282 + ,30 + ,62452 + ,320016 + ,12 + ,0 + ,335 + ,33 + ,64175 + ,320398 + ,23 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,1324 + ,117 + ,67440 + ,291841 + ,29 + ,0 + ,176 + ,28 + ,68136 + ,310670 + ,8 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,249 + ,72 + ,56726 + ,313491 + ,26 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,333 + ,13 + ,70811 + ,331323 + ,9 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,601 + ,6 + ,60720 + ,319210 + ,5 + ,0 + ,30 + ,4 + ,62045 + ,318098 + ,3 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,249 + ,62 + ,54323 + ,292754 + ,13 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,165 + ,24 + ,62841 + ,325176 + ,12 + ,0 + ,453 + ,21 + ,81125 + ,365959 + ,19 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,53 + ,14 + ,59506 + ,302409 + ,10 + ,1 + ,382 + ,21 + ,59365 + ,340968 + ,9 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,9 + ,1 + ,30 + ,4 + ,60798 + ,313164 + ,4 + ,0 + ,290 + ,2 + ,58790 + ,301164 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,366 + ,53 + ,61808 + ,344425 + ,14 + ,0 + ,2 + ,9 + ,60735 + ,315394 + ,12 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,209 + ,13 + ,64016 + ,316647 + ,19 + ,0 + ,384 + ,22 + ,54683 + ,309836 + ,17 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,365 + ,83 + ,87192 + ,346611 + ,32 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,49 + ,8 + ,64107 + ,322031 + ,14 + ,0 + ,3 + ,4 + ,60761 + ,315656 + ,8 + ,0 + ,133 + ,14 + ,65990 + ,339445 + ,4 + ,0 + ,32 + ,1 + ,59988 + ,314964 + ,0 + ,0 + ,368 + ,17 + ,61167 + ,297141 + ,20 + ,0 + ,1 + ,6 + ,60719 + ,315372 + ,5 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,22 + ,2 + ,60722 + ,312502 + ,1 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,96 + ,5 + ,60379 + ,313729 + ,4 + ,0 + ,1 + ,2 + ,60727 + ,315388 + ,1 + ,0 + ,314 + ,5 + ,60720 + ,315371 + ,4 + ,0 + ,844 + ,78 + ,60925 + ,296139 + ,20 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,26 + ,1 + ,60896 + ,313880 + ,1 + ,0 + ,125 + ,13 + ,59734 + ,317698 + ,10 + ,0 + ,304 + ,15 + ,62969 + ,295580 + ,12 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,621 + ,48 + ,60720 + ,308256 + ,13 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,119 + ,6 + ,59118 + ,303677 + ,3 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,1595 + ,17 + ,60720 + ,319369 + ,10 + ,0 + ,312 + ,14 + ,58598 + ,318690 + ,3 + ,0 + ,60 + ,10 + ,61124 + ,314049 + ,7 + ,0 + ,587 + ,12 + ,59595 + ,325699 + ,10 + ,0 + ,135 + ,2 + ,62065 + ,314210 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,514 + ,52 + ,78780 + ,322378 + ,15 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,1 + ,4 + ,60722 + ,315398 + ,4 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,1763 + ,24 + ,61600 + ,308336 + ,28 + ,0 + ,180 + ,11 + ,59635 + ,316386 + ,9 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,218 + ,21 + ,60720 + ,315553 + ,7 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,448 + ,40 + ,59781 + ,323361 + ,7 + ,0 + ,227 + ,9 + ,76644 + ,336639 + ,7 + ,0 + ,174 + ,1 + ,64820 + ,307424 + ,3 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,121 + ,24 + ,56178 + ,295370 + ,11 + ,0 + ,607 + ,11 + ,60436 + ,322340 + ,7 + ,0 + ,2212 + ,14 + ,60720 + ,319864 + ,10 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,530 + ,60 + ,73433 + ,317291 + ,18 + ,0 + ,571 + ,80 + ,41477 + ,280398 + ,14 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,78 + ,16 + ,62700 + ,317330 + ,12 + ,0 + ,2489 + ,40 + ,67804 + ,238125 + ,29 + ,0 + ,131 + ,6 + ,59661 + ,327071 + ,3 + ,0 + ,923 + ,8 + ,58620 + ,309038 + ,6 + ,0 + ,72 + ,3 + ,60398 + ,314210 + ,3 + ,0 + ,572 + ,16 + ,58580 + ,307930 + ,8 + ,0 + ,397 + ,10 + ,62710 + ,322327 + ,10 + ,0 + ,450 + ,8 + ,59325 + ,292136 + ,6 + ,0 + ,622 + ,7 + ,60950 + ,263276 + ,8 + ,0 + ,694 + ,8 + ,68060 + ,367655 + ,6 + ,1 + ,3425 + ,12 + ,83620 + ,283910 + ,9 + ,0 + ,562 + ,13 + ,58456 + ,283587 + ,8 + ,0 + ,4917 + ,42 + ,52811 + ,243650 + ,26 + ,1 + ,1442 + ,118 + ,121173 + ,438493 + ,239 + ,0 + ,529 + ,9 + ,63870 + ,296261 + ,7 + ,1 + ,2126 + ,138 + ,21001 + ,230621 + ,41 + ,0 + ,1061 + ,5 + ,70415 + ,304252 + ,3 + ,0 + ,776 + ,9 + ,64230 + ,333505 + ,8 + ,0 + ,611 + ,8 + ,59190 + ,296919 + ,6 + ,1 + ,1526 + ,25 + ,69351 + ,278990 + ,21 + ,0 + ,592 + ,7 + ,64270 + ,276898 + ,7 + ,0 + ,1182 + ,13 + ,70694 + ,327007 + ,11 + ,0 + ,621 + ,16 + ,68005 + ,317046 + ,11 + ,0 + ,989 + ,11 + ,58930 + ,304555 + ,12 + ,0 + ,438 + ,11 + ,58320 + ,298096 + ,9 + ,0 + ,726 + ,3 + ,69980 + ,231861 + ,3 + ,0 + ,1303 + ,61 + ,69863 + ,309422 + ,57 + ,1 + ,7419 + ,29 + ,63255 + ,286963 + ,21 + ,1 + ,1164 + ,17 + ,57320 + ,269753 + ,15 + ,1 + ,3310 + ,33 + ,75230 + ,448243 + ,32 + ,0 + ,1920 + ,15 + ,79420 + ,165404 + ,11 + ,0 + ,965 + ,3 + ,73490 + ,204325 + ,2 + ,0 + ,3256 + ,66 + ,35250 + ,407159 + ,23 + ,1 + ,1135 + ,17 + ,62285 + ,290476 + ,20 + ,0 + ,1270 + ,26 + ,69206 + ,275311 + ,24 + ,0 + ,661 + ,3 + ,65920 + ,246541 + ,1 + ,0 + ,1013 + ,2 + ,69770 + ,253468 + ,1 + ,0 + ,2844 + ,67 + ,72683 + ,240897 + ,74 + ,1 + ,11528 + ,70 + ,-14545 + ,-83265 + ,68 + ,0 + ,6526 + ,26 + ,55830 + ,-42143 + ,20 + ,0 + ,2264 + ,24 + ,55174 + ,272713 + ,20 + ,1 + ,5109 + ,97 + ,67038 + ,215362 + ,82 + ,0 + ,3999 + ,30 + ,51252 + ,42754 + ,21 + ,0 + ,35624 + ,223 + ,157278 + ,306275 + ,244 + ,0 + ,9252 + ,48 + ,79510 + ,253537 + ,32 + ,0 + ,15236 + ,90 + ,77440 + ,372631 + ,86 + ,0 + ,18073 + ,180 + ,27284 + ,-7170 + ,69) + ,dim=c(6 + ,431) + ,dimnames=list(c('Group' + ,'Costs' + ,'Trades' + ,'Dividends' + ,'Wealth' + ,'Orders ') + ,1:431)) > y <- array(NA,dim=c(6,431),dimnames=list(c('Group','Costs','Trades','Dividends','Wealth','Orders '),1:431)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par3 = 'No Linear Trend' > par2 = 'Do not include Seasonal Dummies' > par1 = '5' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > library(lattice) > library(lmtest) Loading required package: zoo Attaching package: 'zoo' The following object(s) are masked from package:base : as.Date.numeric > n25 <- 25 #minimum number of obs. for Goldfeld-Quandt test > par1 <- as.numeric(par1) > x <- t(y) > k <- length(x[1,]) > n <- length(x[,1]) > x1 <- cbind(x[,par1], x[,1:k!=par1]) > mycolnames <- c(colnames(x)[par1], colnames(x)[1:k!=par1]) > colnames(x1) <- mycolnames #colnames(x)[par1] > x <- x1 > if (par3 == 'First Differences'){ + x2 <- array(0, dim=c(n-1,k), dimnames=list(1:(n-1), paste('(1-B)',colnames(x),sep=''))) + for (i in 1:n-1) { + for (j in 1:k) { + x2[i,j] <- x[i+1,j] - x[i,j] + } + } + x <- x2 + } > if (par2 == 'Include Monthly Dummies'){ + x2 <- array(0, dim=c(n,11), dimnames=list(1:n, paste('M', seq(1:11), sep =''))) + for (i in 1:11){ + x2[seq(i,n,12),i] <- 1 + } + x <- cbind(x, x2) + } > if (par2 == 'Include Quarterly Dummies'){ + x2 <- array(0, dim=c(n,3), dimnames=list(1:n, paste('Q', seq(1:3), sep =''))) + for (i in 1:3){ + x2[seq(i,n,4),i] <- 1 + } + x <- cbind(x, x2) + } > k <- length(x[1,]) > if (par3 == 'Linear Trend'){ + x <- cbind(x, c(1:n)) + colnames(x)[k+1] <- 't' + } > x Wealth Group Costs Trades Dividends Orders\r\r 1 6282154 1 162556 1081 213118 807 2 4321023 1 29790 309 81767 444 3 4111912 1 87550 458 153198 412 4 223193 0 84738 588 -26007 428 5 1491348 1 54660 302 126942 315 6 1629616 1 42634 156 157214 168 7 1398893 0 40949 481 129352 263 8 1926517 1 45187 353 234817 267 9 983660 1 37704 452 60448 228 10 1443586 1 16275 109 47818 129 11 1073089 0 25830 115 245546 104 12 984885 0 12679 110 48020 122 13 1405225 1 18014 239 -1710 393 14 227132 0 43556 247 32648 190 15 929118 1 24811 505 95350 280 16 1071292 0 6575 159 151352 63 17 638830 0 7123 109 288170 102 18 856956 1 21950 519 114337 265 19 992426 1 37597 248 37884 234 20 444477 0 17821 373 122844 277 21 857217 1 12988 119 82340 73 22 711969 1 22330 84 79801 67 23 702380 0 13326 102 165548 103 24 358589 0 16189 295 116384 290 25 297978 0 7146 105 134028 83 26 585715 0 15824 64 63838 56 27 657954 1 27664 282 74996 236 28 209458 0 11920 182 31080 73 29 786690 0 8568 37 32168 34 30 439798 0 14416 361 49857 139 31 688779 1 3369 28 87161 26 32 574339 1 11819 85 106113 70 33 741409 1 6984 45 80570 40 34 597793 1 4519 49 102129 42 35 644190 0 2220 22 301670 12 36 377934 0 18562 155 102313 211 37 640273 0 10327 91 88577 74 38 697458 1 5336 81 112477 80 39 550608 1 2365 79 191778 83 40 207393 0 4069 145 79804 131 41 301607 0 8636 855 128294 203 42 345783 0 13718 61 96448 56 43 501749 0 4525 226 93811 89 44 379983 0 6869 105 117520 88 45 387475 0 4628 62 69159 39 46 377305 1 3689 25 101792 25 47 370837 1 4891 217 210568 49 48 430866 1 7489 322 136996 149 49 469107 0 4901 84 121920 58 50 194493 0 2284 33 76403 41 51 530670 1 3160 108 108094 90 52 518365 1 4150 150 134759 136 53 491303 1 7285 115 188873 97 54 527021 1 1134 162 146216 63 55 233773 1 4658 158 156608 114 56 405972 0 2384 97 61348 77 57 652925 0 3748 9 50350 6 58 446211 0 5371 66 87720 47 59 341340 0 1285 107 99489 51 60 387699 1 9327 101 87419 85 61 493408 1 5565 47 94355 43 62 146494 0 1528 38 60326 32 63 414462 1 3122 34 94670 25 64 364304 1 7561 87 82425 77 65 355178 0 2675 79 59017 54 66 357760 0 13253 947 90829 251 67 261216 0 880 74 80791 15 68 397144 1 2053 53 100423 44 69 374943 0 1424 94 131116 73 70 424898 1 4036 63 100269 85 71 202055 1 3045 58 27330 49 72 378525 0 5119 49 39039 38 73 310768 0 1431 34 106885 35 74 325738 0 554 11 79285 9 75 394510 0 1975 35 118881 34 76 247060 1 1765 20 77623 20 77 368078 0 1012 47 114768 29 78 236761 0 810 43 74015 11 79 312378 0 1280 117 69465 52 80 339836 1 666 171 117869 13 81 347385 0 1380 26 60982 29 82 426280 1 4677 75 90131 66 83 352850 0 876 59 138971 33 84 301881 0 814 18 39625 15 85 377516 0 514 15 102725 15 86 357312 1 5692 72 64239 68 87 458343 0 3642 86 90262 100 88 354228 0 540 14 103960 13 89 308636 0 2099 64 106611 45 90 386212 0 567 11 103345 14 91 393343 0 2001 52 95551 36 92 378509 1 2949 41 82903 40 93 452469 0 2253 99 63593 68 94 364839 1 6533 75 126910 29 95 358649 0 1889 45 37527 43 96 376641 1 3055 43 60247 30 97 429112 0 272 8 112995 9 98 330546 1 1414 198 70184 22 99 403560 0 2564 22 130140 19 100 317892 1 1383 11 73221 9 101 307528 1 1261 33 76114 31 102 235133 0 975 23 90534 19 103 299243 0 3366 80 108479 55 104 314073 0 576 18 113761 8 105 368186 0 1686 40 68696 28 106 269661 0 746 23 71561 29 107 125390 1 3192 60 59831 48 108 510834 1 2045 20 97890 16 109 321896 0 5702 61 101481 47 110 249898 1 1932 36 72954 20 111 408881 0 936 30 67939 22 112 158492 1 3437 47 48022 33 113 292154 0 5131 71 86111 44 114 289513 1 2397 14 74020 13 115 378049 1 1389 9 57530 6 116 343466 0 1503 39 56364 35 117 332743 0 402 26 84990 8 118 442882 0 2239 21 88590 17 119 214215 1 2234 16 77200 11 120 315688 0 837 69 61262 21 121 375195 0 10579 92 110309 92 122 334280 0 875 14 67000 12 123 355864 0 1585 107 93099 112 124 480382 1 1659 29 107577 25 125 353058 1 2647 37 62920 17 126 217193 1 3294 23 75832 23 127 315380 0 0 0 60720 0 128 314533 0 94 7 60793 10 129 318056 0 422 28 57935 23 130 315380 0 0 0 60720 0 131 314353 0 34 8 60630 7 132 369448 0 1558 63 55637 25 133 315380 0 0 0 60720 0 134 312846 0 43 3 60887 20 135 312075 0 645 5 60720 4 136 315009 0 316 9 60505 4 137 318903 0 115 13 60945 10 138 314887 0 5 2 60720 1 139 314913 0 897 5 60720 4 140 315380 0 0 0 60720 0 141 325506 0 389 14 58990 8 142 315380 0 0 0 60720 0 143 298568 0 1002 15 56750 11 144 315834 0 36 3 60894 4 145 329784 0 460 15 63346 15 146 312878 0 309 11 56535 9 147 315380 0 0 0 60720 0 148 314987 0 9 6 60835 7 149 325249 0 271 2 60720 2 150 315877 0 14 1 61016 0 151 291650 0 520 10 58650 7 152 305959 0 1766 73 60438 46 153 315380 0 0 0 60720 5 154 297765 0 458 11 58625 7 155 315245 0 20 3 60938 2 156 315380 0 0 0 60720 0 157 315380 0 0 0 60720 0 158 315236 0 98 2 61490 2 159 336425 0 405 7 60845 5 160 315380 0 0 0 60720 0 161 315380 0 0 0 60720 0 162 315380 0 0 0 60720 0 163 315380 0 0 0 60720 0 164 306268 0 483 27 60830 7 165 302187 0 454 51 63261 24 166 314882 0 47 3 60720 1 167 315380 0 0 0 60720 0 168 382712 0 757 19 45689 18 169 341570 0 4655 393 60720 55 170 315380 0 0 0 60720 0 171 315380 0 0 0 60720 0 172 312412 0 36 4 61564 3 173 315380 0 0 0 60720 0 174 309596 0 203 9 61938 9 175 315380 0 0 0 60720 0 176 315547 0 126 10 60951 8 177 313267 0 400 152 60720 113 178 316176 1 71 1 60745 0 179 315380 0 0 0 60720 0 180 315380 0 0 0 60720 0 181 359335 0 972 34 71642 19 182 330068 1 531 10 71641 11 183 314289 0 2461 57 55792 25 184 297413 1 378 52 71873 16 185 314806 1 23 5 62555 5 186 333210 1 638 14 60370 11 187 352108 1 2300 29 64873 23 188 313332 0 149 5 62041 6 189 291787 0 226 5 65745 5 190 315380 0 0 0 60720 0 191 318745 0 275 4 59500 7 192 315380 0 0 0 60720 0 193 315366 0 141 6 61630 7 194 315380 0 0 0 60720 0 195 315688 0 28 2 60890 3 196 315380 1 0 0 60720 0 197 409642 1 4980 91 113521 89 198 315380 1 0 0 60720 0 199 315380 1 0 0 60720 0 200 269587 1 472 20 80045 19 201 315380 1 0 0 60720 0 202 315380 1 0 0 60720 0 203 315380 1 0 0 60720 0 204 300962 1 203 27 50804 12 205 325479 1 496 17 87390 12 206 316155 1 10 2 61656 5 207 318574 1 63 4 65688 2 208 315380 1 0 0 60720 0 209 343613 1 1136 32 48522 26 210 306948 1 265 31 60720 3 211 315380 0 0 0 60720 0 212 315380 0 0 0 60720 0 213 330059 0 267 32 57640 11 214 288985 0 474 20 61977 10 215 304485 0 534 7 62620 5 216 315380 0 0 0 60720 2 217 315688 0 15 8 60831 6 218 317736 0 397 28 60646 7 219 315380 0 0 0 60720 2 220 322331 0 1866 29 56225 28 221 296656 0 288 4 60510 3 222 315380 0 0 0 60720 0 223 315354 0 3 2 60698 1 224 312161 0 468 21 60720 20 225 315576 0 20 2 60805 1 226 314922 0 278 26 61404 22 227 314551 0 61 14 60720 9 228 315380 0 0 0 60720 0 229 312339 0 192 4 65276 2 230 315380 0 0 0 60720 0 231 298700 0 317 9 63915 7 232 321376 0 738 10 60720 9 233 315380 0 0 0 60720 0 234 303230 0 368 17 61686 13 235 315380 0 0 0 60720 0 236 315487 0 2 1 60743 0 237 315380 0 0 0 60720 0 238 315793 0 53 6 60349 6 239 315380 0 0 0 60720 0 240 315380 0 0 0 60720 0 241 315380 0 0 0 60720 0 242 312887 0 94 3 61360 3 243 315380 0 0 0 60720 0 244 315637 0 24 8 59818 7 245 324385 0 2332 4 72680 2 246 315380 1 0 0 60720 0 247 315380 1 0 0 60720 0 248 308989 0 131 11 61808 15 249 315380 1 0 0 60720 0 250 315380 1 0 0 60720 0 251 296702 0 206 9 53110 9 252 315380 1 0 0 60720 0 253 307322 0 167 2 64245 1 254 304376 0 622 73 73007 38 255 253588 0 2328 94 82732 57 256 315380 0 0 0 60720 0 257 309560 0 365 8 54820 7 258 298466 0 364 35 47705 26 259 315380 1 0 0 60720 0 260 315380 0 0 0 60720 0 261 315380 1 0 0 60720 0 262 315380 1 0 0 60720 0 263 343929 0 226 12 72835 13 264 331955 0 307 15 58856 10 265 315380 1 0 0 60720 0 266 315380 0 0 0 60720 0 267 315380 1 0 0 60720 0 268 381180 0 188 11 77655 9 269 315380 1 0 0 60720 0 270 331420 0 138 6 69817 26 271 315380 1 0 0 60720 0 272 315380 1 0 0 60720 0 273 315380 1 0 0 60720 0 274 310201 0 125 12 60798 19 275 315380 0 0 0 60720 0 276 320016 0 282 30 62452 12 277 320398 0 335 33 64175 23 278 315380 0 0 0 60720 0 279 291841 1 1324 117 67440 29 280 310670 0 176 28 68136 8 281 315380 1 0 0 60720 0 282 315380 0 0 0 60720 0 283 313491 0 249 72 56726 26 284 315380 0 0 0 60720 0 285 331323 0 333 13 70811 9 286 315380 0 0 0 60720 0 287 319210 1 601 6 60720 5 288 318098 0 30 4 62045 3 289 315380 1 0 0 60720 0 290 292754 0 249 62 54323 13 291 315380 1 0 0 60720 0 292 325176 0 165 24 62841 12 293 365959 0 453 21 81125 19 294 315380 0 0 0 60720 0 295 302409 0 53 14 59506 10 296 340968 1 382 21 59365 9 297 315380 0 0 0 60720 0 298 315380 1 0 0 60720 0 299 315380 0 0 0 60720 0 300 315380 0 0 0 60720 9 301 313164 1 30 4 60798 4 302 301164 0 290 2 58790 1 303 315380 0 0 0 60720 1 304 315380 1 0 0 60720 0 305 344425 0 366 53 61808 14 306 315394 0 2 9 60735 12 307 315380 1 0 0 60720 0 308 316647 1 209 13 64016 19 309 309836 0 384 22 54683 17 310 315380 1 0 0 60720 0 311 315380 1 0 0 60720 0 312 346611 0 365 83 87192 32 313 315380 1 0 0 60720 0 314 322031 1 49 8 64107 14 315 315656 0 3 4 60761 8 316 339445 0 133 14 65990 4 317 314964 0 32 1 59988 0 318 297141 0 368 17 61167 20 319 315372 0 1 6 60719 5 320 315380 0 0 0 60720 0 321 315380 0 0 0 60720 0 322 315380 0 0 0 60720 0 323 315380 0 0 0 60720 0 324 315380 1 0 0 60720 0 325 315380 0 0 0 60720 0 326 312502 0 22 2 60722 1 327 315380 1 0 0 60720 0 328 315380 1 0 0 60720 0 329 315380 0 0 0 60720 0 330 315380 0 0 0 60720 0 331 315380 0 0 0 60720 0 332 315380 0 0 0 60720 0 333 315380 0 0 0 60720 0 334 313729 0 96 5 60379 4 335 315388 0 1 2 60727 1 336 315371 0 314 5 60720 4 337 296139 0 844 78 60925 20 338 315380 0 0 0 60720 0 339 313880 0 26 1 60896 1 340 317698 0 125 13 59734 10 341 295580 0 304 15 62969 12 342 315380 0 0 0 60720 0 343 315380 0 0 0 60720 0 344 315380 0 0 0 60720 0 345 308256 0 621 48 60720 13 346 315380 0 0 0 60720 0 347 303677 0 119 6 59118 3 348 315380 0 0 0 60720 0 349 315380 0 0 0 60720 0 350 319369 0 1595 17 60720 10 351 318690 0 312 14 58598 3 352 314049 0 60 10 61124 7 353 325699 0 587 12 59595 10 354 314210 0 135 2 62065 1 355 315380 0 0 0 60720 0 356 315380 0 0 0 60720 0 357 322378 0 514 52 78780 15 358 315380 0 0 0 60720 0 359 315380 0 0 0 60720 0 360 315380 0 0 0 60720 0 361 315398 0 1 4 60722 4 362 315380 0 0 0 60720 0 363 315380 0 0 0 60720 0 364 308336 1 1763 24 61600 28 365 316386 0 180 11 59635 9 366 315380 0 0 0 60720 0 367 315380 0 0 0 60720 0 368 315380 0 0 0 60720 0 369 315380 0 0 0 60720 0 370 315553 0 218 21 60720 7 371 315380 0 0 0 60720 0 372 323361 0 448 40 59781 7 373 336639 0 227 9 76644 7 374 307424 0 174 1 64820 3 375 315380 0 0 0 60720 0 376 315380 0 0 0 60720 0 377 295370 0 121 24 56178 11 378 322340 0 607 11 60436 7 379 319864 0 2212 14 60720 10 380 315380 0 0 0 60720 0 381 315380 0 0 0 60720 0 382 317291 0 530 60 73433 18 383 280398 0 571 80 41477 14 384 315380 0 0 0 60720 0 385 317330 0 78 16 62700 12 386 238125 0 2489 40 67804 29 387 327071 0 131 6 59661 3 388 309038 0 923 8 58620 6 389 314210 0 72 3 60398 3 390 307930 0 572 16 58580 8 391 322327 0 397 10 62710 10 392 292136 0 450 8 59325 6 393 263276 0 622 7 60950 8 394 367655 0 694 8 68060 6 395 283910 1 3425 12 83620 9 396 283587 0 562 13 58456 8 397 243650 0 4917 42 52811 26 398 438493 1 1442 118 121173 239 399 296261 0 529 9 63870 7 400 230621 1 2126 138 21001 41 401 304252 0 1061 5 70415 3 402 333505 0 776 9 64230 8 403 296919 0 611 8 59190 6 404 278990 1 1526 25 69351 21 405 276898 0 592 7 64270 7 406 327007 0 1182 13 70694 11 407 317046 0 621 16 68005 11 408 304555 0 989 11 58930 12 409 298096 0 438 11 58320 9 410 231861 0 726 3 69980 3 411 309422 0 1303 61 69863 57 412 286963 1 7419 29 63255 21 413 269753 1 1164 17 57320 15 414 448243 1 3310 33 75230 32 415 165404 0 1920 15 79420 11 416 204325 0 965 3 73490 2 417 407159 0 3256 66 35250 23 418 290476 1 1135 17 62285 20 419 275311 0 1270 26 69206 24 420 246541 0 661 3 65920 1 421 253468 0 1013 2 69770 1 422 240897 0 2844 67 72683 74 423 -83265 1 11528 70 -14545 68 424 -42143 0 6526 26 55830 20 425 272713 0 2264 24 55174 20 426 215362 1 5109 97 67038 82 427 42754 0 3999 30 51252 21 428 306275 0 35624 223 157278 244 429 253537 0 9252 48 79510 32 430 372631 0 15236 90 77440 86 431 -7170 0 18073 180 27284 69 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Group Costs Trades Dividends 155435.184 37009.233 21.738 -1157.434 1.669 `Orders\r\r` 2499.590 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -2120094 -11235 34448 58624 2592406 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 1.554e+05 2.986e+04 5.205 3.03e-07 *** Group 3.701e+04 2.706e+04 1.368 0.172 Costs 2.174e+01 2.079e+00 10.458 < 2e-16 *** Trades -1.157e+03 2.075e+02 -5.579 4.31e-08 *** Dividends 1.669e+00 3.987e-01 4.185 3.46e-05 *** `Orders\r\r` 2.500e+03 4.380e+02 5.707 2.16e-08 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 244500 on 425 degrees of freedom Multiple R-squared: 0.6813, Adjusted R-squared: 0.6776 F-statistic: 181.7 on 5 and 425 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 9.974460e-44 4.987230e-44 [2,] 1.0000000 1.331465e-56 6.657326e-57 [3,] 1.0000000 3.386240e-58 1.693120e-58 [4,] 1.0000000 2.143603e-69 1.071801e-69 [5,] 1.0000000 3.132541e-95 1.566271e-95 [6,] 1.0000000 6.821950e-97 3.410975e-97 [7,] 1.0000000 4.695350e-100 2.347675e-100 [8,] 1.0000000 5.381015e-112 2.690507e-112 [9,] 1.0000000 2.505376e-119 1.252688e-119 [10,] 1.0000000 7.021203e-122 3.510601e-122 [11,] 1.0000000 9.838501e-127 4.919250e-127 [12,] 1.0000000 1.197743e-128 5.988716e-129 [13,] 1.0000000 2.719007e-134 1.359503e-134 [14,] 1.0000000 9.146647e-136 4.573323e-136 [15,] 1.0000000 3.136624e-136 1.568312e-136 [16,] 1.0000000 9.139798e-140 4.569899e-140 [17,] 1.0000000 2.180360e-139 1.090180e-139 [18,] 1.0000000 8.671937e-144 4.335969e-144 [19,] 1.0000000 7.624892e-149 3.812446e-149 [20,] 1.0000000 2.918407e-151 1.459204e-151 [21,] 1.0000000 3.739420e-167 1.869710e-167 [22,] 1.0000000 2.372397e-170 1.186199e-170 [23,] 1.0000000 2.406106e-174 1.203053e-174 [24,] 1.0000000 6.504338e-176 3.252169e-176 [25,] 1.0000000 8.757654e-184 4.378827e-184 [26,] 1.0000000 1.027542e-185 5.137711e-186 [27,] 1.0000000 6.846050e-185 3.423025e-185 [28,] 1.0000000 5.178241e-186 2.589120e-186 [29,] 1.0000000 2.884381e-192 1.442190e-192 [30,] 1.0000000 1.097182e-197 5.485912e-198 [31,] 1.0000000 3.130504e-198 1.565252e-198 [32,] 1.0000000 1.165698e-198 5.828489e-199 [33,] 1.0000000 9.045465e-201 4.522733e-201 [34,] 1.0000000 4.955653e-200 2.477827e-200 [35,] 1.0000000 1.752112e-201 8.760562e-202 [36,] 1.0000000 2.088923e-200 1.044461e-200 [37,] 1.0000000 1.119034e-200 5.595171e-201 [38,] 1.0000000 8.642654e-200 4.321327e-200 [39,] 1.0000000 8.808475e-201 4.404237e-201 [40,] 1.0000000 2.117767e-200 1.058883e-200 [41,] 1.0000000 4.560586e-200 2.280293e-200 [42,] 1.0000000 1.971335e-200 9.856676e-201 [43,] 1.0000000 6.608758e-201 3.304379e-201 [44,] 1.0000000 3.746165e-201 1.873082e-201 [45,] 1.0000000 2.295774e-201 1.147887e-201 [46,] 1.0000000 7.235020e-201 3.617510e-201 [47,] 1.0000000 5.084592e-205 2.542296e-205 [48,] 1.0000000 5.164669e-205 2.582334e-205 [49,] 1.0000000 3.872729e-219 1.936364e-219 [50,] 1.0000000 9.255818e-220 4.627909e-220 [51,] 1.0000000 5.048187e-219 2.524093e-219 [52,] 1.0000000 7.279062e-219 3.639531e-219 [53,] 1.0000000 2.317114e-220 1.158557e-220 [54,] 1.0000000 6.229618e-222 3.114809e-222 [55,] 1.0000000 2.342625e-221 1.171313e-221 [56,] 1.0000000 7.841021e-221 3.920511e-221 [57,] 1.0000000 2.544665e-220 1.272332e-220 [58,] 1.0000000 3.602077e-220 1.801039e-220 [59,] 1.0000000 5.512893e-220 2.756447e-220 [60,] 1.0000000 5.420487e-219 2.710244e-219 [61,] 1.0000000 2.974699e-218 1.487350e-218 [62,] 1.0000000 1.072594e-217 5.362970e-218 [63,] 1.0000000 7.797233e-217 3.898617e-217 [64,] 1.0000000 3.990194e-218 1.995097e-218 [65,] 1.0000000 1.195943e-217 5.979713e-218 [66,] 1.0000000 6.821109e-217 3.410555e-217 [67,] 1.0000000 6.917766e-216 3.458883e-216 [68,] 1.0000000 1.573713e-215 7.868564e-216 [69,] 1.0000000 1.125573e-214 5.627867e-215 [70,] 1.0000000 1.098403e-214 5.492016e-215 [71,] 1.0000000 7.561085e-214 3.780543e-214 [72,] 1.0000000 2.294503e-213 1.147251e-213 [73,] 1.0000000 1.045789e-212 5.228945e-213 [74,] 1.0000000 1.574856e-212 7.874280e-213 [75,] 1.0000000 1.314312e-212 6.571559e-213 [76,] 1.0000000 5.045541e-212 2.522770e-212 [77,] 1.0000000 3.980773e-211 1.990387e-211 [78,] 1.0000000 8.482941e-211 4.241471e-211 [79,] 1.0000000 6.368253e-211 3.184127e-211 [80,] 1.0000000 3.933560e-210 1.966780e-210 [81,] 1.0000000 7.409978e-210 3.704989e-210 [82,] 1.0000000 6.023619e-209 3.011809e-209 [83,] 1.0000000 4.661104e-208 2.330552e-208 [84,] 1.0000000 2.957481e-207 1.478741e-207 [85,] 1.0000000 8.978947e-209 4.489473e-209 [86,] 1.0000000 7.403393e-208 3.701696e-208 [87,] 1.0000000 2.091849e-208 1.045925e-208 [88,] 1.0000000 2.033519e-208 1.016760e-208 [89,] 1.0000000 1.524491e-207 7.622456e-208 [90,] 1.0000000 9.484618e-207 4.742309e-207 [91,] 1.0000000 7.294933e-206 3.647466e-206 [92,] 1.0000000 7.830817e-205 3.915408e-205 [93,] 1.0000000 6.983084e-204 3.491542e-204 [94,] 1.0000000 4.854181e-205 2.427091e-205 [95,] 1.0000000 3.645952e-205 1.822976e-205 [96,] 1.0000000 4.109047e-206 2.054524e-206 [97,] 1.0000000 1.530656e-205 7.653278e-206 [98,] 1.0000000 5.007308e-205 2.503654e-205 [99,] 1.0000000 2.458930e-207 1.229465e-207 [100,] 1.0000000 1.357079e-208 6.785396e-209 [101,] 1.0000000 9.411013e-208 4.705506e-208 [102,] 1.0000000 1.782732e-207 8.913662e-208 [103,] 1.0000000 1.030704e-207 5.153522e-208 [104,] 1.0000000 5.776787e-208 2.888393e-208 [105,] 1.0000000 4.026076e-207 2.013038e-207 [106,] 1.0000000 3.190264e-206 1.595132e-206 [107,] 1.0000000 2.355389e-206 1.177695e-206 [108,] 1.0000000 7.820737e-206 3.910368e-206 [109,] 1.0000000 4.127436e-205 2.063718e-205 [110,] 1.0000000 3.031649e-205 1.515825e-205 [111,] 1.0000000 2.678520e-206 1.339260e-206 [112,] 1.0000000 1.987369e-205 9.936844e-206 [113,] 1.0000000 8.449207e-205 4.224604e-205 [114,] 1.0000000 6.549712e-204 3.274856e-204 [115,] 1.0000000 5.085971e-203 2.542986e-203 [116,] 1.0000000 9.430968e-203 4.715484e-203 [117,] 1.0000000 2.928238e-202 1.464119e-202 [118,] 1.0000000 7.892139e-203 3.946069e-203 [119,] 1.0000000 5.908090e-202 2.954045e-202 [120,] 1.0000000 4.937940e-201 2.468970e-201 [121,] 1.0000000 4.008111e-200 2.004056e-200 [122,] 1.0000000 3.098253e-199 1.549126e-199 [123,] 1.0000000 2.551252e-198 1.275626e-198 [124,] 1.0000000 1.883034e-198 9.415169e-199 [125,] 1.0000000 1.511286e-197 7.556430e-198 [126,] 1.0000000 1.459092e-196 7.295462e-197 [127,] 1.0000000 1.244129e-195 6.220647e-196 [128,] 1.0000000 1.025973e-194 5.129864e-195 [129,] 1.0000000 8.684065e-194 4.342033e-194 [130,] 1.0000000 7.199991e-193 3.599995e-193 [131,] 1.0000000 6.026371e-192 3.013186e-192 [132,] 1.0000000 5.009055e-191 2.504527e-191 [133,] 1.0000000 3.584103e-190 1.792052e-190 [134,] 1.0000000 3.003007e-189 1.501503e-189 [135,] 1.0000000 2.759688e-188 1.379844e-188 [136,] 1.0000000 2.376722e-187 1.188361e-187 [137,] 1.0000000 1.965088e-186 9.825441e-187 [138,] 1.0000000 1.600984e-185 8.004922e-186 [139,] 1.0000000 1.351842e-184 6.759209e-185 [140,] 1.0000000 1.188672e-183 5.943361e-184 [141,] 1.0000000 9.092506e-183 4.546253e-183 [142,] 1.0000000 7.693173e-182 3.846587e-182 [143,] 1.0000000 6.892251e-181 3.446126e-181 [144,] 1.0000000 6.657374e-180 3.328687e-180 [145,] 1.0000000 5.807407e-179 2.903703e-179 [146,] 1.0000000 5.215368e-178 2.607684e-178 [147,] 1.0000000 4.442202e-177 2.221101e-177 [148,] 1.0000000 3.745329e-176 1.872664e-176 [149,] 1.0000000 3.156819e-175 1.578409e-175 [150,] 1.0000000 2.708358e-174 1.354179e-174 [151,] 1.0000000 1.697324e-173 8.486618e-174 [152,] 1.0000000 1.428335e-172 7.141675e-173 [153,] 1.0000000 1.200777e-171 6.003887e-172 [154,] 1.0000000 1.008314e-170 5.041571e-171 [155,] 1.0000000 8.456060e-170 4.228030e-170 [156,] 1.0000000 7.225614e-169 3.612807e-169 [157,] 1.0000000 6.065334e-168 3.032667e-168 [158,] 1.0000000 5.071152e-167 2.535576e-167 [159,] 1.0000000 4.222341e-166 2.111170e-166 [160,] 1.0000000 4.007606e-167 2.003803e-167 [161,] 1.0000000 4.989211e-167 2.494606e-167 [162,] 1.0000000 4.257436e-166 2.128718e-166 [163,] 1.0000000 3.624777e-165 1.812388e-165 [164,] 1.0000000 3.156103e-164 1.578051e-164 [165,] 1.0000000 2.673803e-163 1.336902e-163 [166,] 1.0000000 2.368311e-162 1.184155e-162 [167,] 1.0000000 1.995139e-161 9.975693e-162 [168,] 1.0000000 1.705029e-160 8.525144e-161 [169,] 1.0000000 1.502878e-159 7.514390e-160 [170,] 1.0000000 1.337343e-158 6.686715e-159 [171,] 1.0000000 1.110310e-157 5.551551e-158 [172,] 1.0000000 9.192717e-157 4.596359e-157 [173,] 1.0000000 5.503871e-156 2.751935e-156 [174,] 1.0000000 4.889033e-155 2.444516e-155 [175,] 1.0000000 2.485759e-154 1.242879e-154 [176,] 1.0000000 1.434550e-153 7.172752e-154 [177,] 1.0000000 1.257698e-152 6.288492e-153 [178,] 1.0000000 8.765372e-152 4.382686e-152 [179,] 1.0000000 3.445473e-151 1.722737e-151 [180,] 1.0000000 2.881730e-150 1.440865e-150 [181,] 1.0000000 1.943279e-149 9.716393e-150 [182,] 1.0000000 1.563303e-148 7.816515e-149 [183,] 1.0000000 1.205935e-147 6.029674e-148 [184,] 1.0000000 9.638910e-147 4.819455e-147 [185,] 1.0000000 7.863879e-146 3.931939e-146 [186,] 1.0000000 6.243712e-145 3.121856e-145 [187,] 1.0000000 4.979794e-144 2.489897e-144 [188,] 1.0000000 4.150601e-143 2.075301e-143 [189,] 1.0000000 2.513182e-142 1.256591e-142 [190,] 1.0000000 2.075130e-141 1.037565e-141 [191,] 1.0000000 1.705968e-140 8.529840e-141 [192,] 1.0000000 1.459537e-140 7.297686e-141 [193,] 1.0000000 1.204166e-139 6.020828e-140 [194,] 1.0000000 9.891965e-139 4.945982e-139 [195,] 1.0000000 8.090767e-138 4.045383e-138 [196,] 1.0000000 6.161996e-137 3.080998e-137 [197,] 1.0000000 2.329648e-136 1.164824e-136 [198,] 1.0000000 1.899446e-135 9.497231e-136 [199,] 1.0000000 1.512641e-134 7.563207e-135 [200,] 1.0000000 1.215832e-133 6.079160e-134 [201,] 1.0000000 2.010215e-133 1.005107e-133 [202,] 1.0000000 1.597931e-132 7.989654e-133 [203,] 1.0000000 1.208301e-131 6.041506e-132 [204,] 1.0000000 9.105597e-131 4.552798e-131 [205,] 1.0000000 4.911331e-130 2.455665e-130 [206,] 1.0000000 3.507629e-129 1.753814e-129 [207,] 1.0000000 2.710332e-128 1.355166e-128 [208,] 1.0000000 2.028478e-127 1.014239e-127 [209,] 1.0000000 1.512727e-126 7.563634e-127 [210,] 1.0000000 1.061814e-125 5.309071e-126 [211,] 1.0000000 7.859702e-125 3.929851e-125 [212,] 1.0000000 3.491157e-124 1.745578e-124 [213,] 1.0000000 2.597173e-123 1.298586e-123 [214,] 1.0000000 1.893315e-122 9.466576e-123 [215,] 1.0000000 1.374892e-121 6.874462e-122 [216,] 1.0000000 1.014911e-120 5.074555e-121 [217,] 1.0000000 7.309882e-120 3.654941e-120 [218,] 1.0000000 5.350771e-119 2.675385e-119 [219,] 1.0000000 3.852648e-118 1.926324e-118 [220,] 1.0000000 2.742208e-117 1.371104e-117 [221,] 1.0000000 1.979850e-116 9.899251e-117 [222,] 1.0000000 1.398480e-115 6.992399e-116 [223,] 1.0000000 9.665062e-115 4.832531e-115 [224,] 1.0000000 6.226271e-114 3.113136e-114 [225,] 1.0000000 4.347513e-113 2.173757e-113 [226,] 1.0000000 3.133984e-112 1.566992e-112 [227,] 1.0000000 2.170497e-111 1.085248e-111 [228,] 1.0000000 1.495137e-110 7.475685e-111 [229,] 1.0000000 1.027087e-109 5.135435e-110 [230,] 1.0000000 7.003932e-109 3.501966e-109 [231,] 1.0000000 4.771286e-108 2.385643e-108 [232,] 1.0000000 3.236830e-107 1.618415e-107 [233,] 1.0000000 2.186676e-106 1.093338e-106 [234,] 1.0000000 1.497211e-105 7.486054e-106 [235,] 1.0000000 1.002908e-104 5.014542e-105 [236,] 1.0000000 6.620763e-104 3.310382e-104 [237,] 1.0000000 4.536778e-103 2.268389e-103 [238,] 1.0000000 3.113333e-102 1.556667e-102 [239,] 1.0000000 2.126388e-101 1.063194e-101 [240,] 1.0000000 1.445377e-100 7.226885e-101 [241,] 1.0000000 9.777797e-100 4.888899e-100 [242,] 1.0000000 6.582893e-99 3.291447e-99 [243,] 1.0000000 4.223707e-98 2.111854e-98 [244,] 1.0000000 2.817123e-97 1.408562e-97 [245,] 1.0000000 1.828819e-96 9.144094e-97 [246,] 1.0000000 9.954439e-96 4.977220e-96 [247,] 1.0000000 7.311556e-96 3.655778e-96 [248,] 1.0000000 4.678196e-95 2.339098e-95 [249,] 1.0000000 2.698555e-94 1.349278e-94 [250,] 1.0000000 1.411163e-93 7.055814e-94 [251,] 1.0000000 9.244191e-93 4.622095e-93 [252,] 1.0000000 5.811616e-92 2.905808e-92 [253,] 1.0000000 3.770383e-91 1.885192e-91 [254,] 1.0000000 2.433789e-90 1.216895e-90 [255,] 1.0000000 1.544276e-89 7.721379e-90 [256,] 1.0000000 7.471609e-89 3.735804e-89 [257,] 1.0000000 4.760559e-88 2.380280e-88 [258,] 1.0000000 2.912610e-87 1.456305e-87 [259,] 1.0000000 1.837527e-86 9.187635e-87 [260,] 1.0000000 9.263582e-86 4.631791e-86 [261,] 1.0000000 5.792544e-85 2.896272e-85 [262,] 1.0000000 3.618307e-84 1.809153e-84 [263,] 1.0000000 2.239391e-83 1.119696e-83 [264,] 1.0000000 1.378628e-82 6.893139e-83 [265,] 1.0000000 8.441771e-82 4.220886e-82 [266,] 1.0000000 5.127502e-81 2.563751e-81 [267,] 1.0000000 3.006248e-80 1.503124e-80 [268,] 1.0000000 1.754672e-79 8.773361e-80 [269,] 1.0000000 1.044374e-78 5.221868e-79 [270,] 1.0000000 6.035695e-78 3.017848e-78 [271,] 1.0000000 2.551655e-77 1.275827e-77 [272,] 1.0000000 1.404023e-76 7.020114e-77 [273,] 1.0000000 8.301931e-76 4.150966e-76 [274,] 1.0000000 4.707267e-75 2.353633e-75 [275,] 1.0000000 2.614894e-74 1.307447e-74 [276,] 1.0000000 1.468868e-73 7.344341e-74 [277,] 1.0000000 8.481277e-73 4.240638e-73 [278,] 1.0000000 4.717298e-72 2.358649e-72 [279,] 1.0000000 2.621144e-71 1.310572e-71 [280,] 1.0000000 1.452902e-70 7.264509e-71 [281,] 1.0000000 8.241609e-70 4.120804e-70 [282,] 1.0000000 4.548347e-69 2.274173e-69 [283,] 1.0000000 2.552170e-68 1.276085e-68 [284,] 1.0000000 1.363193e-67 6.815967e-68 [285,] 1.0000000 7.507074e-67 3.753537e-67 [286,] 1.0000000 4.010445e-66 2.005222e-66 [287,] 1.0000000 2.203106e-65 1.101553e-65 [288,] 1.0000000 9.412203e-65 4.706101e-65 [289,] 1.0000000 4.960013e-64 2.480006e-64 [290,] 1.0000000 2.675802e-63 1.337901e-63 [291,] 1.0000000 1.395108e-62 6.975539e-63 [292,] 1.0000000 7.247511e-62 3.623756e-62 [293,] 1.0000000 3.867825e-61 1.933913e-61 [294,] 1.0000000 2.035083e-60 1.017542e-60 [295,] 1.0000000 1.038378e-59 5.191889e-60 [296,] 1.0000000 5.413444e-59 2.706722e-59 [297,] 1.0000000 2.245641e-58 1.122820e-58 [298,] 1.0000000 1.136489e-57 5.682447e-58 [299,] 1.0000000 5.833760e-57 2.916880e-57 [300,] 1.0000000 3.003672e-56 1.501836e-56 [301,] 1.0000000 1.363325e-55 6.816627e-56 [302,] 1.0000000 6.877006e-55 3.438503e-55 [303,] 1.0000000 3.442900e-54 1.721450e-54 [304,] 1.0000000 1.275543e-53 6.377713e-54 [305,] 1.0000000 6.312187e-53 3.156093e-53 [306,] 1.0000000 3.122885e-52 1.561443e-52 [307,] 1.0000000 1.498254e-51 7.491270e-52 [308,] 1.0000000 6.826651e-51 3.413325e-51 [309,] 1.0000000 3.201368e-50 1.600684e-50 [310,] 1.0000000 1.563448e-49 7.817241e-50 [311,] 1.0000000 7.346504e-49 3.673252e-49 [312,] 1.0000000 3.413845e-48 1.706922e-48 [313,] 1.0000000 1.576562e-47 7.882809e-48 [314,] 1.0000000 7.235355e-47 3.617677e-47 [315,] 1.0000000 3.299644e-46 1.649822e-46 [316,] 1.0000000 1.522859e-45 7.614296e-46 [317,] 1.0000000 6.859128e-45 3.429564e-45 [318,] 1.0000000 3.109828e-44 1.554914e-44 [319,] 1.0000000 1.402465e-43 7.012323e-44 [320,] 1.0000000 6.245412e-43 3.122706e-43 [321,] 1.0000000 2.742599e-42 1.371300e-42 [322,] 1.0000000 1.196384e-41 5.981920e-42 [323,] 1.0000000 5.183925e-41 2.591963e-41 [324,] 1.0000000 2.230996e-40 1.115498e-40 [325,] 1.0000000 9.535908e-40 4.767954e-40 [326,] 1.0000000 4.065385e-39 2.032692e-39 [327,] 1.0000000 1.717402e-38 8.587008e-39 [328,] 1.0000000 7.129937e-38 3.564968e-38 [329,] 1.0000000 2.953427e-37 1.476714e-37 [330,] 1.0000000 1.219889e-36 6.099447e-37 [331,] 1.0000000 5.047328e-36 2.523664e-36 [332,] 1.0000000 2.021038e-35 1.010519e-35 [333,] 1.0000000 8.471949e-35 4.235974e-35 [334,] 1.0000000 3.402743e-34 1.701372e-34 [335,] 1.0000000 1.356422e-33 6.782111e-34 [336,] 1.0000000 5.365894e-33 2.682947e-33 [337,] 1.0000000 2.193437e-32 1.096718e-32 [338,] 1.0000000 8.551125e-32 4.275562e-32 [339,] 1.0000000 3.404075e-31 1.702038e-31 [340,] 1.0000000 1.306558e-30 6.532789e-31 [341,] 1.0000000 4.974531e-30 2.487266e-30 [342,] 1.0000000 1.711571e-29 8.557853e-30 [343,] 1.0000000 6.135520e-29 3.067760e-29 [344,] 1.0000000 2.326670e-28 1.163335e-28 [345,] 1.0000000 7.694758e-28 3.847379e-28 [346,] 1.0000000 2.869614e-27 1.434807e-27 [347,] 1.0000000 1.040071e-26 5.200357e-27 [348,] 1.0000000 3.736351e-26 1.868176e-26 [349,] 1.0000000 1.187715e-25 5.938573e-26 [350,] 1.0000000 4.208751e-25 2.104375e-25 [351,] 1.0000000 1.477710e-24 7.388550e-25 [352,] 1.0000000 5.139890e-24 2.569945e-24 [353,] 1.0000000 1.779826e-23 8.899129e-24 [354,] 1.0000000 6.073455e-23 3.036728e-23 [355,] 1.0000000 2.052138e-22 1.026069e-22 [356,] 1.0000000 6.862306e-22 3.431153e-22 [357,] 1.0000000 2.247030e-21 1.123515e-21 [358,] 1.0000000 7.382186e-21 3.691093e-21 [359,] 1.0000000 2.399845e-20 1.199923e-20 [360,] 1.0000000 7.718034e-20 3.859017e-20 [361,] 1.0000000 2.455004e-19 1.227502e-19 [362,] 1.0000000 8.005635e-19 4.002818e-19 [363,] 1.0000000 2.492720e-18 1.246360e-18 [364,] 1.0000000 7.878511e-18 3.939255e-18 [365,] 1.0000000 2.606939e-17 1.303470e-17 [366,] 1.0000000 8.430119e-17 4.215060e-17 [367,] 1.0000000 2.519399e-16 1.259699e-16 [368,] 1.0000000 7.435345e-16 3.717673e-16 [369,] 1.0000000 2.307569e-15 1.153784e-15 [370,] 1.0000000 6.295799e-15 3.147899e-15 [371,] 1.0000000 1.545826e-14 7.729130e-15 [372,] 1.0000000 4.333658e-14 2.166829e-14 [373,] 1.0000000 1.197359e-13 5.986794e-14 [374,] 1.0000000 3.415083e-13 1.707542e-13 [375,] 1.0000000 9.950943e-13 4.975472e-13 [376,] 1.0000000 2.663265e-12 1.331632e-12 [377,] 1.0000000 7.381387e-12 3.690694e-12 [378,] 1.0000000 1.816267e-11 9.081333e-12 [379,] 1.0000000 4.355462e-11 2.177731e-11 [380,] 1.0000000 1.071530e-10 5.357651e-11 [381,] 1.0000000 2.693241e-10 1.346621e-10 [382,] 1.0000000 6.772349e-10 3.386174e-10 [383,] 1.0000000 1.632293e-09 8.161463e-10 [384,] 1.0000000 4.215518e-09 2.107759e-09 [385,] 1.0000000 1.140337e-08 5.701687e-09 [386,] 1.0000000 2.093450e-08 1.046725e-08 [387,] 1.0000000 4.762370e-08 2.381185e-08 [388,] 0.9999999 1.191779e-07 5.958896e-08 [389,] 0.9999999 2.814458e-07 1.407229e-07 [390,] 0.9999997 5.818837e-07 2.909419e-07 [391,] 0.9999993 1.409663e-06 7.048315e-07 [392,] 0.9999982 3.501477e-06 1.750739e-06 [393,] 0.9999959 8.283890e-06 4.141945e-06 [394,] 0.9999920 1.593713e-05 7.968567e-06 [395,] 0.9999834 3.318217e-05 1.659109e-05 [396,] 0.9999642 7.165457e-05 3.582729e-05 [397,] 0.9999197 1.605918e-04 8.029588e-05 [398,] 0.9998380 3.240952e-04 1.620476e-04 [399,] 0.9996707 6.585617e-04 3.292808e-04 [400,] 0.9994272 1.145610e-03 5.728050e-04 [401,] 0.9989885 2.022990e-03 1.011495e-03 [402,] 0.9979016 4.196783e-03 2.098392e-03 [403,] 0.9958308 8.338443e-03 4.169221e-03 [404,] 0.9922712 1.545765e-02 7.728824e-03 [405,] 0.9852129 2.957411e-02 1.478706e-02 [406,] 0.9839799 3.204016e-02 1.602008e-02 [407,] 0.9756373 4.872541e-02 2.436271e-02 [408,] 0.9573504 8.529912e-02 4.264956e-02 [409,] 0.9837381 3.252381e-02 1.626191e-02 [410,] 0.9667559 6.648817e-02 3.324409e-02 [411,] 0.9333787 1.332427e-01 6.662133e-02 [412,] 0.8719773 2.560454e-01 1.280227e-01 [413,] 0.7708062 4.583876e-01 2.291938e-01 [414,] 0.6147456 7.705087e-01 3.852544e-01 > postscript(file="/var/www/html/freestat/rcomp/tmp/1qe781291224108.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/2in6a1291224108.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/3in6a1291224108.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/4in6a1291224108.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/5bwov1291224108.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 1.434540e+06 2.592406e+06 1.260986e+06 -2.120094e+06 -5.389179e+05 6 7 8 9 10 8.706658e+03 3.681755e+04 1.011745e+05 -1.759896e+05 6.212848e+05 11 12 13 14 15 -1.804099e+05 2.960787e+05 1.183419e+05 -1.118617e+06 -7.714344e+04 16 17 18 19 20 5.469366e+05 -2.810947e+05 -6.509910e+04 -3.783602e+05 -5.639902e+05 21 22 23 24 25 2.003133e+05 -1.692827e+05 -1.583710e+05 -7.263975e+05 -3.223750e+05 26 27 28 29 30 -8.611960e+04 -5.244867e+05 -2.287674e+05 3.491694e+05 -4.180886e+04 31 32 33 34 35 2.450783e+05 -1.286764e+05 2.148073e+05 8.843046e+04 -6.741604e+04 36 37 38 39 40 -6.997293e+05 3.290698e+04 9.512216e+04 -1.292850e+05 -3.292756e+05 41 42 43 44 45 2.265576e+05 -3.381594e+05 1.305302e+05 -2.193006e+05 -9.686861e+03 46 47 48 49 50 -9.873843e+04 -1.506070e+05 -1.527146e+05 -4.405757e+04 -2.023685e+05 51 52 53 54 55 -1.079644e+04 -1.554851e+05 -2.840178e+05 9.597283e+04 -4.233277e+05 56 57 58 59 60 1.614878e+04 3.274203e+05 -1.344067e+04 -1.167402e+04 -2.489276e+05 61 62 63 64 65 -3.053437e+04 -1.788236e+05 -2.695521e+04 -2.218082e+05 -4.244611e+02 66 67 68 69 70 2.313680e+05 -3.961293e+00 -5.613649e+04 -1.039046e+05 -1.621399e+05 71 72 73 74 75 -1.575332e+05 8.402731e+03 -1.022603e+05 1.619668e+04 -4.670361e+04 76 77 78 79 80 -1.401200e+05 -1.895204e+04 -3.751264e+04 1.864700e+04 1.016586e+05 81 82 83 84 85 1.779951e+04 -9.639347e+04 -6.771926e+04 4.597112e+04 1.936319e+04 86 87 88 89 90 -1.526918e+05 -7.729575e+04 -2.709022e+03 -1.087286e+05 2.374239e+04 91 92 93 94 95 5.171042e+03 -6.890441e+04 8.655845e+04 -1.670658e+05 4.413430e+04 96 97 98 99 100 -7.960861e+03 6.597813e+04 1.644330e+05 -4.679715e+04 -3.656027e+04 101 102 103 104 105 -7.862713e+04 -1.134372e+05 -1.552571e+05 -4.287334e+04 3.778051e+04 106 107 108 109 110 -6.726786e+04 -2.868120e+05 9.374730e+04 -1.737003e+05 -1.146023e+05 111 112 113 114 115 9.946503e+04 -2.168834e+05 -1.463099e+05 -9.484045e+04 5.483301e+04 116 117 118 119 120 1.896172e+04 3.684722e+04 7.276375e+04 -1.645875e+05 6.720513e+04 121 122 123 124 125 -3.177472e+05 3.423384e+04 -1.454836e+05 4.344236e+04 -1.585187e+03 126 127 128 129 130 -2.042616e+05 5.862432e+04 3.871832e+04 3.169185e+04 5.862432e+04 131 132 133 134 135 4.877076e+04 9.773557e+04 5.862432e+04 8.357447e+03 3.708742e+04 136 137 138 139 140 5.216156e+04 4.932280e+04 5.783791e+04 3.444756e+04 5.862432e+04 141 142 143 144 145 5.938854e+04 5.862432e+04 1.652187e+04 5.147936e+04 3.851483e+04 146 147 148 149 150 4.662419e+04 5.862432e+04 4.729126e+04 5.991813e+04 5.948051e+04 151 152 153 154 155 2.112212e+04 -1.920308e+04 4.612637e+04 2.978399e+04 5.616392e+04 156 157 158 159 160 5.862432e+04 5.862432e+04 5.238087e+04 6.626112e+04 5.862432e+04 161 162 163 164 165 5.862432e+04 5.862432e+04 5.862432e+04 5.258312e+04 3.036140e+04 166 167 168 169 170 5.807737e+04 5.862432e+04 1.115811e+05 3.010201e+05 5.862432e+04 171 172 173 174 175 5.862432e+04 5.059639e+04 5.862432e+04 3.431578e+04 5.862432e+04 176 177 178 179 180 4.724455e+04 -5.870740e+04 2.198344e+04 5.862432e+04 5.862432e+04 181 182 183 184 185 5.508597e+04 -9.384035e+03 1.574433e+04 -2.986051e+03 1.076837e+04 186 187 188 189 190 1.486915e+04 -2.250814e+04 4.192277e+04 1.502288e+04 5.862432e+04 191 192 193 194 195 4.517986e+04 5.862432e+04 4.347433e+04 5.862432e+04 5.285609e+04 196 197 198 199 200 2.161509e+04 -1.976193e+05 2.161509e+04 2.161509e+04 -9.102825e+04 201 202 203 204 205 2.161509e+04 2.161509e+04 2.161509e+04 2.058634e+04 -3.388936e+04 206 207 208 209 210 1.042777e+04 1.478032e+04 2.161509e+04 1.755700e+04 3.580431e+04 211 212 213 214 215 5.862432e+04 5.862432e+04 8.218123e+04 1.798101e+04 2.855512e+04 216 217 218 219 220 5.362514e+04 5.268297e+04 6.738501e+04 5.362514e+04 -3.909274e+03 221 222 223 224 225 3.112129e+04 5.862432e+04 5.838509e+04 1.954647e+04 5.805901e+04 226 227 228 229 230 2.608423e+04 5.017709e+04 5.862432e+04 4.343789e+04 5.862432e+04 231 232 233 234 235 2.264195e+04 3.765605e+04 5.862432e+04 2.404469e+04 5.862432e+04 236 237 238 239 240 5.980690e+04 5.862432e+04 5.045136e+04 5.862432e+04 5.862432e+04 241 242 243 244 245 5.862432e+04 4.899359e+04 5.862432e+04 5.162708e+04 -3.389133e+03 246 247 248 249 250 2.161509e+04 2.161509e+04 2.280813e+04 2.161509e+04 2.161509e+04 251 252 253 254 255 3.608742e+04 2.161509e+04 4.086943e+04 3.105108e+03 -1.241809e+05 256 257 258 259 260 5.862432e+04 4.647750e+04 3.103620e+04 2.161509e+04 5.862432e+04 261 262 263 264 265 2.161509e+04 2.161509e+04 4.343946e+04 6.400187e+04 2.161509e+04 266 267 268 269 270 5.862432e+04 2.161509e+04 8.231452e+04 2.161509e+04 -1.559911e+03 271 272 273 274 275 2.161509e+04 2.161509e+04 2.161509e+04 1.699497e+04 5.862432e+04 276 277 278 279 280 5.896816e+04 3.129980e+04 5.862432e+04 2.101390e+04 5.012522e+04 281 282 283 284 285 2.161509e+04 5.862432e+04 7.633316e+04 5.862432e+04 4.304069e+04 286 287 288 289 290 5.862432e+04 6.827479e+03 5.561020e+04 2.161509e+04 8.052626e+04 291 292 293 294 295 2.161509e+04 5.907775e+04 4.212129e+04 5.862432e+04 3.773515e+04 296 297 298 299 300 4.297017e+04 5.862432e+04 2.161509e+04 5.862432e+04 3.612801e+04 301 302 303 304 305 1.324818e+04 4.114021e+04 5.612473e+04 2.161509e+04 1.042476e+05 306 307 308 309 310 3.899164e+04 2.161509e+04 -1.960650e+04 3.777727e+04 2.161509e+04 311 312 313 314 315 2.161509e+04 5.382871e+04 2.161509e+04 -4.185561e+03 4.339971e+04 316 317 318 319 320 7.721015e+04 5.989160e+04 1.324597e+03 5.304290e+04 5.862432e+04 321 322 323 324 325 5.862432e+04 5.862432e+04 5.862432e+04 2.161509e+04 5.862432e+04 326 327 328 329 330 5.508003e+04 2.161509e+04 2.161509e+04 5.862432e+04 5.862432e+04 331 332 333 334 335 5.862432e+04 5.862432e+04 5.862432e+04 5.124433e+04 5.841418e+04 336 337 338 339 340 4.757854e+04 6.098280e+04 5.862432e+04 5.492330e+04 4.992116e+04 341 342 343 344 345 1.582974e+04 5.862432e+04 5.862432e+04 5.862432e+04 6.106346e+04 346 347 348 349 350 5.862432e+04 4.645356e+04 5.862432e+04 5.862432e+04 2.262242e+04 351 352 353 354 355 6.739839e+04 4.939214e+04 4.695392e+04 5.209069e+04 5.862432e+04 356 357 358 359 360 5.862432e+04 4.700609e+04 5.862432e+04 5.862432e+04 5.862432e+04 361 362 363 364 365 5.324862e+04 5.862432e+04 5.862432e+04 -6.743071e+04 4.776351e+04 366 367 368 369 370 5.862432e+04 5.862432e+04 5.862432e+04 5.862432e+04 6.086752e+04 371 372 373 374 375 5.862432e+04 8.723399e+04 4.129707e+04 3.370318e+04 5.862432e+04 376 377 378 379 380 5.862432e+04 4.384601e+04 4.809817e+04 6.233059e+03 5.862432e+04 381 382 383 384 385 5.862432e+04 5.225427e+04 1.009405e+05 5.862432e+04 4.409872e+04 386 387 388 389 390 -1.107469e+05 6.868064e+04 2.998467e+04 5.240005e+04 4.083358e+04 391 392 393 394 395 4.019934e+04 2.218813e+04 -1.927890e+04 7.782750e+04 -1.311252e+05 396 397 398 399 400 1.344257e+04 -1.231689e+05 -4.483172e+05 1.566969e+04 1.416191e+04 401 402 403 404 405 6.543617e+03 4.444421e+04 2.369665e+04 -8.590407e+04 -8.045107e+03 406 407 408 409 410 1.546557e+04 2.565864e+04 1.202447e+04 2.605950e+04 -6.015431e+04 411 412 413 414 415 -6.278733e+04 -1.912285e+05 -6.145847e+04 1.652315e+04 -1.744255e+05 416 417 418 419 420 -9.624296e+04 1.410265e+05 -6.088788e+04 -5.310840e+04 -3.228747e+04 421 422 423 424 425 -4.059382e+04 -2.050639e+05 -5.909810e+05 -4.524967e+05 -4.621551e+04 426 427 428 429 430 -2.926978e+05 -3.028997e+05 -1.237773e+06 -2.601184e+05 -3.540133e+05 431 -5.651288e+05 > postscript(file="/var/www/html/freestat/rcomp/tmp/6bwov1291224108.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 1.434540e+06 NA 1 2.592406e+06 1.434540e+06 2 1.260986e+06 2.592406e+06 3 -2.120094e+06 1.260986e+06 4 -5.389179e+05 -2.120094e+06 5 8.706658e+03 -5.389179e+05 6 3.681755e+04 8.706658e+03 7 1.011745e+05 3.681755e+04 8 -1.759896e+05 1.011745e+05 9 6.212848e+05 -1.759896e+05 10 -1.804099e+05 6.212848e+05 11 2.960787e+05 -1.804099e+05 12 1.183419e+05 2.960787e+05 13 -1.118617e+06 1.183419e+05 14 -7.714344e+04 -1.118617e+06 15 5.469366e+05 -7.714344e+04 16 -2.810947e+05 5.469366e+05 17 -6.509910e+04 -2.810947e+05 18 -3.783602e+05 -6.509910e+04 19 -5.639902e+05 -3.783602e+05 20 2.003133e+05 -5.639902e+05 21 -1.692827e+05 2.003133e+05 22 -1.583710e+05 -1.692827e+05 23 -7.263975e+05 -1.583710e+05 24 -3.223750e+05 -7.263975e+05 25 -8.611960e+04 -3.223750e+05 26 -5.244867e+05 -8.611960e+04 27 -2.287674e+05 -5.244867e+05 28 3.491694e+05 -2.287674e+05 29 -4.180886e+04 3.491694e+05 30 2.450783e+05 -4.180886e+04 31 -1.286764e+05 2.450783e+05 32 2.148073e+05 -1.286764e+05 33 8.843046e+04 2.148073e+05 34 -6.741604e+04 8.843046e+04 35 -6.997293e+05 -6.741604e+04 36 3.290698e+04 -6.997293e+05 37 9.512216e+04 3.290698e+04 38 -1.292850e+05 9.512216e+04 39 -3.292756e+05 -1.292850e+05 40 2.265576e+05 -3.292756e+05 41 -3.381594e+05 2.265576e+05 42 1.305302e+05 -3.381594e+05 43 -2.193006e+05 1.305302e+05 44 -9.686861e+03 -2.193006e+05 45 -9.873843e+04 -9.686861e+03 46 -1.506070e+05 -9.873843e+04 47 -1.527146e+05 -1.506070e+05 48 -4.405757e+04 -1.527146e+05 49 -2.023685e+05 -4.405757e+04 50 -1.079644e+04 -2.023685e+05 51 -1.554851e+05 -1.079644e+04 52 -2.840178e+05 -1.554851e+05 53 9.597283e+04 -2.840178e+05 54 -4.233277e+05 9.597283e+04 55 1.614878e+04 -4.233277e+05 56 3.274203e+05 1.614878e+04 57 -1.344067e+04 3.274203e+05 58 -1.167402e+04 -1.344067e+04 59 -2.489276e+05 -1.167402e+04 60 -3.053437e+04 -2.489276e+05 61 -1.788236e+05 -3.053437e+04 62 -2.695521e+04 -1.788236e+05 63 -2.218082e+05 -2.695521e+04 64 -4.244611e+02 -2.218082e+05 65 2.313680e+05 -4.244611e+02 66 -3.961293e+00 2.313680e+05 67 -5.613649e+04 -3.961293e+00 68 -1.039046e+05 -5.613649e+04 69 -1.621399e+05 -1.039046e+05 70 -1.575332e+05 -1.621399e+05 71 8.402731e+03 -1.575332e+05 72 -1.022603e+05 8.402731e+03 73 1.619668e+04 -1.022603e+05 74 -4.670361e+04 1.619668e+04 75 -1.401200e+05 -4.670361e+04 76 -1.895204e+04 -1.401200e+05 77 -3.751264e+04 -1.895204e+04 78 1.864700e+04 -3.751264e+04 79 1.016586e+05 1.864700e+04 80 1.779951e+04 1.016586e+05 81 -9.639347e+04 1.779951e+04 82 -6.771926e+04 -9.639347e+04 83 4.597112e+04 -6.771926e+04 84 1.936319e+04 4.597112e+04 85 -1.526918e+05 1.936319e+04 86 -7.729575e+04 -1.526918e+05 87 -2.709022e+03 -7.729575e+04 88 -1.087286e+05 -2.709022e+03 89 2.374239e+04 -1.087286e+05 90 5.171042e+03 2.374239e+04 91 -6.890441e+04 5.171042e+03 92 8.655845e+04 -6.890441e+04 93 -1.670658e+05 8.655845e+04 94 4.413430e+04 -1.670658e+05 95 -7.960861e+03 4.413430e+04 96 6.597813e+04 -7.960861e+03 97 1.644330e+05 6.597813e+04 98 -4.679715e+04 1.644330e+05 99 -3.656027e+04 -4.679715e+04 100 -7.862713e+04 -3.656027e+04 101 -1.134372e+05 -7.862713e+04 102 -1.552571e+05 -1.134372e+05 103 -4.287334e+04 -1.552571e+05 104 3.778051e+04 -4.287334e+04 105 -6.726786e+04 3.778051e+04 106 -2.868120e+05 -6.726786e+04 107 9.374730e+04 -2.868120e+05 108 -1.737003e+05 9.374730e+04 109 -1.146023e+05 -1.737003e+05 110 9.946503e+04 -1.146023e+05 111 -2.168834e+05 9.946503e+04 112 -1.463099e+05 -2.168834e+05 113 -9.484045e+04 -1.463099e+05 114 5.483301e+04 -9.484045e+04 115 1.896172e+04 5.483301e+04 116 3.684722e+04 1.896172e+04 117 7.276375e+04 3.684722e+04 118 -1.645875e+05 7.276375e+04 119 6.720513e+04 -1.645875e+05 120 -3.177472e+05 6.720513e+04 121 3.423384e+04 -3.177472e+05 122 -1.454836e+05 3.423384e+04 123 4.344236e+04 -1.454836e+05 124 -1.585187e+03 4.344236e+04 125 -2.042616e+05 -1.585187e+03 126 5.862432e+04 -2.042616e+05 127 3.871832e+04 5.862432e+04 128 3.169185e+04 3.871832e+04 129 5.862432e+04 3.169185e+04 130 4.877076e+04 5.862432e+04 131 9.773557e+04 4.877076e+04 132 5.862432e+04 9.773557e+04 133 8.357447e+03 5.862432e+04 134 3.708742e+04 8.357447e+03 135 5.216156e+04 3.708742e+04 136 4.932280e+04 5.216156e+04 137 5.783791e+04 4.932280e+04 138 3.444756e+04 5.783791e+04 139 5.862432e+04 3.444756e+04 140 5.938854e+04 5.862432e+04 141 5.862432e+04 5.938854e+04 142 1.652187e+04 5.862432e+04 143 5.147936e+04 1.652187e+04 144 3.851483e+04 5.147936e+04 145 4.662419e+04 3.851483e+04 146 5.862432e+04 4.662419e+04 147 4.729126e+04 5.862432e+04 148 5.991813e+04 4.729126e+04 149 5.948051e+04 5.991813e+04 150 2.112212e+04 5.948051e+04 151 -1.920308e+04 2.112212e+04 152 4.612637e+04 -1.920308e+04 153 2.978399e+04 4.612637e+04 154 5.616392e+04 2.978399e+04 155 5.862432e+04 5.616392e+04 156 5.862432e+04 5.862432e+04 157 5.238087e+04 5.862432e+04 158 6.626112e+04 5.238087e+04 159 5.862432e+04 6.626112e+04 160 5.862432e+04 5.862432e+04 161 5.862432e+04 5.862432e+04 162 5.862432e+04 5.862432e+04 163 5.258312e+04 5.862432e+04 164 3.036140e+04 5.258312e+04 165 5.807737e+04 3.036140e+04 166 5.862432e+04 5.807737e+04 167 1.115811e+05 5.862432e+04 168 3.010201e+05 1.115811e+05 169 5.862432e+04 3.010201e+05 170 5.862432e+04 5.862432e+04 171 5.059639e+04 5.862432e+04 172 5.862432e+04 5.059639e+04 173 3.431578e+04 5.862432e+04 174 5.862432e+04 3.431578e+04 175 4.724455e+04 5.862432e+04 176 -5.870740e+04 4.724455e+04 177 2.198344e+04 -5.870740e+04 178 5.862432e+04 2.198344e+04 179 5.862432e+04 5.862432e+04 180 5.508597e+04 5.862432e+04 181 -9.384035e+03 5.508597e+04 182 1.574433e+04 -9.384035e+03 183 -2.986051e+03 1.574433e+04 184 1.076837e+04 -2.986051e+03 185 1.486915e+04 1.076837e+04 186 -2.250814e+04 1.486915e+04 187 4.192277e+04 -2.250814e+04 188 1.502288e+04 4.192277e+04 189 5.862432e+04 1.502288e+04 190 4.517986e+04 5.862432e+04 191 5.862432e+04 4.517986e+04 192 4.347433e+04 5.862432e+04 193 5.862432e+04 4.347433e+04 194 5.285609e+04 5.862432e+04 195 2.161509e+04 5.285609e+04 196 -1.976193e+05 2.161509e+04 197 2.161509e+04 -1.976193e+05 198 2.161509e+04 2.161509e+04 199 -9.102825e+04 2.161509e+04 200 2.161509e+04 -9.102825e+04 201 2.161509e+04 2.161509e+04 202 2.161509e+04 2.161509e+04 203 2.058634e+04 2.161509e+04 204 -3.388936e+04 2.058634e+04 205 1.042777e+04 -3.388936e+04 206 1.478032e+04 1.042777e+04 207 2.161509e+04 1.478032e+04 208 1.755700e+04 2.161509e+04 209 3.580431e+04 1.755700e+04 210 5.862432e+04 3.580431e+04 211 5.862432e+04 5.862432e+04 212 8.218123e+04 5.862432e+04 213 1.798101e+04 8.218123e+04 214 2.855512e+04 1.798101e+04 215 5.362514e+04 2.855512e+04 216 5.268297e+04 5.362514e+04 217 6.738501e+04 5.268297e+04 218 5.362514e+04 6.738501e+04 219 -3.909274e+03 5.362514e+04 220 3.112129e+04 -3.909274e+03 221 5.862432e+04 3.112129e+04 222 5.838509e+04 5.862432e+04 223 1.954647e+04 5.838509e+04 224 5.805901e+04 1.954647e+04 225 2.608423e+04 5.805901e+04 226 5.017709e+04 2.608423e+04 227 5.862432e+04 5.017709e+04 228 4.343789e+04 5.862432e+04 229 5.862432e+04 4.343789e+04 230 2.264195e+04 5.862432e+04 231 3.765605e+04 2.264195e+04 232 5.862432e+04 3.765605e+04 233 2.404469e+04 5.862432e+04 234 5.862432e+04 2.404469e+04 235 5.980690e+04 5.862432e+04 236 5.862432e+04 5.980690e+04 237 5.045136e+04 5.862432e+04 238 5.862432e+04 5.045136e+04 239 5.862432e+04 5.862432e+04 240 5.862432e+04 5.862432e+04 241 4.899359e+04 5.862432e+04 242 5.862432e+04 4.899359e+04 243 5.162708e+04 5.862432e+04 244 -3.389133e+03 5.162708e+04 245 2.161509e+04 -3.389133e+03 246 2.161509e+04 2.161509e+04 247 2.280813e+04 2.161509e+04 248 2.161509e+04 2.280813e+04 249 2.161509e+04 2.161509e+04 250 3.608742e+04 2.161509e+04 251 2.161509e+04 3.608742e+04 252 4.086943e+04 2.161509e+04 253 3.105108e+03 4.086943e+04 254 -1.241809e+05 3.105108e+03 255 5.862432e+04 -1.241809e+05 256 4.647750e+04 5.862432e+04 257 3.103620e+04 4.647750e+04 258 2.161509e+04 3.103620e+04 259 5.862432e+04 2.161509e+04 260 2.161509e+04 5.862432e+04 261 2.161509e+04 2.161509e+04 262 4.343946e+04 2.161509e+04 263 6.400187e+04 4.343946e+04 264 2.161509e+04 6.400187e+04 265 5.862432e+04 2.161509e+04 266 2.161509e+04 5.862432e+04 267 8.231452e+04 2.161509e+04 268 2.161509e+04 8.231452e+04 269 -1.559911e+03 2.161509e+04 270 2.161509e+04 -1.559911e+03 271 2.161509e+04 2.161509e+04 272 2.161509e+04 2.161509e+04 273 1.699497e+04 2.161509e+04 274 5.862432e+04 1.699497e+04 275 5.896816e+04 5.862432e+04 276 3.129980e+04 5.896816e+04 277 5.862432e+04 3.129980e+04 278 2.101390e+04 5.862432e+04 279 5.012522e+04 2.101390e+04 280 2.161509e+04 5.012522e+04 281 5.862432e+04 2.161509e+04 282 7.633316e+04 5.862432e+04 283 5.862432e+04 7.633316e+04 284 4.304069e+04 5.862432e+04 285 5.862432e+04 4.304069e+04 286 6.827479e+03 5.862432e+04 287 5.561020e+04 6.827479e+03 288 2.161509e+04 5.561020e+04 289 8.052626e+04 2.161509e+04 290 2.161509e+04 8.052626e+04 291 5.907775e+04 2.161509e+04 292 4.212129e+04 5.907775e+04 293 5.862432e+04 4.212129e+04 294 3.773515e+04 5.862432e+04 295 4.297017e+04 3.773515e+04 296 5.862432e+04 4.297017e+04 297 2.161509e+04 5.862432e+04 298 5.862432e+04 2.161509e+04 299 3.612801e+04 5.862432e+04 300 1.324818e+04 3.612801e+04 301 4.114021e+04 1.324818e+04 302 5.612473e+04 4.114021e+04 303 2.161509e+04 5.612473e+04 304 1.042476e+05 2.161509e+04 305 3.899164e+04 1.042476e+05 306 2.161509e+04 3.899164e+04 307 -1.960650e+04 2.161509e+04 308 3.777727e+04 -1.960650e+04 309 2.161509e+04 3.777727e+04 310 2.161509e+04 2.161509e+04 311 5.382871e+04 2.161509e+04 312 2.161509e+04 5.382871e+04 313 -4.185561e+03 2.161509e+04 314 4.339971e+04 -4.185561e+03 315 7.721015e+04 4.339971e+04 316 5.989160e+04 7.721015e+04 317 1.324597e+03 5.989160e+04 318 5.304290e+04 1.324597e+03 319 5.862432e+04 5.304290e+04 320 5.862432e+04 5.862432e+04 321 5.862432e+04 5.862432e+04 322 5.862432e+04 5.862432e+04 323 2.161509e+04 5.862432e+04 324 5.862432e+04 2.161509e+04 325 5.508003e+04 5.862432e+04 326 2.161509e+04 5.508003e+04 327 2.161509e+04 2.161509e+04 328 5.862432e+04 2.161509e+04 329 5.862432e+04 5.862432e+04 330 5.862432e+04 5.862432e+04 331 5.862432e+04 5.862432e+04 332 5.862432e+04 5.862432e+04 333 5.124433e+04 5.862432e+04 334 5.841418e+04 5.124433e+04 335 4.757854e+04 5.841418e+04 336 6.098280e+04 4.757854e+04 337 5.862432e+04 6.098280e+04 338 5.492330e+04 5.862432e+04 339 4.992116e+04 5.492330e+04 340 1.582974e+04 4.992116e+04 341 5.862432e+04 1.582974e+04 342 5.862432e+04 5.862432e+04 343 5.862432e+04 5.862432e+04 344 6.106346e+04 5.862432e+04 345 5.862432e+04 6.106346e+04 346 4.645356e+04 5.862432e+04 347 5.862432e+04 4.645356e+04 348 5.862432e+04 5.862432e+04 349 2.262242e+04 5.862432e+04 350 6.739839e+04 2.262242e+04 351 4.939214e+04 6.739839e+04 352 4.695392e+04 4.939214e+04 353 5.209069e+04 4.695392e+04 354 5.862432e+04 5.209069e+04 355 5.862432e+04 5.862432e+04 356 4.700609e+04 5.862432e+04 357 5.862432e+04 4.700609e+04 358 5.862432e+04 5.862432e+04 359 5.862432e+04 5.862432e+04 360 5.324862e+04 5.862432e+04 361 5.862432e+04 5.324862e+04 362 5.862432e+04 5.862432e+04 363 -6.743071e+04 5.862432e+04 364 4.776351e+04 -6.743071e+04 365 5.862432e+04 4.776351e+04 366 5.862432e+04 5.862432e+04 367 5.862432e+04 5.862432e+04 368 5.862432e+04 5.862432e+04 369 6.086752e+04 5.862432e+04 370 5.862432e+04 6.086752e+04 371 8.723399e+04 5.862432e+04 372 4.129707e+04 8.723399e+04 373 3.370318e+04 4.129707e+04 374 5.862432e+04 3.370318e+04 375 5.862432e+04 5.862432e+04 376 4.384601e+04 5.862432e+04 377 4.809817e+04 4.384601e+04 378 6.233059e+03 4.809817e+04 379 5.862432e+04 6.233059e+03 380 5.862432e+04 5.862432e+04 381 5.225427e+04 5.862432e+04 382 1.009405e+05 5.225427e+04 383 5.862432e+04 1.009405e+05 384 4.409872e+04 5.862432e+04 385 -1.107469e+05 4.409872e+04 386 6.868064e+04 -1.107469e+05 387 2.998467e+04 6.868064e+04 388 5.240005e+04 2.998467e+04 389 4.083358e+04 5.240005e+04 390 4.019934e+04 4.083358e+04 391 2.218813e+04 4.019934e+04 392 -1.927890e+04 2.218813e+04 393 7.782750e+04 -1.927890e+04 394 -1.311252e+05 7.782750e+04 395 1.344257e+04 -1.311252e+05 396 -1.231689e+05 1.344257e+04 397 -4.483172e+05 -1.231689e+05 398 1.566969e+04 -4.483172e+05 399 1.416191e+04 1.566969e+04 400 6.543617e+03 1.416191e+04 401 4.444421e+04 6.543617e+03 402 2.369665e+04 4.444421e+04 403 -8.590407e+04 2.369665e+04 404 -8.045107e+03 -8.590407e+04 405 1.546557e+04 -8.045107e+03 406 2.565864e+04 1.546557e+04 407 1.202447e+04 2.565864e+04 408 2.605950e+04 1.202447e+04 409 -6.015431e+04 2.605950e+04 410 -6.278733e+04 -6.015431e+04 411 -1.912285e+05 -6.278733e+04 412 -6.145847e+04 -1.912285e+05 413 1.652315e+04 -6.145847e+04 414 -1.744255e+05 1.652315e+04 415 -9.624296e+04 -1.744255e+05 416 1.410265e+05 -9.624296e+04 417 -6.088788e+04 1.410265e+05 418 -5.310840e+04 -6.088788e+04 419 -3.228747e+04 -5.310840e+04 420 -4.059382e+04 -3.228747e+04 421 -2.050639e+05 -4.059382e+04 422 -5.909810e+05 -2.050639e+05 423 -4.524967e+05 -5.909810e+05 424 -4.621551e+04 -4.524967e+05 425 -2.926978e+05 -4.621551e+04 426 -3.028997e+05 -2.926978e+05 427 -1.237773e+06 -3.028997e+05 428 -2.601184e+05 -1.237773e+06 429 -3.540133e+05 -2.601184e+05 430 -5.651288e+05 -3.540133e+05 431 NA -5.651288e+05 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 2.592406e+06 1.434540e+06 [2,] 1.260986e+06 2.592406e+06 [3,] -2.120094e+06 1.260986e+06 [4,] -5.389179e+05 -2.120094e+06 [5,] 8.706658e+03 -5.389179e+05 [6,] 3.681755e+04 8.706658e+03 [7,] 1.011745e+05 3.681755e+04 [8,] -1.759896e+05 1.011745e+05 [9,] 6.212848e+05 -1.759896e+05 [10,] -1.804099e+05 6.212848e+05 [11,] 2.960787e+05 -1.804099e+05 [12,] 1.183419e+05 2.960787e+05 [13,] -1.118617e+06 1.183419e+05 [14,] -7.714344e+04 -1.118617e+06 [15,] 5.469366e+05 -7.714344e+04 [16,] -2.810947e+05 5.469366e+05 [17,] -6.509910e+04 -2.810947e+05 [18,] -3.783602e+05 -6.509910e+04 [19,] -5.639902e+05 -3.783602e+05 [20,] 2.003133e+05 -5.639902e+05 [21,] -1.692827e+05 2.003133e+05 [22,] -1.583710e+05 -1.692827e+05 [23,] -7.263975e+05 -1.583710e+05 [24,] -3.223750e+05 -7.263975e+05 [25,] -8.611960e+04 -3.223750e+05 [26,] -5.244867e+05 -8.611960e+04 [27,] -2.287674e+05 -5.244867e+05 [28,] 3.491694e+05 -2.287674e+05 [29,] -4.180886e+04 3.491694e+05 [30,] 2.450783e+05 -4.180886e+04 [31,] -1.286764e+05 2.450783e+05 [32,] 2.148073e+05 -1.286764e+05 [33,] 8.843046e+04 2.148073e+05 [34,] -6.741604e+04 8.843046e+04 [35,] -6.997293e+05 -6.741604e+04 [36,] 3.290698e+04 -6.997293e+05 [37,] 9.512216e+04 3.290698e+04 [38,] -1.292850e+05 9.512216e+04 [39,] -3.292756e+05 -1.292850e+05 [40,] 2.265576e+05 -3.292756e+05 [41,] -3.381594e+05 2.265576e+05 [42,] 1.305302e+05 -3.381594e+05 [43,] -2.193006e+05 1.305302e+05 [44,] -9.686861e+03 -2.193006e+05 [45,] -9.873843e+04 -9.686861e+03 [46,] -1.506070e+05 -9.873843e+04 [47,] -1.527146e+05 -1.506070e+05 [48,] -4.405757e+04 -1.527146e+05 [49,] -2.023685e+05 -4.405757e+04 [50,] -1.079644e+04 -2.023685e+05 [51,] -1.554851e+05 -1.079644e+04 [52,] -2.840178e+05 -1.554851e+05 [53,] 9.597283e+04 -2.840178e+05 [54,] -4.233277e+05 9.597283e+04 [55,] 1.614878e+04 -4.233277e+05 [56,] 3.274203e+05 1.614878e+04 [57,] -1.344067e+04 3.274203e+05 [58,] -1.167402e+04 -1.344067e+04 [59,] -2.489276e+05 -1.167402e+04 [60,] -3.053437e+04 -2.489276e+05 [61,] -1.788236e+05 -3.053437e+04 [62,] -2.695521e+04 -1.788236e+05 [63,] -2.218082e+05 -2.695521e+04 [64,] -4.244611e+02 -2.218082e+05 [65,] 2.313680e+05 -4.244611e+02 [66,] -3.961293e+00 2.313680e+05 [67,] -5.613649e+04 -3.961293e+00 [68,] -1.039046e+05 -5.613649e+04 [69,] -1.621399e+05 -1.039046e+05 [70,] -1.575332e+05 -1.621399e+05 [71,] 8.402731e+03 -1.575332e+05 [72,] -1.022603e+05 8.402731e+03 [73,] 1.619668e+04 -1.022603e+05 [74,] -4.670361e+04 1.619668e+04 [75,] -1.401200e+05 -4.670361e+04 [76,] -1.895204e+04 -1.401200e+05 [77,] -3.751264e+04 -1.895204e+04 [78,] 1.864700e+04 -3.751264e+04 [79,] 1.016586e+05 1.864700e+04 [80,] 1.779951e+04 1.016586e+05 [81,] -9.639347e+04 1.779951e+04 [82,] -6.771926e+04 -9.639347e+04 [83,] 4.597112e+04 -6.771926e+04 [84,] 1.936319e+04 4.597112e+04 [85,] -1.526918e+05 1.936319e+04 [86,] -7.729575e+04 -1.526918e+05 [87,] -2.709022e+03 -7.729575e+04 [88,] -1.087286e+05 -2.709022e+03 [89,] 2.374239e+04 -1.087286e+05 [90,] 5.171042e+03 2.374239e+04 [91,] -6.890441e+04 5.171042e+03 [92,] 8.655845e+04 -6.890441e+04 [93,] -1.670658e+05 8.655845e+04 [94,] 4.413430e+04 -1.670658e+05 [95,] -7.960861e+03 4.413430e+04 [96,] 6.597813e+04 -7.960861e+03 [97,] 1.644330e+05 6.597813e+04 [98,] -4.679715e+04 1.644330e+05 [99,] -3.656027e+04 -4.679715e+04 [100,] -7.862713e+04 -3.656027e+04 [101,] -1.134372e+05 -7.862713e+04 [102,] -1.552571e+05 -1.134372e+05 [103,] -4.287334e+04 -1.552571e+05 [104,] 3.778051e+04 -4.287334e+04 [105,] -6.726786e+04 3.778051e+04 [106,] -2.868120e+05 -6.726786e+04 [107,] 9.374730e+04 -2.868120e+05 [108,] -1.737003e+05 9.374730e+04 [109,] -1.146023e+05 -1.737003e+05 [110,] 9.946503e+04 -1.146023e+05 [111,] -2.168834e+05 9.946503e+04 [112,] -1.463099e+05 -2.168834e+05 [113,] -9.484045e+04 -1.463099e+05 [114,] 5.483301e+04 -9.484045e+04 [115,] 1.896172e+04 5.483301e+04 [116,] 3.684722e+04 1.896172e+04 [117,] 7.276375e+04 3.684722e+04 [118,] -1.645875e+05 7.276375e+04 [119,] 6.720513e+04 -1.645875e+05 [120,] -3.177472e+05 6.720513e+04 [121,] 3.423384e+04 -3.177472e+05 [122,] -1.454836e+05 3.423384e+04 [123,] 4.344236e+04 -1.454836e+05 [124,] -1.585187e+03 4.344236e+04 [125,] -2.042616e+05 -1.585187e+03 [126,] 5.862432e+04 -2.042616e+05 [127,] 3.871832e+04 5.862432e+04 [128,] 3.169185e+04 3.871832e+04 [129,] 5.862432e+04 3.169185e+04 [130,] 4.877076e+04 5.862432e+04 [131,] 9.773557e+04 4.877076e+04 [132,] 5.862432e+04 9.773557e+04 [133,] 8.357447e+03 5.862432e+04 [134,] 3.708742e+04 8.357447e+03 [135,] 5.216156e+04 3.708742e+04 [136,] 4.932280e+04 5.216156e+04 [137,] 5.783791e+04 4.932280e+04 [138,] 3.444756e+04 5.783791e+04 [139,] 5.862432e+04 3.444756e+04 [140,] 5.938854e+04 5.862432e+04 [141,] 5.862432e+04 5.938854e+04 [142,] 1.652187e+04 5.862432e+04 [143,] 5.147936e+04 1.652187e+04 [144,] 3.851483e+04 5.147936e+04 [145,] 4.662419e+04 3.851483e+04 [146,] 5.862432e+04 4.662419e+04 [147,] 4.729126e+04 5.862432e+04 [148,] 5.991813e+04 4.729126e+04 [149,] 5.948051e+04 5.991813e+04 [150,] 2.112212e+04 5.948051e+04 [151,] -1.920308e+04 2.112212e+04 [152,] 4.612637e+04 -1.920308e+04 [153,] 2.978399e+04 4.612637e+04 [154,] 5.616392e+04 2.978399e+04 [155,] 5.862432e+04 5.616392e+04 [156,] 5.862432e+04 5.862432e+04 [157,] 5.238087e+04 5.862432e+04 [158,] 6.626112e+04 5.238087e+04 [159,] 5.862432e+04 6.626112e+04 [160,] 5.862432e+04 5.862432e+04 [161,] 5.862432e+04 5.862432e+04 [162,] 5.862432e+04 5.862432e+04 [163,] 5.258312e+04 5.862432e+04 [164,] 3.036140e+04 5.258312e+04 [165,] 5.807737e+04 3.036140e+04 [166,] 5.862432e+04 5.807737e+04 [167,] 1.115811e+05 5.862432e+04 [168,] 3.010201e+05 1.115811e+05 [169,] 5.862432e+04 3.010201e+05 [170,] 5.862432e+04 5.862432e+04 [171,] 5.059639e+04 5.862432e+04 [172,] 5.862432e+04 5.059639e+04 [173,] 3.431578e+04 5.862432e+04 [174,] 5.862432e+04 3.431578e+04 [175,] 4.724455e+04 5.862432e+04 [176,] -5.870740e+04 4.724455e+04 [177,] 2.198344e+04 -5.870740e+04 [178,] 5.862432e+04 2.198344e+04 [179,] 5.862432e+04 5.862432e+04 [180,] 5.508597e+04 5.862432e+04 [181,] -9.384035e+03 5.508597e+04 [182,] 1.574433e+04 -9.384035e+03 [183,] -2.986051e+03 1.574433e+04 [184,] 1.076837e+04 -2.986051e+03 [185,] 1.486915e+04 1.076837e+04 [186,] -2.250814e+04 1.486915e+04 [187,] 4.192277e+04 -2.250814e+04 [188,] 1.502288e+04 4.192277e+04 [189,] 5.862432e+04 1.502288e+04 [190,] 4.517986e+04 5.862432e+04 [191,] 5.862432e+04 4.517986e+04 [192,] 4.347433e+04 5.862432e+04 [193,] 5.862432e+04 4.347433e+04 [194,] 5.285609e+04 5.862432e+04 [195,] 2.161509e+04 5.285609e+04 [196,] -1.976193e+05 2.161509e+04 [197,] 2.161509e+04 -1.976193e+05 [198,] 2.161509e+04 2.161509e+04 [199,] -9.102825e+04 2.161509e+04 [200,] 2.161509e+04 -9.102825e+04 [201,] 2.161509e+04 2.161509e+04 [202,] 2.161509e+04 2.161509e+04 [203,] 2.058634e+04 2.161509e+04 [204,] -3.388936e+04 2.058634e+04 [205,] 1.042777e+04 -3.388936e+04 [206,] 1.478032e+04 1.042777e+04 [207,] 2.161509e+04 1.478032e+04 [208,] 1.755700e+04 2.161509e+04 [209,] 3.580431e+04 1.755700e+04 [210,] 5.862432e+04 3.580431e+04 [211,] 5.862432e+04 5.862432e+04 [212,] 8.218123e+04 5.862432e+04 [213,] 1.798101e+04 8.218123e+04 [214,] 2.855512e+04 1.798101e+04 [215,] 5.362514e+04 2.855512e+04 [216,] 5.268297e+04 5.362514e+04 [217,] 6.738501e+04 5.268297e+04 [218,] 5.362514e+04 6.738501e+04 [219,] -3.909274e+03 5.362514e+04 [220,] 3.112129e+04 -3.909274e+03 [221,] 5.862432e+04 3.112129e+04 [222,] 5.838509e+04 5.862432e+04 [223,] 1.954647e+04 5.838509e+04 [224,] 5.805901e+04 1.954647e+04 [225,] 2.608423e+04 5.805901e+04 [226,] 5.017709e+04 2.608423e+04 [227,] 5.862432e+04 5.017709e+04 [228,] 4.343789e+04 5.862432e+04 [229,] 5.862432e+04 4.343789e+04 [230,] 2.264195e+04 5.862432e+04 [231,] 3.765605e+04 2.264195e+04 [232,] 5.862432e+04 3.765605e+04 [233,] 2.404469e+04 5.862432e+04 [234,] 5.862432e+04 2.404469e+04 [235,] 5.980690e+04 5.862432e+04 [236,] 5.862432e+04 5.980690e+04 [237,] 5.045136e+04 5.862432e+04 [238,] 5.862432e+04 5.045136e+04 [239,] 5.862432e+04 5.862432e+04 [240,] 5.862432e+04 5.862432e+04 [241,] 4.899359e+04 5.862432e+04 [242,] 5.862432e+04 4.899359e+04 [243,] 5.162708e+04 5.862432e+04 [244,] -3.389133e+03 5.162708e+04 [245,] 2.161509e+04 -3.389133e+03 [246,] 2.161509e+04 2.161509e+04 [247,] 2.280813e+04 2.161509e+04 [248,] 2.161509e+04 2.280813e+04 [249,] 2.161509e+04 2.161509e+04 [250,] 3.608742e+04 2.161509e+04 [251,] 2.161509e+04 3.608742e+04 [252,] 4.086943e+04 2.161509e+04 [253,] 3.105108e+03 4.086943e+04 [254,] -1.241809e+05 3.105108e+03 [255,] 5.862432e+04 -1.241809e+05 [256,] 4.647750e+04 5.862432e+04 [257,] 3.103620e+04 4.647750e+04 [258,] 2.161509e+04 3.103620e+04 [259,] 5.862432e+04 2.161509e+04 [260,] 2.161509e+04 5.862432e+04 [261,] 2.161509e+04 2.161509e+04 [262,] 4.343946e+04 2.161509e+04 [263,] 6.400187e+04 4.343946e+04 [264,] 2.161509e+04 6.400187e+04 [265,] 5.862432e+04 2.161509e+04 [266,] 2.161509e+04 5.862432e+04 [267,] 8.231452e+04 2.161509e+04 [268,] 2.161509e+04 8.231452e+04 [269,] -1.559911e+03 2.161509e+04 [270,] 2.161509e+04 -1.559911e+03 [271,] 2.161509e+04 2.161509e+04 [272,] 2.161509e+04 2.161509e+04 [273,] 1.699497e+04 2.161509e+04 [274,] 5.862432e+04 1.699497e+04 [275,] 5.896816e+04 5.862432e+04 [276,] 3.129980e+04 5.896816e+04 [277,] 5.862432e+04 3.129980e+04 [278,] 2.101390e+04 5.862432e+04 [279,] 5.012522e+04 2.101390e+04 [280,] 2.161509e+04 5.012522e+04 [281,] 5.862432e+04 2.161509e+04 [282,] 7.633316e+04 5.862432e+04 [283,] 5.862432e+04 7.633316e+04 [284,] 4.304069e+04 5.862432e+04 [285,] 5.862432e+04 4.304069e+04 [286,] 6.827479e+03 5.862432e+04 [287,] 5.561020e+04 6.827479e+03 [288,] 2.161509e+04 5.561020e+04 [289,] 8.052626e+04 2.161509e+04 [290,] 2.161509e+04 8.052626e+04 [291,] 5.907775e+04 2.161509e+04 [292,] 4.212129e+04 5.907775e+04 [293,] 5.862432e+04 4.212129e+04 [294,] 3.773515e+04 5.862432e+04 [295,] 4.297017e+04 3.773515e+04 [296,] 5.862432e+04 4.297017e+04 [297,] 2.161509e+04 5.862432e+04 [298,] 5.862432e+04 2.161509e+04 [299,] 3.612801e+04 5.862432e+04 [300,] 1.324818e+04 3.612801e+04 [301,] 4.114021e+04 1.324818e+04 [302,] 5.612473e+04 4.114021e+04 [303,] 2.161509e+04 5.612473e+04 [304,] 1.042476e+05 2.161509e+04 [305,] 3.899164e+04 1.042476e+05 [306,] 2.161509e+04 3.899164e+04 [307,] -1.960650e+04 2.161509e+04 [308,] 3.777727e+04 -1.960650e+04 [309,] 2.161509e+04 3.777727e+04 [310,] 2.161509e+04 2.161509e+04 [311,] 5.382871e+04 2.161509e+04 [312,] 2.161509e+04 5.382871e+04 [313,] -4.185561e+03 2.161509e+04 [314,] 4.339971e+04 -4.185561e+03 [315,] 7.721015e+04 4.339971e+04 [316,] 5.989160e+04 7.721015e+04 [317,] 1.324597e+03 5.989160e+04 [318,] 5.304290e+04 1.324597e+03 [319,] 5.862432e+04 5.304290e+04 [320,] 5.862432e+04 5.862432e+04 [321,] 5.862432e+04 5.862432e+04 [322,] 5.862432e+04 5.862432e+04 [323,] 2.161509e+04 5.862432e+04 [324,] 5.862432e+04 2.161509e+04 [325,] 5.508003e+04 5.862432e+04 [326,] 2.161509e+04 5.508003e+04 [327,] 2.161509e+04 2.161509e+04 [328,] 5.862432e+04 2.161509e+04 [329,] 5.862432e+04 5.862432e+04 [330,] 5.862432e+04 5.862432e+04 [331,] 5.862432e+04 5.862432e+04 [332,] 5.862432e+04 5.862432e+04 [333,] 5.124433e+04 5.862432e+04 [334,] 5.841418e+04 5.124433e+04 [335,] 4.757854e+04 5.841418e+04 [336,] 6.098280e+04 4.757854e+04 [337,] 5.862432e+04 6.098280e+04 [338,] 5.492330e+04 5.862432e+04 [339,] 4.992116e+04 5.492330e+04 [340,] 1.582974e+04 4.992116e+04 [341,] 5.862432e+04 1.582974e+04 [342,] 5.862432e+04 5.862432e+04 [343,] 5.862432e+04 5.862432e+04 [344,] 6.106346e+04 5.862432e+04 [345,] 5.862432e+04 6.106346e+04 [346,] 4.645356e+04 5.862432e+04 [347,] 5.862432e+04 4.645356e+04 [348,] 5.862432e+04 5.862432e+04 [349,] 2.262242e+04 5.862432e+04 [350,] 6.739839e+04 2.262242e+04 [351,] 4.939214e+04 6.739839e+04 [352,] 4.695392e+04 4.939214e+04 [353,] 5.209069e+04 4.695392e+04 [354,] 5.862432e+04 5.209069e+04 [355,] 5.862432e+04 5.862432e+04 [356,] 4.700609e+04 5.862432e+04 [357,] 5.862432e+04 4.700609e+04 [358,] 5.862432e+04 5.862432e+04 [359,] 5.862432e+04 5.862432e+04 [360,] 5.324862e+04 5.862432e+04 [361,] 5.862432e+04 5.324862e+04 [362,] 5.862432e+04 5.862432e+04 [363,] -6.743071e+04 5.862432e+04 [364,] 4.776351e+04 -6.743071e+04 [365,] 5.862432e+04 4.776351e+04 [366,] 5.862432e+04 5.862432e+04 [367,] 5.862432e+04 5.862432e+04 [368,] 5.862432e+04 5.862432e+04 [369,] 6.086752e+04 5.862432e+04 [370,] 5.862432e+04 6.086752e+04 [371,] 8.723399e+04 5.862432e+04 [372,] 4.129707e+04 8.723399e+04 [373,] 3.370318e+04 4.129707e+04 [374,] 5.862432e+04 3.370318e+04 [375,] 5.862432e+04 5.862432e+04 [376,] 4.384601e+04 5.862432e+04 [377,] 4.809817e+04 4.384601e+04 [378,] 6.233059e+03 4.809817e+04 [379,] 5.862432e+04 6.233059e+03 [380,] 5.862432e+04 5.862432e+04 [381,] 5.225427e+04 5.862432e+04 [382,] 1.009405e+05 5.225427e+04 [383,] 5.862432e+04 1.009405e+05 [384,] 4.409872e+04 5.862432e+04 [385,] -1.107469e+05 4.409872e+04 [386,] 6.868064e+04 -1.107469e+05 [387,] 2.998467e+04 6.868064e+04 [388,] 5.240005e+04 2.998467e+04 [389,] 4.083358e+04 5.240005e+04 [390,] 4.019934e+04 4.083358e+04 [391,] 2.218813e+04 4.019934e+04 [392,] -1.927890e+04 2.218813e+04 [393,] 7.782750e+04 -1.927890e+04 [394,] -1.311252e+05 7.782750e+04 [395,] 1.344257e+04 -1.311252e+05 [396,] -1.231689e+05 1.344257e+04 [397,] -4.483172e+05 -1.231689e+05 [398,] 1.566969e+04 -4.483172e+05 [399,] 1.416191e+04 1.566969e+04 [400,] 6.543617e+03 1.416191e+04 [401,] 4.444421e+04 6.543617e+03 [402,] 2.369665e+04 4.444421e+04 [403,] -8.590407e+04 2.369665e+04 [404,] -8.045107e+03 -8.590407e+04 [405,] 1.546557e+04 -8.045107e+03 [406,] 2.565864e+04 1.546557e+04 [407,] 1.202447e+04 2.565864e+04 [408,] 2.605950e+04 1.202447e+04 [409,] -6.015431e+04 2.605950e+04 [410,] -6.278733e+04 -6.015431e+04 [411,] -1.912285e+05 -6.278733e+04 [412,] -6.145847e+04 -1.912285e+05 [413,] 1.652315e+04 -6.145847e+04 [414,] -1.744255e+05 1.652315e+04 [415,] -9.624296e+04 -1.744255e+05 [416,] 1.410265e+05 -9.624296e+04 [417,] -6.088788e+04 1.410265e+05 [418,] -5.310840e+04 -6.088788e+04 [419,] -3.228747e+04 -5.310840e+04 [420,] -4.059382e+04 -3.228747e+04 [421,] -2.050639e+05 -4.059382e+04 [422,] -5.909810e+05 -2.050639e+05 [423,] -4.524967e+05 -5.909810e+05 [424,] -4.621551e+04 -4.524967e+05 [425,] -2.926978e+05 -4.621551e+04 [426,] -3.028997e+05 -2.926978e+05 [427,] -1.237773e+06 -3.028997e+05 [428,] -2.601184e+05 -1.237773e+06 [429,] -3.540133e+05 -2.601184e+05 [430,] -5.651288e+05 -3.540133e+05 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 2.592406e+06 1.434540e+06 2 1.260986e+06 2.592406e+06 3 -2.120094e+06 1.260986e+06 4 -5.389179e+05 -2.120094e+06 5 8.706658e+03 -5.389179e+05 6 3.681755e+04 8.706658e+03 7 1.011745e+05 3.681755e+04 8 -1.759896e+05 1.011745e+05 9 6.212848e+05 -1.759896e+05 10 -1.804099e+05 6.212848e+05 11 2.960787e+05 -1.804099e+05 12 1.183419e+05 2.960787e+05 13 -1.118617e+06 1.183419e+05 14 -7.714344e+04 -1.118617e+06 15 5.469366e+05 -7.714344e+04 16 -2.810947e+05 5.469366e+05 17 -6.509910e+04 -2.810947e+05 18 -3.783602e+05 -6.509910e+04 19 -5.639902e+05 -3.783602e+05 20 2.003133e+05 -5.639902e+05 21 -1.692827e+05 2.003133e+05 22 -1.583710e+05 -1.692827e+05 23 -7.263975e+05 -1.583710e+05 24 -3.223750e+05 -7.263975e+05 25 -8.611960e+04 -3.223750e+05 26 -5.244867e+05 -8.611960e+04 27 -2.287674e+05 -5.244867e+05 28 3.491694e+05 -2.287674e+05 29 -4.180886e+04 3.491694e+05 30 2.450783e+05 -4.180886e+04 31 -1.286764e+05 2.450783e+05 32 2.148073e+05 -1.286764e+05 33 8.843046e+04 2.148073e+05 34 -6.741604e+04 8.843046e+04 35 -6.997293e+05 -6.741604e+04 36 3.290698e+04 -6.997293e+05 37 9.512216e+04 3.290698e+04 38 -1.292850e+05 9.512216e+04 39 -3.292756e+05 -1.292850e+05 40 2.265576e+05 -3.292756e+05 41 -3.381594e+05 2.265576e+05 42 1.305302e+05 -3.381594e+05 43 -2.193006e+05 1.305302e+05 44 -9.686861e+03 -2.193006e+05 45 -9.873843e+04 -9.686861e+03 46 -1.506070e+05 -9.873843e+04 47 -1.527146e+05 -1.506070e+05 48 -4.405757e+04 -1.527146e+05 49 -2.023685e+05 -4.405757e+04 50 -1.079644e+04 -2.023685e+05 51 -1.554851e+05 -1.079644e+04 52 -2.840178e+05 -1.554851e+05 53 9.597283e+04 -2.840178e+05 54 -4.233277e+05 9.597283e+04 55 1.614878e+04 -4.233277e+05 56 3.274203e+05 1.614878e+04 57 -1.344067e+04 3.274203e+05 58 -1.167402e+04 -1.344067e+04 59 -2.489276e+05 -1.167402e+04 60 -3.053437e+04 -2.489276e+05 61 -1.788236e+05 -3.053437e+04 62 -2.695521e+04 -1.788236e+05 63 -2.218082e+05 -2.695521e+04 64 -4.244611e+02 -2.218082e+05 65 2.313680e+05 -4.244611e+02 66 -3.961293e+00 2.313680e+05 67 -5.613649e+04 -3.961293e+00 68 -1.039046e+05 -5.613649e+04 69 -1.621399e+05 -1.039046e+05 70 -1.575332e+05 -1.621399e+05 71 8.402731e+03 -1.575332e+05 72 -1.022603e+05 8.402731e+03 73 1.619668e+04 -1.022603e+05 74 -4.670361e+04 1.619668e+04 75 -1.401200e+05 -4.670361e+04 76 -1.895204e+04 -1.401200e+05 77 -3.751264e+04 -1.895204e+04 78 1.864700e+04 -3.751264e+04 79 1.016586e+05 1.864700e+04 80 1.779951e+04 1.016586e+05 81 -9.639347e+04 1.779951e+04 82 -6.771926e+04 -9.639347e+04 83 4.597112e+04 -6.771926e+04 84 1.936319e+04 4.597112e+04 85 -1.526918e+05 1.936319e+04 86 -7.729575e+04 -1.526918e+05 87 -2.709022e+03 -7.729575e+04 88 -1.087286e+05 -2.709022e+03 89 2.374239e+04 -1.087286e+05 90 5.171042e+03 2.374239e+04 91 -6.890441e+04 5.171042e+03 92 8.655845e+04 -6.890441e+04 93 -1.670658e+05 8.655845e+04 94 4.413430e+04 -1.670658e+05 95 -7.960861e+03 4.413430e+04 96 6.597813e+04 -7.960861e+03 97 1.644330e+05 6.597813e+04 98 -4.679715e+04 1.644330e+05 99 -3.656027e+04 -4.679715e+04 100 -7.862713e+04 -3.656027e+04 101 -1.134372e+05 -7.862713e+04 102 -1.552571e+05 -1.134372e+05 103 -4.287334e+04 -1.552571e+05 104 3.778051e+04 -4.287334e+04 105 -6.726786e+04 3.778051e+04 106 -2.868120e+05 -6.726786e+04 107 9.374730e+04 -2.868120e+05 108 -1.737003e+05 9.374730e+04 109 -1.146023e+05 -1.737003e+05 110 9.946503e+04 -1.146023e+05 111 -2.168834e+05 9.946503e+04 112 -1.463099e+05 -2.168834e+05 113 -9.484045e+04 -1.463099e+05 114 5.483301e+04 -9.484045e+04 115 1.896172e+04 5.483301e+04 116 3.684722e+04 1.896172e+04 117 7.276375e+04 3.684722e+04 118 -1.645875e+05 7.276375e+04 119 6.720513e+04 -1.645875e+05 120 -3.177472e+05 6.720513e+04 121 3.423384e+04 -3.177472e+05 122 -1.454836e+05 3.423384e+04 123 4.344236e+04 -1.454836e+05 124 -1.585187e+03 4.344236e+04 125 -2.042616e+05 -1.585187e+03 126 5.862432e+04 -2.042616e+05 127 3.871832e+04 5.862432e+04 128 3.169185e+04 3.871832e+04 129 5.862432e+04 3.169185e+04 130 4.877076e+04 5.862432e+04 131 9.773557e+04 4.877076e+04 132 5.862432e+04 9.773557e+04 133 8.357447e+03 5.862432e+04 134 3.708742e+04 8.357447e+03 135 5.216156e+04 3.708742e+04 136 4.932280e+04 5.216156e+04 137 5.783791e+04 4.932280e+04 138 3.444756e+04 5.783791e+04 139 5.862432e+04 3.444756e+04 140 5.938854e+04 5.862432e+04 141 5.862432e+04 5.938854e+04 142 1.652187e+04 5.862432e+04 143 5.147936e+04 1.652187e+04 144 3.851483e+04 5.147936e+04 145 4.662419e+04 3.851483e+04 146 5.862432e+04 4.662419e+04 147 4.729126e+04 5.862432e+04 148 5.991813e+04 4.729126e+04 149 5.948051e+04 5.991813e+04 150 2.112212e+04 5.948051e+04 151 -1.920308e+04 2.112212e+04 152 4.612637e+04 -1.920308e+04 153 2.978399e+04 4.612637e+04 154 5.616392e+04 2.978399e+04 155 5.862432e+04 5.616392e+04 156 5.862432e+04 5.862432e+04 157 5.238087e+04 5.862432e+04 158 6.626112e+04 5.238087e+04 159 5.862432e+04 6.626112e+04 160 5.862432e+04 5.862432e+04 161 5.862432e+04 5.862432e+04 162 5.862432e+04 5.862432e+04 163 5.258312e+04 5.862432e+04 164 3.036140e+04 5.258312e+04 165 5.807737e+04 3.036140e+04 166 5.862432e+04 5.807737e+04 167 1.115811e+05 5.862432e+04 168 3.010201e+05 1.115811e+05 169 5.862432e+04 3.010201e+05 170 5.862432e+04 5.862432e+04 171 5.059639e+04 5.862432e+04 172 5.862432e+04 5.059639e+04 173 3.431578e+04 5.862432e+04 174 5.862432e+04 3.431578e+04 175 4.724455e+04 5.862432e+04 176 -5.870740e+04 4.724455e+04 177 2.198344e+04 -5.870740e+04 178 5.862432e+04 2.198344e+04 179 5.862432e+04 5.862432e+04 180 5.508597e+04 5.862432e+04 181 -9.384035e+03 5.508597e+04 182 1.574433e+04 -9.384035e+03 183 -2.986051e+03 1.574433e+04 184 1.076837e+04 -2.986051e+03 185 1.486915e+04 1.076837e+04 186 -2.250814e+04 1.486915e+04 187 4.192277e+04 -2.250814e+04 188 1.502288e+04 4.192277e+04 189 5.862432e+04 1.502288e+04 190 4.517986e+04 5.862432e+04 191 5.862432e+04 4.517986e+04 192 4.347433e+04 5.862432e+04 193 5.862432e+04 4.347433e+04 194 5.285609e+04 5.862432e+04 195 2.161509e+04 5.285609e+04 196 -1.976193e+05 2.161509e+04 197 2.161509e+04 -1.976193e+05 198 2.161509e+04 2.161509e+04 199 -9.102825e+04 2.161509e+04 200 2.161509e+04 -9.102825e+04 201 2.161509e+04 2.161509e+04 202 2.161509e+04 2.161509e+04 203 2.058634e+04 2.161509e+04 204 -3.388936e+04 2.058634e+04 205 1.042777e+04 -3.388936e+04 206 1.478032e+04 1.042777e+04 207 2.161509e+04 1.478032e+04 208 1.755700e+04 2.161509e+04 209 3.580431e+04 1.755700e+04 210 5.862432e+04 3.580431e+04 211 5.862432e+04 5.862432e+04 212 8.218123e+04 5.862432e+04 213 1.798101e+04 8.218123e+04 214 2.855512e+04 1.798101e+04 215 5.362514e+04 2.855512e+04 216 5.268297e+04 5.362514e+04 217 6.738501e+04 5.268297e+04 218 5.362514e+04 6.738501e+04 219 -3.909274e+03 5.362514e+04 220 3.112129e+04 -3.909274e+03 221 5.862432e+04 3.112129e+04 222 5.838509e+04 5.862432e+04 223 1.954647e+04 5.838509e+04 224 5.805901e+04 1.954647e+04 225 2.608423e+04 5.805901e+04 226 5.017709e+04 2.608423e+04 227 5.862432e+04 5.017709e+04 228 4.343789e+04 5.862432e+04 229 5.862432e+04 4.343789e+04 230 2.264195e+04 5.862432e+04 231 3.765605e+04 2.264195e+04 232 5.862432e+04 3.765605e+04 233 2.404469e+04 5.862432e+04 234 5.862432e+04 2.404469e+04 235 5.980690e+04 5.862432e+04 236 5.862432e+04 5.980690e+04 237 5.045136e+04 5.862432e+04 238 5.862432e+04 5.045136e+04 239 5.862432e+04 5.862432e+04 240 5.862432e+04 5.862432e+04 241 4.899359e+04 5.862432e+04 242 5.862432e+04 4.899359e+04 243 5.162708e+04 5.862432e+04 244 -3.389133e+03 5.162708e+04 245 2.161509e+04 -3.389133e+03 246 2.161509e+04 2.161509e+04 247 2.280813e+04 2.161509e+04 248 2.161509e+04 2.280813e+04 249 2.161509e+04 2.161509e+04 250 3.608742e+04 2.161509e+04 251 2.161509e+04 3.608742e+04 252 4.086943e+04 2.161509e+04 253 3.105108e+03 4.086943e+04 254 -1.241809e+05 3.105108e+03 255 5.862432e+04 -1.241809e+05 256 4.647750e+04 5.862432e+04 257 3.103620e+04 4.647750e+04 258 2.161509e+04 3.103620e+04 259 5.862432e+04 2.161509e+04 260 2.161509e+04 5.862432e+04 261 2.161509e+04 2.161509e+04 262 4.343946e+04 2.161509e+04 263 6.400187e+04 4.343946e+04 264 2.161509e+04 6.400187e+04 265 5.862432e+04 2.161509e+04 266 2.161509e+04 5.862432e+04 267 8.231452e+04 2.161509e+04 268 2.161509e+04 8.231452e+04 269 -1.559911e+03 2.161509e+04 270 2.161509e+04 -1.559911e+03 271 2.161509e+04 2.161509e+04 272 2.161509e+04 2.161509e+04 273 1.699497e+04 2.161509e+04 274 5.862432e+04 1.699497e+04 275 5.896816e+04 5.862432e+04 276 3.129980e+04 5.896816e+04 277 5.862432e+04 3.129980e+04 278 2.101390e+04 5.862432e+04 279 5.012522e+04 2.101390e+04 280 2.161509e+04 5.012522e+04 281 5.862432e+04 2.161509e+04 282 7.633316e+04 5.862432e+04 283 5.862432e+04 7.633316e+04 284 4.304069e+04 5.862432e+04 285 5.862432e+04 4.304069e+04 286 6.827479e+03 5.862432e+04 287 5.561020e+04 6.827479e+03 288 2.161509e+04 5.561020e+04 289 8.052626e+04 2.161509e+04 290 2.161509e+04 8.052626e+04 291 5.907775e+04 2.161509e+04 292 4.212129e+04 5.907775e+04 293 5.862432e+04 4.212129e+04 294 3.773515e+04 5.862432e+04 295 4.297017e+04 3.773515e+04 296 5.862432e+04 4.297017e+04 297 2.161509e+04 5.862432e+04 298 5.862432e+04 2.161509e+04 299 3.612801e+04 5.862432e+04 300 1.324818e+04 3.612801e+04 301 4.114021e+04 1.324818e+04 302 5.612473e+04 4.114021e+04 303 2.161509e+04 5.612473e+04 304 1.042476e+05 2.161509e+04 305 3.899164e+04 1.042476e+05 306 2.161509e+04 3.899164e+04 307 -1.960650e+04 2.161509e+04 308 3.777727e+04 -1.960650e+04 309 2.161509e+04 3.777727e+04 310 2.161509e+04 2.161509e+04 311 5.382871e+04 2.161509e+04 312 2.161509e+04 5.382871e+04 313 -4.185561e+03 2.161509e+04 314 4.339971e+04 -4.185561e+03 315 7.721015e+04 4.339971e+04 316 5.989160e+04 7.721015e+04 317 1.324597e+03 5.989160e+04 318 5.304290e+04 1.324597e+03 319 5.862432e+04 5.304290e+04 320 5.862432e+04 5.862432e+04 321 5.862432e+04 5.862432e+04 322 5.862432e+04 5.862432e+04 323 2.161509e+04 5.862432e+04 324 5.862432e+04 2.161509e+04 325 5.508003e+04 5.862432e+04 326 2.161509e+04 5.508003e+04 327 2.161509e+04 2.161509e+04 328 5.862432e+04 2.161509e+04 329 5.862432e+04 5.862432e+04 330 5.862432e+04 5.862432e+04 331 5.862432e+04 5.862432e+04 332 5.862432e+04 5.862432e+04 333 5.124433e+04 5.862432e+04 334 5.841418e+04 5.124433e+04 335 4.757854e+04 5.841418e+04 336 6.098280e+04 4.757854e+04 337 5.862432e+04 6.098280e+04 338 5.492330e+04 5.862432e+04 339 4.992116e+04 5.492330e+04 340 1.582974e+04 4.992116e+04 341 5.862432e+04 1.582974e+04 342 5.862432e+04 5.862432e+04 343 5.862432e+04 5.862432e+04 344 6.106346e+04 5.862432e+04 345 5.862432e+04 6.106346e+04 346 4.645356e+04 5.862432e+04 347 5.862432e+04 4.645356e+04 348 5.862432e+04 5.862432e+04 349 2.262242e+04 5.862432e+04 350 6.739839e+04 2.262242e+04 351 4.939214e+04 6.739839e+04 352 4.695392e+04 4.939214e+04 353 5.209069e+04 4.695392e+04 354 5.862432e+04 5.209069e+04 355 5.862432e+04 5.862432e+04 356 4.700609e+04 5.862432e+04 357 5.862432e+04 4.700609e+04 358 5.862432e+04 5.862432e+04 359 5.862432e+04 5.862432e+04 360 5.324862e+04 5.862432e+04 361 5.862432e+04 5.324862e+04 362 5.862432e+04 5.862432e+04 363 -6.743071e+04 5.862432e+04 364 4.776351e+04 -6.743071e+04 365 5.862432e+04 4.776351e+04 366 5.862432e+04 5.862432e+04 367 5.862432e+04 5.862432e+04 368 5.862432e+04 5.862432e+04 369 6.086752e+04 5.862432e+04 370 5.862432e+04 6.086752e+04 371 8.723399e+04 5.862432e+04 372 4.129707e+04 8.723399e+04 373 3.370318e+04 4.129707e+04 374 5.862432e+04 3.370318e+04 375 5.862432e+04 5.862432e+04 376 4.384601e+04 5.862432e+04 377 4.809817e+04 4.384601e+04 378 6.233059e+03 4.809817e+04 379 5.862432e+04 6.233059e+03 380 5.862432e+04 5.862432e+04 381 5.225427e+04 5.862432e+04 382 1.009405e+05 5.225427e+04 383 5.862432e+04 1.009405e+05 384 4.409872e+04 5.862432e+04 385 -1.107469e+05 4.409872e+04 386 6.868064e+04 -1.107469e+05 387 2.998467e+04 6.868064e+04 388 5.240005e+04 2.998467e+04 389 4.083358e+04 5.240005e+04 390 4.019934e+04 4.083358e+04 391 2.218813e+04 4.019934e+04 392 -1.927890e+04 2.218813e+04 393 7.782750e+04 -1.927890e+04 394 -1.311252e+05 7.782750e+04 395 1.344257e+04 -1.311252e+05 396 -1.231689e+05 1.344257e+04 397 -4.483172e+05 -1.231689e+05 398 1.566969e+04 -4.483172e+05 399 1.416191e+04 1.566969e+04 400 6.543617e+03 1.416191e+04 401 4.444421e+04 6.543617e+03 402 2.369665e+04 4.444421e+04 403 -8.590407e+04 2.369665e+04 404 -8.045107e+03 -8.590407e+04 405 1.546557e+04 -8.045107e+03 406 2.565864e+04 1.546557e+04 407 1.202447e+04 2.565864e+04 408 2.605950e+04 1.202447e+04 409 -6.015431e+04 2.605950e+04 410 -6.278733e+04 -6.015431e+04 411 -1.912285e+05 -6.278733e+04 412 -6.145847e+04 -1.912285e+05 413 1.652315e+04 -6.145847e+04 414 -1.744255e+05 1.652315e+04 415 -9.624296e+04 -1.744255e+05 416 1.410265e+05 -9.624296e+04 417 -6.088788e+04 1.410265e+05 418 -5.310840e+04 -6.088788e+04 419 -3.228747e+04 -5.310840e+04 420 -4.059382e+04 -3.228747e+04 421 -2.050639e+05 -4.059382e+04 422 -5.909810e+05 -2.050639e+05 423 -4.524967e+05 -5.909810e+05 424 -4.621551e+04 -4.524967e+05 425 -2.926978e+05 -4.621551e+04 426 -3.028997e+05 -2.926978e+05 427 -1.237773e+06 -3.028997e+05 428 -2.601184e+05 -1.237773e+06 429 -3.540133e+05 -2.601184e+05 430 -5.651288e+05 -3.540133e+05 > 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/74ong1291224108.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/8wx411291224108.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/9wx411291224108.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/10wx411291224108.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) + plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint') + grid() + dev.off() + } null device 1 > > #Note: the /var/www/html/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/freestat/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE) > a<-table.row.end(a) > myeq <- colnames(x)[1] > myeq <- paste(myeq, '[t] = ', sep='') > for (i in 1:k){ + if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '') + myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ') + if (rownames(mysum$coefficients)[i] != '(Intercept)') { + myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='') + if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='') + } + } > myeq <- paste(myeq, ' + e[t]') > a<-table.row.start(a) > a<-table.element(a, myeq) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/freestat/rcomp/tmp/11bp2a1291224108.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Variable',header=TRUE) > a<-table.element(a,'Parameter',header=TRUE) > a<-table.element(a,'S.D.',header=TRUE) > a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE) > a<-table.element(a,'2-tail p-value',header=TRUE) > a<-table.element(a,'1-tail p-value',header=TRUE) > a<-table.row.end(a) > for (i in 1:k){ + a<-table.row.start(a) + a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE) + a<-table.element(a,mysum$coefficients[i,1]) + a<-table.element(a, round(mysum$coefficients[i,2],6)) + a<-table.element(a, round(mysum$coefficients[i,3],4)) + a<-table.element(a, round(mysum$coefficients[i,4],6)) + a<-table.element(a, round(mysum$coefficients[i,4]/2,6)) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/freestat/rcomp/tmp/123y1d1291224108.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple R',1,TRUE) > a<-table.element(a, sqrt(mysum$r.squared)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'R-squared',1,TRUE) > a<-table.element(a, mysum$r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Adjusted R-squared',1,TRUE) > a<-table.element(a, mysum$adj.r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (value)',1,TRUE) > a<-table.element(a, mysum$fstatistic[1]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[2]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[3]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'p-value',1,TRUE) > a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3])) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Residual Standard Deviation',1,TRUE) > a<-table.element(a, mysum$sigma) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Sum Squared Residuals',1,TRUE) > a<-table.element(a, sum(myerror*myerror)) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/freestat/rcomp/tmp/13shy71291224108.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Time or Index', 1, TRUE) > a<-table.element(a, 'Actuals', 1, TRUE) > a<-table.element(a, 'Interpolation
Forecast', 1, TRUE) > a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,i, 1, TRUE) + a<-table.element(a,x[i]) + a<-table.element(a,x[i]-mysum$resid[i]) + a<-table.element(a,mysum$resid[i]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/freestat/rcomp/tmp/143qfr1291224108.tab") > if (n > n25) { + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'p-values',header=TRUE) + a<-table.element(a,'Alternative Hypothesis',3,header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'breakpoint index',header=TRUE) + a<-table.element(a,'greater',header=TRUE) + a<-table.element(a,'2-sided',header=TRUE) + a<-table.element(a,'less',header=TRUE) + a<-table.row.end(a) + for (mypoint in kp3:nmkm3) { + a<-table.row.start(a) + a<-table.element(a,mypoint,header=TRUE) + a<-table.element(a,gqarr[mypoint-kp3+1,1]) + a<-table.element(a,gqarr[mypoint-kp3+1,2]) + a<-table.element(a,gqarr[mypoint-kp3+1,3]) + a<-table.row.end(a) + } + a<-table.end(a) + table.save(a,file="/var/www/html/freestat/rcomp/tmp/1569wx1291224108.tab") + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'Description',header=TRUE) + a<-table.element(a,'# significant tests',header=TRUE) + a<-table.element(a,'% significant tests',header=TRUE) + a<-table.element(a,'OK/NOK',header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'1% type I error level',header=TRUE) + a<-table.element(a,numsignificant1) + a<-table.element(a,numsignificant1/numgqtests) + if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'5% type I error level',header=TRUE) + a<-table.element(a,numsignificant5) + a<-table.element(a,numsignificant5/numgqtests) + if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'10% type I error level',header=TRUE) + a<-table.element(a,numsignificant10) + a<-table.element(a,numsignificant10/numgqtests) + if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.end(a) + table.save(a,file="/var/www/html/freestat/rcomp/tmp/16k1u61291224108.tab") + } > > try(system("convert tmp/1qe781291224108.ps tmp/1qe781291224108.png",intern=TRUE)) character(0) > try(system("convert tmp/2in6a1291224108.ps tmp/2in6a1291224108.png",intern=TRUE)) character(0) > try(system("convert tmp/3in6a1291224108.ps tmp/3in6a1291224108.png",intern=TRUE)) character(0) > try(system("convert tmp/4in6a1291224108.ps tmp/4in6a1291224108.png",intern=TRUE)) character(0) > try(system("convert tmp/5bwov1291224108.ps tmp/5bwov1291224108.png",intern=TRUE)) character(0) > try(system("convert tmp/6bwov1291224108.ps tmp/6bwov1291224108.png",intern=TRUE)) character(0) > try(system("convert tmp/74ong1291224108.ps tmp/74ong1291224108.png",intern=TRUE)) character(0) > try(system("convert tmp/8wx411291224108.ps tmp/8wx411291224108.png",intern=TRUE)) character(0) > try(system("convert tmp/9wx411291224108.ps tmp/9wx411291224108.png",intern=TRUE)) character(0) > try(system("convert tmp/10wx411291224108.ps tmp/10wx411291224108.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 16.410 3.224 17.664