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(0 + ,162556 + ,1081 + ,807 + ,213118 + ,6282154 + ,1 + ,5626 + ,37 + ,56289 + ,0 + ,29790 + ,309 + ,444 + ,81767 + ,4321023 + ,2 + ,13337 + ,138 + ,28328 + ,0 + ,87550 + ,458 + ,412 + ,153198 + ,4111912 + ,3 + ,8541 + ,45 + ,17936 + ,1 + ,84738 + ,588 + ,428 + ,-26007 + ,223193 + ,415 + ,39 + ,0 + ,4145 + ,0 + ,54660 + ,302 + ,315 + ,126942 + ,1491348 + ,6 + ,4276 + ,24 + ,3040 + ,0 + ,42634 + ,156 + ,168 + ,157214 + ,1629616 + ,5 + ,9164 + ,34 + ,2964 + ,1 + ,40949 + ,481 + ,263 + ,129352 + ,1398893 + ,9 + ,2493 + ,29 + ,2865 + ,0 + ,45187 + ,353 + ,267 + ,234817 + ,1926517 + ,4 + ,4891 + ,38 + ,2854 + ,0 + ,37704 + ,452 + ,228 + ,60448 + ,983660 + ,14 + ,1734 + ,21 + ,2167 + ,0 + ,16275 + ,109 + ,129 + ,47818 + ,1443586 + ,7 + ,11409 + ,76 + ,1974 + ,1 + ,25830 + ,115 + ,104 + ,245546 + ,1073089 + ,10 + ,7592 + ,34 + ,1910 + ,1 + ,12679 + ,110 + ,122 + ,48020 + ,984885 + ,13 + ,7135 + ,62 + ,1871 + ,0 + ,18014 + ,239 + ,393 + ,-1710 + ,1405225 + ,8 + ,5043 + ,67 + ,943 + ,1 + ,43556 + ,247 + ,190 + ,32648 + ,227132 + ,414 + ,110 + ,1 + ,929 + ,0 + ,24811 + ,505 + ,280 + ,95350 + ,929118 + ,15 + ,1444 + ,29 + ,822 + ,1 + ,6575 + ,159 + ,63 + ,151352 + ,1071292 + ,11 + ,5480 + ,133 + ,819 + ,1 + ,7123 + ,109 + ,102 + ,288170 + ,638830 + ,28 + ,4026 + ,62 + ,769 + ,0 + ,21950 + ,519 + ,265 + ,114337 + ,856956 + ,17 + ,1266 + ,30 + ,745 + ,0 + ,37597 + ,248 + ,234 + ,37884 + ,992426 + ,12 + ,3195 + ,21 + ,652 + ,1 + ,17821 + ,373 + ,277 + ,122844 + ,444477 + ,46 + ,655 + ,14 + ,643 + ,0 + ,12988 + ,119 + ,73 + ,82340 + ,857217 + ,16 + ,5523 + ,51 + ,601 + ,0 + ,22330 + ,84 + ,67 + ,79801 + ,711969 + ,20 + ,6095 + ,23 + ,446 + ,1 + ,13326 + ,102 + ,103 + ,165548 + ,702380 + ,21 + ,4925 + ,38 + ,436 + ,1 + ,16189 + ,295 + ,290 + ,116384 + ,358589 + ,89 + ,538 + ,10 + ,379 + ,1 + ,7146 + ,105 + ,83 + ,134028 + ,297978 + ,367 + ,933 + ,14 + ,305 + ,1 + ,15824 + ,64 + ,56 + ,63838 + ,585715 + ,30 + ,6027 + ,24 + ,284 + ,0 + ,27664 + ,282 + ,236 + ,74996 + ,657954 + ,24 + ,1624 + ,17 + ,247 + ,1 + ,11920 + ,182 + ,73 + ,31080 + ,209458 + ,419 + ,52 + ,1 + ,238 + ,1 + ,8568 + ,37 + ,34 + ,32168 + ,786690 + ,18 + ,15856 + ,68 + ,223 + ,1 + ,14416 + ,361 + ,139 + ,49857 + ,439798 + ,48 + ,664 + ,17 + ,220 + ,0 + ,3369 + ,28 + ,26 + ,87161 + ,688779 + ,23 + ,17456 + ,145 + ,217 + ,0 + ,11819 + ,85 + ,70 + ,106113 + ,574339 + ,31 + ,4404 + ,32 + ,199 + ,0 + ,6984 + ,45 + ,40 + ,80570 + ,741409 + ,19 + ,12031 + ,78 + ,195 + ,0 + ,4519 + ,49 + ,42 + ,102129 + ,597793 + ,29 + ,8118 + ,88 + ,163 + ,1 + ,2220 + ,22 + ,12 + ,301670 + ,644190 + ,26 + ,20190 + ,200 + ,154 + ,1 + ,18562 + ,155 + ,211 + ,102313 + ,377934 + ,72 + ,1148 + ,10 + ,143 + ,1 + ,10327 + ,91 + ,74 + ,88577 + ,640273 + ,27 + ,4838 + ,43 + ,135 + ,0 + ,5336 + ,81 + ,80 + ,112477 + ,697458 + ,22 + ,6141 + ,93 + ,132 + ,0 + ,2365 + ,79 + ,83 + ,191778 + ,550608 + ,32 + ,4438 + ,148 + ,120 + ,1 + ,4069 + ,145 + ,131 + ,79804 + ,207393 + ,420 + ,51 + ,2 + ,119 + ,1 + ,8636 + ,855 + ,203 + ,128294 + ,301607 + ,359 + ,119 + ,12 + ,108 + ,1 + ,13718 + ,61 + ,56 + ,96448 + ,345783 + ,100 + ,2390 + ,11 + ,104 + ,1 + ,4525 + ,226 + ,89 + ,93811 + ,501749 + ,37 + ,1335 + ,67 + ,101 + ,1 + ,6869 + ,105 + ,88 + ,117520 + ,379983 + ,68 + ,1714 + ,26 + ,90 + ,1 + ,4628 + ,62 + ,39 + ,69159 + ,387475 + ,64 + ,3024 + ,41 + ,85 + ,0 + ,3689 + ,25 + ,25 + ,101792 + ,377305 + ,74 + ,7092 + ,48 + ,74 + ,0 + ,4891 + ,217 + ,49 + ,210568 + ,370837 + ,79 + ,787 + ,35 + ,73 + ,0 + ,7489 + ,322 + ,149 + ,136996 + ,430866 + ,50 + ,717 + ,31 + ,70 + ,1 + ,4901 + ,84 + ,58 + ,121920 + ,469107 + ,41 + ,3204 + ,55 + ,67 + ,1 + ,2284 + ,33 + ,41 + ,76403 + ,194493 + ,423 + ,-167 + ,-2 + ,66 + ,0 + ,3160 + ,108 + ,90 + ,108094 + ,530670 + ,33 + ,3062 + ,105 + ,66 + ,0 + ,4150 + ,150 + ,136 + ,134759 + ,518365 + ,35 + ,2122 + ,77 + ,66 + ,0 + ,7285 + ,115 + ,97 + ,188873 + ,491303 + ,39 + ,2533 + ,40 + ,59 + ,0 + ,1134 + ,162 + ,63 + ,146216 + ,527021 + ,34 + ,2019 + ,288 + ,58 + ,0 + ,4658 + ,158 + ,114 + ,156608 + ,233773 + ,411 + ,214 + ,7 + ,58 + ,1 + ,2384 + ,97 + ,77 + ,61348 + ,405972 + ,58 + ,2123 + ,86 + ,54 + ,1 + ,3748 + ,9 + ,6 + ,50350 + ,652925 + ,25 + ,50325 + ,121 + ,53 + ,1 + ,5371 + ,66 + ,47 + ,87720 + ,446211 + ,45 + ,3730 + ,46 + ,49 + ,1 + ,1285 + ,107 + ,51 + ,99489 + ,341340 + ,106 + ,1321 + ,110 + ,49 + ,0 + ,9327 + ,101 + ,85 + ,87419 + ,387699 + ,63 + ,1858 + ,20 + ,44 + ,0 + ,5565 + ,47 + ,43 + ,94355 + ,493408 + ,38 + ,6243 + ,53 + ,39 + ,1 + ,1528 + ,38 + ,32 + ,60326 + ,146494 + ,426 + ,-1408 + ,-35 + ,38 + ,0 + ,3122 + ,34 + ,25 + ,94670 + ,414462 + ,54 + ,6308 + ,69 + ,38 + ,0 + ,7561 + ,87 + ,77 + ,82425 + ,364304 + ,86 + ,1889 + ,22 + ,37 + ,1 + ,2675 + ,79 + ,54 + ,59017 + ,355178 + ,93 + ,1964 + ,58 + ,36 + ,1 + ,13253 + ,947 + ,251 + ,90829 + ,357760 + ,90 + ,167 + ,12 + ,35 + ,1 + ,880 + ,74 + ,15 + ,80791 + ,261216 + ,399 + ,827 + ,70 + ,34 + ,0 + ,2053 + ,53 + ,44 + ,100423 + ,397144 + ,60 + ,3720 + ,96 + ,33 + ,1 + ,1424 + ,94 + ,73 + ,131116 + ,374943 + ,77 + ,1861 + ,123 + ,32 + ,0 + ,4036 + ,63 + ,85 + ,100269 + ,424898 + ,53 + ,3570 + ,56 + ,32 + ,0 + ,3045 + ,58 + ,49 + ,27330 + ,202055 + ,422 + ,35 + ,1 + ,31 + ,1 + ,5119 + ,49 + ,38 + ,39039 + ,378525 + ,69 + ,3643 + ,35 + ,31 + ,1 + ,1431 + ,34 + ,35 + ,106885 + ,310768 + ,330 + ,3258 + ,77 + ,30 + ,1 + ,554 + ,11 + ,9 + ,79285 + ,325738 + ,125 + ,11431 + ,227 + ,30 + ,1 + ,1975 + ,35 + ,34 + ,118881 + ,394510 + ,61 + ,5557 + ,98 + ,30 + ,0 + ,1765 + ,20 + ,20 + ,77623 + ,247060 + ,404 + ,2353 + ,27 + ,28 + ,1 + ,1012 + ,47 + ,29 + ,114768 + ,368078 + ,82 + ,3576 + ,166 + ,26 + ,1 + ,810 + ,43 + ,11 + ,74015 + ,236761 + ,409 + ,855 + ,45 + ,26 + ,1 + ,1280 + ,117 + ,52 + ,69465 + ,312378 + ,326 + ,960 + ,88 + ,26 + ,0 + ,666 + ,171 + ,13 + ,117869 + ,339836 + ,108 + ,818 + ,210 + ,25 + ,1 + ,1380 + ,26 + ,29 + ,60982 + ,347385 + ,98 + ,5669 + ,107 + ,25 + ,0 + ,4677 + ,75 + ,66 + ,90131 + ,426280 + ,52 + ,3017 + ,48 + ,24 + ,1 + ,876 + ,59 + ,33 + ,138971 + ,352850 + ,96 + ,2591 + ,174 + ,24 + ,1 + ,814 + ,18 + ,15 + ,39625 + ,301881 + ,358 + ,5660 + ,125 + ,23 + ,1 + ,514 + ,15 + ,15 + ,102725 + ,377516 + ,73 + ,11834 + ,346 + ,23 + ,0 + ,5692 + ,72 + ,68 + ,64239 + ,357312 + ,91 + ,2185 + ,28 + ,22 + ,1 + ,3642 + ,86 + ,100 + ,90262 + ,458343 + ,42 + ,3004 + ,71 + ,22 + ,1 + ,540 + ,14 + ,13 + ,103960 + ,354228 + ,94 + ,11016 + ,286 + ,22 + ,1 + ,2099 + ,64 + ,45 + ,106611 + ,308636 + ,339 + ,1697 + ,52 + ,20 + ,1 + ,567 + ,11 + ,14 + ,103345 + ,386212 + ,65 + ,16928 + ,329 + ,20 + ,1 + ,2001 + ,52 + ,36 + ,95551 + ,393343 + ,62 + ,3718 + ,97 + ,19 + ,0 + ,2949 + ,41 + ,40 + ,82903 + ,378509 + ,70 + ,4354 + ,61 + ,18 + ,1 + ,2253 + ,99 + ,68 + ,63593 + ,452469 + ,43 + ,2550 + ,112 + ,17 + ,0 + ,6533 + ,75 + ,29 + ,126910 + ,364839 + ,85 + ,2198 + ,25 + ,16 + ,1 + ,1889 + ,45 + ,43 + ,37527 + ,358649 + ,88 + ,3526 + ,84 + ,14 + ,0 + ,3055 + ,43 + ,30 + ,60247 + ,376641 + ,75 + ,4108 + ,58 + ,14 + ,1 + ,272 + ,8 + ,9 + ,112995 + ,429112 + ,51 + ,28639 + ,844 + ,14 + ,0 + ,1414 + ,198 + ,22 + ,70184 + ,330546 + ,119 + ,659 + ,92 + ,13 + ,1 + ,2564 + ,22 + ,19 + ,130140 + ,403560 + ,59 + ,9253 + ,79 + ,13 + ,0 + ,1383 + ,11 + ,9 + ,73221 + ,317892 + ,151 + ,10717 + ,85 + ,13 + ,0 + ,1261 + ,33 + ,31 + ,76114 + ,307528 + ,343 + ,3258 + ,85 + ,13 + ,1 + ,975 + ,23 + ,19 + ,90534 + ,235133 + ,410 + ,1528 + ,36 + ,13 + ,1 + ,3366 + ,80 + ,55 + ,108479 + ,299243 + ,362 + ,1241 + ,29 + ,12 + ,1 + ,576 + ,18 + ,8 + ,113761 + ,314073 + ,313 + ,6337 + ,198 + ,12 + ,1 + ,1686 + ,40 + ,28 + ,68696 + ,368186 + ,81 + ,4205 + ,100 + ,11 + ,1 + ,746 + ,23 + ,29 + ,71561 + ,269661 + ,396 + ,3029 + ,93 + ,10 + ,0 + ,3192 + ,60 + ,48 + ,59831 + ,125390 + ,427 + ,-1243 + ,-23 + ,10 + ,0 + ,2045 + ,20 + ,16 + ,97890 + ,510834 + ,36 + ,15542 + ,152 + ,10 + ,1 + ,5702 + ,61 + ,47 + ,101481 + ,321896 + ,138 + ,1998 + ,21 + ,10 + ,0 + ,1932 + ,36 + ,20 + ,72954 + ,249898 + ,403 + ,1386 + ,26 + ,10 + ,1 + ,936 + ,30 + ,22 + ,67939 + ,408881 + ,56 + ,6963 + ,223 + ,9 + ,0 + ,3437 + ,47 + ,33 + ,48022 + ,158492 + ,425 + ,-883 + ,-12 + ,9 + ,1 + ,5131 + ,71 + ,44 + ,86111 + ,292154 + ,379 + ,1298 + ,18 + ,9 + ,0 + ,2397 + ,14 + ,13 + ,74020 + ,289513 + ,385 + ,6394 + ,37 + ,6 + ,0 + ,1389 + ,9 + ,6 + ,57530 + ,378049 + ,71 + ,19783 + ,128 + ,6 + ,1 + ,1503 + ,39 + ,35 + ,56364 + ,343466 + ,104 + ,3679 + ,95 + ,6 + ,1 + ,402 + ,26 + ,8 + ,84990 + ,332743 + ,115 + ,5106 + ,330 + ,5 + ,1 + ,2239 + ,21 + ,17 + ,88590 + ,442882 + ,47 + ,11566 + ,108 + ,5 + ,0 + ,2234 + ,16 + ,11 + ,77200 + ,214215 + ,418 + ,888 + ,6 + ,5 + ,1 + ,837 + ,69 + ,21 + ,61262 + ,315688 + ,165 + ,1677 + ,138 + ,4 + ,1 + ,10579 + ,92 + ,92 + ,110309 + ,375195 + ,76 + ,1904 + ,17 + ,4 + ,1 + ,875 + ,14 + ,12 + ,67000 + ,334280 + ,112 + ,9591 + ,153 + ,3 + ,1 + ,1585 + ,107 + ,112 + ,93099 + ,355864 + ,92 + ,1457 + ,98 + ,2 + ,0 + ,1659 + ,29 + ,25 + ,107577 + ,480382 + ,40 + ,9668 + ,169 + ,2 + ,0 + ,2647 + ,37 + ,17 + ,62920 + ,353058 + ,95 + ,4137 + ,58 + ,1 + ,0 + ,3294 + ,23 + ,23 + ,75832 + ,217193 + ,416 + ,748 + ,5 + ,1 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,191 + ,0 + ,1 + ,94 + ,7 + ,10 + ,60793 + ,314533 + ,308 + ,16362 + ,1212 + ,0 + ,1 + ,422 + ,28 + ,23 + ,57935 + ,318056 + ,150 + ,4216 + ,280 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,197 + ,0 + ,1 + ,34 + ,8 + ,7 + ,60630 + ,314353 + ,309 + ,14294 + ,3340 + ,0 + ,1 + ,1558 + ,63 + ,25 + ,55637 + ,369448 + ,80 + ,2690 + ,109 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,196 + ,0 + ,1 + ,43 + ,3 + ,20 + ,60887 + ,312846 + ,323 + ,37615 + ,2613 + ,0 + ,1 + ,645 + ,5 + ,4 + ,60720 + ,312075 + ,329 + ,22415 + ,174 + ,0 + ,1 + ,316 + ,9 + ,4 + ,60505 + ,315009 + ,299 + ,12779 + ,364 + ,0 + ,1 + ,115 + ,13 + ,10 + ,60945 + ,318903 + ,145 + ,9146 + ,1037 + ,0 + ,1 + ,5 + ,2 + ,1 + ,60720 + ,314887 + ,304 + ,57443 + ,24814 + ,0 + ,1 + ,897 + ,5 + ,4 + ,60720 + ,314913 + ,303 + ,22983 + ,128 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,188 + ,0 + ,1 + ,389 + ,14 + ,8 + ,58990 + ,325506 + ,127 + ,8965 + ,322 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,220 + ,0 + ,1 + ,1002 + ,15 + ,11 + ,56750 + ,298568 + ,364 + ,6571 + ,98 + ,0 + ,1 + ,36 + ,3 + ,4 + ,60894 + ,315834 + ,162 + ,38611 + ,3215 + ,0 + ,1 + ,460 + ,15 + ,15 + ,63346 + ,329784 + ,122 + ,8652 + ,282 + ,0 + ,1 + ,309 + ,11 + ,9 + ,56535 + ,312878 + ,322 + ,10262 + ,365 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,203 + ,0 + ,1 + ,9 + ,6 + ,7 + ,60835 + ,314987 + ,300 + ,19165 + ,12458 + ,0 + ,1 + ,271 + ,2 + ,2 + ,60720 + ,325249 + ,129 + ,62625 + ,463 + ,0 + ,1 + ,14 + ,1 + ,0 + ,61016 + ,315877 + ,161 + ,115877 + ,8041 + ,0 + ,1 + ,520 + ,10 + ,7 + ,58650 + ,291650 + ,383 + ,9165 + ,176 + ,0 + ,1 + ,1766 + ,73 + ,46 + ,60438 + ,305959 + ,349 + ,1451 + ,60 + ,0 + ,1 + ,0 + ,0 + ,5 + ,60720 + ,315380 + ,208 + ,0 + ,1 + ,458 + ,11 + ,7 + ,58625 + ,297765 + ,368 + ,8888 + ,213 + ,0 + ,1 + ,20 + ,3 + ,2 + ,60938 + ,315245 + ,297 + ,38415 + ,5674 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,217 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,219 + ,0 + ,1 + ,98 + ,2 + ,2 + ,61490 + ,315236 + ,298 + ,57618 + ,1179 + ,0 + ,1 + ,405 + ,7 + ,5 + ,60845 + ,336425 + ,111 + ,19489 + ,337 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,213 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,212 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,207 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,218 + ,0 + ,1 + ,483 + ,27 + ,7 + ,60830 + ,306268 + ,348 + ,3936 + ,220 + ,0 + ,1 + ,454 + ,51 + ,24 + ,63261 + ,302187 + ,357 + ,2004 + ,225 + ,0 + ,1 + ,47 + ,3 + ,1 + ,60720 + ,314882 + ,305 + ,38294 + ,2432 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,193 + ,0 + ,1 + ,757 + ,19 + ,18 + ,45689 + ,382712 + ,66 + ,9616 + ,242 + ,0 + ,1 + ,4655 + ,393 + ,55 + ,60720 + ,341570 + ,105 + ,360 + ,30 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,192 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,211 + ,0 + ,1 + ,36 + ,4 + ,3 + ,61564 + ,312412 + ,325 + ,28103 + ,3146 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,189 + ,0 + ,1 + ,203 + ,9 + ,9 + ,61938 + ,309596 + ,334 + ,12177 + ,541 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,205 + ,0 + ,1 + ,126 + ,10 + ,8 + ,60951 + ,315547 + ,171 + ,11555 + ,921 + ,0 + ,1 + ,400 + ,152 + ,113 + ,60720 + ,313267 + ,319 + ,745 + ,283 + ,0 + ,0 + ,71 + ,1 + ,0 + ,60745 + ,316176 + ,159 + ,116176 + ,1646 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,209 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,190 + ,0 + ,1 + ,972 + ,34 + ,19 + ,71642 + ,359335 + ,87 + ,4686 + ,164 + ,0 + ,0 + ,531 + ,10 + ,11 + ,71641 + ,330068 + ,120 + ,13007 + ,245 + ,0 + ,1 + ,2461 + ,57 + ,25 + ,55792 + ,314289 + ,310 + ,2005 + ,46 + ,0 + ,0 + ,378 + ,52 + ,16 + ,71873 + ,297413 + ,369 + ,1873 + ,258 + ,0 + ,0 + ,23 + ,5 + ,5 + ,62555 + ,314806 + ,306 + ,22961 + ,4974 + ,0 + ,0 + ,638 + ,14 + ,11 + ,60370 + ,333210 + ,114 + ,9515 + ,209 + ,0 + ,0 + ,2300 + ,29 + ,23 + ,64873 + ,352108 + ,97 + ,5245 + ,66 + ,0 + ,1 + ,149 + ,5 + ,6 + ,62041 + ,313332 + ,318 + ,22666 + ,761 + ,0 + ,1 + ,226 + ,5 + ,5 + ,65745 + ,291787 + ,382 + ,18357 + ,406 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,184 + ,0 + ,1 + ,275 + ,4 + ,7 + ,59500 + ,318745 + ,146 + ,29686 + ,431 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,179 + ,0 + ,1 + ,141 + ,6 + ,7 + ,61630 + ,315366 + ,295 + ,19228 + ,819 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,186 + ,0 + ,1 + ,28 + ,2 + ,3 + ,60890 + ,315688 + ,166 + ,57844 + ,4148 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,263 + ,0 + ,0 + ,4980 + ,91 + ,89 + ,113521 + ,409642 + ,55 + ,2304 + ,42 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,259 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,261 + ,0 + ,0 + ,472 + ,20 + ,19 + ,80045 + ,269587 + ,397 + ,3479 + ,148 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,243 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,250 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,247 + ,0 + ,0 + ,203 + ,27 + ,12 + ,50804 + ,300962 + ,361 + ,3739 + ,498 + ,0 + ,0 + ,496 + ,17 + ,12 + ,87390 + ,325479 + ,128 + ,7381 + ,253 + ,0 + ,0 + ,10 + ,2 + ,5 + ,61656 + ,316155 + ,160 + ,58077 + ,11639 + ,0 + ,0 + ,63 + ,4 + ,2 + ,65688 + ,318574 + ,148 + ,29643 + ,1891 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,228 + ,0 + ,0 + ,1136 + ,32 + ,26 + ,48522 + ,343613 + ,103 + ,4488 + ,126 + ,0 + ,0 + ,265 + ,31 + ,3 + ,60720 + ,306948 + ,346 + ,3450 + ,404 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,210 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,206 + ,0 + ,1 + ,267 + ,32 + ,11 + ,57640 + ,330059 + ,121 + ,4064 + ,487 + ,0 + ,1 + ,474 + ,20 + ,10 + ,61977 + ,288985 + ,386 + ,4449 + ,188 + ,0 + ,1 + ,534 + ,7 + ,5 + ,62620 + ,304485 + ,351 + ,14926 + ,196 + ,0 + ,1 + ,0 + ,0 + ,2 + ,60720 + ,315380 + ,201 + ,0 + ,1 + ,15 + ,8 + ,6 + ,60831 + ,315688 + ,164 + ,14461 + ,7687 + ,0 + ,1 + ,397 + ,28 + ,7 + ,60646 + ,317736 + ,152 + ,4205 + ,297 + ,0 + ,1 + ,0 + ,0 + ,2 + ,60720 + ,315380 + ,199 + ,0 + ,1 + ,1866 + ,29 + ,28 + ,56225 + ,322331 + ,135 + ,4218 + ,66 + ,0 + ,1 + ,288 + ,4 + ,3 + ,60510 + ,296656 + ,373 + ,24164 + ,335 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,198 + ,0 + ,1 + ,3 + ,2 + ,1 + ,60698 + ,315354 + ,296 + ,57677 + ,36505 + ,0 + ,1 + ,468 + ,21 + ,20 + ,60720 + ,312161 + ,328 + ,5341 + ,239 + ,0 + ,1 + ,20 + ,2 + ,1 + ,60805 + ,315576 + ,169 + ,57788 + ,5924 + ,0 + ,1 + ,278 + ,26 + ,22 + ,61404 + ,314922 + ,302 + ,4420 + ,413 + ,0 + ,1 + ,61 + ,14 + ,9 + ,60720 + ,314551 + ,307 + ,8182 + ,1874 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,215 + ,0 + ,1 + ,192 + ,4 + ,2 + ,65276 + ,312339 + ,327 + ,28085 + ,586 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,214 + ,0 + ,1 + ,317 + ,9 + ,7 + ,63915 + ,298700 + ,363 + ,10967 + ,311 + ,0 + ,1 + ,738 + ,10 + ,9 + ,60720 + ,321376 + ,139 + ,12138 + ,164 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,216 + ,0 + ,1 + ,368 + ,17 + ,13 + ,61686 + ,303230 + ,355 + ,6072 + ,281 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,204 + ,0 + ,1 + ,2 + ,1 + ,0 + ,60743 + ,315487 + ,172 + ,115487 + ,52257 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,202 + ,0 + ,1 + ,53 + ,6 + ,6 + ,60349 + ,315793 + ,163 + ,19299 + ,2197 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,195 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,194 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,200 + ,0 + ,1 + ,94 + ,3 + ,3 + ,61360 + ,312887 + ,321 + ,37629 + ,1200 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,237 + ,0 + ,1 + ,24 + ,8 + ,7 + ,59818 + ,315637 + ,168 + ,14455 + ,4772 + ,0 + ,1 + ,2332 + ,4 + ,2 + ,72680 + ,324385 + ,131 + ,31096 + ,53 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,236 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,233 + ,0 + ,1 + ,131 + ,11 + ,15 + ,61808 + ,308989 + ,338 + ,9908 + ,834 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,260 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,234 + ,0 + ,1 + ,206 + ,9 + ,9 + ,53110 + ,296702 + ,372 + ,10745 + ,469 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,240 + ,0 + ,1 + ,167 + ,2 + ,1 + ,64245 + ,307322 + ,345 + ,53661 + ,642 + ,0 + ,1 + ,622 + ,73 + ,38 + ,73007 + ,304376 + ,352 + ,1430 + ,168 + ,0 + ,1 + ,2328 + ,94 + ,57 + ,82732 + ,253588 + ,400 + ,570 + ,23 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,238 + ,0 + ,1 + ,365 + ,8 + ,7 + ,54820 + ,309560 + ,335 + ,13695 + ,300 + ,0 + ,1 + ,364 + ,35 + ,26 + ,47705 + ,298466 + ,365 + ,2813 + ,271 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,239 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,235 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,262 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,283 + ,0 + ,1 + ,226 + ,12 + ,13 + ,72835 + ,343929 + ,102 + ,11994 + ,636 + ,0 + ,1 + ,307 + ,15 + ,10 + ,58856 + ,331955 + ,116 + ,8797 + ,430 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,288 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,281 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,282 + ,0 + ,1 + ,188 + ,11 + ,9 + ,77655 + ,381180 + ,67 + ,16471 + ,964 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,221 + ,0 + ,1 + ,138 + ,6 + ,26 + ,69817 + ,331420 + ,117 + ,21903 + ,954 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,286 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,285 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,265 + ,0 + ,1 + ,125 + ,12 + ,19 + ,60798 + ,310201 + ,332 + ,9183 + ,883 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,284 + ,0 + ,1 + ,282 + ,30 + ,12 + ,62452 + ,320016 + ,141 + ,4001 + ,425 + ,0 + ,1 + ,335 + ,33 + ,23 + ,64175 + ,320398 + ,140 + ,3648 + ,359 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,232 + ,0 + ,0 + ,1324 + ,117 + ,29 + ,67440 + ,291841 + ,381 + ,785 + ,69 + ,0 + ,1 + ,176 + ,28 + ,8 + ,68136 + ,310670 + ,331 + ,3952 + ,629 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,256 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,255 + ,0 + ,1 + ,249 + ,72 + ,26 + ,56726 + ,313491 + ,317 + ,1576 + ,455 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,242 + ,0 + ,1 + ,333 + ,13 + ,9 + ,70811 + ,331323 + ,118 + ,10102 + ,394 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,241 + ,0 + ,0 + ,601 + ,6 + ,5 + ,60720 + ,319210 + ,144 + ,19868 + ,198 + ,0 + ,1 + ,30 + ,4 + ,3 + ,62045 + ,318098 + ,149 + ,29525 + ,3912 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,254 + ,0 + ,1 + ,249 + ,62 + ,13 + ,54323 + ,292754 + ,378 + ,1496 + ,373 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,258 + ,0 + ,1 + ,165 + ,24 + ,12 + ,62841 + ,325176 + ,130 + ,5216 + ,759 + ,0 + ,1 + ,453 + ,21 + ,19 + ,81125 + ,365959 + ,84 + ,7903 + ,366 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,257 + ,0 + ,1 + ,53 + ,14 + ,10 + ,59506 + ,302409 + ,356 + ,7315 + ,1932 + ,0 + ,0 + ,382 + ,21 + ,9 + ,59365 + ,340968 + ,107 + ,6713 + ,369 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,229 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,230 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,224 + ,0 + ,1 + ,0 + ,0 + ,9 + ,60720 + ,315380 + ,223 + ,0 + ,0 + ,30 + ,4 + ,4 + ,60798 + ,313164 + ,320 + ,28291 + ,3808 + ,0 + ,1 + ,290 + ,2 + ,1 + ,58790 + ,301164 + ,360 + ,50582 + ,348 + ,0 + ,1 + ,0 + ,0 + ,1 + ,60720 + ,315380 + ,225 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,231 + ,0 + ,1 + ,366 + ,53 + ,14 + ,61808 + ,344425 + ,101 + ,2725 + ,394 + ,0 + ,1 + ,2 + ,9 + ,12 + ,60735 + ,315394 + ,174 + ,12822 + ,63057 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,222 + ,0 + ,0 + ,209 + ,13 + ,19 + ,64016 + ,316647 + ,157 + ,8973 + ,558 + ,0 + ,1 + ,384 + ,22 + ,17 + ,54683 + ,309836 + ,333 + ,4993 + ,286 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,227 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,244 + ,0 + ,1 + ,365 + ,83 + ,32 + ,87192 + ,346611 + ,99 + ,1766 + ,402 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,226 + ,0 + ,0 + ,49 + ,8 + ,14 + ,64107 + ,322031 + ,137 + ,15254 + ,2500 + ,0 + ,1 + ,3 + ,4 + ,8 + ,60761 + ,315656 + ,167 + ,28914 + ,37920 + ,0 + ,1 + ,133 + ,14 + ,4 + ,65990 + ,339445 + ,109 + ,9960 + ,1052 + ,0 + ,1 + ,32 + ,1 + ,0 + ,59988 + ,314964 + ,301 + ,114964 + ,3603 + ,0 + ,1 + ,368 + ,17 + ,20 + ,61167 + ,297141 + ,370 + ,5714 + ,264 + ,0 + ,1 + ,1 + ,6 + ,5 + ,60719 + ,315372 + ,293 + ,19229 + ,84213 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,274 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,246 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,177 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,185 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,245 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,273 + ,0 + ,1 + ,22 + ,2 + ,1 + ,60722 + ,312502 + ,324 + ,56251 + ,5154 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,252 + ,0 + ,0 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,249 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,253 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,277 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,278 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,275 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,276 + ,0 + ,1 + ,96 + ,5 + ,4 + ,60379 + ,313729 + ,316 + ,22746 + ,1188 + ,0 + ,1 + ,1 + ,2 + ,1 + ,60727 + ,315388 + ,175 + ,57694 + ,117743 + ,0 + ,1 + ,314 + ,5 + ,4 + ,60720 + ,315371 + ,294 + ,23074 + ,367 + ,0 + ,1 + ,844 + ,78 + ,20 + ,60925 + ,296139 + ,375 + ,1233 + ,114 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,180 + ,0 + ,1 + ,26 + ,1 + ,1 + ,60896 + ,313880 + ,315 + ,113880 + ,4328 + ,0 + ,1 + ,125 + ,13 + ,10 + ,59734 + ,317698 + ,153 + ,9054 + ,944 + ,0 + ,1 + ,304 + ,15 + ,12 + ,62969 + ,295580 + ,376 + ,6372 + ,314 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,187 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,183 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,182 + ,0 + ,1 + ,621 + ,48 + ,13 + ,60720 + ,308256 + ,341 + ,2255 + ,174 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,181 + ,0 + ,1 + ,119 + ,6 + ,3 + ,59118 + ,303677 + ,354 + ,17280 + ,872 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,178 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,176 + ,0 + ,1 + ,1595 + ,17 + ,10 + ,60720 + ,319369 + ,143 + ,7022 + ,75 + ,0 + ,1 + ,312 + ,14 + ,3 + ,58598 + ,318690 + ,147 + ,8478 + ,381 + ,0 + ,1 + ,60 + ,10 + ,7 + ,61124 + ,314049 + ,314 + ,11405 + ,1916 + ,0 + ,1 + ,587 + ,12 + ,10 + ,59595 + ,325699 + ,126 + ,10475 + ,214 + ,0 + ,1 + ,135 + ,2 + ,1 + ,62065 + ,314210 + ,311 + ,57105 + ,847 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,272 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,269 + ,0 + ,1 + ,514 + ,52 + ,15 + ,78780 + ,322378 + ,133 + ,2353 + ,238 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,279 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,251 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,270 + ,0 + ,1 + ,1 + ,4 + ,4 + ,60722 + ,315398 + ,173 + ,28849 + ,172235 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,268 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,289 + ,0 + ,0 + ,1763 + ,24 + ,28 + ,61600 + ,308336 + ,340 + ,4514 + ,61 + ,0 + ,1 + ,180 + ,11 + ,9 + ,59635 + ,316386 + ,158 + ,10581 + ,647 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,287 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,248 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,291 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,292 + ,0 + ,1 + ,218 + ,21 + ,7 + ,60720 + ,315553 + ,170 + ,5503 + ,531 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,290 + ,0 + ,1 + ,448 + ,40 + ,7 + ,59781 + ,323361 + ,132 + ,3084 + ,275 + ,0 + ,1 + ,227 + ,9 + ,7 + ,76644 + ,336639 + ,110 + ,15182 + ,603 + ,0 + ,1 + ,174 + ,1 + ,3 + ,64820 + ,307424 + ,344 + ,107424 + ,618 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,266 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,280 + ,0 + ,1 + ,121 + ,24 + ,11 + ,56178 + ,295370 + ,377 + ,3974 + ,788 + ,0 + ,1 + ,607 + ,11 + ,7 + ,60436 + ,322340 + ,134 + ,11122 + ,201 + ,0 + ,1 + ,2212 + ,14 + ,10 + ,60720 + ,319864 + ,142 + ,8562 + ,54 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,271 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,264 + ,0 + ,1 + ,530 + ,60 + ,18 + ,73433 + ,317291 + ,155 + ,1955 + ,221 + ,0 + ,1 + ,571 + ,80 + ,14 + ,41477 + ,280398 + ,390 + ,1005 + ,141 + ,0 + ,1 + ,0 + ,0 + ,0 + ,60720 + ,315380 + ,267 + ,0 + ,1 + ,78 + ,16 + ,12 + ,62700 + ,317330 + ,154 + ,7333 + ,1511 + ,0 + ,1 + ,2489 + ,40 + ,29 + ,67804 + ,238125 + ,408 + ,953 + ,15 + ,0 + ,1 + ,131 + ,6 + ,3 + ,59661 + ,327071 + ,123 + ,21179 + ,969 + ,0 + ,1 + ,923 + ,8 + ,6 + ,58620 + ,309038 + ,337 + ,13630 + ,118 + ,0 + ,1 + ,72 + ,3 + ,3 + ,60398 + ,314210 + ,312 + ,38070 + ,1591 + ,0 + ,1 + ,572 + ,16 + ,8 + ,58580 + ,307930 + ,342 + ,6746 + ,189 + ,-1 + ,1 + ,397 + ,10 + ,10 + ,62710 + ,322327 + ,136 + ,12233 + ,308 + ,-1 + ,1 + ,450 + ,8 + ,6 + ,59325 + ,292136 + ,380 + ,11517 + ,205 + ,-1 + ,1 + ,622 + ,7 + ,8 + ,60950 + ,263276 + ,398 + ,9039 + ,102 + ,-2 + ,1 + ,694 + ,8 + ,6 + ,68060 + ,367655 + ,83 + ,20957 + ,241 + ,-2 + ,0 + ,3425 + ,12 + ,9 + ,83620 + ,283910 + ,388 + ,6993 + ,25 + ,-3 + ,1 + ,562 + ,13 + ,8 + ,58456 + ,283587 + ,389 + ,6430 + ,149 + ,-3 + ,1 + ,4917 + ,42 + ,26 + ,52811 + ,243650 + ,406 + ,1039 + ,9 + ,-3 + ,0 + ,1442 + ,118 + ,239 + ,121173 + ,438493 + ,49 + ,2021 + ,165 + ,-3 + ,1 + ,529 + ,9 + ,7 + ,63870 + ,296261 + ,374 + ,10696 + ,182 + ,-3 + ,0 + ,2126 + ,138 + ,41 + ,21001 + ,230621 + ,413 + ,222 + ,14 + ,-3 + ,1 + ,1061 + ,5 + ,3 + ,70415 + ,304252 + ,353 + ,20850 + ,98 + ,-4 + ,1 + ,776 + ,9 + ,8 + ,64230 + ,333505 + ,113 + ,14834 + ,172 + ,-4 + ,1 + ,611 + ,8 + ,6 + ,59190 + ,296919 + ,371 + ,12115 + ,159 + ,-4 + ,0 + ,1526 + ,25 + ,21 + ,69351 + ,278990 + ,391 + ,3160 + ,52 + ,-4 + ,1 + ,592 + ,7 + ,7 + ,64270 + ,276898 + ,392 + ,10985 + ,130 + ,-4 + ,1 + ,1182 + ,13 + ,11 + ,70694 + ,327007 + ,124 + ,9770 + ,107 + ,-5 + ,1 + ,621 + ,16 + ,11 + ,68005 + ,317046 + ,156 + ,7315 + ,188 + ,-5 + ,1 + ,989 + ,11 + ,12 + ,58930 + ,304555 + ,350 + ,9505 + ,106 + ,-5 + ,1 + ,438 + ,11 + ,9 + ,58320 + ,298096 + ,366 + ,8918 + ,224 + ,-5 + ,1 + ,726 + ,3 + ,3 + ,69980 + ,231861 + ,412 + ,10620 + ,44 + ,-5 + ,1 + ,1303 + ,61 + ,57 + ,69863 + ,309422 + ,336 + ,1794 + ,84 + ,-5 + ,0 + ,7419 + ,29 + ,21 + ,63255 + ,286963 + ,387 + ,2999 + ,12 + ,-5 + ,0 + ,1164 + ,17 + ,15 + ,57320 + ,269753 + ,395 + ,4103 + ,60 + ,-6 + ,0 + ,3310 + ,33 + ,32 + ,75230 + ,448243 + ,44 + ,7523 + ,75 + ,-6 + ,1 + ,1920 + ,15 + ,11 + ,79420 + ,165404 + ,424 + ,-2306 + ,-18 + ,-6 + ,1 + ,965 + ,3 + ,2 + ,73490 + ,204325 + ,421 + ,1442 + ,4 + ,-8 + ,1 + ,3256 + ,66 + ,23 + ,35250 + ,407159 + ,57 + ,3139 + ,64 + ,-8 + ,0 + ,1135 + ,17 + ,20 + ,62285 + ,290476 + ,384 + ,5322 + ,80 + ,-9 + ,1 + ,1270 + ,26 + ,24 + ,69206 + ,275311 + ,393 + ,2897 + ,59 + ,-10 + ,1 + ,661 + ,3 + ,1 + ,65920 + ,246541 + ,405 + ,15514 + ,70 + ,-12 + ,1 + ,1013 + ,2 + ,1 + ,69770 + ,253468 + ,402 + ,26734 + ,53 + ,-14 + ,1 + ,2844 + ,67 + ,74 + ,72683 + ,240897 + ,407 + ,610 + ,14 + ,-15 + ,0 + ,11528 + ,70 + ,68 + ,-14545 + ,-83265 + ,431 + ,-4047 + ,-25 + ,-25 + ,1 + ,6526 + ,26 + ,20 + ,55830 + ,-42143 + ,430 + ,-9313 + ,-37 + ,-28 + ,1 + ,2264 + ,24 + ,20 + ,55174 + ,272713 + ,394 + ,3030 + ,32 + ,-31 + ,0 + ,5109 + ,97 + ,82 + ,67038 + ,215362 + ,417 + ,158 + ,3 + ,-31 + ,1 + ,3999 + ,30 + ,21 + ,51252 + ,42754 + ,428 + ,-5242 + ,-39 + ,-58 + ,1 + ,35624 + ,223 + ,244 + ,157278 + ,306275 + ,347 + ,477 + ,3 + ,-70 + ,1 + ,9252 + ,48 + ,32 + ,79510 + ,253537 + ,401 + ,1115 + ,6 + ,-129 + ,1 + ,15236 + ,90 + ,86 + ,77440 + ,372631 + ,78 + ,1918 + ,11 + ,-158 + ,1 + ,18073 + ,180 + ,69 + ,27284 + ,-7170 + ,429 + ,-1151 + ,-11 + ,-366) + ,dim=c(10 + ,431) + ,dimnames=list(c('groep' + ,'Costs' + ,'Trades' + ,'Orders' + ,'Dividends' + ,'Wealth' + ,'WRank' + ,'Profit/Trades' + ,'Profit/Cost' + ,'CCScore ') + ,1:431)) > y <- array(NA,dim=c(10,431),dimnames=list(c('groep','Costs','Trades','Orders','Dividends','Wealth','WRank','Profit/Trades','Profit/Cost','CCScore '),1:431)) > for (i in 1:dim(x)[1]) + { + for (j in 1:dim(x)[2]) + { + y[i,j] <- as.numeric(x[i,j]) + } + } > par20 = '' > par19 = '' > par18 = '' > par17 = '' > par16 = '' > par15 = '' > par14 = '' > par13 = '' > par12 = '' > par11 = '' > par10 = '' > par9 = '' > par8 = '' > par7 = '' > par6 = '' > par5 = '' > par4 = '' > par3 = 'No Linear Trend' > par2 = 'Do not include Seasonal Dummies' > par1 = '10' > ylab = '' > xlab = '' > main = '' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > library(lattice) > library(lmtest) Loading required package: zoo 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 CCScore\r groep Costs Trades Orders Dividends Wealth WRank 1 56289 0 162556 1081 807 213118 6282154 1 2 28328 0 29790 309 444 81767 4321023 2 3 17936 0 87550 458 412 153198 4111912 3 4 4145 1 84738 588 428 -26007 223193 415 5 3040 0 54660 302 315 126942 1491348 6 6 2964 0 42634 156 168 157214 1629616 5 7 2865 1 40949 481 263 129352 1398893 9 8 2854 0 45187 353 267 234817 1926517 4 9 2167 0 37704 452 228 60448 983660 14 10 1974 0 16275 109 129 47818 1443586 7 11 1910 1 25830 115 104 245546 1073089 10 12 1871 1 12679 110 122 48020 984885 13 13 943 0 18014 239 393 -1710 1405225 8 14 929 1 43556 247 190 32648 227132 414 15 822 0 24811 505 280 95350 929118 15 16 819 1 6575 159 63 151352 1071292 11 17 769 1 7123 109 102 288170 638830 28 18 745 0 21950 519 265 114337 856956 17 19 652 0 37597 248 234 37884 992426 12 20 643 1 17821 373 277 122844 444477 46 21 601 0 12988 119 73 82340 857217 16 22 446 0 22330 84 67 79801 711969 20 23 436 1 13326 102 103 165548 702380 21 24 379 1 16189 295 290 116384 358589 89 25 305 1 7146 105 83 134028 297978 367 26 284 1 15824 64 56 63838 585715 30 27 247 0 27664 282 236 74996 657954 24 28 238 1 11920 182 73 31080 209458 419 29 223 1 8568 37 34 32168 786690 18 30 220 1 14416 361 139 49857 439798 48 31 217 0 3369 28 26 87161 688779 23 32 199 0 11819 85 70 106113 574339 31 33 195 0 6984 45 40 80570 741409 19 34 163 0 4519 49 42 102129 597793 29 35 154 1 2220 22 12 301670 644190 26 36 143 1 18562 155 211 102313 377934 72 37 135 1 10327 91 74 88577 640273 27 38 132 0 5336 81 80 112477 697458 22 39 120 0 2365 79 83 191778 550608 32 40 119 1 4069 145 131 79804 207393 420 41 108 1 8636 855 203 128294 301607 359 42 104 1 13718 61 56 96448 345783 100 43 101 1 4525 226 89 93811 501749 37 44 90 1 6869 105 88 117520 379983 68 45 85 1 4628 62 39 69159 387475 64 46 74 0 3689 25 25 101792 377305 74 47 73 0 4891 217 49 210568 370837 79 48 70 0 7489 322 149 136996 430866 50 49 67 1 4901 84 58 121920 469107 41 50 66 1 2284 33 41 76403 194493 423 51 66 0 3160 108 90 108094 530670 33 52 66 0 4150 150 136 134759 518365 35 53 59 0 7285 115 97 188873 491303 39 54 58 0 1134 162 63 146216 527021 34 55 58 0 4658 158 114 156608 233773 411 56 54 1 2384 97 77 61348 405972 58 57 53 1 3748 9 6 50350 652925 25 58 49 1 5371 66 47 87720 446211 45 59 49 1 1285 107 51 99489 341340 106 60 44 0 9327 101 85 87419 387699 63 61 39 0 5565 47 43 94355 493408 38 62 38 1 1528 38 32 60326 146494 426 63 38 0 3122 34 25 94670 414462 54 64 37 0 7561 87 77 82425 364304 86 65 36 1 2675 79 54 59017 355178 93 66 35 1 13253 947 251 90829 357760 90 67 34 1 880 74 15 80791 261216 399 68 33 0 2053 53 44 100423 397144 60 69 32 1 1424 94 73 131116 374943 77 70 32 0 4036 63 85 100269 424898 53 71 31 0 3045 58 49 27330 202055 422 72 31 1 5119 49 38 39039 378525 69 73 30 1 1431 34 35 106885 310768 330 74 30 1 554 11 9 79285 325738 125 75 30 1 1975 35 34 118881 394510 61 76 28 0 1765 20 20 77623 247060 404 77 26 1 1012 47 29 114768 368078 82 78 26 1 810 43 11 74015 236761 409 79 26 1 1280 117 52 69465 312378 326 80 25 0 666 171 13 117869 339836 108 81 25 1 1380 26 29 60982 347385 98 82 24 0 4677 75 66 90131 426280 52 83 24 1 876 59 33 138971 352850 96 84 23 1 814 18 15 39625 301881 358 85 23 1 514 15 15 102725 377516 73 86 22 0 5692 72 68 64239 357312 91 87 22 1 3642 86 100 90262 458343 42 88 22 1 540 14 13 103960 354228 94 89 20 1 2099 64 45 106611 308636 339 90 20 1 567 11 14 103345 386212 65 91 19 1 2001 52 36 95551 393343 62 92 18 0 2949 41 40 82903 378509 70 93 17 1 2253 99 68 63593 452469 43 94 16 0 6533 75 29 126910 364839 85 95 14 1 1889 45 43 37527 358649 88 96 14 0 3055 43 30 60247 376641 75 97 14 1 272 8 9 112995 429112 51 98 13 0 1414 198 22 70184 330546 119 99 13 1 2564 22 19 130140 403560 59 100 13 0 1383 11 9 73221 317892 151 101 13 0 1261 33 31 76114 307528 343 102 13 1 975 23 19 90534 235133 410 103 12 1 3366 80 55 108479 299243 362 104 12 1 576 18 8 113761 314073 313 105 11 1 1686 40 28 68696 368186 81 106 10 1 746 23 29 71561 269661 396 107 10 0 3192 60 48 59831 125390 427 108 10 0 2045 20 16 97890 510834 36 109 10 1 5702 61 47 101481 321896 138 110 10 0 1932 36 20 72954 249898 403 111 9 1 936 30 22 67939 408881 56 112 9 0 3437 47 33 48022 158492 425 113 9 1 5131 71 44 86111 292154 379 114 6 0 2397 14 13 74020 289513 385 115 6 0 1389 9 6 57530 378049 71 116 6 1 1503 39 35 56364 343466 104 117 5 1 402 26 8 84990 332743 115 118 5 1 2239 21 17 88590 442882 47 119 5 0 2234 16 11 77200 214215 418 120 4 1 837 69 21 61262 315688 165 121 4 1 10579 92 92 110309 375195 76 122 3 1 875 14 12 67000 334280 112 123 2 1 1585 107 112 93099 355864 92 124 2 0 1659 29 25 107577 480382 40 125 1 0 2647 37 17 62920 353058 95 126 1 0 3294 23 23 75832 217193 416 127 94 1 0 0 0 60720 315380 191 128 422 7 10 60793 314533 308 16362 1212 129 0 28 23 57935 318056 150 4216 280 130 7 0 0 60720 315380 197 0 1 131 25 60630 314353 309 14294 3340 0 1 132 0 55637 369448 80 2690 109 0 1 133 312846 60720 315380 196 0 1 43 3 134 312075 323 37615 2613 0 1 645 5 135 315009 329 22415 174 0 1 316 9 136 318903 299 12779 364 0 1 115 13 137 314887 145 9146 1037 0 1 5 2 138 314913 304 57443 24814 0 1 897 5 139 315380 303 22983 128 0 1 0 0 140 8965 188 0 1 389 14 8 58990 141 0 322 0 1 0 0 0 60720 142 0 1 1002 15 11 56750 298568 364 143 0 1 36 3 4 60894 315834 162 144 0 1 460 15 15 63346 329784 122 145 0 1 309 11 9 56535 312878 322 146 9 1 0 0 0 60720 315380 203 147 271 6 7 60835 314987 300 19165 12458 148 14 2 2 60720 325249 129 62625 463 149 520 1 0 61016 315877 161 115877 8041 150 1766 10 7 58650 291650 383 9165 176 151 0 73 46 60438 305959 349 1451 60 152 7 0 5 60720 315380 208 0 1 153 2 58625 297765 368 8888 213 0 1 154 0 60938 315245 297 38415 5674 0 1 155 315380 60720 315380 217 0 1 0 0 156 57618 219 0 1 98 2 2 61490 157 19489 1179 0 1 405 7 5 60845 158 0 337 0 1 0 0 0 60720 159 0 1 0 0 0 60720 315380 212 160 0 0 0 60720 315380 207 0 1 161 306268 60720 315380 218 0 1 483 27 162 302187 348 3936 220 0 1 454 51 163 314882 357 2004 225 0 1 47 3 164 315380 305 38294 2432 0 1 0 0 165 9616 193 0 1 757 19 18 45689 166 360 242 0 1 4655 393 55 60720 167 0 30 0 1 0 0 0 60720 168 36 1 0 0 0 60720 315380 211 169 0 4 3 61564 312412 325 28103 3146 170 9 0 0 60720 315380 189 0 1 171 0 61938 309596 334 12177 541 0 1 172 315547 60720 315380 205 0 1 126 10 173 313267 171 11555 921 0 1 400 152 174 316176 319 745 283 0 0 71 1 175 315380 159 116176 1646 0 1 0 0 176 0 209 0 1 0 0 0 60720 177 0 1 972 34 19 71642 359335 87 178 0 0 531 10 11 71641 330068 120 179 0 1 2461 57 25 55792 314289 310 180 0 0 378 52 16 71873 297413 369 181 0 0 23 5 5 62555 314806 306 182 0 0 638 14 11 60370 333210 114 183 0 0 2300 29 23 64873 352108 97 184 0 1 149 5 6 62041 313332 318 185 0 1 226 5 5 65745 291787 382 186 275 1 0 0 0 60720 315380 184 187 0 4 7 59500 318745 146 29686 431 188 7 0 0 60720 315380 179 0 1 189 0 61630 315366 295 19228 819 0 1 190 315688 60720 315380 186 0 1 28 2 191 315380 166 57844 4148 0 0 0 0 192 2304 263 0 0 4980 91 89 113521 193 0 42 0 0 0 0 0 60720 194 472 0 0 0 0 60720 315380 261 195 0 20 19 80045 269587 397 3479 148 196 0 0 0 60720 315380 243 0 0 197 315380 60720 315380 250 0 0 0 0 198 3739 247 0 0 203 27 12 50804 199 7381 498 0 0 496 17 12 87390 200 58077 253 0 0 10 2 5 61656 201 29643 11639 0 0 63 4 2 65688 202 0 1891 0 0 0 0 0 60720 203 0 0 1136 32 26 48522 343613 103 204 1 0 265 31 3 60720 306948 346 205 0 0 0 0 60720 315380 210 0 206 57640 0 60720 315380 206 0 1 267 207 61977 330059 121 4064 487 0 1 474 208 62620 288985 386 4449 188 0 1 534 209 60720 304485 351 14926 196 0 1 0 210 164 315380 201 0 1 15 8 6 211 152 14461 7687 0 1 397 28 7 212 199 4205 297 0 1 0 0 2 213 66 0 1 1866 29 28 56225 322331 214 335 0 1 288 4 3 60510 296656 215 1 0 1 0 0 0 60720 315380 216 1 3 2 1 60698 315354 296 57677 217 1 468 21 20 60720 312161 328 5341 218 1 20 2 1 60805 315576 169 57788 219 1 278 26 22 61404 314922 302 4420 220 1 61 14 9 60720 314551 307 8182 221 4 0 0 0 60720 315380 215 0 222 0 2 65276 312339 327 28085 586 0 223 63915 0 60720 315380 214 0 1 317 224 60720 298700 363 10967 311 0 1 738 225 60720 321376 139 12138 164 0 1 0 226 355 315380 216 0 1 368 17 13 227 204 6072 281 0 1 0 0 0 228 52257 0 1 2 1 0 60743 315487 229 1 0 1 0 0 0 60720 315380 230 1 53 6 6 60349 315793 163 19299 231 0 0 0 0 60720 315380 195 0 232 60720 0 60720 315380 194 0 1 0 233 321 315380 200 0 1 94 3 3 234 237 37629 1200 0 1 0 0 0 235 4772 0 1 24 8 7 59818 315637 236 53 0 1 2332 4 2 72680 324385 237 0 0 0 0 0 0 60720 315380 238 11 0 0 0 60720 315380 233 0 239 0 15 61808 308989 338 9908 834 0 240 60720 0 60720 315380 260 0 0 0 241 372 315380 234 0 1 206 9 9 242 240 10745 469 0 0 0 0 0 243 642 0 1 167 2 1 64245 307322 244 168 0 1 622 73 38 73007 304376 245 23 0 1 2328 94 57 82732 253588 246 1 0 1 0 0 0 60720 315380 247 1 365 8 7 54820 309560 335 13695 248 0 364 35 26 47705 298466 365 2813 249 0 0 0 0 60720 315380 239 0 250 60720 0 60720 315380 235 0 0 0 251 283 315380 262 0 0 0 0 0 252 636 0 1 226 12 13 72835 343929 253 430 0 1 307 15 10 58856 331955 254 1 0 0 0 0 0 60720 315380 255 0 0 0 0 60720 315380 281 0 256 77655 0 60720 315380 282 0 1 188 257 60720 381180 67 16471 964 0 0 0 258 117 315380 221 0 1 138 6 26 259 286 21903 954 0 0 0 0 0 260 0 0 0 0 0 0 60720 315380 261 12 0 0 0 60720 315380 265 0 262 0 19 60798 310201 332 9183 883 0 263 62452 0 60720 315380 284 0 1 282 264 64175 320016 141 4001 425 0 1 335 265 60720 320398 140 3648 359 0 1 0 266 381 315380 232 0 0 1324 117 29 267 331 785 69 0 1 176 28 8 268 256 3952 629 0 0 0 0 0 269 1 0 1 0 0 0 60720 315380 270 1 249 72 26 56726 313491 317 1576 271 13 0 0 0 60720 315380 242 0 272 0 9 70811 331323 118 10102 394 0 273 60720 0 60720 315380 241 0 0 601 274 62045 319210 144 19868 198 0 1 30 275 60720 318098 149 29525 3912 0 0 0 276 378 315380 254 0 1 249 62 13 277 258 1496 373 0 0 0 0 0 278 759 0 1 165 24 12 62841 325176 279 366 0 1 453 21 19 81125 365959 280 1 0 1 0 0 0 60720 315380 281 0 53 14 10 59506 302409 356 7315 282 1 382 21 9 59365 340968 107 6713 283 0 0 0 0 60720 315380 229 0 284 60720 0 60720 315380 230 0 1 0 285 223 315380 224 0 1 0 0 9 286 3808 0 0 30 4 4 60798 313164 287 348 0 1 290 2 1 58790 301164 288 0 0 1 0 0 1 60720 315380 289 53 0 0 0 60720 315380 231 0 290 9 14 61808 344425 101 2725 394 0 291 0 12 60735 315394 174 12822 63057 0 292 64016 0 60720 315380 222 0 0 209 293 54683 316647 157 8973 558 0 1 384 294 60720 309836 333 4993 286 0 0 0 295 244 315380 227 0 0 0 0 0 296 402 0 1 365 83 32 87192 346611 297 0 0 0 0 0 0 60720 315380 298 1 49 8 14 64107 322031 137 15254 299 1 3 4 8 60761 315656 167 28914 300 1 133 14 4 65990 339445 109 9960 301 1 32 1 0 59988 314964 301 114964 302 1 368 17 20 61167 297141 370 5714 303 1 1 6 5 60719 315372 293 19229 304 0 0 0 0 60720 315380 274 0 305 60720 0 60720 315380 246 0 1 0 306 185 315380 177 0 1 0 0 0 307 1 0 0 0 0 0 60720 315380 308 2 0 0 0 60720 315380 273 0 309 0 1 60722 312502 324 56251 5154 0 310 60720 0 60720 315380 252 0 0 0 311 253 315380 249 0 1 0 0 0 312 1 0 1 0 0 0 60720 315380 313 0 0 0 0 60720 315380 278 0 314 60720 0 60720 315380 275 0 1 0 315 316 315380 276 0 1 96 5 4 316 175 22746 1188 0 1 1 2 1 317 294 57694 117743 0 1 314 5 4 318 375 23074 367 0 1 844 78 20 319 180 1233 114 0 1 0 0 0 320 4328 0 1 26 1 1 60896 313880 321 944 0 1 125 13 10 59734 317698 322 314 0 1 304 15 12 62969 295580 323 1 0 1 0 0 0 60720 315380 324 0 0 0 0 60720 315380 183 0 325 60720 0 60720 315380 182 0 1 621 326 60720 308256 341 2255 174 0 1 0 327 354 315380 181 0 1 119 6 3 328 178 17280 872 0 1 0 0 0 329 1 0 1 0 0 0 60720 315380 330 1 1595 17 10 60720 319369 143 7022 331 1 312 14 3 58598 318690 147 8478 332 1 60 10 7 61124 314049 314 11405 333 1 587 12 10 59595 325699 126 10475 334 1 135 2 1 62065 314210 311 57105 335 0 0 0 0 60720 315380 272 0 336 78780 0 60720 315380 269 0 1 514 337 60720 322378 133 2353 238 0 1 0 338 251 315380 279 0 1 0 0 0 339 1 0 1 0 0 0 60720 315380 340 1 1 4 4 60722 315398 173 28849 341 0 0 0 0 60720 315380 268 0 342 61600 0 60720 315380 289 0 0 1763 343 59635 308336 340 4514 61 0 1 180 344 60720 316386 158 10581 647 0 1 0 345 248 315380 287 0 1 0 0 0 346 1 0 1 0 0 0 60720 315380 347 21 0 0 0 60720 315380 292 0 348 0 7 60720 315553 170 5503 531 0 349 59781 0 60720 315380 290 0 1 448 350 76644 323361 132 3084 275 0 1 227 351 64820 336639 110 15182 603 0 1 174 352 60720 307424 344 107424 618 0 1 0 353 280 315380 266 0 1 0 0 0 354 788 0 1 121 24 11 56178 295370 355 201 0 1 607 11 7 60436 322340 356 54 0 1 2212 14 10 60720 319864 357 1 0 1 0 0 0 60720 315380 358 60 0 0 0 60720 315380 264 0 359 80 18 73433 317291 155 1955 221 0 360 0 14 41477 280398 390 1005 141 0 361 62700 0 60720 315380 267 0 1 78 362 67804 317330 154 7333 1511 0 1 2489 363 59661 238125 408 953 15 0 1 131 364 58620 327071 123 21179 969 0 1 923 365 60398 309038 337 13630 118 0 1 72 366 58580 314210 312 38070 1591 0 1 572 367 62710 307930 342 6746 189 -1 1 397 368 59325 322327 136 12233 308 -1 1 450 369 60950 292136 380 11517 205 -1 1 622 370 68060 263276 398 9039 102 -2 1 694 371 83620 367655 83 20957 241 -2 0 3425 372 58456 283910 388 6993 25 -3 1 562 373 52811 283587 389 6430 149 -3 1 4917 374 121173 243650 406 1039 9 -3 0 1442 375 63870 438493 49 2021 165 -3 1 529 376 21001 296261 374 10696 182 -3 0 2126 377 70415 230621 413 222 14 -3 1 1061 378 64230 304252 353 20850 98 -4 1 776 379 59190 333505 113 14834 172 -4 1 611 380 69351 296919 371 12115 159 -4 0 1526 381 64270 278990 391 3160 52 -4 1 592 382 70694 276898 392 10985 130 -4 1 1182 383 68005 327007 124 9770 107 -5 1 621 384 58930 317046 156 7315 188 -5 1 989 385 58320 304555 350 9505 106 -5 1 438 386 69980 298096 366 8918 224 -5 1 726 387 69863 231861 412 10620 44 -5 1 1303 388 63255 309422 336 1794 84 -5 0 7419 389 57320 286963 387 2999 12 -5 0 1164 390 75230 269753 395 4103 60 -6 0 3310 391 79420 448243 44 7523 75 -6 1 1920 392 73490 165404 424 -2306 -18 -6 1 965 393 35250 204325 421 1442 4 -8 1 3256 394 62285 407159 57 3139 64 -8 0 1135 395 69206 290476 384 5322 80 -9 1 1270 396 65920 275311 393 2897 59 -10 1 661 397 69770 246541 405 15514 70 -12 1 1013 398 72683 253468 402 26734 53 -14 1 2844 399 -14545 240897 407 610 14 -15 0 11528 400 55830 -83265 431 -4047 -25 -25 1 6526 401 55174 -42143 430 -9313 -37 -28 1 2264 402 67038 272713 394 3030 32 -31 0 5109 403 51252 215362 417 158 3 -31 1 3999 404 157278 42754 428 -5242 -39 -58 1 35624 405 79510 306275 347 477 3 -70 1 9252 406 77440 253537 401 1115 6 -129 1 15236 407 27284 372631 78 1918 11 -158 1 18073 408 213118 -7170 429 -1151 -11 -366 0 162556 409 81767 6282154 1 5626 37 56289 0 29790 410 153198 4321023 2 13337 138 28328 0 87550 411 -26007 4111912 3 8541 45 17936 1 84738 412 126942 223193 415 39 0 4145 0 54660 413 157214 1491348 6 4276 24 3040 0 42634 414 129352 1629616 5 9164 34 2964 1 40949 415 234817 1398893 9 2493 29 2865 0 45187 416 60448 1926517 4 4891 38 2854 0 37704 417 47818 983660 14 1734 21 2167 0 16275 418 245546 1443586 7 11409 76 1974 1 25830 419 48020 1073089 10 7592 34 1910 1 12679 420 -1710 984885 13 7135 62 1871 0 18014 421 32648 1405225 8 5043 67 943 1 43556 422 95350 227132 414 110 1 929 0 24811 423 151352 929118 15 1444 29 822 1 6575 424 288170 1071292 11 5480 133 819 1 7123 425 114337 638830 28 4026 62 769 0 21950 426 37884 856956 17 1266 30 745 0 37597 427 122844 992426 12 3195 21 652 1 17821 428 82340 444477 46 655 14 643 0 12988 429 79801 857217 16 5523 51 601 0 22330 430 165548 711969 20 6095 23 446 1 13326 431 116384 702380 21 4925 38 436 1 16189 Profit/Trades Profit/Cost 1 5626 37 2 13337 138 3 8541 45 4 39 0 5 4276 24 6 9164 34 7 2493 29 8 4891 38 9 1734 21 10 11409 76 11 7592 34 12 7135 62 13 5043 67 14 110 1 15 1444 29 16 5480 133 17 4026 62 18 1266 30 19 3195 21 20 655 14 21 5523 51 22 6095 23 23 4925 38 24 538 10 25 933 14 26 6027 24 27 1624 17 28 52 1 29 15856 68 30 664 17 31 17456 145 32 4404 32 33 12031 78 34 8118 88 35 20190 200 36 1148 10 37 4838 43 38 6141 93 39 4438 148 40 51 2 41 119 12 42 2390 11 43 1335 67 44 1714 26 45 3024 41 46 7092 48 47 787 35 48 717 31 49 3204 55 50 -167 -2 51 3062 105 52 2122 77 53 2533 40 54 2019 288 55 214 7 56 2123 86 57 50325 121 58 3730 46 59 1321 110 60 1858 20 61 6243 53 62 -1408 -35 63 6308 69 64 1889 22 65 1964 58 66 167 12 67 827 70 68 3720 96 69 1861 123 70 3570 56 71 35 1 72 3643 35 73 3258 77 74 11431 227 75 5557 98 76 2353 27 77 3576 166 78 855 45 79 960 88 80 818 210 81 5669 107 82 3017 48 83 2591 174 84 5660 125 85 11834 346 86 2185 28 87 3004 71 88 11016 286 89 1697 52 90 16928 329 91 3718 97 92 4354 61 93 2550 112 94 2198 25 95 3526 84 96 4108 58 97 28639 844 98 659 92 99 9253 79 100 10717 85 101 3258 85 102 1528 36 103 1241 29 104 6337 198 105 4205 100 106 3029 93 107 -1243 -23 108 15542 152 109 1998 21 110 1386 26 111 6963 223 112 -883 -12 113 1298 18 114 6394 37 115 19783 128 116 3679 95 117 5106 330 118 11566 108 119 888 6 120 1677 138 121 1904 17 122 9591 153 123 1457 98 124 9668 169 125 4137 58 126 748 5 127 0 1 128 0 1 129 0 1 130 34 8 131 1558 63 132 0 0 133 20 60887 134 4 60720 135 4 60505 136 10 60945 137 1 60720 138 4 60720 139 0 60720 140 325506 127 141 315380 220 142 6571 98 143 38611 3215 144 8652 282 145 10262 365 146 0 1 147 0 1 148 0 1 149 0 1 150 0 1 151 0 1 152 458 11 153 20 3 154 0 0 155 0 60720 156 315236 298 157 336425 111 158 315380 213 159 0 1 160 0 0 161 7 60830 162 24 63261 163 1 60720 164 0 60720 165 382712 66 166 341570 105 167 315380 192 168 0 1 169 0 1 170 203 9 171 0 0 172 8 60951 173 113 60720 174 0 60745 175 0 60720 176 315380 190 177 4686 164 178 13007 245 179 2005 46 180 1873 258 181 22961 4974 182 9515 209 183 5245 66 184 22666 761 185 18357 406 186 0 1 187 0 1 188 141 6 189 0 0 190 3 60890 191 0 60720 192 409642 55 193 315380 259 194 0 0 195 0 0 196 0 0 197 0 60720 198 300962 361 199 325479 128 200 316155 160 201 318574 148 202 315380 228 203 4488 126 204 3450 404 205 1 0 206 32 11 207 20 10 208 7 5 209 0 2 210 60831 315688 211 60646 317736 212 60720 315380 213 135 4218 214 373 24164 215 198 0 216 36505 0 217 239 0 218 5924 0 219 413 0 220 1874 0 221 1 192 222 1 0 223 9 7 224 10 9 225 0 0 226 61686 303230 227 60720 315380 228 172 115487 229 202 0 230 2197 0 231 1 0 232 0 0 233 61360 312887 234 60720 315380 235 168 14455 236 131 31096 237 236 0 238 1 131 239 0 0 240 0 0 241 53110 296702 242 60720 315380 243 345 53661 244 352 1430 245 400 570 246 238 0 247 300 0 248 271 0 249 1 0 250 0 0 251 60720 315380 252 102 11994 253 116 8797 254 288 0 255 0 0 256 11 9 257 0 0 258 69817 331420 259 60720 315380 260 285 0 261 1 125 262 1 0 263 30 12 264 33 23 265 0 0 266 67440 291841 267 68136 310670 268 60720 315380 269 255 0 270 455 0 271 1 333 272 1 0 273 6 5 274 4 3 275 0 0 276 54323 292754 277 60720 315380 278 130 5216 279 84 7903 280 257 0 281 1932 0 282 369 0 283 0 0 284 0 0 285 60720 315380 286 320 28291 287 360 50582 288 225 0 289 1 366 290 1 2 291 0 0 292 13 19 293 22 17 294 0 0 295 60720 315380 296 99 1766 297 226 0 298 2500 0 299 37920 0 300 1052 0 301 3603 0 302 264 0 303 84213 0 304 1 0 305 0 0 306 60720 315380 307 245 0 308 1 22 309 0 0 310 0 0 311 60720 315380 312 277 0 313 1 0 314 0 0 315 60379 313729 316 60727 315388 317 60720 315371 318 60925 296139 319 60720 315380 320 315 113880 321 153 9054 322 376 6372 323 187 0 324 1 0 325 48 13 326 0 0 327 59118 303677 328 60720 315380 329 176 0 330 75 0 331 381 0 332 1916 0 333 214 0 334 847 0 335 1 0 336 52 15 337 0 0 338 60720 315380 339 270 0 340 172235 0 341 1 0 342 24 28 343 11 9 344 0 0 345 60720 315380 346 291 0 347 1 218 348 1 0 349 40 7 350 9 7 351 1 3 352 0 0 353 60720 315380 354 377 3974 355 134 11122 356 142 8562 357 271 0 358 1 530 359 1 571 360 1 0 361 16 12 362 40 29 363 6 3 364 8 6 365 3 3 366 16 8 367 10 10 368 8 6 369 7 8 370 8 6 371 12 9 372 13 8 373 42 26 374 118 239 375 9 7 376 138 41 377 5 3 378 9 8 379 8 6 380 25 21 381 7 7 382 13 11 383 16 11 384 11 12 385 11 9 386 3 3 387 61 57 388 29 21 389 17 15 390 33 32 391 15 11 392 3 2 393 66 23 394 17 20 395 26 24 396 3 1 397 2 1 398 67 74 399 70 68 400 26 20 401 24 20 402 97 82 403 30 21 404 223 244 405 48 32 406 90 86 407 180 69 408 1081 807 409 309 444 410 458 412 411 588 428 412 302 315 413 156 168 414 481 263 415 353 267 416 452 228 417 109 129 418 115 104 419 110 122 420 239 393 421 247 190 422 505 280 423 159 63 424 109 102 425 519 265 426 248 234 427 373 277 428 119 73 429 84 67 430 102 103 431 295 290 > k <- length(x[1,]) > df <- as.data.frame(x) > (mylm <- lm(df)) Call: lm(formula = df) Coefficients: (Intercept) groep Costs Trades 5.858e+04 1.693e-02 3.853e-01 -1.255e-01 Orders Dividends Wealth WRank -1.389e-01 -2.128e-01 -2.612e-02 -1.516e-01 `Profit/Trades` `Profit/Cost` -1.405e-01 -9.062e-02 > (mysum <- summary(mylm)) Call: lm(formula = df) Residuals: Min 1Q Median 3Q Max -201474 -29102 -9139 15042 262849 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 5.858e+04 5.623e+03 10.417 < 2e-16 *** groep 1.693e-02 6.514e-03 2.598 0.00969 ** Costs 3.853e-01 5.711e-02 6.748 4.99e-11 *** Trades -1.255e-01 4.158e-02 -3.019 0.00269 ** Orders -1.389e-01 5.116e-02 -2.714 0.00692 ** Dividends -2.128e-01 3.508e-02 -6.064 2.95e-09 *** Wealth -2.612e-02 6.503e-03 -4.016 7.02e-05 *** WRank -1.516e-01 3.564e-02 -4.253 2.60e-05 *** `Profit/Trades` -1.405e-01 4.855e-02 -2.893 0.00401 ** `Profit/Cost` -9.062e-02 3.927e-02 -2.308 0.02148 * --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Residual standard error: 61930 on 421 degrees of freedom Multiple R-squared: 0.2737, Adjusted R-squared: 0.2582 F-statistic: 17.63 on 9 and 421 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,] 7.342694e-04 1.468539e-03 9.992657e-01 [2,] 1.671036e-04 3.342072e-04 9.998329e-01 [3,] 1.413869e-05 2.827738e-05 9.999859e-01 [4,] 1.132602e-06 2.265204e-06 9.999989e-01 [5,] 1.189625e-07 2.379249e-07 9.999999e-01 [6,] 8.473423e-09 1.694685e-08 1.000000e+00 [7,] 8.827562e-10 1.765512e-09 1.000000e+00 [8,] 6.080367e-11 1.216073e-10 1.000000e+00 [9,] 9.120068e-12 1.824014e-11 1.000000e+00 [10,] 1.414215e-12 2.828431e-12 1.000000e+00 [11,] 1.108327e-13 2.216654e-13 1.000000e+00 [12,] 8.083261e-15 1.616652e-14 1.000000e+00 [13,] 5.243550e-16 1.048710e-15 1.000000e+00 [14,] 3.724208e-17 7.448416e-17 1.000000e+00 [15,] 2.897221e-18 5.794441e-18 1.000000e+00 [16,] 1.818472e-19 3.636943e-19 1.000000e+00 [17,] 1.139150e-20 2.278300e-20 1.000000e+00 [18,] 6.552352e-22 1.310470e-21 1.000000e+00 [19,] 6.980413e-23 1.396083e-22 1.000000e+00 [20,] 9.759528e-24 1.951906e-23 1.000000e+00 [21,] 7.293052e-25 1.458610e-24 1.000000e+00 [22,] 6.696854e-26 1.339371e-25 1.000000e+00 [23,] 4.277184e-27 8.554368e-27 1.000000e+00 [24,] 3.342157e-28 6.684314e-28 1.000000e+00 [25,] 1.970605e-29 3.941210e-29 1.000000e+00 [26,] 1.391374e-30 2.782747e-30 1.000000e+00 [27,] 9.087902e-32 1.817580e-31 1.000000e+00 [28,] 4.953342e-33 9.906683e-33 1.000000e+00 [29,] 3.501977e-34 7.003955e-34 1.000000e+00 [30,] 2.386415e-35 4.772830e-35 1.000000e+00 [31,] 1.246882e-36 2.493765e-36 1.000000e+00 [32,] 7.881900e-38 1.576380e-37 1.000000e+00 [33,] 4.534801e-39 9.069603e-39 1.000000e+00 [34,] 6.174878e-40 1.234976e-39 1.000000e+00 [35,] 5.263815e-41 1.052763e-40 1.000000e+00 [36,] 3.686819e-42 7.373638e-42 1.000000e+00 [37,] 1.950022e-43 3.900044e-43 1.000000e+00 [38,] 9.911705e-45 1.982341e-44 1.000000e+00 [39,] 5.129447e-46 1.025889e-45 1.000000e+00 [40,] 2.822262e-47 5.644525e-47 1.000000e+00 [41,] 1.872533e-48 3.745067e-48 1.000000e+00 [42,] 1.442526e-49 2.885053e-49 1.000000e+00 [43,] 6.992315e-51 1.398463e-50 1.000000e+00 [44,] 3.431522e-52 6.863044e-52 1.000000e+00 [45,] 3.195138e-53 6.390275e-53 1.000000e+00 [46,] 1.685090e-54 3.370179e-54 1.000000e+00 [47,] 7.988731e-56 1.597746e-55 1.000000e+00 [48,] 6.051215e-57 1.210243e-56 1.000000e+00 [49,] 3.837768e-58 7.675537e-58 1.000000e+00 [50,] 1.863146e-59 3.726293e-59 1.000000e+00 [51,] 1.222453e-60 2.444907e-60 1.000000e+00 [52,] 8.009968e-62 1.601994e-61 1.000000e+00 [53,] 3.910441e-63 7.820881e-63 1.000000e+00 [54,] 1.737013e-64 3.474026e-64 1.000000e+00 [55,] 8.723242e-66 1.744648e-65 1.000000e+00 [56,] 4.878229e-67 9.756458e-67 1.000000e+00 [57,] 2.182721e-68 4.365441e-68 1.000000e+00 [58,] 1.282330e-69 2.564659e-69 1.000000e+00 [59,] 5.849079e-71 1.169816e-70 1.000000e+00 [60,] 2.823939e-72 5.647878e-72 1.000000e+00 [61,] 1.250837e-73 2.501674e-73 1.000000e+00 [62,] 5.219231e-75 1.043846e-74 1.000000e+00 [63,] 2.348134e-76 4.696269e-76 1.000000e+00 [64,] 9.809470e-78 1.961894e-77 1.000000e+00 [65,] 4.015775e-79 8.031550e-79 1.000000e+00 [66,] 1.772331e-80 3.544662e-80 1.000000e+00 [67,] 7.867738e-82 1.573548e-81 1.000000e+00 [68,] 3.137371e-83 6.274742e-83 1.000000e+00 [69,] 1.315035e-84 2.630071e-84 1.000000e+00 [70,] 6.915642e-86 1.383128e-85 1.000000e+00 [71,] 2.714751e-87 5.429503e-87 1.000000e+00 [72,] 1.216563e-88 2.433125e-88 1.000000e+00 [73,] 4.940744e-90 9.881487e-90 1.000000e+00 [74,] 2.644442e-91 5.288884e-91 1.000000e+00 [75,] 1.083441e-92 2.166881e-92 1.000000e+00 [76,] 4.110435e-94 8.220871e-94 1.000000e+00 [77,] 1.592398e-95 3.184796e-95 1.000000e+00 [78,] 5.995667e-97 1.199133e-96 1.000000e+00 [79,] 2.366887e-98 4.733774e-98 1.000000e+00 [80,] 1.204595e-99 2.409189e-99 1.000000e+00 [81,] 4.502275e-101 9.004551e-101 1.000000e+00 [82,] 2.196001e-102 4.392001e-102 1.000000e+00 [83,] 8.451217e-104 1.690243e-103 1.000000e+00 [84,] 3.895505e-105 7.791010e-105 1.000000e+00 [85,] 2.237688e-106 4.475376e-106 1.000000e+00 [86,] 8.529629e-108 1.705926e-107 1.000000e+00 [87,] 3.592146e-109 7.184292e-109 1.000000e+00 [88,] 1.688979e-110 3.377958e-110 1.000000e+00 [89,] 5.952125e-112 1.190425e-111 1.000000e+00 [90,] 2.198982e-113 4.397964e-113 1.000000e+00 [91,] 7.978659e-115 1.595732e-114 1.000000e+00 [92,] 2.968050e-116 5.936101e-116 1.000000e+00 [93,] 1.069912e-117 2.139824e-117 1.000000e+00 [94,] 3.955853e-119 7.911707e-119 1.000000e+00 [95,] 1.622237e-120 3.244474e-120 1.000000e+00 [96,] 6.998540e-122 1.399708e-121 1.000000e+00 [97,] 2.526290e-123 5.052579e-123 1.000000e+00 [98,] 8.573506e-125 1.714701e-124 1.000000e+00 [99,] 2.892748e-126 5.785495e-126 1.000000e+00 [100,] 1.078990e-127 2.157979e-127 1.000000e+00 [101,] 3.807350e-129 7.614700e-129 1.000000e+00 [102,] 1.243757e-130 2.487514e-130 1.000000e+00 [103,] 6.399957e-132 1.279991e-131 1.000000e+00 [104,] 2.182655e-133 4.365310e-133 1.000000e+00 [105,] 7.075193e-135 1.415039e-134 1.000000e+00 [106,] 2.551251e-136 5.102502e-136 1.000000e+00 [107,] 8.399734e-138 1.679947e-137 1.000000e+00 [108,] 2.653058e-139 5.306117e-139 1.000000e+00 [109,] 9.318868e-141 1.863774e-140 1.000000e+00 [110,] 3.096816e-142 6.193633e-142 1.000000e+00 [111,] 1.068004e-143 2.136008e-143 1.000000e+00 [112,] 4.072451e-145 8.144903e-145 1.000000e+00 [113,] 1.512083e-146 3.024165e-146 1.000000e+00 [114,] 4.774809e-148 9.549618e-148 1.000000e+00 [115,] 1.485975e-149 2.971949e-149 1.000000e+00 [116,] 1.034168e-150 2.068337e-150 1.000000e+00 [117,] 3.210147e-152 6.420295e-152 1.000000e+00 [118,] 9.399032e-152 1.879806e-151 1.000000e+00 [119,] 1.553888e-152 3.107775e-152 1.000000e+00 [120,] 7.990306e-153 1.598061e-152 1.000000e+00 [121,] 4.739132e-154 9.478264e-154 1.000000e+00 [122,] 1.818213e-154 3.636426e-154 1.000000e+00 [123,] 2.260372e-154 4.520744e-154 1.000000e+00 [124,] 6.397645e-154 1.279529e-153 1.000000e+00 [125,] 6.354008e-154 1.270802e-153 1.000000e+00 [126,] 3.761040e-153 7.522079e-153 1.000000e+00 [127,] 5.873642e-153 1.174728e-152 1.000000e+00 [128,] 1.237845e-151 2.475691e-151 1.000000e+00 [129,] 6.008162e-150 1.201632e-149 1.000000e+00 [130,] 2.582766e-151 5.165531e-151 1.000000e+00 [131,] 9.064090e-148 1.812818e-147 1.000000e+00 [132,] 4.001539e-149 8.003079e-149 1.000000e+00 [133,] 1.736612e-150 3.473223e-150 1.000000e+00 [134,] 7.645654e-152 1.529131e-151 1.000000e+00 [135,] 5.925764e-151 1.185153e-150 1.000000e+00 [136,] 2.978189e-152 5.956377e-152 1.000000e+00 [137,] 2.938185e-153 5.876369e-153 1.000000e+00 [138,] 4.101842e-154 8.203684e-154 1.000000e+00 [139,] 2.225561e-155 4.451122e-155 1.000000e+00 [140,] 1.116614e-156 2.233228e-156 1.000000e+00 [141,] 2.594041e-156 5.188083e-156 1.000000e+00 [142,] 1.322628e-155 2.645255e-155 1.000000e+00 [143,] 1.164702e-156 2.329405e-156 1.000000e+00 [144,] 1.134212e-123 2.268424e-123 1.000000e+00 [145,] 9.106685e-125 1.821337e-124 1.000000e+00 [146,] 1.762328e-121 3.524656e-121 1.000000e+00 [147,] 1.498869e-122 2.997739e-122 1.000000e+00 [148,] 1.350807e-123 2.701614e-123 1.000000e+00 [149,] 1.240879e-123 2.481758e-123 1.000000e+00 [150,] 1.265193e-117 2.530386e-117 1.000000e+00 [151,] 1.158204e-115 2.316408e-115 1.000000e+00 [152,] 7.427545e-114 1.485509e-113 1.000000e+00 [153,] 1.575370e-114 3.150741e-114 1.000000e+00 [154,] 1.260643e-113 2.521286e-113 1.000000e+00 [155,] 7.340168e-113 1.468034e-112 1.000000e+00 [156,] 7.924513e-114 1.584903e-113 1.000000e+00 [157,] 8.472328e-115 1.694466e-114 1.000000e+00 [158,] 8.954579e-116 1.790916e-115 1.000000e+00 [159,] 2.538928e-113 5.077856e-113 1.000000e+00 [160,] 5.954263e-114 1.190853e-113 1.000000e+00 [161,] 5.516338e-111 1.103268e-110 1.000000e+00 [162,] 7.853296e-107 1.570659e-106 1.000000e+00 [163,] 1.890279e-104 3.780557e-104 1.000000e+00 [164,] 4.885430e-104 9.770860e-104 1.000000e+00 [165,] 6.484921e-105 1.296984e-104 1.000000e+00 [166,] 8.444318e-106 1.688864e-105 1.000000e+00 [167,] 1.091504e-106 2.183008e-106 1.000000e+00 [168,] 1.389103e-107 2.778206e-107 1.000000e+00 [169,] 4.487009e-104 8.974018e-104 1.000000e+00 [170,] 6.095636e-105 1.219127e-104 1.000000e+00 [171,] 8.487882e-106 1.697576e-105 1.000000e+00 [172,] 1.208207e-106 2.416414e-106 1.000000e+00 [173,] 1.604372e-107 3.208743e-107 1.000000e+00 [174,] 2.232679e-108 4.465357e-108 1.000000e+00 [175,] 3.014372e-109 6.028744e-109 1.000000e+00 [176,] 4.004255e-110 8.008510e-110 1.000000e+00 [177,] 5.612622e-103 1.122524e-102 1.000000e+00 [178,] 8.572358e-104 1.714472e-103 1.000000e+00 [179,] 1.200102e-98 2.400205e-98 1.000000e+00 [180,] 2.595177e-97 5.190354e-97 1.000000e+00 [181,] 1.597128e-97 3.194257e-97 1.000000e+00 [182,] 2.837916e-98 5.675832e-98 1.000000e+00 [183,] 4.565324e-99 9.130648e-99 1.000000e+00 [184,] 1.534657e-99 3.069314e-99 1.000000e+00 [185,] 3.472320e-99 6.944639e-99 1.000000e+00 [186,] 8.946582e-100 1.789316e-99 1.000000e+00 [187,] 1.619234e-100 3.238468e-100 1.000000e+00 [188,] 1.479917e-86 2.959834e-86 1.000000e+00 [189,] 2.859177e-85 5.718354e-85 1.000000e+00 [190,] 3.100290e-85 6.200580e-85 1.000000e+00 [191,] 7.280909e-86 1.456182e-85 1.000000e+00 [192,] 2.453082e-84 4.906165e-84 1.000000e+00 [193,] 6.455723e-85 1.291145e-84 1.000000e+00 [194,] 1.899283e-85 3.798566e-85 1.000000e+00 [195,] 1.722217e-85 3.444434e-85 1.000000e+00 [196,] 6.433824e-86 1.286765e-85 1.000000e+00 [197,] 1.361017e-86 2.722035e-86 1.000000e+00 [198,] 1.000000e+00 2.872050e-18 1.436025e-18 [199,] 1.000000e+00 2.483596e-24 1.241798e-24 [200,] 1.000000e+00 2.886336e-26 1.443168e-26 [201,] 1.000000e+00 8.422094e-27 4.211047e-27 [202,] 1.000000e+00 1.717982e-26 8.589909e-27 [203,] 1.000000e+00 3.606516e-26 1.803258e-26 [204,] 1.000000e+00 3.555737e-26 1.777868e-26 [205,] 1.000000e+00 5.051883e-26 2.525942e-26 [206,] 1.000000e+00 7.232400e-26 3.616200e-26 [207,] 1.000000e+00 1.208167e-25 6.040836e-26 [208,] 1.000000e+00 2.110087e-25 1.055043e-25 [209,] 1.000000e+00 3.838929e-25 1.919465e-25 [210,] 1.000000e+00 5.475544e-25 2.737772e-25 [211,] 1.000000e+00 1.072212e-24 5.361061e-25 [212,] 1.000000e+00 1.702378e-24 8.511888e-25 [213,] 1.000000e+00 2.891772e-24 1.445886e-24 [214,] 1.000000e+00 1.462230e-24 7.311152e-25 [215,] 1.000000e+00 5.252751e-25 2.626375e-25 [216,] 1.000000e+00 3.256883e-25 1.628442e-25 [217,] 1.000000e+00 6.802858e-25 3.401429e-25 [218,] 1.000000e+00 1.244132e-24 6.220662e-25 [219,] 1.000000e+00 2.359053e-24 1.179526e-24 [220,] 1.000000e+00 4.765485e-24 2.382743e-24 [221,] 1.000000e+00 5.974073e-24 2.987036e-24 [222,] 1.000000e+00 5.203212e-24 2.601606e-24 [223,] 1.000000e+00 1.099459e-23 5.497294e-24 [224,] 1.000000e+00 2.186432e-23 1.093216e-23 [225,] 1.000000e+00 4.488388e-23 2.244194e-23 [226,] 1.000000e+00 8.506042e-23 4.253021e-23 [227,] 1.000000e+00 1.025648e-22 5.128241e-23 [228,] 1.000000e+00 2.034321e-22 1.017160e-22 [229,] 1.000000e+00 3.405531e-22 1.702765e-22 [230,] 1.000000e+00 3.860711e-22 1.930355e-22 [231,] 1.000000e+00 7.750780e-22 3.875390e-22 [232,] 1.000000e+00 1.561917e-21 7.809587e-22 [233,] 1.000000e+00 2.434959e-21 1.217480e-21 [234,] 1.000000e+00 4.906950e-21 2.453475e-21 [235,] 1.000000e+00 9.184709e-21 4.592354e-21 [236,] 1.000000e+00 1.692756e-20 8.463780e-21 [237,] 1.000000e+00 3.212635e-20 1.606317e-20 [238,] 1.000000e+00 6.195787e-20 3.097894e-20 [239,] 1.000000e+00 1.102990e-19 5.514951e-20 [240,] 1.000000e+00 2.152313e-19 1.076157e-19 [241,] 1.000000e+00 4.159148e-19 2.079574e-19 [242,] 1.000000e+00 8.094882e-19 4.047441e-19 [243,] 1.000000e+00 1.510116e-18 7.550578e-19 [244,] 1.000000e+00 2.306479e-18 1.153239e-18 [245,] 1.000000e+00 3.554633e-18 1.777316e-18 [246,] 1.000000e+00 6.197187e-18 3.098593e-18 [247,] 1.000000e+00 8.868997e-18 4.434499e-18 [248,] 1.000000e+00 1.689403e-17 8.447014e-18 [249,] 1.000000e+00 3.112156e-17 1.556078e-17 [250,] 1.000000e+00 3.771515e-17 1.885758e-17 [251,] 1.000000e+00 6.869922e-17 3.434961e-17 [252,] 1.000000e+00 1.159316e-16 5.796579e-17 [253,] 1.000000e+00 1.974091e-16 9.870456e-17 [254,] 1.000000e+00 3.506268e-16 1.753134e-16 [255,] 1.000000e+00 5.426555e-16 2.713277e-16 [256,] 1.000000e+00 8.610821e-16 4.305411e-16 [257,] 1.000000e+00 1.586604e-15 7.933022e-16 [258,] 1.000000e+00 2.832254e-15 1.416127e-15 [259,] 1.000000e+00 5.089670e-15 2.544835e-15 [260,] 1.000000e+00 6.169040e-15 3.084520e-15 [261,] 1.000000e+00 1.105601e-14 5.528003e-15 [262,] 1.000000e+00 1.890923e-14 9.454617e-15 [263,] 1.000000e+00 3.115656e-14 1.557828e-14 [264,] 1.000000e+00 5.409336e-14 2.704668e-14 [265,] 1.000000e+00 8.760751e-14 4.380376e-14 [266,] 1.000000e+00 1.580633e-13 7.903165e-14 [267,] 1.000000e+00 2.650410e-13 1.325205e-13 [268,] 1.000000e+00 4.681801e-13 2.340900e-13 [269,] 1.000000e+00 8.114773e-13 4.057386e-13 [270,] 1.000000e+00 1.292039e-12 6.460197e-13 [271,] 1.000000e+00 2.235785e-12 1.117892e-12 [272,] 1.000000e+00 3.833883e-12 1.916941e-12 [273,] 1.000000e+00 6.515085e-12 3.257542e-12 [274,] 1.000000e+00 1.140828e-11 5.704140e-12 [275,] 1.000000e+00 1.982738e-11 9.913692e-12 [276,] 1.000000e+00 3.366661e-11 1.683330e-11 [277,] 1.000000e+00 5.668184e-11 2.834092e-11 [278,] 1.000000e+00 6.434433e-11 3.217216e-11 [279,] 1.000000e+00 1.393868e-11 6.969341e-12 [280,] 1.000000e+00 2.288294e-11 1.144147e-11 [281,] 1.000000e+00 3.671452e-11 1.835726e-11 [282,] 1.000000e+00 6.051392e-11 3.025696e-11 [283,] 1.000000e+00 1.014109e-10 5.070545e-11 [284,] 1.000000e+00 6.004096e-11 3.002048e-11 [285,] 1.000000e+00 1.038719e-10 5.193595e-11 [286,] 1.000000e+00 1.739296e-10 8.696480e-11 [287,] 1.000000e+00 2.864878e-10 1.432439e-10 [288,] 1.000000e+00 4.724611e-10 2.362306e-10 [289,] 1.000000e+00 3.463082e-10 1.731541e-10 [290,] 1.000000e+00 5.784774e-10 2.892387e-10 [291,] 1.000000e+00 9.532554e-10 4.766277e-10 [292,] 1.000000e+00 1.587571e-09 7.937857e-10 [293,] 1.000000e+00 2.554539e-09 1.277270e-09 [294,] 1.000000e+00 4.166889e-09 2.083445e-09 [295,] 1.000000e+00 6.911764e-09 3.455882e-09 [296,] 1.000000e+00 1.126637e-08 5.633183e-09 [297,] 1.000000e+00 1.768985e-08 8.844925e-09 [298,] 1.000000e+00 2.748604e-08 1.374302e-08 [299,] 1.000000e+00 4.382361e-08 2.191180e-08 [300,] 1.000000e+00 7.071430e-08 3.535715e-08 [301,] 9.999999e-01 1.122734e-07 5.613672e-08 [302,] 9.999999e-01 1.695239e-07 8.476197e-08 [303,] 9.999999e-01 2.651428e-07 1.325714e-07 [304,] 9.999998e-01 4.002681e-07 2.001341e-07 [305,] 9.999997e-01 5.325398e-07 2.662699e-07 [306,] 9.999996e-01 8.081179e-07 4.040589e-07 [307,] 9.999994e-01 1.222764e-06 6.113821e-07 [308,] 9.999991e-01 1.718039e-06 8.590193e-07 [309,] 9.999987e-01 2.641022e-06 1.320511e-06 [310,] 9.999982e-01 3.667798e-06 1.833899e-06 [311,] 9.999972e-01 5.596657e-06 2.798328e-06 [312,] 9.999958e-01 8.398375e-06 4.199187e-06 [313,] 9.999940e-01 1.195680e-05 5.978401e-06 [314,] 9.999913e-01 1.749095e-05 8.745474e-06 [315,] 9.999871e-01 2.586047e-05 1.293023e-05 [316,] 9.999812e-01 3.752039e-05 1.876019e-05 [317,] 9.999724e-01 5.520940e-05 2.760470e-05 [318,] 9.999600e-01 8.009134e-05 4.004567e-05 [319,] 9.999426e-01 1.147351e-04 5.736755e-05 [320,] 9.999177e-01 1.646470e-04 8.232349e-05 [321,] 9.998849e-01 2.301897e-04 1.150948e-04 [322,] 9.998698e-01 2.604956e-04 1.302478e-04 [323,] 9.998152e-01 3.696313e-04 1.848156e-04 [324,] 9.997856e-01 4.288529e-04 2.144265e-04 [325,] 9.997032e-01 5.936397e-04 2.968198e-04 [326,] 9.995859e-01 8.282803e-04 4.141402e-04 [327,] 9.994234e-01 1.153129e-03 5.765643e-04 [328,] 9.992253e-01 1.549375e-03 7.746877e-04 [329,] 9.989378e-01 2.124414e-03 1.062207e-03 [330,] 9.986652e-01 2.669662e-03 1.334831e-03 [331,] 9.982147e-01 3.570654e-03 1.785327e-03 [332,] 9.976335e-01 4.733004e-03 2.366502e-03 [333,] 9.968420e-01 6.315955e-03 3.157977e-03 [334,] 9.957991e-01 8.401827e-03 4.200913e-03 [335,] 9.944619e-01 1.107627e-02 5.538136e-03 [336,] 9.930724e-01 1.385518e-02 6.927588e-03 [337,] 9.916411e-01 1.671781e-02 8.358907e-03 [338,] 9.893242e-01 2.135170e-02 1.067585e-02 [339,] 9.863721e-01 2.725590e-02 1.362795e-02 [340,] 9.827538e-01 3.449247e-02 1.724623e-02 [341,] 9.787515e-01 4.249696e-02 2.124848e-02 [342,] 9.733271e-01 5.334576e-02 2.667288e-02 [343,] 9.665944e-01 6.681128e-02 3.340564e-02 [344,] 9.586726e-01 8.265474e-02 4.132737e-02 [345,] 9.585874e-01 8.282526e-02 4.141263e-02 [346,] 9.518916e-01 9.621682e-02 4.810841e-02 [347,] 9.418464e-01 1.163071e-01 5.815357e-02 [348,] 9.337118e-01 1.325765e-01 6.628823e-02 [349,] 9.491276e-01 1.017447e-01 5.087235e-02 [350,] 9.376856e-01 1.246288e-01 6.231438e-02 [351,] 9.239844e-01 1.520312e-01 7.601560e-02 [352,] 9.096983e-01 1.806035e-01 9.030173e-02 [353,] 8.917917e-01 2.164166e-01 1.082083e-01 [354,] 8.714156e-01 2.571688e-01 1.285844e-01 [355,] 8.479915e-01 3.040170e-01 1.520085e-01 [356,] 8.270036e-01 3.459929e-01 1.729964e-01 [357,] 7.984924e-01 4.030153e-01 2.015076e-01 [358,] 7.672686e-01 4.654628e-01 2.327314e-01 [359,] 7.343527e-01 5.312946e-01 2.656473e-01 [360,] 6.976266e-01 6.047468e-01 3.023734e-01 [361,] 6.609061e-01 6.781879e-01 3.390939e-01 [362,] 6.676885e-01 6.646230e-01 3.323115e-01 [363,] 6.443026e-01 7.113948e-01 3.556974e-01 [364,] 6.183036e-01 7.633929e-01 3.816964e-01 [365,] 5.755672e-01 8.488656e-01 4.244328e-01 [366,] 5.307762e-01 9.384477e-01 4.692238e-01 [367,] 5.152458e-01 9.695084e-01 4.847542e-01 [368,] 4.697219e-01 9.394438e-01 5.302781e-01 [369,] 4.234798e-01 8.469595e-01 5.765202e-01 [370,] 3.783802e-01 7.567605e-01 6.216198e-01 [371,] 3.524361e-01 7.048722e-01 6.475639e-01 [372,] 3.600814e-01 7.201627e-01 6.399186e-01 [373,] 3.252119e-01 6.504238e-01 6.747881e-01 [374,] 3.455599e-01 6.911197e-01 6.544401e-01 [375,] 3.010644e-01 6.021289e-01 6.989356e-01 [376,] 2.669735e-01 5.339470e-01 7.330265e-01 [377,] 2.329101e-01 4.658202e-01 7.670899e-01 [378,] 1.973429e-01 3.946859e-01 8.026571e-01 [379,] 1.970760e-01 3.941520e-01 8.029240e-01 [380,] 1.680193e-01 3.360386e-01 8.319807e-01 [381,] 1.409200e-01 2.818401e-01 8.590800e-01 [382,] 1.309024e-01 2.618048e-01 8.690976e-01 [383,] 1.166711e-01 2.333422e-01 8.833289e-01 [384,] 1.026449e-01 2.052898e-01 8.973551e-01 [385,] 8.602611e-02 1.720522e-01 9.139739e-01 [386,] 6.587762e-02 1.317552e-01 9.341224e-01 [387,] 6.963542e-02 1.392708e-01 9.303646e-01 [388,] 5.493476e-02 1.098695e-01 9.450652e-01 [389,] 4.243365e-02 8.486731e-02 9.575663e-01 [390,] 3.060331e-02 6.120662e-02 9.693967e-01 [391,] 2.540188e-02 5.080376e-02 9.745981e-01 [392,] 2.809360e-02 5.618720e-02 9.719064e-01 [393,] 1.947636e-02 3.895273e-02 9.805236e-01 [394,] 1.391882e-02 2.783764e-02 9.860812e-01 [395,] 3.699661e-02 7.399321e-02 9.630034e-01 [396,] 4.086751e-02 8.173502e-02 9.591325e-01 [397,] 3.302705e-02 6.605411e-02 9.669729e-01 [398,] 2.565638e-02 5.131275e-02 9.743436e-01 [399,] 6.563556e-02 1.312711e-01 9.343644e-01 [400,] 6.732975e-02 1.346595e-01 9.326703e-01 [401,] 5.616855e-02 1.123371e-01 9.438315e-01 [402,] 4.745379e-02 9.490759e-02 9.525462e-01 [403,] 4.261369e-01 8.522738e-01 5.738631e-01 [404,] 3.127712e-01 6.255424e-01 6.872288e-01 [405,] 2.419994e-01 4.839989e-01 7.580006e-01 [406,] 8.007104e-01 3.985792e-01 1.992896e-01 > postscript(file="/var/www/html/rcomp/tmp/1j9zp1291404105.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/2cizs1291404105.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/3cizs1291404105.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/4cizs1291404105.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/55ayv1291404105.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > qqnorm(mysum$resid, main='Residual Normal Q-Q Plot') > qqline(mysum$resid) > grid() > dev.off() null device 1 > (myerror <- as.ts(mysum$resid)) Time Series: Start = 1 End = 431 Frequency = 1 1 2 3 4 5 6 145522.7289 90502.1320 66923.0051 -86588.0519 -9957.1434 5301.5072 7 8 9 10 11 12 -6984.8505 27911.8102 -32052.4725 -13358.0981 14749.1775 -24613.2783 13 14 15 16 17 18 -27441.0986 -61418.0626 -22453.5278 701.2505 18048.2723 -19297.6555 19 20 21 22 23 24 -37918.1849 -26870.2544 -22266.7383 -30281.8283 -8983.6467 -30141.1847 25 26 27 28 29 30 -24513.9376 -34642.5331 -35548.8347 -50745.5872 -32021.1376 -41651.8651 31 32 33 34 35 36 -20649.5177 -24708.8810 -22857.0338 -21649.0627 24594.2029 -33725.7474 37 38 39 40 41 42 -26144.5943 -17459.3018 -3520.9934 -37522.5937 -26415.2579 -33840.2753 43 44 45 46 47 48 -26916.1847 -25926.8864 -34990.0435 -27394.4139 -5741.5214 -20819.5735 49 50 51 52 53 54 -21726.7947 -38005.3901 -22400.1202 -17551.9752 -7914.5792 -13738.1365 55 56 57 58 59 60 -20758.5474 -35448.9761 -25118.4803 -29730.7906 -28708.5831 -33105.1848 61 62 63 64 65 66 -26821.5868 -42593.8217 -27867.0400 -34100.0918 -37426.6184 -34788.2364 67 68 69 70 71 72 -34676.1753 -27043.4415 -21098.1599 -27135.3218 -48544.7856 -41790.3067 73 74 75 76 77 78 -27716.9248 -31736.5474 -22903.1157 -35861.8158 -24369.5589 -36738.8726 79 80 81 82 83 84 -35892.2641 -24681.0853 -36208.5432 -29591.0053 -19700.5470 -41688.1648 85 86 87 88 89 90 -25327.6803 -37407.2420 -28324.1366 -25801.6609 -28313.3152 -24279.9608 91 92 93 94 95 96 -28174.4394 -31532.8118 -33684.4625 -24210.4439 -41412.3816 -36482.2282 97 98 99 100 101 102 -19310.2338 -35396.7236 -20001.7129 -33677.3486 -34298.2341 -33250.7074 103 104 105 106 107 108 -28716.6305 -25421.0106 -34363.1246 -36085.4182 -43890.0703 -22979.2732 109 110 111 112 113 114 -30448.0100 -36997.5478 -32781.7528 -45586.0873 -34337.2549 -35221.3719 115 116 117 118 119 120 -34189.7069 -37637.1046 -31185.0300 -27372.9765 -37220.6192 -37331.9891 121 122 123 124 125 126 -29075.8186 -34544.4095 -29827.1351 -22393.2265 -36380.6076 -37863.8889 127 128 129 130 131 132 -37298.4386 -6173.0898 -6962.1681 -7104.7631 -177750.9378 -201473.6058 133 134 135 136 137 138 137260.7680 244846.8956 253304.2055 260971.7333 258416.5632 242838.6944 139 140 141 142 143 144 253460.3053 5117.4458 -5058.6583 -38100.9019 -31646.2323 -35400.5289 145 146 147 148 149 150 -36970.3993 -37381.6194 -4478.2257 -4042.4752 -2253.6210 -8603.6966 151 152 153 154 155 156 -8400.5903 -7044.5213 -172978.1529 -174504.7953 139777.8833 52678.7444 157 158 159 160 161 162 17438.7996 -5059.5465 -37389.2550 -7115.1361 130693.6677 247870.8106 163 164 165 166 167 168 261059.3415 247849.6687 11834.1492 -297.5484 -5056.2529 -37353.4066 169 170 171 172 173 174 -6186.6401 -7080.6365 -177075.9749 139970.2411 255902.0738 262848.6397 175 176 177 178 179 180 217742.7744 -5059.4642 -33631.0470 -33048.6284 -39103.4120 -35312.1309 181 182 183 184 185 186 -33331.2475 -35900.1317 -35700.6713 -33948.9541 -34381.0366 -37118.4998 187 188 189 190 191 192 -5976.2294 -7093.7448 -178260.7668 140098.8532 240533.9467 19189.3050 193 194 195 196 197 198 -5050.5098 -36909.9004 -10901.9479 -7107.6279 139781.8136 -4799.9257 199 200 201 202 203 204 7845.3632 53266.6459 25597.5517 -5084.6181 -39052.2320 -37164.8190 205 206 207 208 209 210 16964.3181 15335.4581 -1580.6319 -330.8313 -1245.4368 -26671.5726 211 212 213 214 215 216 -24232.8131 -21453.6527 -7533.1168 -9410.3867 -9152.1501 30829.4584 217 218 219 220 221 222 17113.0973 26609.2914 17681.2257 18289.6287 16985.8485 -38479.9818 223 224 225 226 227 228 21615.5558 -1519.2978 -1804.3255 -27418.8472 -21474.3942 53583.3154 229 230 231 232 233 234 -9151.5882 20232.1426 16963.9264 18367.8233 -26677.4970 -22329.6993 235 236 237 238 239 240 -3054.3746 -4320.2660 -9147.4272 16987.7905 -41424.1582 18376.9622 241 242 243 244 245 246 -29240.2684 -21590.0785 -4735.5861 -10084.9169 -17525.3641 -9146.5316 247 248 249 250 251 252 17020.7058 12010.2079 16965.0754 18373.4906 -26624.0306 -2766.5122 253 254 255 256 257 258 -5430.6314 -9139.1231 16966.0318 35345.9045 -2133.8446 -24009.2340 259 260 261 262 263 264 -21919.8434 -9140.5445 16989.0825 -41036.5446 20160.3731 745.0774 265 266 267 268 269 270 -2826.9926 -27414.6001 -20521.9849 -21520.7434 -9144.1437 16285.5589 271 272 273 274 275 276 17008.3316 -42089.8303 18466.7302 536.0135 950.6720 -29418.2383 277 278 279 280 281 282 -21378.5236 -6363.9565 177.4511 -9143.8628 15414.8125 23274.8412 283 284 285 286 287 288 16964.6738 18372.8224 -26667.8846 -3092.5135 -6367.0712 -9149.1448 289 290 291 292 293 294 17051.0349 -38538.5964 -37984.2586 21702.9173 -8048.0187 -2563.8728 295 296 297 298 299 300 -26649.5438 -3114.9548 -9148.8318 21510.2352 26737.5365 24465.2693 301 302 303 304 305 306 34710.8968 14043.7096 31707.6553 16965.9895 18375.0442 -26689.1382 307 308 309 310 311 312 -9145.1630 16969.9571 -30592.9841 18375.8513 -26648.8823 -9141.0535 313 314 315 316 317 318 16966.0939 18379.0713 -26772.6413 -22133.0181 -67453.7499 -23154.5891 319 320 321 322 323 324 -21352.1309 5292.3718 -7049.6788 -11136.7996 -9153.6952 16963.6130 325 326 327 328 329 330 18468.2198 -2899.4911 -27781.3429 -21917.8511 -9155.2402 18854.9420 331 332 333 334 335 336 18701.6097 18734.9081 20607.1118 25678.7033 16965.9373 36524.8230 337 338 339 340 341 342 -3037.2008 -26662.4424 -9142.0368 45533.2814 16965.8328 19534.1656 343 344 345 346 347 348 -3687.8868 -1855.5881 -26668.5251 -9139.0871 17007.2156 -41149.0371 349 350 351 352 353 354 17516.3231 13003.7708 2518.4085 10379.0069 -26628.4330 -11111.1012 355 356 357 358 359 360 -6826.2858 -7369.5384 -9141.8963 17073.7591 -46463.1364 -39084.5284 361 362 363 364 365 366 20373.1199 5311.7414 -2961.7350 -2606.0918 -1801.2296 -346.0918 367 368 369 370 371 372 -276.1497 -3112.6701 -1148.7137 6128.1256 21973.0515 -4108.2105 373 374 375 376 377 378 -9135.6784 58702.6752 -1790.7359 -41022.5544 7965.6545 3116.4803 379 380 381 382 383 384 -3096.5854 7383.5677 1313.8368 8856.7682 5182.2422 -3978.2731 385 386 387 388 389 390 -4271.8223 7476.0284 8751.7002 676.9163 -5706.7481 12965.6691 391 392 393 394 395 396 14485.8632 11802.9822 -26263.4451 -2628.9920 6439.3599 3001.3730 397 398 399 400 401 402 8972.1360 13468.2362 -75518.0470 -1026.0610 -3687.6104 4866.2730 403 404 405 406 407 408 -10505.5715 102590.0857 17072.0931 16859.9296 -34651.0483 179140.9522 409 410 411 412 413 414 -65864.0831 42570.7665 -136327.6304 73670.7658 81076.8645 51272.4613 415 416 417 418 419 420 160407.0204 -23715.4828 -24239.7462 168333.9309 -25413.0401 -72861.4773 421 422 423 424 425 426 -42221.3985 36837.5340 78426.2163 213439.2120 49037.2086 -29128.7692 427 428 429 430 431 50784.6062 18434.2770 10938.3181 97818.2382 49147.3151 > postscript(file="/var/www/html/rcomp/tmp/65ayv1291404105.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 145522.7289 NA 1 90502.1320 145522.7289 2 66923.0051 90502.1320 3 -86588.0519 66923.0051 4 -9957.1434 -86588.0519 5 5301.5072 -9957.1434 6 -6984.8505 5301.5072 7 27911.8102 -6984.8505 8 -32052.4725 27911.8102 9 -13358.0981 -32052.4725 10 14749.1775 -13358.0981 11 -24613.2783 14749.1775 12 -27441.0986 -24613.2783 13 -61418.0626 -27441.0986 14 -22453.5278 -61418.0626 15 701.2505 -22453.5278 16 18048.2723 701.2505 17 -19297.6555 18048.2723 18 -37918.1849 -19297.6555 19 -26870.2544 -37918.1849 20 -22266.7383 -26870.2544 21 -30281.8283 -22266.7383 22 -8983.6467 -30281.8283 23 -30141.1847 -8983.6467 24 -24513.9376 -30141.1847 25 -34642.5331 -24513.9376 26 -35548.8347 -34642.5331 27 -50745.5872 -35548.8347 28 -32021.1376 -50745.5872 29 -41651.8651 -32021.1376 30 -20649.5177 -41651.8651 31 -24708.8810 -20649.5177 32 -22857.0338 -24708.8810 33 -21649.0627 -22857.0338 34 24594.2029 -21649.0627 35 -33725.7474 24594.2029 36 -26144.5943 -33725.7474 37 -17459.3018 -26144.5943 38 -3520.9934 -17459.3018 39 -37522.5937 -3520.9934 40 -26415.2579 -37522.5937 41 -33840.2753 -26415.2579 42 -26916.1847 -33840.2753 43 -25926.8864 -26916.1847 44 -34990.0435 -25926.8864 45 -27394.4139 -34990.0435 46 -5741.5214 -27394.4139 47 -20819.5735 -5741.5214 48 -21726.7947 -20819.5735 49 -38005.3901 -21726.7947 50 -22400.1202 -38005.3901 51 -17551.9752 -22400.1202 52 -7914.5792 -17551.9752 53 -13738.1365 -7914.5792 54 -20758.5474 -13738.1365 55 -35448.9761 -20758.5474 56 -25118.4803 -35448.9761 57 -29730.7906 -25118.4803 58 -28708.5831 -29730.7906 59 -33105.1848 -28708.5831 60 -26821.5868 -33105.1848 61 -42593.8217 -26821.5868 62 -27867.0400 -42593.8217 63 -34100.0918 -27867.0400 64 -37426.6184 -34100.0918 65 -34788.2364 -37426.6184 66 -34676.1753 -34788.2364 67 -27043.4415 -34676.1753 68 -21098.1599 -27043.4415 69 -27135.3218 -21098.1599 70 -48544.7856 -27135.3218 71 -41790.3067 -48544.7856 72 -27716.9248 -41790.3067 73 -31736.5474 -27716.9248 74 -22903.1157 -31736.5474 75 -35861.8158 -22903.1157 76 -24369.5589 -35861.8158 77 -36738.8726 -24369.5589 78 -35892.2641 -36738.8726 79 -24681.0853 -35892.2641 80 -36208.5432 -24681.0853 81 -29591.0053 -36208.5432 82 -19700.5470 -29591.0053 83 -41688.1648 -19700.5470 84 -25327.6803 -41688.1648 85 -37407.2420 -25327.6803 86 -28324.1366 -37407.2420 87 -25801.6609 -28324.1366 88 -28313.3152 -25801.6609 89 -24279.9608 -28313.3152 90 -28174.4394 -24279.9608 91 -31532.8118 -28174.4394 92 -33684.4625 -31532.8118 93 -24210.4439 -33684.4625 94 -41412.3816 -24210.4439 95 -36482.2282 -41412.3816 96 -19310.2338 -36482.2282 97 -35396.7236 -19310.2338 98 -20001.7129 -35396.7236 99 -33677.3486 -20001.7129 100 -34298.2341 -33677.3486 101 -33250.7074 -34298.2341 102 -28716.6305 -33250.7074 103 -25421.0106 -28716.6305 104 -34363.1246 -25421.0106 105 -36085.4182 -34363.1246 106 -43890.0703 -36085.4182 107 -22979.2732 -43890.0703 108 -30448.0100 -22979.2732 109 -36997.5478 -30448.0100 110 -32781.7528 -36997.5478 111 -45586.0873 -32781.7528 112 -34337.2549 -45586.0873 113 -35221.3719 -34337.2549 114 -34189.7069 -35221.3719 115 -37637.1046 -34189.7069 116 -31185.0300 -37637.1046 117 -27372.9765 -31185.0300 118 -37220.6192 -27372.9765 119 -37331.9891 -37220.6192 120 -29075.8186 -37331.9891 121 -34544.4095 -29075.8186 122 -29827.1351 -34544.4095 123 -22393.2265 -29827.1351 124 -36380.6076 -22393.2265 125 -37863.8889 -36380.6076 126 -37298.4386 -37863.8889 127 -6173.0898 -37298.4386 128 -6962.1681 -6173.0898 129 -7104.7631 -6962.1681 130 -177750.9378 -7104.7631 131 -201473.6058 -177750.9378 132 137260.7680 -201473.6058 133 244846.8956 137260.7680 134 253304.2055 244846.8956 135 260971.7333 253304.2055 136 258416.5632 260971.7333 137 242838.6944 258416.5632 138 253460.3053 242838.6944 139 5117.4458 253460.3053 140 -5058.6583 5117.4458 141 -38100.9019 -5058.6583 142 -31646.2323 -38100.9019 143 -35400.5289 -31646.2323 144 -36970.3993 -35400.5289 145 -37381.6194 -36970.3993 146 -4478.2257 -37381.6194 147 -4042.4752 -4478.2257 148 -2253.6210 -4042.4752 149 -8603.6966 -2253.6210 150 -8400.5903 -8603.6966 151 -7044.5213 -8400.5903 152 -172978.1529 -7044.5213 153 -174504.7953 -172978.1529 154 139777.8833 -174504.7953 155 52678.7444 139777.8833 156 17438.7996 52678.7444 157 -5059.5465 17438.7996 158 -37389.2550 -5059.5465 159 -7115.1361 -37389.2550 160 130693.6677 -7115.1361 161 247870.8106 130693.6677 162 261059.3415 247870.8106 163 247849.6687 261059.3415 164 11834.1492 247849.6687 165 -297.5484 11834.1492 166 -5056.2529 -297.5484 167 -37353.4066 -5056.2529 168 -6186.6401 -37353.4066 169 -7080.6365 -6186.6401 170 -177075.9749 -7080.6365 171 139970.2411 -177075.9749 172 255902.0738 139970.2411 173 262848.6397 255902.0738 174 217742.7744 262848.6397 175 -5059.4642 217742.7744 176 -33631.0470 -5059.4642 177 -33048.6284 -33631.0470 178 -39103.4120 -33048.6284 179 -35312.1309 -39103.4120 180 -33331.2475 -35312.1309 181 -35900.1317 -33331.2475 182 -35700.6713 -35900.1317 183 -33948.9541 -35700.6713 184 -34381.0366 -33948.9541 185 -37118.4998 -34381.0366 186 -5976.2294 -37118.4998 187 -7093.7448 -5976.2294 188 -178260.7668 -7093.7448 189 140098.8532 -178260.7668 190 240533.9467 140098.8532 191 19189.3050 240533.9467 192 -5050.5098 19189.3050 193 -36909.9004 -5050.5098 194 -10901.9479 -36909.9004 195 -7107.6279 -10901.9479 196 139781.8136 -7107.6279 197 -4799.9257 139781.8136 198 7845.3632 -4799.9257 199 53266.6459 7845.3632 200 25597.5517 53266.6459 201 -5084.6181 25597.5517 202 -39052.2320 -5084.6181 203 -37164.8190 -39052.2320 204 16964.3181 -37164.8190 205 15335.4581 16964.3181 206 -1580.6319 15335.4581 207 -330.8313 -1580.6319 208 -1245.4368 -330.8313 209 -26671.5726 -1245.4368 210 -24232.8131 -26671.5726 211 -21453.6527 -24232.8131 212 -7533.1168 -21453.6527 213 -9410.3867 -7533.1168 214 -9152.1501 -9410.3867 215 30829.4584 -9152.1501 216 17113.0973 30829.4584 217 26609.2914 17113.0973 218 17681.2257 26609.2914 219 18289.6287 17681.2257 220 16985.8485 18289.6287 221 -38479.9818 16985.8485 222 21615.5558 -38479.9818 223 -1519.2978 21615.5558 224 -1804.3255 -1519.2978 225 -27418.8472 -1804.3255 226 -21474.3942 -27418.8472 227 53583.3154 -21474.3942 228 -9151.5882 53583.3154 229 20232.1426 -9151.5882 230 16963.9264 20232.1426 231 18367.8233 16963.9264 232 -26677.4970 18367.8233 233 -22329.6993 -26677.4970 234 -3054.3746 -22329.6993 235 -4320.2660 -3054.3746 236 -9147.4272 -4320.2660 237 16987.7905 -9147.4272 238 -41424.1582 16987.7905 239 18376.9622 -41424.1582 240 -29240.2684 18376.9622 241 -21590.0785 -29240.2684 242 -4735.5861 -21590.0785 243 -10084.9169 -4735.5861 244 -17525.3641 -10084.9169 245 -9146.5316 -17525.3641 246 17020.7058 -9146.5316 247 12010.2079 17020.7058 248 16965.0754 12010.2079 249 18373.4906 16965.0754 250 -26624.0306 18373.4906 251 -2766.5122 -26624.0306 252 -5430.6314 -2766.5122 253 -9139.1231 -5430.6314 254 16966.0318 -9139.1231 255 35345.9045 16966.0318 256 -2133.8446 35345.9045 257 -24009.2340 -2133.8446 258 -21919.8434 -24009.2340 259 -9140.5445 -21919.8434 260 16989.0825 -9140.5445 261 -41036.5446 16989.0825 262 20160.3731 -41036.5446 263 745.0774 20160.3731 264 -2826.9926 745.0774 265 -27414.6001 -2826.9926 266 -20521.9849 -27414.6001 267 -21520.7434 -20521.9849 268 -9144.1437 -21520.7434 269 16285.5589 -9144.1437 270 17008.3316 16285.5589 271 -42089.8303 17008.3316 272 18466.7302 -42089.8303 273 536.0135 18466.7302 274 950.6720 536.0135 275 -29418.2383 950.6720 276 -21378.5236 -29418.2383 277 -6363.9565 -21378.5236 278 177.4511 -6363.9565 279 -9143.8628 177.4511 280 15414.8125 -9143.8628 281 23274.8412 15414.8125 282 16964.6738 23274.8412 283 18372.8224 16964.6738 284 -26667.8846 18372.8224 285 -3092.5135 -26667.8846 286 -6367.0712 -3092.5135 287 -9149.1448 -6367.0712 288 17051.0349 -9149.1448 289 -38538.5964 17051.0349 290 -37984.2586 -38538.5964 291 21702.9173 -37984.2586 292 -8048.0187 21702.9173 293 -2563.8728 -8048.0187 294 -26649.5438 -2563.8728 295 -3114.9548 -26649.5438 296 -9148.8318 -3114.9548 297 21510.2352 -9148.8318 298 26737.5365 21510.2352 299 24465.2693 26737.5365 300 34710.8968 24465.2693 301 14043.7096 34710.8968 302 31707.6553 14043.7096 303 16965.9895 31707.6553 304 18375.0442 16965.9895 305 -26689.1382 18375.0442 306 -9145.1630 -26689.1382 307 16969.9571 -9145.1630 308 -30592.9841 16969.9571 309 18375.8513 -30592.9841 310 -26648.8823 18375.8513 311 -9141.0535 -26648.8823 312 16966.0939 -9141.0535 313 18379.0713 16966.0939 314 -26772.6413 18379.0713 315 -22133.0181 -26772.6413 316 -67453.7499 -22133.0181 317 -23154.5891 -67453.7499 318 -21352.1309 -23154.5891 319 5292.3718 -21352.1309 320 -7049.6788 5292.3718 321 -11136.7996 -7049.6788 322 -9153.6952 -11136.7996 323 16963.6130 -9153.6952 324 18468.2198 16963.6130 325 -2899.4911 18468.2198 326 -27781.3429 -2899.4911 327 -21917.8511 -27781.3429 328 -9155.2402 -21917.8511 329 18854.9420 -9155.2402 330 18701.6097 18854.9420 331 18734.9081 18701.6097 332 20607.1118 18734.9081 333 25678.7033 20607.1118 334 16965.9373 25678.7033 335 36524.8230 16965.9373 336 -3037.2008 36524.8230 337 -26662.4424 -3037.2008 338 -9142.0368 -26662.4424 339 45533.2814 -9142.0368 340 16965.8328 45533.2814 341 19534.1656 16965.8328 342 -3687.8868 19534.1656 343 -1855.5881 -3687.8868 344 -26668.5251 -1855.5881 345 -9139.0871 -26668.5251 346 17007.2156 -9139.0871 347 -41149.0371 17007.2156 348 17516.3231 -41149.0371 349 13003.7708 17516.3231 350 2518.4085 13003.7708 351 10379.0069 2518.4085 352 -26628.4330 10379.0069 353 -11111.1012 -26628.4330 354 -6826.2858 -11111.1012 355 -7369.5384 -6826.2858 356 -9141.8963 -7369.5384 357 17073.7591 -9141.8963 358 -46463.1364 17073.7591 359 -39084.5284 -46463.1364 360 20373.1199 -39084.5284 361 5311.7414 20373.1199 362 -2961.7350 5311.7414 363 -2606.0918 -2961.7350 364 -1801.2296 -2606.0918 365 -346.0918 -1801.2296 366 -276.1497 -346.0918 367 -3112.6701 -276.1497 368 -1148.7137 -3112.6701 369 6128.1256 -1148.7137 370 21973.0515 6128.1256 371 -4108.2105 21973.0515 372 -9135.6784 -4108.2105 373 58702.6752 -9135.6784 374 -1790.7359 58702.6752 375 -41022.5544 -1790.7359 376 7965.6545 -41022.5544 377 3116.4803 7965.6545 378 -3096.5854 3116.4803 379 7383.5677 -3096.5854 380 1313.8368 7383.5677 381 8856.7682 1313.8368 382 5182.2422 8856.7682 383 -3978.2731 5182.2422 384 -4271.8223 -3978.2731 385 7476.0284 -4271.8223 386 8751.7002 7476.0284 387 676.9163 8751.7002 388 -5706.7481 676.9163 389 12965.6691 -5706.7481 390 14485.8632 12965.6691 391 11802.9822 14485.8632 392 -26263.4451 11802.9822 393 -2628.9920 -26263.4451 394 6439.3599 -2628.9920 395 3001.3730 6439.3599 396 8972.1360 3001.3730 397 13468.2362 8972.1360 398 -75518.0470 13468.2362 399 -1026.0610 -75518.0470 400 -3687.6104 -1026.0610 401 4866.2730 -3687.6104 402 -10505.5715 4866.2730 403 102590.0857 -10505.5715 404 17072.0931 102590.0857 405 16859.9296 17072.0931 406 -34651.0483 16859.9296 407 179140.9522 -34651.0483 408 -65864.0831 179140.9522 409 42570.7665 -65864.0831 410 -136327.6304 42570.7665 411 73670.7658 -136327.6304 412 81076.8645 73670.7658 413 51272.4613 81076.8645 414 160407.0204 51272.4613 415 -23715.4828 160407.0204 416 -24239.7462 -23715.4828 417 168333.9309 -24239.7462 418 -25413.0401 168333.9309 419 -72861.4773 -25413.0401 420 -42221.3985 -72861.4773 421 36837.5340 -42221.3985 422 78426.2163 36837.5340 423 213439.2120 78426.2163 424 49037.2086 213439.2120 425 -29128.7692 49037.2086 426 50784.6062 -29128.7692 427 18434.2770 50784.6062 428 10938.3181 18434.2770 429 97818.2382 10938.3181 430 49147.3151 97818.2382 431 NA 49147.3151 > dum1 <- dum[2:length(myerror),] > dum1 lag(myerror, k = 1) myerror [1,] 90502.1320 145522.7289 [2,] 66923.0051 90502.1320 [3,] -86588.0519 66923.0051 [4,] -9957.1434 -86588.0519 [5,] 5301.5072 -9957.1434 [6,] -6984.8505 5301.5072 [7,] 27911.8102 -6984.8505 [8,] -32052.4725 27911.8102 [9,] -13358.0981 -32052.4725 [10,] 14749.1775 -13358.0981 [11,] -24613.2783 14749.1775 [12,] -27441.0986 -24613.2783 [13,] -61418.0626 -27441.0986 [14,] -22453.5278 -61418.0626 [15,] 701.2505 -22453.5278 [16,] 18048.2723 701.2505 [17,] -19297.6555 18048.2723 [18,] -37918.1849 -19297.6555 [19,] -26870.2544 -37918.1849 [20,] -22266.7383 -26870.2544 [21,] -30281.8283 -22266.7383 [22,] -8983.6467 -30281.8283 [23,] -30141.1847 -8983.6467 [24,] -24513.9376 -30141.1847 [25,] -34642.5331 -24513.9376 [26,] -35548.8347 -34642.5331 [27,] -50745.5872 -35548.8347 [28,] -32021.1376 -50745.5872 [29,] -41651.8651 -32021.1376 [30,] -20649.5177 -41651.8651 [31,] -24708.8810 -20649.5177 [32,] -22857.0338 -24708.8810 [33,] -21649.0627 -22857.0338 [34,] 24594.2029 -21649.0627 [35,] -33725.7474 24594.2029 [36,] -26144.5943 -33725.7474 [37,] -17459.3018 -26144.5943 [38,] -3520.9934 -17459.3018 [39,] -37522.5937 -3520.9934 [40,] -26415.2579 -37522.5937 [41,] -33840.2753 -26415.2579 [42,] -26916.1847 -33840.2753 [43,] -25926.8864 -26916.1847 [44,] -34990.0435 -25926.8864 [45,] -27394.4139 -34990.0435 [46,] -5741.5214 -27394.4139 [47,] -20819.5735 -5741.5214 [48,] -21726.7947 -20819.5735 [49,] -38005.3901 -21726.7947 [50,] -22400.1202 -38005.3901 [51,] -17551.9752 -22400.1202 [52,] -7914.5792 -17551.9752 [53,] -13738.1365 -7914.5792 [54,] -20758.5474 -13738.1365 [55,] -35448.9761 -20758.5474 [56,] -25118.4803 -35448.9761 [57,] -29730.7906 -25118.4803 [58,] -28708.5831 -29730.7906 [59,] -33105.1848 -28708.5831 [60,] -26821.5868 -33105.1848 [61,] -42593.8217 -26821.5868 [62,] -27867.0400 -42593.8217 [63,] -34100.0918 -27867.0400 [64,] -37426.6184 -34100.0918 [65,] -34788.2364 -37426.6184 [66,] -34676.1753 -34788.2364 [67,] -27043.4415 -34676.1753 [68,] -21098.1599 -27043.4415 [69,] -27135.3218 -21098.1599 [70,] -48544.7856 -27135.3218 [71,] -41790.3067 -48544.7856 [72,] -27716.9248 -41790.3067 [73,] -31736.5474 -27716.9248 [74,] -22903.1157 -31736.5474 [75,] -35861.8158 -22903.1157 [76,] -24369.5589 -35861.8158 [77,] -36738.8726 -24369.5589 [78,] -35892.2641 -36738.8726 [79,] -24681.0853 -35892.2641 [80,] -36208.5432 -24681.0853 [81,] -29591.0053 -36208.5432 [82,] -19700.5470 -29591.0053 [83,] -41688.1648 -19700.5470 [84,] -25327.6803 -41688.1648 [85,] -37407.2420 -25327.6803 [86,] -28324.1366 -37407.2420 [87,] -25801.6609 -28324.1366 [88,] -28313.3152 -25801.6609 [89,] -24279.9608 -28313.3152 [90,] -28174.4394 -24279.9608 [91,] -31532.8118 -28174.4394 [92,] -33684.4625 -31532.8118 [93,] -24210.4439 -33684.4625 [94,] -41412.3816 -24210.4439 [95,] -36482.2282 -41412.3816 [96,] -19310.2338 -36482.2282 [97,] -35396.7236 -19310.2338 [98,] -20001.7129 -35396.7236 [99,] -33677.3486 -20001.7129 [100,] -34298.2341 -33677.3486 [101,] -33250.7074 -34298.2341 [102,] -28716.6305 -33250.7074 [103,] -25421.0106 -28716.6305 [104,] -34363.1246 -25421.0106 [105,] -36085.4182 -34363.1246 [106,] -43890.0703 -36085.4182 [107,] -22979.2732 -43890.0703 [108,] -30448.0100 -22979.2732 [109,] -36997.5478 -30448.0100 [110,] -32781.7528 -36997.5478 [111,] -45586.0873 -32781.7528 [112,] -34337.2549 -45586.0873 [113,] -35221.3719 -34337.2549 [114,] -34189.7069 -35221.3719 [115,] -37637.1046 -34189.7069 [116,] -31185.0300 -37637.1046 [117,] -27372.9765 -31185.0300 [118,] -37220.6192 -27372.9765 [119,] -37331.9891 -37220.6192 [120,] -29075.8186 -37331.9891 [121,] -34544.4095 -29075.8186 [122,] -29827.1351 -34544.4095 [123,] -22393.2265 -29827.1351 [124,] -36380.6076 -22393.2265 [125,] -37863.8889 -36380.6076 [126,] -37298.4386 -37863.8889 [127,] -6173.0898 -37298.4386 [128,] -6962.1681 -6173.0898 [129,] -7104.7631 -6962.1681 [130,] -177750.9378 -7104.7631 [131,] -201473.6058 -177750.9378 [132,] 137260.7680 -201473.6058 [133,] 244846.8956 137260.7680 [134,] 253304.2055 244846.8956 [135,] 260971.7333 253304.2055 [136,] 258416.5632 260971.7333 [137,] 242838.6944 258416.5632 [138,] 253460.3053 242838.6944 [139,] 5117.4458 253460.3053 [140,] -5058.6583 5117.4458 [141,] -38100.9019 -5058.6583 [142,] -31646.2323 -38100.9019 [143,] -35400.5289 -31646.2323 [144,] -36970.3993 -35400.5289 [145,] -37381.6194 -36970.3993 [146,] -4478.2257 -37381.6194 [147,] -4042.4752 -4478.2257 [148,] -2253.6210 -4042.4752 [149,] -8603.6966 -2253.6210 [150,] -8400.5903 -8603.6966 [151,] -7044.5213 -8400.5903 [152,] -172978.1529 -7044.5213 [153,] -174504.7953 -172978.1529 [154,] 139777.8833 -174504.7953 [155,] 52678.7444 139777.8833 [156,] 17438.7996 52678.7444 [157,] -5059.5465 17438.7996 [158,] -37389.2550 -5059.5465 [159,] -7115.1361 -37389.2550 [160,] 130693.6677 -7115.1361 [161,] 247870.8106 130693.6677 [162,] 261059.3415 247870.8106 [163,] 247849.6687 261059.3415 [164,] 11834.1492 247849.6687 [165,] -297.5484 11834.1492 [166,] -5056.2529 -297.5484 [167,] -37353.4066 -5056.2529 [168,] -6186.6401 -37353.4066 [169,] -7080.6365 -6186.6401 [170,] -177075.9749 -7080.6365 [171,] 139970.2411 -177075.9749 [172,] 255902.0738 139970.2411 [173,] 262848.6397 255902.0738 [174,] 217742.7744 262848.6397 [175,] -5059.4642 217742.7744 [176,] -33631.0470 -5059.4642 [177,] -33048.6284 -33631.0470 [178,] -39103.4120 -33048.6284 [179,] -35312.1309 -39103.4120 [180,] -33331.2475 -35312.1309 [181,] -35900.1317 -33331.2475 [182,] -35700.6713 -35900.1317 [183,] -33948.9541 -35700.6713 [184,] -34381.0366 -33948.9541 [185,] -37118.4998 -34381.0366 [186,] -5976.2294 -37118.4998 [187,] -7093.7448 -5976.2294 [188,] -178260.7668 -7093.7448 [189,] 140098.8532 -178260.7668 [190,] 240533.9467 140098.8532 [191,] 19189.3050 240533.9467 [192,] -5050.5098 19189.3050 [193,] -36909.9004 -5050.5098 [194,] -10901.9479 -36909.9004 [195,] -7107.6279 -10901.9479 [196,] 139781.8136 -7107.6279 [197,] -4799.9257 139781.8136 [198,] 7845.3632 -4799.9257 [199,] 53266.6459 7845.3632 [200,] 25597.5517 53266.6459 [201,] -5084.6181 25597.5517 [202,] -39052.2320 -5084.6181 [203,] -37164.8190 -39052.2320 [204,] 16964.3181 -37164.8190 [205,] 15335.4581 16964.3181 [206,] -1580.6319 15335.4581 [207,] -330.8313 -1580.6319 [208,] -1245.4368 -330.8313 [209,] -26671.5726 -1245.4368 [210,] -24232.8131 -26671.5726 [211,] -21453.6527 -24232.8131 [212,] -7533.1168 -21453.6527 [213,] -9410.3867 -7533.1168 [214,] -9152.1501 -9410.3867 [215,] 30829.4584 -9152.1501 [216,] 17113.0973 30829.4584 [217,] 26609.2914 17113.0973 [218,] 17681.2257 26609.2914 [219,] 18289.6287 17681.2257 [220,] 16985.8485 18289.6287 [221,] -38479.9818 16985.8485 [222,] 21615.5558 -38479.9818 [223,] -1519.2978 21615.5558 [224,] -1804.3255 -1519.2978 [225,] -27418.8472 -1804.3255 [226,] -21474.3942 -27418.8472 [227,] 53583.3154 -21474.3942 [228,] -9151.5882 53583.3154 [229,] 20232.1426 -9151.5882 [230,] 16963.9264 20232.1426 [231,] 18367.8233 16963.9264 [232,] -26677.4970 18367.8233 [233,] -22329.6993 -26677.4970 [234,] -3054.3746 -22329.6993 [235,] -4320.2660 -3054.3746 [236,] -9147.4272 -4320.2660 [237,] 16987.7905 -9147.4272 [238,] -41424.1582 16987.7905 [239,] 18376.9622 -41424.1582 [240,] -29240.2684 18376.9622 [241,] -21590.0785 -29240.2684 [242,] -4735.5861 -21590.0785 [243,] -10084.9169 -4735.5861 [244,] -17525.3641 -10084.9169 [245,] -9146.5316 -17525.3641 [246,] 17020.7058 -9146.5316 [247,] 12010.2079 17020.7058 [248,] 16965.0754 12010.2079 [249,] 18373.4906 16965.0754 [250,] -26624.0306 18373.4906 [251,] -2766.5122 -26624.0306 [252,] -5430.6314 -2766.5122 [253,] -9139.1231 -5430.6314 [254,] 16966.0318 -9139.1231 [255,] 35345.9045 16966.0318 [256,] -2133.8446 35345.9045 [257,] -24009.2340 -2133.8446 [258,] -21919.8434 -24009.2340 [259,] -9140.5445 -21919.8434 [260,] 16989.0825 -9140.5445 [261,] -41036.5446 16989.0825 [262,] 20160.3731 -41036.5446 [263,] 745.0774 20160.3731 [264,] -2826.9926 745.0774 [265,] -27414.6001 -2826.9926 [266,] -20521.9849 -27414.6001 [267,] -21520.7434 -20521.9849 [268,] -9144.1437 -21520.7434 [269,] 16285.5589 -9144.1437 [270,] 17008.3316 16285.5589 [271,] -42089.8303 17008.3316 [272,] 18466.7302 -42089.8303 [273,] 536.0135 18466.7302 [274,] 950.6720 536.0135 [275,] -29418.2383 950.6720 [276,] -21378.5236 -29418.2383 [277,] -6363.9565 -21378.5236 [278,] 177.4511 -6363.9565 [279,] -9143.8628 177.4511 [280,] 15414.8125 -9143.8628 [281,] 23274.8412 15414.8125 [282,] 16964.6738 23274.8412 [283,] 18372.8224 16964.6738 [284,] -26667.8846 18372.8224 [285,] -3092.5135 -26667.8846 [286,] -6367.0712 -3092.5135 [287,] -9149.1448 -6367.0712 [288,] 17051.0349 -9149.1448 [289,] -38538.5964 17051.0349 [290,] -37984.2586 -38538.5964 [291,] 21702.9173 -37984.2586 [292,] -8048.0187 21702.9173 [293,] -2563.8728 -8048.0187 [294,] -26649.5438 -2563.8728 [295,] -3114.9548 -26649.5438 [296,] -9148.8318 -3114.9548 [297,] 21510.2352 -9148.8318 [298,] 26737.5365 21510.2352 [299,] 24465.2693 26737.5365 [300,] 34710.8968 24465.2693 [301,] 14043.7096 34710.8968 [302,] 31707.6553 14043.7096 [303,] 16965.9895 31707.6553 [304,] 18375.0442 16965.9895 [305,] -26689.1382 18375.0442 [306,] -9145.1630 -26689.1382 [307,] 16969.9571 -9145.1630 [308,] -30592.9841 16969.9571 [309,] 18375.8513 -30592.9841 [310,] -26648.8823 18375.8513 [311,] -9141.0535 -26648.8823 [312,] 16966.0939 -9141.0535 [313,] 18379.0713 16966.0939 [314,] -26772.6413 18379.0713 [315,] -22133.0181 -26772.6413 [316,] -67453.7499 -22133.0181 [317,] -23154.5891 -67453.7499 [318,] -21352.1309 -23154.5891 [319,] 5292.3718 -21352.1309 [320,] -7049.6788 5292.3718 [321,] -11136.7996 -7049.6788 [322,] -9153.6952 -11136.7996 [323,] 16963.6130 -9153.6952 [324,] 18468.2198 16963.6130 [325,] -2899.4911 18468.2198 [326,] -27781.3429 -2899.4911 [327,] -21917.8511 -27781.3429 [328,] -9155.2402 -21917.8511 [329,] 18854.9420 -9155.2402 [330,] 18701.6097 18854.9420 [331,] 18734.9081 18701.6097 [332,] 20607.1118 18734.9081 [333,] 25678.7033 20607.1118 [334,] 16965.9373 25678.7033 [335,] 36524.8230 16965.9373 [336,] -3037.2008 36524.8230 [337,] -26662.4424 -3037.2008 [338,] -9142.0368 -26662.4424 [339,] 45533.2814 -9142.0368 [340,] 16965.8328 45533.2814 [341,] 19534.1656 16965.8328 [342,] -3687.8868 19534.1656 [343,] -1855.5881 -3687.8868 [344,] -26668.5251 -1855.5881 [345,] -9139.0871 -26668.5251 [346,] 17007.2156 -9139.0871 [347,] -41149.0371 17007.2156 [348,] 17516.3231 -41149.0371 [349,] 13003.7708 17516.3231 [350,] 2518.4085 13003.7708 [351,] 10379.0069 2518.4085 [352,] -26628.4330 10379.0069 [353,] -11111.1012 -26628.4330 [354,] -6826.2858 -11111.1012 [355,] -7369.5384 -6826.2858 [356,] -9141.8963 -7369.5384 [357,] 17073.7591 -9141.8963 [358,] -46463.1364 17073.7591 [359,] -39084.5284 -46463.1364 [360,] 20373.1199 -39084.5284 [361,] 5311.7414 20373.1199 [362,] -2961.7350 5311.7414 [363,] -2606.0918 -2961.7350 [364,] -1801.2296 -2606.0918 [365,] -346.0918 -1801.2296 [366,] -276.1497 -346.0918 [367,] -3112.6701 -276.1497 [368,] -1148.7137 -3112.6701 [369,] 6128.1256 -1148.7137 [370,] 21973.0515 6128.1256 [371,] -4108.2105 21973.0515 [372,] -9135.6784 -4108.2105 [373,] 58702.6752 -9135.6784 [374,] -1790.7359 58702.6752 [375,] -41022.5544 -1790.7359 [376,] 7965.6545 -41022.5544 [377,] 3116.4803 7965.6545 [378,] -3096.5854 3116.4803 [379,] 7383.5677 -3096.5854 [380,] 1313.8368 7383.5677 [381,] 8856.7682 1313.8368 [382,] 5182.2422 8856.7682 [383,] -3978.2731 5182.2422 [384,] -4271.8223 -3978.2731 [385,] 7476.0284 -4271.8223 [386,] 8751.7002 7476.0284 [387,] 676.9163 8751.7002 [388,] -5706.7481 676.9163 [389,] 12965.6691 -5706.7481 [390,] 14485.8632 12965.6691 [391,] 11802.9822 14485.8632 [392,] -26263.4451 11802.9822 [393,] -2628.9920 -26263.4451 [394,] 6439.3599 -2628.9920 [395,] 3001.3730 6439.3599 [396,] 8972.1360 3001.3730 [397,] 13468.2362 8972.1360 [398,] -75518.0470 13468.2362 [399,] -1026.0610 -75518.0470 [400,] -3687.6104 -1026.0610 [401,] 4866.2730 -3687.6104 [402,] -10505.5715 4866.2730 [403,] 102590.0857 -10505.5715 [404,] 17072.0931 102590.0857 [405,] 16859.9296 17072.0931 [406,] -34651.0483 16859.9296 [407,] 179140.9522 -34651.0483 [408,] -65864.0831 179140.9522 [409,] 42570.7665 -65864.0831 [410,] -136327.6304 42570.7665 [411,] 73670.7658 -136327.6304 [412,] 81076.8645 73670.7658 [413,] 51272.4613 81076.8645 [414,] 160407.0204 51272.4613 [415,] -23715.4828 160407.0204 [416,] -24239.7462 -23715.4828 [417,] 168333.9309 -24239.7462 [418,] -25413.0401 168333.9309 [419,] -72861.4773 -25413.0401 [420,] -42221.3985 -72861.4773 [421,] 36837.5340 -42221.3985 [422,] 78426.2163 36837.5340 [423,] 213439.2120 78426.2163 [424,] 49037.2086 213439.2120 [425,] -29128.7692 49037.2086 [426,] 50784.6062 -29128.7692 [427,] 18434.2770 50784.6062 [428,] 10938.3181 18434.2770 [429,] 97818.2382 10938.3181 [430,] 49147.3151 97818.2382 > z <- as.data.frame(dum1) > z lag(myerror, k = 1) myerror 1 90502.1320 145522.7289 2 66923.0051 90502.1320 3 -86588.0519 66923.0051 4 -9957.1434 -86588.0519 5 5301.5072 -9957.1434 6 -6984.8505 5301.5072 7 27911.8102 -6984.8505 8 -32052.4725 27911.8102 9 -13358.0981 -32052.4725 10 14749.1775 -13358.0981 11 -24613.2783 14749.1775 12 -27441.0986 -24613.2783 13 -61418.0626 -27441.0986 14 -22453.5278 -61418.0626 15 701.2505 -22453.5278 16 18048.2723 701.2505 17 -19297.6555 18048.2723 18 -37918.1849 -19297.6555 19 -26870.2544 -37918.1849 20 -22266.7383 -26870.2544 21 -30281.8283 -22266.7383 22 -8983.6467 -30281.8283 23 -30141.1847 -8983.6467 24 -24513.9376 -30141.1847 25 -34642.5331 -24513.9376 26 -35548.8347 -34642.5331 27 -50745.5872 -35548.8347 28 -32021.1376 -50745.5872 29 -41651.8651 -32021.1376 30 -20649.5177 -41651.8651 31 -24708.8810 -20649.5177 32 -22857.0338 -24708.8810 33 -21649.0627 -22857.0338 34 24594.2029 -21649.0627 35 -33725.7474 24594.2029 36 -26144.5943 -33725.7474 37 -17459.3018 -26144.5943 38 -3520.9934 -17459.3018 39 -37522.5937 -3520.9934 40 -26415.2579 -37522.5937 41 -33840.2753 -26415.2579 42 -26916.1847 -33840.2753 43 -25926.8864 -26916.1847 44 -34990.0435 -25926.8864 45 -27394.4139 -34990.0435 46 -5741.5214 -27394.4139 47 -20819.5735 -5741.5214 48 -21726.7947 -20819.5735 49 -38005.3901 -21726.7947 50 -22400.1202 -38005.3901 51 -17551.9752 -22400.1202 52 -7914.5792 -17551.9752 53 -13738.1365 -7914.5792 54 -20758.5474 -13738.1365 55 -35448.9761 -20758.5474 56 -25118.4803 -35448.9761 57 -29730.7906 -25118.4803 58 -28708.5831 -29730.7906 59 -33105.1848 -28708.5831 60 -26821.5868 -33105.1848 61 -42593.8217 -26821.5868 62 -27867.0400 -42593.8217 63 -34100.0918 -27867.0400 64 -37426.6184 -34100.0918 65 -34788.2364 -37426.6184 66 -34676.1753 -34788.2364 67 -27043.4415 -34676.1753 68 -21098.1599 -27043.4415 69 -27135.3218 -21098.1599 70 -48544.7856 -27135.3218 71 -41790.3067 -48544.7856 72 -27716.9248 -41790.3067 73 -31736.5474 -27716.9248 74 -22903.1157 -31736.5474 75 -35861.8158 -22903.1157 76 -24369.5589 -35861.8158 77 -36738.8726 -24369.5589 78 -35892.2641 -36738.8726 79 -24681.0853 -35892.2641 80 -36208.5432 -24681.0853 81 -29591.0053 -36208.5432 82 -19700.5470 -29591.0053 83 -41688.1648 -19700.5470 84 -25327.6803 -41688.1648 85 -37407.2420 -25327.6803 86 -28324.1366 -37407.2420 87 -25801.6609 -28324.1366 88 -28313.3152 -25801.6609 89 -24279.9608 -28313.3152 90 -28174.4394 -24279.9608 91 -31532.8118 -28174.4394 92 -33684.4625 -31532.8118 93 -24210.4439 -33684.4625 94 -41412.3816 -24210.4439 95 -36482.2282 -41412.3816 96 -19310.2338 -36482.2282 97 -35396.7236 -19310.2338 98 -20001.7129 -35396.7236 99 -33677.3486 -20001.7129 100 -34298.2341 -33677.3486 101 -33250.7074 -34298.2341 102 -28716.6305 -33250.7074 103 -25421.0106 -28716.6305 104 -34363.1246 -25421.0106 105 -36085.4182 -34363.1246 106 -43890.0703 -36085.4182 107 -22979.2732 -43890.0703 108 -30448.0100 -22979.2732 109 -36997.5478 -30448.0100 110 -32781.7528 -36997.5478 111 -45586.0873 -32781.7528 112 -34337.2549 -45586.0873 113 -35221.3719 -34337.2549 114 -34189.7069 -35221.3719 115 -37637.1046 -34189.7069 116 -31185.0300 -37637.1046 117 -27372.9765 -31185.0300 118 -37220.6192 -27372.9765 119 -37331.9891 -37220.6192 120 -29075.8186 -37331.9891 121 -34544.4095 -29075.8186 122 -29827.1351 -34544.4095 123 -22393.2265 -29827.1351 124 -36380.6076 -22393.2265 125 -37863.8889 -36380.6076 126 -37298.4386 -37863.8889 127 -6173.0898 -37298.4386 128 -6962.1681 -6173.0898 129 -7104.7631 -6962.1681 130 -177750.9378 -7104.7631 131 -201473.6058 -177750.9378 132 137260.7680 -201473.6058 133 244846.8956 137260.7680 134 253304.2055 244846.8956 135 260971.7333 253304.2055 136 258416.5632 260971.7333 137 242838.6944 258416.5632 138 253460.3053 242838.6944 139 5117.4458 253460.3053 140 -5058.6583 5117.4458 141 -38100.9019 -5058.6583 142 -31646.2323 -38100.9019 143 -35400.5289 -31646.2323 144 -36970.3993 -35400.5289 145 -37381.6194 -36970.3993 146 -4478.2257 -37381.6194 147 -4042.4752 -4478.2257 148 -2253.6210 -4042.4752 149 -8603.6966 -2253.6210 150 -8400.5903 -8603.6966 151 -7044.5213 -8400.5903 152 -172978.1529 -7044.5213 153 -174504.7953 -172978.1529 154 139777.8833 -174504.7953 155 52678.7444 139777.8833 156 17438.7996 52678.7444 157 -5059.5465 17438.7996 158 -37389.2550 -5059.5465 159 -7115.1361 -37389.2550 160 130693.6677 -7115.1361 161 247870.8106 130693.6677 162 261059.3415 247870.8106 163 247849.6687 261059.3415 164 11834.1492 247849.6687 165 -297.5484 11834.1492 166 -5056.2529 -297.5484 167 -37353.4066 -5056.2529 168 -6186.6401 -37353.4066 169 -7080.6365 -6186.6401 170 -177075.9749 -7080.6365 171 139970.2411 -177075.9749 172 255902.0738 139970.2411 173 262848.6397 255902.0738 174 217742.7744 262848.6397 175 -5059.4642 217742.7744 176 -33631.0470 -5059.4642 177 -33048.6284 -33631.0470 178 -39103.4120 -33048.6284 179 -35312.1309 -39103.4120 180 -33331.2475 -35312.1309 181 -35900.1317 -33331.2475 182 -35700.6713 -35900.1317 183 -33948.9541 -35700.6713 184 -34381.0366 -33948.9541 185 -37118.4998 -34381.0366 186 -5976.2294 -37118.4998 187 -7093.7448 -5976.2294 188 -178260.7668 -7093.7448 189 140098.8532 -178260.7668 190 240533.9467 140098.8532 191 19189.3050 240533.9467 192 -5050.5098 19189.3050 193 -36909.9004 -5050.5098 194 -10901.9479 -36909.9004 195 -7107.6279 -10901.9479 196 139781.8136 -7107.6279 197 -4799.9257 139781.8136 198 7845.3632 -4799.9257 199 53266.6459 7845.3632 200 25597.5517 53266.6459 201 -5084.6181 25597.5517 202 -39052.2320 -5084.6181 203 -37164.8190 -39052.2320 204 16964.3181 -37164.8190 205 15335.4581 16964.3181 206 -1580.6319 15335.4581 207 -330.8313 -1580.6319 208 -1245.4368 -330.8313 209 -26671.5726 -1245.4368 210 -24232.8131 -26671.5726 211 -21453.6527 -24232.8131 212 -7533.1168 -21453.6527 213 -9410.3867 -7533.1168 214 -9152.1501 -9410.3867 215 30829.4584 -9152.1501 216 17113.0973 30829.4584 217 26609.2914 17113.0973 218 17681.2257 26609.2914 219 18289.6287 17681.2257 220 16985.8485 18289.6287 221 -38479.9818 16985.8485 222 21615.5558 -38479.9818 223 -1519.2978 21615.5558 224 -1804.3255 -1519.2978 225 -27418.8472 -1804.3255 226 -21474.3942 -27418.8472 227 53583.3154 -21474.3942 228 -9151.5882 53583.3154 229 20232.1426 -9151.5882 230 16963.9264 20232.1426 231 18367.8233 16963.9264 232 -26677.4970 18367.8233 233 -22329.6993 -26677.4970 234 -3054.3746 -22329.6993 235 -4320.2660 -3054.3746 236 -9147.4272 -4320.2660 237 16987.7905 -9147.4272 238 -41424.1582 16987.7905 239 18376.9622 -41424.1582 240 -29240.2684 18376.9622 241 -21590.0785 -29240.2684 242 -4735.5861 -21590.0785 243 -10084.9169 -4735.5861 244 -17525.3641 -10084.9169 245 -9146.5316 -17525.3641 246 17020.7058 -9146.5316 247 12010.2079 17020.7058 248 16965.0754 12010.2079 249 18373.4906 16965.0754 250 -26624.0306 18373.4906 251 -2766.5122 -26624.0306 252 -5430.6314 -2766.5122 253 -9139.1231 -5430.6314 254 16966.0318 -9139.1231 255 35345.9045 16966.0318 256 -2133.8446 35345.9045 257 -24009.2340 -2133.8446 258 -21919.8434 -24009.2340 259 -9140.5445 -21919.8434 260 16989.0825 -9140.5445 261 -41036.5446 16989.0825 262 20160.3731 -41036.5446 263 745.0774 20160.3731 264 -2826.9926 745.0774 265 -27414.6001 -2826.9926 266 -20521.9849 -27414.6001 267 -21520.7434 -20521.9849 268 -9144.1437 -21520.7434 269 16285.5589 -9144.1437 270 17008.3316 16285.5589 271 -42089.8303 17008.3316 272 18466.7302 -42089.8303 273 536.0135 18466.7302 274 950.6720 536.0135 275 -29418.2383 950.6720 276 -21378.5236 -29418.2383 277 -6363.9565 -21378.5236 278 177.4511 -6363.9565 279 -9143.8628 177.4511 280 15414.8125 -9143.8628 281 23274.8412 15414.8125 282 16964.6738 23274.8412 283 18372.8224 16964.6738 284 -26667.8846 18372.8224 285 -3092.5135 -26667.8846 286 -6367.0712 -3092.5135 287 -9149.1448 -6367.0712 288 17051.0349 -9149.1448 289 -38538.5964 17051.0349 290 -37984.2586 -38538.5964 291 21702.9173 -37984.2586 292 -8048.0187 21702.9173 293 -2563.8728 -8048.0187 294 -26649.5438 -2563.8728 295 -3114.9548 -26649.5438 296 -9148.8318 -3114.9548 297 21510.2352 -9148.8318 298 26737.5365 21510.2352 299 24465.2693 26737.5365 300 34710.8968 24465.2693 301 14043.7096 34710.8968 302 31707.6553 14043.7096 303 16965.9895 31707.6553 304 18375.0442 16965.9895 305 -26689.1382 18375.0442 306 -9145.1630 -26689.1382 307 16969.9571 -9145.1630 308 -30592.9841 16969.9571 309 18375.8513 -30592.9841 310 -26648.8823 18375.8513 311 -9141.0535 -26648.8823 312 16966.0939 -9141.0535 313 18379.0713 16966.0939 314 -26772.6413 18379.0713 315 -22133.0181 -26772.6413 316 -67453.7499 -22133.0181 317 -23154.5891 -67453.7499 318 -21352.1309 -23154.5891 319 5292.3718 -21352.1309 320 -7049.6788 5292.3718 321 -11136.7996 -7049.6788 322 -9153.6952 -11136.7996 323 16963.6130 -9153.6952 324 18468.2198 16963.6130 325 -2899.4911 18468.2198 326 -27781.3429 -2899.4911 327 -21917.8511 -27781.3429 328 -9155.2402 -21917.8511 329 18854.9420 -9155.2402 330 18701.6097 18854.9420 331 18734.9081 18701.6097 332 20607.1118 18734.9081 333 25678.7033 20607.1118 334 16965.9373 25678.7033 335 36524.8230 16965.9373 336 -3037.2008 36524.8230 337 -26662.4424 -3037.2008 338 -9142.0368 -26662.4424 339 45533.2814 -9142.0368 340 16965.8328 45533.2814 341 19534.1656 16965.8328 342 -3687.8868 19534.1656 343 -1855.5881 -3687.8868 344 -26668.5251 -1855.5881 345 -9139.0871 -26668.5251 346 17007.2156 -9139.0871 347 -41149.0371 17007.2156 348 17516.3231 -41149.0371 349 13003.7708 17516.3231 350 2518.4085 13003.7708 351 10379.0069 2518.4085 352 -26628.4330 10379.0069 353 -11111.1012 -26628.4330 354 -6826.2858 -11111.1012 355 -7369.5384 -6826.2858 356 -9141.8963 -7369.5384 357 17073.7591 -9141.8963 358 -46463.1364 17073.7591 359 -39084.5284 -46463.1364 360 20373.1199 -39084.5284 361 5311.7414 20373.1199 362 -2961.7350 5311.7414 363 -2606.0918 -2961.7350 364 -1801.2296 -2606.0918 365 -346.0918 -1801.2296 366 -276.1497 -346.0918 367 -3112.6701 -276.1497 368 -1148.7137 -3112.6701 369 6128.1256 -1148.7137 370 21973.0515 6128.1256 371 -4108.2105 21973.0515 372 -9135.6784 -4108.2105 373 58702.6752 -9135.6784 374 -1790.7359 58702.6752 375 -41022.5544 -1790.7359 376 7965.6545 -41022.5544 377 3116.4803 7965.6545 378 -3096.5854 3116.4803 379 7383.5677 -3096.5854 380 1313.8368 7383.5677 381 8856.7682 1313.8368 382 5182.2422 8856.7682 383 -3978.2731 5182.2422 384 -4271.8223 -3978.2731 385 7476.0284 -4271.8223 386 8751.7002 7476.0284 387 676.9163 8751.7002 388 -5706.7481 676.9163 389 12965.6691 -5706.7481 390 14485.8632 12965.6691 391 11802.9822 14485.8632 392 -26263.4451 11802.9822 393 -2628.9920 -26263.4451 394 6439.3599 -2628.9920 395 3001.3730 6439.3599 396 8972.1360 3001.3730 397 13468.2362 8972.1360 398 -75518.0470 13468.2362 399 -1026.0610 -75518.0470 400 -3687.6104 -1026.0610 401 4866.2730 -3687.6104 402 -10505.5715 4866.2730 403 102590.0857 -10505.5715 404 17072.0931 102590.0857 405 16859.9296 17072.0931 406 -34651.0483 16859.9296 407 179140.9522 -34651.0483 408 -65864.0831 179140.9522 409 42570.7665 -65864.0831 410 -136327.6304 42570.7665 411 73670.7658 -136327.6304 412 81076.8645 73670.7658 413 51272.4613 81076.8645 414 160407.0204 51272.4613 415 -23715.4828 160407.0204 416 -24239.7462 -23715.4828 417 168333.9309 -24239.7462 418 -25413.0401 168333.9309 419 -72861.4773 -25413.0401 420 -42221.3985 -72861.4773 421 36837.5340 -42221.3985 422 78426.2163 36837.5340 423 213439.2120 78426.2163 424 49037.2086 213439.2120 425 -29128.7692 49037.2086 426 50784.6062 -29128.7692 427 18434.2770 50784.6062 428 10938.3181 18434.2770 429 97818.2382 10938.3181 430 49147.3151 97818.2382 > 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/7fjxg1291404105.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/8fjxg1291404105.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/98swj1291404105.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/108swj1291404105.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/11tbd71291404105.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/12fbbu1291404105.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/13lcro1291404105.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/14elqr1291404105.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/15hm6f1291404105.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/16ee461291404105.tab") + } > > try(system("convert tmp/1j9zp1291404105.ps tmp/1j9zp1291404105.png",intern=TRUE)) character(0) > try(system("convert tmp/2cizs1291404105.ps tmp/2cizs1291404105.png",intern=TRUE)) character(0) > try(system("convert tmp/3cizs1291404105.ps tmp/3cizs1291404105.png",intern=TRUE)) character(0) > try(system("convert tmp/4cizs1291404105.ps tmp/4cizs1291404105.png",intern=TRUE)) character(0) > try(system("convert tmp/55ayv1291404105.ps tmp/55ayv1291404105.png",intern=TRUE)) character(0) > try(system("convert tmp/65ayv1291404105.ps tmp/65ayv1291404105.png",intern=TRUE)) character(0) > try(system("convert tmp/7fjxg1291404105.ps tmp/7fjxg1291404105.png",intern=TRUE)) character(0) > try(system("convert tmp/8fjxg1291404105.ps tmp/8fjxg1291404105.png",intern=TRUE)) character(0) > try(system("convert tmp/98swj1291404105.ps tmp/98swj1291404105.png",intern=TRUE)) character(0) > try(system("convert tmp/108swj1291404105.ps tmp/108swj1291404105.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 15.034 2.148 31.956