R version 2.9.0 (2009-04-17) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > x <- array(list(1 + ,162556 + ,1081 + ,807 + ,213118 + ,6282154 + ,1 + ,29790 + ,309 + ,444 + ,81767 + ,4321023 + ,1 + ,87550 + ,458 + ,412 + ,153198 + ,4111912 + ,0 + ,84738 + ,588 + ,428 + ,-26007 + ,223193 + ,1 + ,54660 + ,302 + ,315 + ,126942 + ,1491348 + ,1 + ,42634 + ,156 + ,168 + ,157214 + ,1629616 + ,0 + ,40949 + ,481 + ,263 + ,129352 + ,1398893 + ,1 + ,45187 + ,353 + ,267 + ,234817 + ,1926517 + ,1 + ,37704 + ,452 + ,228 + ,60448 + ,983660 + ,1 + ,16275 + ,109 + ,129 + ,47818 + ,1443586 + ,0 + ,25830 + ,115 + ,104 + ,245546 + ,1073089 + ,0 + ,12679 + ,110 + ,122 + ,48020 + ,984885 + ,1 + ,18014 + ,239 + ,393 + ,-1710 + ,1405225 + ,0 + ,43556 + ,247 + ,190 + ,32648 + ,227132 + ,1 + ,24811 + ,505 + ,280 + ,95350 + ,929118 + ,0 + ,6575 + ,159 + ,63 + ,151352 + ,1071292 + ,0 + ,7123 + ,109 + ,102 + ,288170 + ,638830 + ,1 + ,21950 + ,519 + ,265 + ,114337 + ,856956 + ,1 + ,37597 + ,248 + ,234 + ,37884 + ,992426 + ,0 + ,17821 + ,373 + ,277 + ,122844 + ,444477 + ,1 + ,12988 + ,119 + ,73 + ,82340 + ,857217 + ,1 + ,22330 + ,84 + ,67 + ,79801 + ,711969 + ,0 + ,13326 + ,102 + ,103 + ,165548 + ,702380 + ,0 + ,16189 + ,295 + ,290 + ,116384 + ,358589 + ,0 + ,7146 + ,105 + ,83 + ,134028 + ,297978 + ,0 + ,15824 + ,64 + ,56 + ,63838 + ,585715 + ,1 + ,27664 + ,282 + ,236 + ,74996 + ,657954 + ,0 + ,11920 + ,182 + ,73 + ,31080 + ,209458 + ,0 + ,8568 + ,37 + ,34 + ,32168 + ,786690 + ,0 + ,14416 + ,361 + ,139 + ,49857 + ,439798 + ,1 + ,3369 + ,28 + ,26 + ,87161 + ,688779 + ,1 + ,11819 + ,85 + ,70 + ,106113 + ,574339 + ,1 + ,6984 + ,45 + ,40 + ,80570 + ,741409 + ,1 + ,4519 + ,49 + ,42 + ,102129 + ,597793 + ,0 + ,2220 + ,22 + ,12 + ,301670 + ,644190 + ,0 + ,18562 + ,155 + ,211 + ,102313 + ,377934 + ,0 + ,10327 + ,91 + ,74 + ,88577 + ,640273 + ,1 + ,5336 + ,81 + ,80 + ,112477 + ,697458 + ,1 + ,2365 + ,79 + ,83 + ,191778 + ,550608 + ,0 + ,4069 + ,145 + ,131 + ,79804 + ,207393 + ,0 + ,8636 + ,855 + ,203 + ,128294 + ,301607 + ,0 + ,13718 + ,61 + ,56 + ,96448 + ,345783 + ,0 + ,4525 + ,226 + ,89 + ,93811 + ,501749 + ,0 + ,6869 + ,105 + ,88 + ,117520 + ,379983 + ,0 + ,4628 + ,62 + ,39 + ,69159 + ,387475 + ,1 + ,3689 + ,25 + ,25 + ,101792 + ,377305 + ,1 + ,4891 + ,217 + ,49 + ,210568 + ,370837 + ,1 + ,7489 + ,322 + ,149 + ,136996 + ,430866 + ,0 + ,4901 + ,84 + ,58 + ,121920 + ,469107 + ,0 + ,2284 + ,33 + ,41 + ,76403 + ,194493 + ,1 + ,3160 + ,108 + ,90 + ,108094 + ,530670 + ,1 + ,4150 + ,150 + ,136 + ,134759 + ,518365 + ,1 + ,7285 + ,115 + ,97 + ,188873 + ,491303 + ,1 + ,1134 + ,162 + ,63 + ,146216 + ,527021 + ,1 + ,4658 + ,158 + ,114 + ,156608 + ,233773 + ,0 + ,2384 + ,97 + ,77 + ,61348 + ,405972 + ,0 + ,3748 + ,9 + ,6 + ,50350 + ,652925 + ,0 + ,5371 + ,66 + ,47 + ,87720 + ,446211 + ,0 + ,1285 + ,107 + ,51 + ,99489 + ,341340 + ,1 + ,9327 + ,101 + ,85 + ,87419 + ,387699 + ,1 + ,5565 + ,47 + ,43 + ,94355 + ,493408 + ,0 + ,1528 + ,38 + ,32 + ,60326 + ,146494 + ,1 + ,3122 + ,34 + ,25 + ,94670 + ,414462 + ,1 + ,7561 + ,87 + ,77 + ,82425 + ,364304 + ,0 + ,2675 + ,79 + ,54 + ,59017 + ,355178 + ,0 + ,13253 + ,947 + ,251 + ,90829 + ,357760 + ,0 + ,880 + ,74 + ,15 + ,80791 + ,261216 + ,1 + ,2053 + ,53 + ,44 + ,100423 + ,397144 + ,0 + ,1424 + ,94 + ,73 + ,131116 + ,374943 + ,1 + ,4036 + ,63 + ,85 + ,100269 + ,424898 + ,1 + ,3045 + ,58 + ,49 + ,27330 + ,202055 + ,0 + ,5119 + ,49 + ,38 + ,39039 + ,378525 + ,0 + ,1431 + ,34 + ,35 + ,106885 + ,310768 + ,0 + ,554 + ,11 + ,9 + ,79285 + ,325738 + ,0 + ,1975 + ,35 + ,34 + ,118881 + ,394510 + ,1 + ,1765 + ,20 + ,20 + ,77623 + ,247060 + ,0 + ,1012 + ,47 + ,29 + ,114768 + ,368078 + ,0 + ,810 + ,43 + ,11 + ,74015 + ,236761 + ,0 + ,1280 + ,117 + ,52 + ,69465 + ,312378 + ,1 + ,666 + ,171 + ,13 + ,117869 + ,339836 + ,0 + ,1380 + ,26 + ,29 + ,60982 + ,347385 + ,1 + ,4677 + ,75 + ,66 + ,90131 + ,426280 + ,0 + ,876 + ,59 + ,33 + ,138971 + ,352850 + ,0 + ,814 + ,18 + ,15 + ,39625 + ,301881 + ,0 + ,514 + ,15 + ,15 + ,102725 + ,377516 + ,1 + ,5692 + ,72 + ,68 + ,64239 + ,357312 + ,0 + ,3642 + ,86 + ,100 + ,90262 + ,458343 + ,0 + ,540 + ,14 + ,13 + ,103960 + ,354228 + ,0 + ,2099 + ,64 + ,45 + ,106611 + ,308636 + ,0 + ,567 + ,11 + ,14 + ,103345 + ,386212 + ,0 + ,2001 + ,52 + ,36 + ,95551 + ,393343 + ,1 + ,2949 + ,41 + ,40 + ,82903 + ,378509 + ,0 + ,2253 + ,99 + ,68 + ,63593 + ,452469 + ,1 + ,6533 + ,75 + ,29 + ,126910 + ,364839 + ,0 + ,1889 + ,45 + ,43 + ,37527 + ,358649 + ,1 + ,3055 + ,43 + ,30 + ,60247 + ,376641 + ,0 + ,272 + ,8 + ,9 + ,112995 + ,429112 + ,1 + ,1414 + ,198 + ,22 + ,70184 + ,330546 + ,0 + ,2564 + ,22 + ,19 + ,130140 + ,403560 + ,1 + ,1383 + ,11 + ,9 + ,73221 + ,317892 + ,1 + ,1261 + ,33 + ,31 + ,76114 + ,307528 + ,0 + ,975 + ,23 + ,19 + ,90534 + ,235133 + ,0 + ,3366 + ,80 + ,55 + ,108479 + ,299243 + ,0 + ,576 + ,18 + ,8 + ,113761 + ,314073 + ,0 + ,1686 + ,40 + ,28 + ,68696 + ,368186 + ,0 + ,746 + ,23 + ,29 + ,71561 + ,269661 + ,1 + ,3192 + ,60 + ,48 + ,59831 + ,125390 + ,1 + ,2045 + ,20 + ,16 + ,97890 + ,510834 + ,0 + ,5702 + ,61 + ,47 + ,101481 + ,321896 + ,1 + ,1932 + ,36 + ,20 + ,72954 + ,249898 + ,0 + ,936 + ,30 + ,22 + ,67939 + ,408881 + ,1 + ,3437 + ,47 + ,33 + ,48022 + ,158492 + ,0 + ,5131 + ,71 + ,44 + ,86111 + ,292154 + ,1 + ,2397 + ,14 + ,13 + ,74020 + ,289513 + ,1 + ,1389 + ,9 + ,6 + ,57530 + ,378049 + ,0 + ,1503 + ,39 + ,35 + ,56364 + ,343466 + ,0 + ,402 + ,26 + ,8 + ,84990 + ,332743 + ,0 + ,2239 + ,21 + ,17 + ,88590 + ,442882 + ,1 + ,2234 + ,16 + ,11 + ,77200 + ,214215 + ,0 + ,837 + ,69 + ,21 + ,61262 + ,315688 + ,0 + ,10579 + ,92 + ,92 + ,110309 + ,375195 + ,0 + ,875 + ,14 + ,12 + ,67000 + ,334280 + ,0 + ,1585 + ,107 + ,112 + ,93099 + ,355864 + ,1 + ,1659 + ,29 + ,25 + ,107577 + ,480382 + ,1 + ,2647 + ,37 + ,17 + ,62920 + ,353058 + ,1 + ,3294 + ,23 + ,23 + ,75832 + ,217193 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,94 + ,7 + ,10 + ,60793 + ,314533 + ,0 + ,422 + ,28 + ,23 + ,57935 + ,318056 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,34 + ,8 + ,7 + ,60630 + ,314353 + ,0 + ,1558 + ,63 + ,25 + ,55637 + ,369448 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,43 + ,3 + ,20 + ,60887 + ,312846 + ,0 + ,645 + ,5 + ,4 + ,60720 + ,312075 + ,0 + ,316 + ,9 + ,4 + ,60505 + ,315009 + ,0 + ,115 + ,13 + ,10 + ,60945 + ,318903 + ,0 + ,5 + ,2 + ,1 + ,60720 + ,314887 + ,0 + ,897 + ,5 + ,4 + ,60720 + ,314913 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,389 + ,14 + ,8 + ,58990 + ,325506 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1002 + ,15 + ,11 + ,56750 + ,298568 + ,0 + ,36 + ,3 + ,4 + ,60894 + ,315834 + ,0 + ,460 + ,15 + ,15 + ,63346 + ,329784 + ,0 + ,309 + ,11 + ,9 + ,56535 + ,312878 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,9 + ,6 + ,7 + ,60835 + ,314987 + ,0 + ,271 + ,2 + ,2 + ,60720 + ,325249 + ,0 + ,14 + ,1 + ,0 + ,61016 + ,315877 + ,0 + ,520 + ,10 + ,7 + ,58650 + ,291650 + ,0 + ,1766 + ,73 + ,46 + ,60438 + ,305959 + ,0 + ,0 + ,0 + ,5 + ,60720 + ,315380 + ,0 + ,458 + ,11 + ,7 + ,58625 + ,297765 + ,0 + ,20 + ,3 + ,2 + ,60938 + ,315245 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,98 + ,2 + ,2 + ,61490 + ,315236 + ,0 + ,405 + ,7 + ,5 + ,60845 + ,336425 + ,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 + ,483 + ,27 + ,7 + ,60830 + ,306268 + ,0 + ,454 + ,51 + ,24 + ,63261 + ,302187 + ,0 + ,47 + ,3 + ,1 + ,60720 + ,314882 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,757 + ,19 + ,18 + ,45689 + ,382712 + ,0 + ,4655 + ,393 + ,55 + ,60720 + ,341570 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,36 + ,4 + ,3 + ,61564 + ,312412 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,203 + ,9 + ,9 + ,61938 + ,309596 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,126 + ,10 + ,8 + ,60951 + ,315547 + ,0 + ,400 + ,152 + ,113 + ,60720 + ,313267 + ,1 + ,71 + ,1 + ,0 + ,60745 + ,316176 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,972 + ,34 + ,19 + ,71642 + ,359335 + ,1 + ,531 + ,10 + ,11 + ,71641 + ,330068 + ,0 + ,2461 + ,57 + ,25 + ,55792 + ,314289 + ,1 + ,378 + ,52 + ,16 + ,71873 + ,297413 + ,1 + ,23 + ,5 + ,5 + ,62555 + ,314806 + ,1 + ,638 + ,14 + ,11 + ,60370 + ,333210 + ,1 + ,2300 + ,29 + ,23 + ,64873 + ,352108 + ,0 + ,149 + ,5 + ,6 + ,62041 + ,313332 + ,0 + ,226 + ,5 + ,5 + ,65745 + ,291787 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,275 + ,4 + ,7 + ,59500 + ,318745 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,141 + ,6 + ,7 + ,61630 + ,315366 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,28 + ,2 + ,3 + ,60890 + ,315688 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,4980 + ,91 + ,89 + ,113521 + ,409642 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,472 + ,20 + ,19 + ,80045 + ,269587 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,203 + ,27 + ,12 + ,50804 + ,300962 + ,1 + ,496 + ,17 + ,12 + ,87390 + ,325479 + ,1 + ,10 + ,2 + ,5 + ,61656 + ,316155 + ,1 + ,63 + ,4 + ,2 + ,65688 + ,318574 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,1136 + ,32 + ,26 + ,48522 + ,343613 + ,1 + ,265 + ,31 + ,3 + ,60720 + ,306948 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,267 + ,32 + ,11 + ,57640 + ,330059 + ,0 + ,474 + ,20 + ,10 + ,61977 + ,288985 + ,0 + ,534 + ,7 + ,5 + ,62620 + ,304485 + ,0 + ,0 + ,0 + ,2 + ,60720 + ,315380 + ,0 + ,15 + ,8 + ,6 + ,60831 + ,315688 + ,0 + ,397 + ,28 + ,7 + ,60646 + ,317736 + ,0 + ,0 + ,0 + ,2 + ,60720 + ,315380 + ,0 + ,1866 + ,29 + ,28 + ,56225 + ,322331 + ,0 + ,288 + ,4 + ,3 + ,60510 + ,296656 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,3 + ,2 + ,1 + ,60698 + ,315354 + ,0 + ,468 + ,21 + ,20 + ,60720 + ,312161 + ,0 + ,20 + ,2 + ,1 + ,60805 + ,315576 + ,0 + ,278 + ,26 + ,22 + ,61404 + ,314922 + ,0 + ,61 + ,14 + ,9 + ,60720 + ,314551 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,192 + ,4 + ,2 + ,65276 + ,312339 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,317 + ,9 + ,7 + ,63915 + ,298700 + ,0 + ,738 + ,10 + ,9 + ,60720 + ,321376 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,368 + ,17 + ,13 + ,61686 + ,303230 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,2 + ,1 + ,0 + ,60743 + ,315487 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,53 + ,6 + ,6 + ,60349 + ,315793 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,94 + ,3 + ,3 + ,61360 + ,312887 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,24 + ,8 + ,7 + ,59818 + ,315637 + ,0 + ,2332 + ,4 + ,2 + ,72680 + ,324385 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,131 + ,11 + ,15 + ,61808 + ,308989 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,206 + ,9 + ,9 + ,53110 + ,296702 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,167 + ,2 + ,1 + ,64245 + ,307322 + ,0 + ,622 + ,73 + ,38 + ,73007 + ,304376 + ,0 + ,2328 + ,94 + ,57 + ,82732 + ,253588 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,365 + ,8 + ,7 + ,54820 + ,309560 + ,0 + ,364 + ,35 + ,26 + ,47705 + ,298466 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,226 + ,12 + ,13 + ,72835 + ,343929 + ,0 + ,307 + ,15 + ,10 + ,58856 + ,331955 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,188 + ,11 + ,9 + ,77655 + ,381180 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,138 + ,6 + ,26 + ,69817 + ,331420 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,125 + ,12 + ,19 + ,60798 + ,310201 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,282 + ,30 + ,12 + ,62452 + ,320016 + ,0 + ,335 + ,33 + ,23 + ,64175 + ,320398 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,1324 + ,117 + ,29 + ,67440 + ,291841 + ,0 + ,176 + ,28 + ,8 + ,68136 + ,310670 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,249 + ,72 + ,26 + ,56726 + ,313491 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,333 + ,13 + ,9 + ,70811 + ,331323 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,601 + ,6 + ,5 + ,60720 + ,319210 + ,0 + ,30 + ,4 + ,3 + ,62045 + ,318098 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,249 + ,62 + ,13 + ,54323 + ,292754 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,165 + ,24 + ,12 + ,62841 + ,325176 + ,0 + ,453 + ,21 + ,19 + ,81125 + ,365959 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,53 + ,14 + ,10 + ,59506 + ,302409 + ,1 + ,382 + ,21 + ,9 + ,59365 + ,340968 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,9 + ,60720 + ,315380 + ,1 + ,30 + ,4 + ,4 + ,60798 + ,313164 + ,0 + ,290 + ,2 + ,1 + ,58790 + ,301164 + ,0 + ,0 + ,0 + ,1 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,366 + ,53 + ,14 + ,61808 + ,344425 + ,0 + ,2 + ,9 + ,12 + ,60735 + ,315394 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,209 + ,13 + ,19 + ,64016 + ,316647 + ,0 + ,384 + ,22 + ,17 + ,54683 + ,309836 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,365 + ,83 + ,32 + ,87192 + ,346611 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,49 + ,8 + ,14 + ,64107 + ,322031 + ,0 + ,3 + ,4 + ,8 + ,60761 + ,315656 + ,0 + ,133 + ,14 + ,4 + ,65990 + ,339445 + ,0 + ,32 + ,1 + ,0 + ,59988 + ,314964 + ,0 + ,368 + ,17 + ,20 + ,61167 + ,297141 + ,0 + ,1 + ,6 + ,5 + ,60719 + ,315372 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,22 + ,2 + ,1 + ,60722 + ,312502 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,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 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,96 + ,5 + ,4 + ,60379 + ,313729 + ,0 + ,1 + ,2 + ,1 + ,60727 + ,315388 + ,0 + ,314 + ,5 + ,4 + ,60720 + ,315371 + ,0 + ,844 + ,78 + ,20 + ,60925 + ,296139 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,26 + ,1 + ,1 + ,60896 + ,313880 + ,0 + ,125 + ,13 + ,10 + ,59734 + ,317698 + ,0 + ,304 + ,15 + ,12 + ,62969 + ,295580 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,621 + ,48 + ,13 + ,60720 + ,308256 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,119 + ,6 + ,3 + ,59118 + ,303677 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1595 + ,17 + ,10 + ,60720 + ,319369 + ,0 + ,312 + ,14 + ,3 + ,58598 + ,318690 + ,0 + ,60 + ,10 + ,7 + ,61124 + ,314049 + ,0 + ,587 + ,12 + ,10 + ,59595 + ,325699 + ,0 + ,135 + ,2 + ,1 + ,62065 + ,314210 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,514 + ,52 + ,15 + ,78780 + ,322378 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,1 + ,4 + ,4 + ,60722 + ,315398 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,1 + ,1763 + ,24 + ,28 + ,61600 + ,308336 + ,0 + ,180 + ,11 + ,9 + ,59635 + ,316386 + ,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 + ,218 + ,21 + ,7 + ,60720 + ,315553 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,448 + ,40 + ,7 + ,59781 + ,323361 + ,0 + ,227 + ,9 + ,7 + ,76644 + ,336639 + ,0 + ,174 + ,1 + ,3 + ,64820 + ,307424 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,121 + ,24 + ,11 + ,56178 + ,295370 + ,0 + ,607 + ,11 + ,7 + ,60436 + ,322340 + ,0 + ,2212 + ,14 + ,10 + ,60720 + ,319864 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,530 + ,60 + ,18 + ,73433 + ,317291 + ,0 + ,571 + ,80 + ,14 + ,41477 + ,280398 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,0 + ,78 + ,16 + ,12 + ,62700 + ,317330 + ,0 + ,2489 + ,40 + ,29 + ,67804 + ,238125 + ,0 + ,131 + ,6 + ,3 + ,59661 + ,327071 + ,0 + ,923 + ,8 + ,6 + ,58620 + ,309038 + ,0 + ,72 + ,3 + ,3 + ,60398 + ,314210 + ,0 + ,572 + ,16 + ,8 + ,58580 + ,307930 + ,0 + ,397 + ,10 + ,10 + ,62710 + ,322327 + ,0 + ,450 + ,8 + ,6 + ,59325 + ,292136 + ,0 + ,622 + ,7 + ,8 + ,60950 + ,263276 + ,0 + ,694 + ,8 + ,6 + ,68060 + ,367655 + ,1 + ,3425 + ,12 + ,9 + ,83620 + ,283910 + ,0 + ,562 + ,13 + ,8 + ,58456 + ,283587 + ,0 + ,4917 + ,42 + ,26 + ,52811 + ,243650 + ,1 + ,1442 + ,118 + ,239 + ,121173 + ,438493 + ,0 + ,529 + ,9 + ,7 + ,63870 + ,296261 + ,1 + ,2126 + ,138 + ,41 + ,21001 + ,230621 + ,0 + ,1061 + ,5 + ,3 + ,70415 + ,304252 + ,0 + ,776 + ,9 + ,8 + ,64230 + ,333505 + ,0 + ,611 + ,8 + ,6 + ,59190 + ,296919 + ,1 + ,1526 + ,25 + ,21 + ,69351 + ,278990 + ,0 + ,592 + ,7 + ,7 + ,64270 + ,276898 + ,0 + ,1182 + ,13 + ,11 + ,70694 + ,327007 + ,0 + ,621 + ,16 + ,11 + ,68005 + ,317046 + ,0 + ,989 + ,11 + ,12 + ,58930 + ,304555 + ,0 + ,438 + ,11 + ,9 + ,58320 + ,298096 + ,0 + ,726 + ,3 + ,3 + ,69980 + ,231861 + ,0 + ,1303 + ,61 + ,57 + ,69863 + ,309422 + ,1 + ,7419 + ,29 + ,21 + ,63255 + ,286963 + ,1 + ,1164 + ,17 + ,15 + ,57320 + ,269753 + ,1 + ,3310 + ,33 + ,32 + ,75230 + ,448243 + ,0 + ,1920 + ,15 + ,11 + ,79420 + ,165404 + ,0 + ,965 + ,3 + ,2 + ,73490 + ,204325 + ,0 + ,3256 + ,66 + ,23 + ,35250 + ,407159 + ,1 + ,1135 + ,17 + ,20 + ,62285 + ,290476 + ,0 + ,1270 + ,26 + ,24 + ,69206 + ,275311 + ,0 + ,661 + ,3 + ,1 + ,65920 + ,246541 + ,0 + ,1013 + ,2 + ,1 + ,69770 + ,253468 + ,0 + ,2844 + ,67 + ,74 + ,72683 + ,240897 + ,1 + ,11528 + ,70 + ,68 + ,-14545 + ,-83265 + ,0 + ,6526 + ,26 + ,20 + ,55830 + ,-42143 + ,0 + ,2264 + ,24 + ,20 + ,55174 + ,272713 + ,1 + ,5109 + ,97 + ,82 + ,67038 + ,215362 + ,0 + ,3999 + ,30 + ,21 + ,51252 + ,42754 + ,0 + ,35624 + ,223 + ,244 + ,157278 + ,306275 + ,0 + ,9252 + ,48 + ,32 + ,79510 + ,253537 + ,0 + ,15236 + ,90 + ,86 + ,77440 + ,372631 + ,0 + ,18073 + ,180 + ,69 + ,27284 + ,-7170) + ,dim=c(6 + ,431) + ,dimnames=list(c('Group' + ,'Costs' + ,'Trades' + ,'Orders' + ,'Dividends' + ,'Wealth') + ,1:431)) > y <- array(NA,dim=c(6,431),dimnames=list(c('Group','Costs','Trades','Orders','Dividends','Wealth'),1:431)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par3 = 'No Linear Trend' > par2 = 'Do not include Seasonal Dummies' > par1 = '6' > #'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 Orders Dividends 1 6282154 1 162556 1081 807 213118 2 4321023 1 29790 309 444 81767 3 4111912 1 87550 458 412 153198 4 223193 0 84738 588 428 -26007 5 1491348 1 54660 302 315 126942 6 1629616 1 42634 156 168 157214 7 1398893 0 40949 481 263 129352 8 1926517 1 45187 353 267 234817 9 983660 1 37704 452 228 60448 10 1443586 1 16275 109 129 47818 11 1073089 0 25830 115 104 245546 12 984885 0 12679 110 122 48020 13 1405225 1 18014 239 393 -1710 14 227132 0 43556 247 190 32648 15 929118 1 24811 505 280 95350 16 1071292 0 6575 159 63 151352 17 638830 0 7123 109 102 288170 18 856956 1 21950 519 265 114337 19 992426 1 37597 248 234 37884 20 444477 0 17821 373 277 122844 21 857217 1 12988 119 73 82340 22 711969 1 22330 84 67 79801 23 702380 0 13326 102 103 165548 24 358589 0 16189 295 290 116384 25 297978 0 7146 105 83 134028 26 585715 0 15824 64 56 63838 27 657954 1 27664 282 236 74996 28 209458 0 11920 182 73 31080 29 786690 0 8568 37 34 32168 30 439798 0 14416 361 139 49857 31 688779 1 3369 28 26 87161 32 574339 1 11819 85 70 106113 33 741409 1 6984 45 40 80570 34 597793 1 4519 49 42 102129 35 644190 0 2220 22 12 301670 36 377934 0 18562 155 211 102313 37 640273 0 10327 91 74 88577 38 697458 1 5336 81 80 112477 39 550608 1 2365 79 83 191778 40 207393 0 4069 145 131 79804 41 301607 0 8636 855 203 128294 42 345783 0 13718 61 56 96448 43 501749 0 4525 226 89 93811 44 379983 0 6869 105 88 117520 45 387475 0 4628 62 39 69159 46 377305 1 3689 25 25 101792 47 370837 1 4891 217 49 210568 48 430866 1 7489 322 149 136996 49 469107 0 4901 84 58 121920 50 194493 0 2284 33 41 76403 51 530670 1 3160 108 90 108094 52 518365 1 4150 150 136 134759 53 491303 1 7285 115 97 188873 54 527021 1 1134 162 63 146216 55 233773 1 4658 158 114 156608 56 405972 0 2384 97 77 61348 57 652925 0 3748 9 6 50350 58 446211 0 5371 66 47 87720 59 341340 0 1285 107 51 99489 60 387699 1 9327 101 85 87419 61 493408 1 5565 47 43 94355 62 146494 0 1528 38 32 60326 63 414462 1 3122 34 25 94670 64 364304 1 7561 87 77 82425 65 355178 0 2675 79 54 59017 66 357760 0 13253 947 251 90829 67 261216 0 880 74 15 80791 68 397144 1 2053 53 44 100423 69 374943 0 1424 94 73 131116 70 424898 1 4036 63 85 100269 71 202055 1 3045 58 49 27330 72 378525 0 5119 49 38 39039 73 310768 0 1431 34 35 106885 74 325738 0 554 11 9 79285 75 394510 0 1975 35 34 118881 76 247060 1 1765 20 20 77623 77 368078 0 1012 47 29 114768 78 236761 0 810 43 11 74015 79 312378 0 1280 117 52 69465 80 339836 1 666 171 13 117869 81 347385 0 1380 26 29 60982 82 426280 1 4677 75 66 90131 83 352850 0 876 59 33 138971 84 301881 0 814 18 15 39625 85 377516 0 514 15 15 102725 86 357312 1 5692 72 68 64239 87 458343 0 3642 86 100 90262 88 354228 0 540 14 13 103960 89 308636 0 2099 64 45 106611 90 386212 0 567 11 14 103345 91 393343 0 2001 52 36 95551 92 378509 1 2949 41 40 82903 93 452469 0 2253 99 68 63593 94 364839 1 6533 75 29 126910 95 358649 0 1889 45 43 37527 96 376641 1 3055 43 30 60247 97 429112 0 272 8 9 112995 98 330546 1 1414 198 22 70184 99 403560 0 2564 22 19 130140 100 317892 1 1383 11 9 73221 101 307528 1 1261 33 31 76114 102 235133 0 975 23 19 90534 103 299243 0 3366 80 55 108479 104 314073 0 576 18 8 113761 105 368186 0 1686 40 28 68696 106 269661 0 746 23 29 71561 107 125390 1 3192 60 48 59831 108 510834 1 2045 20 16 97890 109 321896 0 5702 61 47 101481 110 249898 1 1932 36 20 72954 111 408881 0 936 30 22 67939 112 158492 1 3437 47 33 48022 113 292154 0 5131 71 44 86111 114 289513 1 2397 14 13 74020 115 378049 1 1389 9 6 57530 116 343466 0 1503 39 35 56364 117 332743 0 402 26 8 84990 118 442882 0 2239 21 17 88590 119 214215 1 2234 16 11 77200 120 315688 0 837 69 21 61262 121 375195 0 10579 92 92 110309 122 334280 0 875 14 12 67000 123 355864 0 1585 107 112 93099 124 480382 1 1659 29 25 107577 125 353058 1 2647 37 17 62920 126 217193 1 3294 23 23 75832 127 315380 0 0 0 0 60720 128 314533 0 94 7 10 60793 129 318056 0 422 28 23 57935 130 315380 0 0 0 0 60720 131 314353 0 34 8 7 60630 132 369448 0 1558 63 25 55637 133 315380 0 0 0 0 60720 134 312846 0 43 3 20 60887 135 312075 0 645 5 4 60720 136 315009 0 316 9 4 60505 137 318903 0 115 13 10 60945 138 314887 0 5 2 1 60720 139 314913 0 897 5 4 60720 140 315380 0 0 0 0 60720 141 325506 0 389 14 8 58990 142 315380 0 0 0 0 60720 143 298568 0 1002 15 11 56750 144 315834 0 36 3 4 60894 145 329784 0 460 15 15 63346 146 312878 0 309 11 9 56535 147 315380 0 0 0 0 60720 148 314987 0 9 6 7 60835 149 325249 0 271 2 2 60720 150 315877 0 14 1 0 61016 151 291650 0 520 10 7 58650 152 305959 0 1766 73 46 60438 153 315380 0 0 0 5 60720 154 297765 0 458 11 7 58625 155 315245 0 20 3 2 60938 156 315380 0 0 0 0 60720 157 315380 0 0 0 0 60720 158 315236 0 98 2 2 61490 159 336425 0 405 7 5 60845 160 315380 0 0 0 0 60720 161 315380 0 0 0 0 60720 162 315380 0 0 0 0 60720 163 315380 0 0 0 0 60720 164 306268 0 483 27 7 60830 165 302187 0 454 51 24 63261 166 314882 0 47 3 1 60720 167 315380 0 0 0 0 60720 168 382712 0 757 19 18 45689 169 341570 0 4655 393 55 60720 170 315380 0 0 0 0 60720 171 315380 0 0 0 0 60720 172 312412 0 36 4 3 61564 173 315380 0 0 0 0 60720 174 309596 0 203 9 9 61938 175 315380 0 0 0 0 60720 176 315547 0 126 10 8 60951 177 313267 0 400 152 113 60720 178 316176 1 71 1 0 60745 179 315380 0 0 0 0 60720 180 315380 0 0 0 0 60720 181 359335 0 972 34 19 71642 182 330068 1 531 10 11 71641 183 314289 0 2461 57 25 55792 184 297413 1 378 52 16 71873 185 314806 1 23 5 5 62555 186 333210 1 638 14 11 60370 187 352108 1 2300 29 23 64873 188 313332 0 149 5 6 62041 189 291787 0 226 5 5 65745 190 315380 0 0 0 0 60720 191 318745 0 275 4 7 59500 192 315380 0 0 0 0 60720 193 315366 0 141 6 7 61630 194 315380 0 0 0 0 60720 195 315688 0 28 2 3 60890 196 315380 1 0 0 0 60720 197 409642 1 4980 91 89 113521 198 315380 1 0 0 0 60720 199 315380 1 0 0 0 60720 200 269587 1 472 20 19 80045 201 315380 1 0 0 0 60720 202 315380 1 0 0 0 60720 203 315380 1 0 0 0 60720 204 300962 1 203 27 12 50804 205 325479 1 496 17 12 87390 206 316155 1 10 2 5 61656 207 318574 1 63 4 2 65688 208 315380 1 0 0 0 60720 209 343613 1 1136 32 26 48522 210 306948 1 265 31 3 60720 211 315380 0 0 0 0 60720 212 315380 0 0 0 0 60720 213 330059 0 267 32 11 57640 214 288985 0 474 20 10 61977 215 304485 0 534 7 5 62620 216 315380 0 0 0 2 60720 217 315688 0 15 8 6 60831 218 317736 0 397 28 7 60646 219 315380 0 0 0 2 60720 220 322331 0 1866 29 28 56225 221 296656 0 288 4 3 60510 222 315380 0 0 0 0 60720 223 315354 0 3 2 1 60698 224 312161 0 468 21 20 60720 225 315576 0 20 2 1 60805 226 314922 0 278 26 22 61404 227 314551 0 61 14 9 60720 228 315380 0 0 0 0 60720 229 312339 0 192 4 2 65276 230 315380 0 0 0 0 60720 231 298700 0 317 9 7 63915 232 321376 0 738 10 9 60720 233 315380 0 0 0 0 60720 234 303230 0 368 17 13 61686 235 315380 0 0 0 0 60720 236 315487 0 2 1 0 60743 237 315380 0 0 0 0 60720 238 315793 0 53 6 6 60349 239 315380 0 0 0 0 60720 240 315380 0 0 0 0 60720 241 315380 0 0 0 0 60720 242 312887 0 94 3 3 61360 243 315380 0 0 0 0 60720 244 315637 0 24 8 7 59818 245 324385 0 2332 4 2 72680 246 315380 1 0 0 0 60720 247 315380 1 0 0 0 60720 248 308989 0 131 11 15 61808 249 315380 1 0 0 0 60720 250 315380 1 0 0 0 60720 251 296702 0 206 9 9 53110 252 315380 1 0 0 0 60720 253 307322 0 167 2 1 64245 254 304376 0 622 73 38 73007 255 253588 0 2328 94 57 82732 256 315380 0 0 0 0 60720 257 309560 0 365 8 7 54820 258 298466 0 364 35 26 47705 259 315380 1 0 0 0 60720 260 315380 0 0 0 0 60720 261 315380 1 0 0 0 60720 262 315380 1 0 0 0 60720 263 343929 0 226 12 13 72835 264 331955 0 307 15 10 58856 265 315380 1 0 0 0 60720 266 315380 0 0 0 0 60720 267 315380 1 0 0 0 60720 268 381180 0 188 11 9 77655 269 315380 1 0 0 0 60720 270 331420 0 138 6 26 69817 271 315380 1 0 0 0 60720 272 315380 1 0 0 0 60720 273 315380 1 0 0 0 60720 274 310201 0 125 12 19 60798 275 315380 0 0 0 0 60720 276 320016 0 282 30 12 62452 277 320398 0 335 33 23 64175 278 315380 0 0 0 0 60720 279 291841 1 1324 117 29 67440 280 310670 0 176 28 8 68136 281 315380 1 0 0 0 60720 282 315380 0 0 0 0 60720 283 313491 0 249 72 26 56726 284 315380 0 0 0 0 60720 285 331323 0 333 13 9 70811 286 315380 0 0 0 0 60720 287 319210 1 601 6 5 60720 288 318098 0 30 4 3 62045 289 315380 1 0 0 0 60720 290 292754 0 249 62 13 54323 291 315380 1 0 0 0 60720 292 325176 0 165 24 12 62841 293 365959 0 453 21 19 81125 294 315380 0 0 0 0 60720 295 302409 0 53 14 10 59506 296 340968 1 382 21 9 59365 297 315380 0 0 0 0 60720 298 315380 1 0 0 0 60720 299 315380 0 0 0 0 60720 300 315380 0 0 0 9 60720 301 313164 1 30 4 4 60798 302 301164 0 290 2 1 58790 303 315380 0 0 0 1 60720 304 315380 1 0 0 0 60720 305 344425 0 366 53 14 61808 306 315394 0 2 9 12 60735 307 315380 1 0 0 0 60720 308 316647 1 209 13 19 64016 309 309836 0 384 22 17 54683 310 315380 1 0 0 0 60720 311 315380 1 0 0 0 60720 312 346611 0 365 83 32 87192 313 315380 1 0 0 0 60720 314 322031 1 49 8 14 64107 315 315656 0 3 4 8 60761 316 339445 0 133 14 4 65990 317 314964 0 32 1 0 59988 318 297141 0 368 17 20 61167 319 315372 0 1 6 5 60719 320 315380 0 0 0 0 60720 321 315380 0 0 0 0 60720 322 315380 0 0 0 0 60720 323 315380 0 0 0 0 60720 324 315380 1 0 0 0 60720 325 315380 0 0 0 0 60720 326 312502 0 22 2 1 60722 327 315380 1 0 0 0 60720 328 315380 1 0 0 0 60720 329 315380 0 0 0 0 60720 330 315380 0 0 0 0 60720 331 315380 0 0 0 0 60720 332 315380 0 0 0 0 60720 333 315380 0 0 0 0 60720 334 313729 0 96 5 4 60379 335 315388 0 1 2 1 60727 336 315371 0 314 5 4 60720 337 296139 0 844 78 20 60925 338 315380 0 0 0 0 60720 339 313880 0 26 1 1 60896 340 317698 0 125 13 10 59734 341 295580 0 304 15 12 62969 342 315380 0 0 0 0 60720 343 315380 0 0 0 0 60720 344 315380 0 0 0 0 60720 345 308256 0 621 48 13 60720 346 315380 0 0 0 0 60720 347 303677 0 119 6 3 59118 348 315380 0 0 0 0 60720 349 315380 0 0 0 0 60720 350 319369 0 1595 17 10 60720 351 318690 0 312 14 3 58598 352 314049 0 60 10 7 61124 353 325699 0 587 12 10 59595 354 314210 0 135 2 1 62065 355 315380 0 0 0 0 60720 356 315380 0 0 0 0 60720 357 322378 0 514 52 15 78780 358 315380 0 0 0 0 60720 359 315380 0 0 0 0 60720 360 315380 0 0 0 0 60720 361 315398 0 1 4 4 60722 362 315380 0 0 0 0 60720 363 315380 0 0 0 0 60720 364 308336 1 1763 24 28 61600 365 316386 0 180 11 9 59635 366 315380 0 0 0 0 60720 367 315380 0 0 0 0 60720 368 315380 0 0 0 0 60720 369 315380 0 0 0 0 60720 370 315553 0 218 21 7 60720 371 315380 0 0 0 0 60720 372 323361 0 448 40 7 59781 373 336639 0 227 9 7 76644 374 307424 0 174 1 3 64820 375 315380 0 0 0 0 60720 376 315380 0 0 0 0 60720 377 295370 0 121 24 11 56178 378 322340 0 607 11 7 60436 379 319864 0 2212 14 10 60720 380 315380 0 0 0 0 60720 381 315380 0 0 0 0 60720 382 317291 0 530 60 18 73433 383 280398 0 571 80 14 41477 384 315380 0 0 0 0 60720 385 317330 0 78 16 12 62700 386 238125 0 2489 40 29 67804 387 327071 0 131 6 3 59661 388 309038 0 923 8 6 58620 389 314210 0 72 3 3 60398 390 307930 0 572 16 8 58580 391 322327 0 397 10 10 62710 392 292136 0 450 8 6 59325 393 263276 0 622 7 8 60950 394 367655 0 694 8 6 68060 395 283910 1 3425 12 9 83620 396 283587 0 562 13 8 58456 397 243650 0 4917 42 26 52811 398 438493 1 1442 118 239 121173 399 296261 0 529 9 7 63870 400 230621 1 2126 138 41 21001 401 304252 0 1061 5 3 70415 402 333505 0 776 9 8 64230 403 296919 0 611 8 6 59190 404 278990 1 1526 25 21 69351 405 276898 0 592 7 7 64270 406 327007 0 1182 13 11 70694 407 317046 0 621 16 11 68005 408 304555 0 989 11 12 58930 409 298096 0 438 11 9 58320 410 231861 0 726 3 3 69980 411 309422 0 1303 61 57 69863 412 286963 1 7419 29 21 63255 413 269753 1 1164 17 15 57320 414 448243 1 3310 33 32 75230 415 165404 0 1920 15 11 79420 416 204325 0 965 3 2 73490 417 407159 0 3256 66 23 35250 418 290476 1 1135 17 20 62285 419 275311 0 1270 26 24 69206 420 246541 0 661 3 1 65920 421 253468 0 1013 2 1 69770 422 240897 0 2844 67 74 72683 423 -83265 1 11528 70 68 -14545 424 -42143 0 6526 26 20 55830 425 272713 0 2264 24 20 55174 426 215362 1 5109 97 82 67038 427 42754 0 3999 30 21 51252 428 306275 0 35624 223 244 157278 429 253537 0 9252 48 32 79510 430 372631 0 15236 90 86 77440 431 -7170 0 18073 180 69 27284 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) Group Costs Trades Orders Dividends 155435.184 37009.233 21.738 -1157.434 2499.590 1.669 > (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 *** Orders 2.500e+03 4.380e+02 5.707 2.16e-08 *** Dividends 1.669e+00 3.987e-01 4.185 3.46e-05 *** --- 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/rcomp/tmp/186hp1291392609.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(x[,1], type='l', main='Actuals and Interpolation', ylab='value of Actuals and Interpolation (dots)', xlab='time or index') > points(x[,1]-mysum$resid) > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/2jfgs1291392609.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(mysum$resid, type='b', pch=19, main='Residuals', ylab='value of Residuals', xlab='time or index') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/3jfgs1291392609.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > hist(mysum$resid, main='Residual Histogram', xlab='values of Residuals') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/4jfgs1291392609.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > densityplot(~mysum$resid,col='black',main='Residual Density Plot', xlab='values of Residuals') > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5coxd1291392609.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/rcomp/tmp/6coxd1291392609.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/rcomp/tmp/7mgwg1291392609.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > acf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/8x7d11291392609.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > pacf(mysum$resid, lag.max=length(mysum$resid)/2, main='Residual Partial Autocorrelation Function') > grid() > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/9x7d11291392609.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > opar <- par(mfrow = c(2,2), oma = c(0, 0, 1.1, 0)) > plot(mylm, las = 1, sub='Residual Diagnostics') > par(opar) > dev.off() null device 1 > if (n > n25) { + postscript(file="/var/www/html/rcomp/tmp/10x7d11291392609.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) + plot(kp3:nmkm3,gqarr[,2], main='Goldfeld-Quandt test',ylab='2-sided p-value',xlab='breakpoint') + grid() + dev.off() + } null device 1 > > #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Estimated Regression Equation', 1, TRUE) > a<-table.row.end(a) > myeq <- colnames(x)[1] > myeq <- paste(myeq, '[t] = ', sep='') > for (i in 1:k){ + if (mysum$coefficients[i,1] > 0) myeq <- paste(myeq, '+', '') + myeq <- paste(myeq, mysum$coefficients[i,1], sep=' ') + if (rownames(mysum$coefficients)[i] != '(Intercept)') { + myeq <- paste(myeq, rownames(mysum$coefficients)[i], sep='') + if (rownames(mysum$coefficients)[i] != 't') myeq <- paste(myeq, '[t]', sep='') + } + } > myeq <- paste(myeq, ' + e[t]') > a<-table.row.start(a) > a<-table.element(a, myeq) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/111qcp1291392609.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/ols1.htm','Multiple Linear Regression - Ordinary Least Squares',''), 6, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a,'Variable',header=TRUE) > a<-table.element(a,'Parameter',header=TRUE) > a<-table.element(a,'S.D.',header=TRUE) > a<-table.element(a,'T-STAT
H0: parameter = 0',header=TRUE) > a<-table.element(a,'2-tail p-value',header=TRUE) > a<-table.element(a,'1-tail p-value',header=TRUE) > a<-table.row.end(a) > for (i in 1:k){ + a<-table.row.start(a) + a<-table.element(a,rownames(mysum$coefficients)[i],header=TRUE) + a<-table.element(a,mysum$coefficients[i,1]) + a<-table.element(a, round(mysum$coefficients[i,2],6)) + a<-table.element(a, round(mysum$coefficients[i,3],4)) + a<-table.element(a, round(mysum$coefficients[i,4],6)) + a<-table.element(a, round(mysum$coefficients[i,4]/2,6)) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/12m8tv1291392609.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Regression Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple R',1,TRUE) > a<-table.element(a, sqrt(mysum$r.squared)) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'R-squared',1,TRUE) > a<-table.element(a, mysum$r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Adjusted R-squared',1,TRUE) > a<-table.element(a, mysum$adj.r.squared) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (value)',1,TRUE) > a<-table.element(a, mysum$fstatistic[1]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF numerator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[2]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'F-TEST (DF denominator)',1,TRUE) > a<-table.element(a, mysum$fstatistic[3]) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'p-value',1,TRUE) > a<-table.element(a, 1-pf(mysum$fstatistic[1],mysum$fstatistic[2],mysum$fstatistic[3])) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Residual Statistics', 2, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Residual Standard Deviation',1,TRUE) > a<-table.element(a, mysum$sigma) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Sum Squared Residuals',1,TRUE) > a<-table.element(a, sum(myerror*myerror)) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/13t98o1291392609.tab") > a<-table.start() > a<-table.row.start(a) > a<-table.element(a, 'Multiple Linear Regression - Actuals, Interpolation, and Residuals', 4, TRUE) > a<-table.row.end(a) > a<-table.row.start(a) > a<-table.element(a, 'Time or Index', 1, TRUE) > a<-table.element(a, 'Actuals', 1, TRUE) > a<-table.element(a, 'Interpolation
Forecast', 1, TRUE) > a<-table.element(a, 'Residuals
Prediction Error', 1, TRUE) > a<-table.row.end(a) > for (i in 1:n) { + a<-table.row.start(a) + a<-table.element(a,i, 1, TRUE) + a<-table.element(a,x[i]) + a<-table.element(a,x[i]-mysum$resid[i]) + a<-table.element(a,mysum$resid[i]) + a<-table.row.end(a) + } > a<-table.end(a) > table.save(a,file="/var/www/html/rcomp/tmp/1430p91291392609.tab") > if (n > n25) { + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'p-values',header=TRUE) + a<-table.element(a,'Alternative Hypothesis',3,header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'breakpoint index',header=TRUE) + a<-table.element(a,'greater',header=TRUE) + a<-table.element(a,'2-sided',header=TRUE) + a<-table.element(a,'less',header=TRUE) + a<-table.row.end(a) + for (mypoint in kp3:nmkm3) { + a<-table.row.start(a) + a<-table.element(a,mypoint,header=TRUE) + a<-table.element(a,gqarr[mypoint-kp3+1,1]) + a<-table.element(a,gqarr[mypoint-kp3+1,2]) + a<-table.element(a,gqarr[mypoint-kp3+1,3]) + a<-table.row.end(a) + } + a<-table.end(a) + table.save(a,file="/var/www/html/rcomp/tmp/15pj5f1291392609.tab") + a<-table.start() + a<-table.row.start(a) + a<-table.element(a,'Meta Analysis of Goldfeld-Quandt test for Heteroskedasticity',4,TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'Description',header=TRUE) + a<-table.element(a,'# significant tests',header=TRUE) + a<-table.element(a,'% significant tests',header=TRUE) + a<-table.element(a,'OK/NOK',header=TRUE) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'1% type I error level',header=TRUE) + a<-table.element(a,numsignificant1) + a<-table.element(a,numsignificant1/numgqtests) + if (numsignificant1/numgqtests < 0.01) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'5% type I error level',header=TRUE) + a<-table.element(a,numsignificant5) + a<-table.element(a,numsignificant5/numgqtests) + if (numsignificant5/numgqtests < 0.05) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.row.start(a) + a<-table.element(a,'10% type I error level',header=TRUE) + a<-table.element(a,numsignificant10) + a<-table.element(a,numsignificant10/numgqtests) + if (numsignificant10/numgqtests < 0.1) dum <- 'OK' else dum <- 'NOK' + a<-table.element(a,dum) + a<-table.row.end(a) + a<-table.end(a) + table.save(a,file="/var/www/html/rcomp/tmp/163t361291392609.tab") + } > > try(system("convert tmp/186hp1291392609.ps tmp/186hp1291392609.png",intern=TRUE)) character(0) > try(system("convert tmp/2jfgs1291392609.ps tmp/2jfgs1291392609.png",intern=TRUE)) character(0) > try(system("convert tmp/3jfgs1291392609.ps tmp/3jfgs1291392609.png",intern=TRUE)) character(0) > try(system("convert tmp/4jfgs1291392609.ps tmp/4jfgs1291392609.png",intern=TRUE)) character(0) > try(system("convert tmp/5coxd1291392609.ps tmp/5coxd1291392609.png",intern=TRUE)) character(0) > try(system("convert tmp/6coxd1291392609.ps tmp/6coxd1291392609.png",intern=TRUE)) character(0) > try(system("convert tmp/7mgwg1291392609.ps tmp/7mgwg1291392609.png",intern=TRUE)) character(0) > try(system("convert tmp/8x7d11291392609.ps tmp/8x7d11291392609.png",intern=TRUE)) character(0) > try(system("convert tmp/9x7d11291392609.ps tmp/9x7d11291392609.png",intern=TRUE)) character(0) > try(system("convert tmp/10x7d11291392609.ps tmp/10x7d11291392609.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 12.522 2.144 28.291